Source author record

Joop Schaye

Joop Schaye 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 astro-ph.IM astro-ph.SR

Catalog footprint

What is connected

183works

6topics

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.

preprint2026arXiv

Cosmological back-reaction of baryons on dark matter in the CAMELS simulations

Baryonic processes such as radiative cooling and feedback from massive stars and active galactic nuclei (AGN) directly redistribute baryons in the Universe but also indirectly redistribute dark matter due to changes in the gravitational potential. In this work, we investigate this "back-reaction" of baryons on dark matter using thousands of cosmological hydrodynamic simulations from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project, including parameter variations in the SIMBA, IllustrisTNG, ASTRID, and Swift-EAGLE galaxy formation models. Matching haloes to corresponding N-body (dark matter-only) simulations, we find that virial masses decrease owing to the ejection of baryons by feedback. Relative to N-body simulations, halo profiles show an increased dark matter density in the center (due to radiative cooling) and a decrease in density farther out (due to feedback), with both effects being strongest in SIMBA (> 450% increase at r < 0.01 Rvir). The clustering of dark matter strongly responds to changes in baryonic physics, with dark matter power spectra in some simulations from each model showing as much as 20% suppression or increase in power at k ~ 10 h/Mpc relative to N-body simulations. We find that the dark matter back-reaction depends intrinsically on cosmology (Omega_m and sigma_8) at fixed baryonic physics, and varies strongly with the details of the feedback implementation. These results emphasize the need for marginalizing over uncertainties in baryonic physics to extract cosmological information from weak lensing surveys as well as their potential to constrain feedback models in galaxy evolution.

preprint2026arXiv

Luminosity-Dependent Assembly Bias of Central Galaxies from Weak Lensing and Clustering

Assembly bias, which is the variation in halo clustering at fixed mass driven by formation history, has long been predicted by numerical simulations but remains difficult to confirm observationally. Previous studies have reported evidence for halo assembly bias by dividing samples according to galaxy stellar mass using various methods. In this work, we present observational measurements of halo assembly bias based on the luminosity of spectroscopically confirmed brightest cluster galaxies (BCGs). Using cluster catalogs and shear measurements from the DESI Legacy Imaging Surveys, we employ a mass-dependent halo-bias model to disentangle halo bias from its underlying mass dependence in galaxy-galaxy lensing and clustering measurements. We confirm that brighter BCGs are less strongly clustered on large scales, with a relative bias ratio deviating from unity at the $\sim3σ$ level, suggesting the presence of assembly bias. Similar qualitative trends are also found in the FLAMINGO and MillenniumTNG hydrodynamical simulations, strengthening the connection between galaxy luminosity and halo formation history.

preprint2026arXiv

Modelling the evolution and influence of dust in cosmological simulations that include the cold phase of the interstellar medium

While marginal in mass terms, dust grains play an outsized role in both the physics and observation of the interstellar medium (ISM). However, explicit modelling of this ISM constituent remains uncommon in large cosmological simulations. In this work, we present a model for the life-cycle of dust in the ISM that couples to the forthcoming COLIBRE galaxy formation model, which explicitly simulates the cold ISM. We follow 6 distinct grain types: 3 chemical species, including carbon and two silicate grains, with 2 size bins each. Our dust model accounts for seeding of grains from stellar ejecta, self-consistent element-by-element metal yields and growth by accretion, grain size transfer (shattering and coagulation) and destruction of dust by thermal sputtering in the ISM. We detail the calibration of this model, particularly the use of a clumping factor, to account for unresolved gas clouds in which dust readily evolves. We present a fiducial run in a 25$^3$~cMpc$^3$ cosmological volume that displays good agreement with observations of the cosmic evolution of dust density, as well as the $z=0$ galaxy dust mass function and dust scaling relations. We highlight known tensions between observational datasets of the dust-to-gas ratio as a function of metallicity depending on which metallicity calibrator is used; our model favours higher-normalisation metallicity calibrators, which agree with the observations within 0.1~dex for stellar masses $>10^9 \; {\rm M_\odot}$. We compare the grain size distribution to observations of local galaxies, and find that our simulation suggests a higher concentration of small grains, associated with more diffuse ISM and the warm-neutral medium (WNM), which both play a key role in boosting H$_2$ content. Putting these results and modelling approaches in context, we set the stage for upcoming insights into the dusty ISM of galaxies using the COLIBRE simulations.

preprint2026arXiv

MUSEQuBES: Investigating the Physical and Chemical Properties of Circumgalactic Gas Around Ly$α$ Emitters at $z \approx 3.3$

We present a detailed study of the physical properties of absorbers associated with $z\approx3.3$ \lya\ emitters (LAEs) from the MUSEQuBES survey. Using \HI\ and ionic column densities derived from Voigt profile fitting, we determine the density and metallicity of 75 absorbers associated with 59 LAEs through a custom-built Bayesian framework. The overall sample shows a median density $\log\,(n_{\rm H}/{\rm cm^{-3}})=-2.7\pm0.7$ and metallicity \met$=-1.6^{+1.2}_{-0.9}$ with $\approx15\%$ ($\approx 40\%$ including upper limits) of systems showing metallicities consistent with the IGM at these redshifts (\met$\lesssim -2.8$). Intriguingly, absorbers with $\log\, N(\rm HI)/{\rm cm^{-2}} \gtrsim 16.5$, corresponding to an overdensity of $\sim 100$ at $z=3$, exhibit a bimodal metallicity distribution with peaks at \met$=-3.8\pm 0.2$ and $-1.8\pm 0.6$. The latter are observed at large impact parameters ($\gtrsim150$~pkpc) and often exhibit low-ionization species (e.g., \SiII, \AlII). We interpret the former as pristine inflowing gas from the cosmic web, while the latter likely traces metal-enriched CGM associated with outflows from faint galaxies not detected in Ly$α$ emission. We find no significant correlations between absorber metallicity or density and host galaxy properties, including redshift, impact parameter, SFR, Ly$α$ luminosity, or environment. Absorbers separated by $\lesssim500$~\kms\ show $\sim$1~dex metallicity and $\sim$0.5~dex density variations, indicating physical and chemical inhomogeneity of the medium around these LAEs.

preprint2022arXiv

Galaxy cluster aperture masses are more robust to baryonic effects than 3D halo masses

Systematic uncertainties in the mass measurement of galaxy clusters limit the cosmological constraining power of future surveys that will detect more than $10^5$ clusters. Previously, we argued that aperture masses can be inferred more accurately and precisely than 3D masses without loss of cosmological constraining power. Here, we use the Baryons and Haloes of Massive Systems (BAHAMAS) cosmological, hydrodynamical simulations to show that aperture masses are also less sensitive to changes in mass caused by galaxy formation processes. For haloes with $m_\mathrm{200m,dmo} > 10^{14} \, h^{-1} \, \mathrm{M}_\odot$, binned by their 3D halo mass, baryonic physics affects aperture masses and 3D halo masses similarly when measured within apertures similar to the halo virial radius, reaching a maximum reduction of $\approx 3 \, \%$. For lower-mass haloes, $10^{13.5} < m_\mathrm{200m,dmo} / (h^{-1} \, \mathrm{M}_\odot) < 10^{14}$, and aperture sizes $\sim 1 \, h^{-1} \, \mathrm{cMpc}$, representative of weak lensing observations, the aperture mass is consistently reduced less ($\lesssim 5 \, \%$) than the 3D halo mass ($\lesssim 10 \, \%$ for $m_\mathrm{200m}$). The halo mass reduction evolves only slightly, by up to $2$ percentage points, between redshift 0.25 and 1 for both the aperture mass and $m_\mathrm{200m}$. Varying the strength of the simulated feedback so the mean simulated hot gas fraction covers the observed scatter inferred from X-ray observations, we find that the aperture mass is consistently less biased than the 3D halo mass, by up to $2 \, $ percentage points at $m_\mathrm{200m,dmo} = 10^{14} \, h^{-1} \, \mathrm{M}_\odot$. Therefore, cluster aperture mass calibrations provide a fruitful path forward for future cluster surveys to reduce their sensitivity to systematic uncertainties.

preprint2022arXiv

Galaxy cluster photons alter the ionisation state of the nearby warm-hot intergalactic medium

The physical properties of the faint and extremely tenuous plasma in the far outskirts of galaxy clusters, the circumgalactic media of normal galaxies, and filaments of the cosmic web, remain one of the biggest unknowns in our story of large-scale structure evolution. Modelling the spectral features due to emission and absorption from this very diffuse plasma poses a challenge, as both collisional and photo-ionisation processes must be accounted for. In this paper, we study the ionisation by photons emitted by the intra-cluster medium in addition to the photo-ionisation by the cosmic UV/X-ray background on gas in the vicinity of galaxy clusters. For near massive clusters such as A2029, the ionisation parameter can no longer describe the ionisation balance uniquely. The ionisation fractions (in particular of C IV, C V, C VI, N VII, O VI, O VII, O VIII, Ne VIII, Ne IX, and Fe XVII) obtained by taking into account the photoionisation by the cosmic background are either an upper or lower limit to the ionisation fraction calculated as a function of distance from the emission from the cluster. Using a toy model of a cosmic web filament, we predict how the cluster illumination changes the column densities for two different orientations of the line of sight. For lines of sight passing close to the cluster outskirts, O VI can be suppressed by a factor of up to $4.5$, O VII by a factor of $2.2$, C V by a factor of $3$, and Ne VIII can be boosted by a factor of $2$, for low density gas.

preprint2022arXiv

How baryons can significantly bias cluster count cosmology

We quantify two main pathways through which baryonic physics biases cluster count cosmology. We create mock cluster samples that reproduce the baryon content inferred from X-ray observations. We link clusters to their counterparts in a dark matter-only universe, whose abundances can be predicted robustly, by assuming the dark matter density profile is not significantly affected by baryons. We derive weak lensing halo masses and infer the best-fitting cosmological parameters $Ω_\mathrm{m}$, $S_8=σ_8(Ω_\mathrm{m}/0.3)^{0.2}$, and $w_0$ from the mock cluster sample. We find that because of the need to accommodate the change in the density profile due to the ejection of baryons, weak lensing mass calibrations are only unbiased if the concentration is left free when fitting the reduced shear with NFW profiles. However, even unbiased total mass estimates give rise to biased cosmological parameters if the measured mass functions are compared with predictions from dark matter-only simulations. This bias dominates for haloes with $m_\mathrm{500c} < 10^{14.5} \, h^{-1} \, \mathrm{M_\odot}$. For a stage IV-like cluster survey without mass estimation uncertainties, an area $\approx 15000 \, \mathrm{deg^2}$ and a constant mass cut of $m_\mathrm{200m,min} = 10^{14} \, h^{-1} \, \mathrm{M_\odot}$, the biases are $-11 \pm 1 \, \%$ in $Ω_\mathrm{m}$, $-3.29 \pm 0.04 \, \%$ in $S_8$, and $9 \pm 1.5 \, \%$ in $w_0$. The statistical significance of the baryonic bias depends on how accurately the actual uncertainty on individual cluster mass estimates is known. We suggest that rather than the total halo mass, the (re-scaled) dark matter mass inferred from the combination of weak lensing and observations of the hot gas, should be used for cluster count cosmology.

preprint2022arXiv

Spin-driven jet feedback in idealised simulations of galaxy groups and clusters

We implement a black hole spin evolution and jet feedback model into SWIFT, a smoothed particle hydrodynamics code. The jet power is determined self-consistently assuming Bondi accretion, using a realistic, spin-dependant efficiency. The jets are launched along the spin axis of the black hole, resulting in natural reorientation and precession. We apply the model to idealised simulations of galaxy groups and clusters, finding that jet feedback successfully quenches gas cooling and star formation in all systems. Our group-size halo ($M_\mathrm{200}=10^{13}$ $\mathrm{M}_\odot$) is quenched by a strong jet episode triggered by a cooling flow, and it is kept quenched by a low-power jet fed from hot halo accretion. In more massive systems ($M_\mathrm{200}\geq 10^{14}$ $\mathrm{M}_\odot$), hot halo accretion is insufficient to quench the galaxies, or to keep them quenched after the first cooling episode. These galaxies experience multiple episodes of gas cooling, star formation and jet feedback. In the most massive galaxy cluster that we simulate ($M_\mathrm{200}=10^{15}$ $\mathrm{M}_\odot$), we find peak cold gas masses of $10^{10}$ $\mathrm{M}_\odot$ and peak star formation rates of a few times $100$ $\mathrm{M}_\odot\mathrm{yr}^{-1}$. These values are achieved during strong cooling flows, which also trigger the strongest jets with peak powers of $10^{47}$ $\mathrm{erg}\hspace{0.3mm}\mathrm{s}^{-1}$. These jets subsequently shut off the cooling flows and any associated star formation. Jet-inflated bubbles draw out low-entropy gas that subsequently forms dense cooling filaments in their wakes, as seen in observations.

preprint2022arXiv

The effects of local stellar radiation and dust depletion on non-equilibrium interstellar chemistry

Interstellar chemistry is important for galaxy formation, as it determines the rate at which gas can cool, and enables us to make predictions for observable spectroscopic lines from ions and molecules. We explore two central aspects of modelling the chemistry of the interstellar medium (ISM): (1) the effects of local stellar radiation, which ionises and heats the gas, and (2) the depletion of metals onto dust grains, which reduces the abundance of metals in the gas phase. We run high-resolution (400 M$_\odot$ per baryonic particle) simulations of isolated disc galaxies, from dwarfs to Milky Way-mass, using the FIRE galaxy formation models together with the CHIMES non-equilibrium chemistry and cooling module. In our fiducial model, we couple the chemistry to the stellar fluxes calculated from star particles using an approximate radiative transfer scheme, and we implement an empirical density-dependent prescription for metal depletion. For comparison, we also run simulations with a spatially uniform radiation field, and without metal depletion. Our fiducial model broadly reproduces observed trends in HI and H2 mass with stellar mass, and in line luminosity versus star formation rate for [CII] 158$μ$m, [OI] 63$μ$m, [OIII] 88$μ$m, [NII] 122$μ$m and H$α$ 6563A. Our simulations with a uniform radiation field predict fainter luminosities, by up to an order of magnitude for [OIII] 88$μ$m and H$α$ 6563A, while ignoring metal depletion increases the luminosity of carbon and oxygen lines by a factor $\approx$2. However, the overall evolution of the galaxy is not strongly affected by local stellar fluxes or metal depletion, except in dwarf galaxies where the inclusion of local fluxes leads to weaker outflows and hence higher gas fractions.

preprint2022arXiv

The importance of black hole repositioning for galaxy formation simulations

Active galactic nucleus (AGN) feedback from accreting supermassive black holes (SMBHs) is an essential ingredient of galaxy formation simulations. The orbital evolution of SMBHs is affected by dynamical friction that cannot be predicted self-consistently by contemporary simulations of galaxy formation in representative volumes. Instead, such simulations typically use a simple "repositioning" of SMBHs, but the effects of this approach on SMBH and galaxy properties have not yet been investigated systematically. Based on a suite of smoothed particle hydrodynamics simulations with the SWIFT code and a Bondi-Hoyle-Lyttleton subgrid gas accretion model, we investigate the impact of repositioning on SMBH growth and on other baryonic components through AGN feedback. Across at least a factor ~1000 in mass resolution, SMBH repositioning (or an equivalent approach) is a necessary prerequisite for AGN feedback; without it, black hole growth is negligible. Limiting the effective repositioning speed to $\lesssim$ 10 km/s delays the onset of AGN feedback and severely limits its impact on stellar mass growth in the centre of massive galaxies. Repositioning has three direct physical consequences. It promotes SMBH mergers and thus accelerates their initial growth. In addition, it raises the peak density of the ambient gas and reduces the SMBH velocity relative to it, giving a combined boost to the accretion rate that can reach many orders of magnitude. Our results suggest that a more sophisticated and/or better calibrated treatment of SMBH repositioning is a critical step towards more predictive galaxy formation simulations.

preprint2022arXiv

The importance of the way in which supernova energy is distributed around young stellar populations in simulations of galaxies

Supernova (SN) feedback plays a crucial role in simulations of galaxy formation. Because blastwaves from individual SNe occur on scales that remain unresolved in modern cosmological simulations, SN feedback must be implemented as a subgrid model. Differences in the manner in which SN energy is coupled to the local interstellar medium and in which excessive radiative losses are prevented have resulted in a zoo of models used by different groups. However, the importance of the selection of resolution elements around young stellar particles for SN feedback has largely been overlooked. In this work, we examine various selection methods using the smoothed particle hydrodynamics code SWIFT. We run a suite of isolated disk galaxy simulations of a Milky Way-mass galaxy and small cosmological volumes, all with the thermal stochastic SN feedback model used in the EAGLE simulations. We complement the original mass-weighted neighbour selection with a novel algorithm guaranteeing that the SN energy distribution is as close to isotropic as possible. Additionally, we consider algorithms where the energy is injected into the closest, least dense, or most dense neighbour. We show that different neighbour-selection strategies cause significant variations in star formation rates, gas densities, wind mass loading factors, and galaxy morphology. The isotropic method results in more efficient feedback than the conventional mass-weighted selection. We conclude that the manner in which the feedback energy is distributed among the resolution elements surrounding a feedback event is as important as changing the amount of energy by factors of a few.

preprint2022arXiv

The interplay between AGN feedback and precipitation of the intracluster medium in simulations of galaxy groups and clusters

Using high-resolution hydrodynamical simulations of galaxy clusters, we study the interaction between the brightest cluster galaxy, its supermassive black hole (BH) and the intracluster medium (ICM). We create initial conditions for which the ICM is in hydrostatic equilibrium within the gravitational potential from the galaxy and an NFW dark matter halo. Two free parameters associated with the thermodynamic profiles determine the cluster gas fraction and the central temperature, where the latter can be used to create cool-core or non-cool-core systems. Our simulations include radiative cooling, star formation, BH accretion, and stellar and active galactic nucleus (AGN) feedback. Even though the energy of AGN feedback is injected thermally and isotropically, it leads to anisotropic outflows and buoyantly rising bubbles. We find that the BH accretion rate (BHAR) is highly variable and only correlates strongly with the star formation rate (SFR) and the ICM when it is averaged over more than $1~\rm Myr$. We generally find good agreement with the theoretical precipitation framework. In $10^{13}~\rm M_\odot$ haloes, AGN feedback quenches the central galaxy and converts cool-core systems into non-cool-core systems. In contrast, higher-mass, cool-core clusters evolve cyclically. Episodes of high BHAR raise the entropy of the ICM out to the radius where the ratio of the cooling time and the local dynamical time $t_{\rm cool}/t_{\rm dyn} > 10$, thus suppressing condensation and, after a delay, the BHAR. The corresponding reduction in AGN feedback allows the ICM to cool and become unstable to precipitation, thus initiating a new episode of high SFR and BHAR.

preprint2022arXiv

The MUSE Extremely Deep Field: a first panoramic view of an Mg II emitting intragroup medium

Using the exquisite MUSE eXtremely Deep Field data, we report the discovery of an MgII emission nebula with an area above a 2$σ$ significance level of 1000 proper kpc$^2$, providing the first panoramic view of the spatial distribution of magnesium in the intragroup medium of a low mass group of five star-forming galaxies at z=1.31. The galaxy group members are separated by less than 50 physical kpc in projection and $\approx$120 km/s in velocity space. The most massive galaxy has a stellar mass of 10$^{9.35}$ M$_\odot$ and shows an MgII P-Cygni line profile indicating the presence of an outflow, which is consistent with the spatially resolved spectral analysis showing $\approx+$120 km/s shift of the MgII emission lines with respect to the systemic redshift. The other galaxies are less massive and only show MgII in emission. The detected MgII nebula has a maximal projected extent of $\approx$70 kpc including a low surface brightness (2 $\times$ 10$^{-19}$ erg/s/cm$^{2}$/arcsec$^{2}$) gaseous bridge between two subgroups of galaxies. The presence of absorption features in the spectrum of a background galaxy located at an impact parameter of 19 kpc from the closest galaxy of the group indicates the presence of gas enriched in magnesium even beyond the detected nebula seen in emission, suggesting that we are observing the tip of a larger intragroup medium. The observed MgII velocity gradient suggests an overall rotation of the structure along the major axis of the most massive galaxy. Our MUSE data also reveal extended Fe II* emission in the vicinity of the most massive galaxy, aligned with its minor axis. Extended [OII] emission is found around the galaxy group members and at the location of the MgII bridge. Our results suggest that both tidal stripping effects from galaxy interactions and outflows are enriching the intragroup medium of this system.

preprint2022arXiv

The MUSE eXtremely Deep Field: Individual detections of Lyα haloes around rest-frame UV-selected galaxies at z~2.9-4.4

Hydrogen Ly$α$ haloes (LAHs) are commonly used as a tracer of the circumgalactic medium (CGM) at high redshifts. In this work, we aim to explore the existence of Ly$α$ haloes around individual UV-selected galaxies, rather than around Ly$α$ emitters (LAEs), at high redshifts. Our sample was continuum-selected with F775W<=27.5, and spectroscopic redshifts were assigned or constrained for all the sources thanks to the deepest (100- to 140-hour) existing Very Large Telescope (VLT)/Multi-Unit Spectroscopic Explorer (MUSE) data with adaptive optics. The final sample includes 21 galaxies that are purely F775W-magnitude selected within the redshift range z=2.9-4.4 and within a UV magnitude range -20<=M1500<= -18, thus avoiding any bias toward LAEs. We tested whether galaxy's Ly$α$ emission is significantly more extended than the MUSE PSF-convolved continuum component. We find 17 LAHs and four non-LAHs. We report the first individual detections of extended Ly$α$ emission around non-LAEs. The Ly$α$ halo fraction is thus as high as $81.0^{+10.3}_{-11.2}$%, which is close to that for LAEs at z=3-6 in the literature. This implies that UV-selected galaxies generally have a large amount of hydrogen in their CGM. We derived the mean surface brightness (SB) profile for our LAHs with cosmic dimming corrections and find that Ly$α$ emission extends to 5.4 arcsec (~40 physical kpc at the midpoint redshift z=3.6) above the typical 1$σ$ SB limit. The incidence rate of surrounding gas detected in Ly$α$ per one-dimensional line of sight per unit redshift, dn/dz, is estimated to be $0.76^{+0.09}_{-0.09}$ for galaxies with M1500<= -18 mag at z~3.7. Assuming that Ly$α$ emission and absorption arise in the same gas, this suggests, based on abundance matching, that LAHs trace the same gas as damped Ly$α$ systems (DLAs) and sub-DLAs.

preprint2022arXiv

The warm-hot circumgalactic medium around EAGLE-simulation galaxies and its detection prospects with X-ray line emission

We estimate the detectability of X-ray metal-line emission from the circumgalactic medium (CGM) of galaxies over a large halo mass range ($\mathrm{M}_{\mathrm{200c}} = 10^{11.5}$-$10^{14.5}\,\mathrm{M}_{\odot}$) using the EAGLE simulations. With the XRISM Resolve instrument, a few bright (K-$α$ or Fe L-shell) lines from $\mathrm{M}_{\mathrm{200c}} \gtrsim 10^{13}\,\mathrm{M}_{\odot}$ haloes should be detectable. Using the Athena X-IFU or the Lynx Main Array, emission lines (especially from O VIII and O VII) from the inner CGM of $\mathrm{M}_{\mathrm{200c}} \gtrsim10^{12.5}\,\mathrm{M}_{\odot}$ haloes become detectable, and intragroup and intracluster gas will be detectable out to the virial radius. With the Lynx Ultra-high Resolution Array, the inner CGM of haloes hosting $\mathrm{L}_{*}$ galaxies is accessible. These estimates do assume long exposure times ($\sim 1 \,$Ms) and large spatial bins ($\sim1$-$10\,\mathrm{arcmin}^{2}$). This emission is dominated by collisionally ionized (CI) gas, and tends to come from halo centres. The emission is biased towards temperatures close to the maximum emissivity temperature for CI gas ($\mathrm{T}_\mathrm{peak}$), and high densities and metallicities. However, for the K-$α$ lines, emission can come from hotter gas in haloes where the virialized, volume-filling gas is hotter than $\mathrm{T}_\mathrm{peak}$. Trends of emission with halo mass can largely be explained by differences in virial temperature. Differences in the mass trends of K-$α$, He-$α$-like, and Fe L-shell lines mirror differences in their emissivities as a function of temperature. We conclude that upcoming X-ray missions will open up a new window onto the hot CGM.

preprint2022arXiv

Why are we still using 3D masses for cluster cosmology?

The abundance of clusters of galaxies is highly sensitive to the late-time evolution of the matter distribution, since clusters form at the highest density peaks. However, the 3D cluster mass cannot be inferred without deprojecting the observations, introducing model-dependent biases and uncertainties due to the mismatch between the assumed and the true cluster density profile and the neglected matter along the sightline. Since projected aperture masses can be measured directly in simulations and observationally through weak lensing, we argue that they are better suited for cluster cosmology. Using the Mira-Titan suite of gravity-only simulations, we show that aperture masses correlate strongly with 3D halo masses, albeit with large intrinsic scatter due to the varying matter distribution along the sightline. Nonetheless, aperture masses can be measured $\approx 2-3$ times more precisely from observations, since they do not require assumptions about the density profile and are only affected by the shape noise in the weak lensing measurements. We emulate the cosmology dependence of the aperture mass function directly with a Gaussian process. Comparing the cosmology sensitivity of the aperture mass function and the 3D halo mass function for a fixed survey solid angle and redshift interval, we find the aperture mass sensitivity is higher for $Ω_\mathrm{m}$ and $w_a$, similar for $σ_8$, $n_\mathrm{s}$, and $w_0$, and slightly lower for $h$. With a carefully calibrated aperture mass function emulator, cluster cosmology analyses can use cluster aperture masses directly, reducing the sensitivity to model-dependent mass calibration biases and uncertainties.

preprint2021arXiv

Baryonic mass budgets for haloes in the EAGLE simulation, including ejected and prevented gas

Feedback processes are expected to shape galaxy evolution by ejecting gas from galaxies and their associated dark matter haloes, and also by preventing diffuse gas from ever being accreted. We present predictions from the EAGLE simulation project for the mass budgets associated with "ejected" and "prevented" gas, as well as for ejected metals. We find that most of the baryons that are associated with haloes of mass $10^{11} < M_{200} \, /\mathrm{M_\odot} < 10^{13}$ at $z=0$ have been ejected beyond the virial radius after having been accreted. When the gas ejected from satellites (and their progenitors) is accounted for, the combined ejected mass represents half of the total baryon budget even in the most massive simulated galaxy clusters ($M_{200} \approx 10^{14.5} \, \mathrm{M_\odot}$), with the consequence that the total baryon budget exceeds the cosmic average if ejected gas is included. We find that gas is only prevented from being accreted onto haloes for $M_{200} < 10^{12} \, \mathrm{M_\odot}$, and that this component accounts for about half the total baryon budget for $M_{200} < 10^{11} \, \mathrm{M_\odot}$, with ejected gas making up most of the remaining half. For metals, most of the mass that is not locked into stars has been ejected beyond the virial radius, at least for $M_{200} < 10^{13} \, \mathrm{M_\odot}$. Finally, within the virial radius we find that most of the mass in the circum-galactic medium (CGM) has not passed through the ISM of a progenitor galaxy, for all halo masses and redshifts. About half of the CGM within half the virial radius has passed through the ISM in the past, however.

preprint2021arXiv

How gas flows shape the stellar-halo mass relation in the EAGLE simulation

The difference in shape between the observed galaxy stellar mass function and the predicted dark matter halo mass function is generally explained primarily by feedback processes. Feedback can shape the stellar-halo mass (SHM) relation by driving gas out of galaxies, by modulating the first-time infall of gas onto galaxies (i.e., preventative feedback), and by instigating fountain flows of recycled wind material. We present and apply a method to disentangle these effects for hydrodynamical simulations of galaxy formation. We build a model of linear coupled differential equations that by construction reproduces the flows of gas onto and out of galaxies and haloes in the EAGLE cosmological simulation. By varying individual terms in this model, we isolate the relative effects of star formation, ejection via outflow, first-time inflow and wind recycling on the SHM relation. We find that for halo masses $M_{200} < 10^{12} \, \mathrm{M_\odot}$ the SHM relation is shaped primarily by a combination of ejection from galaxies and haloes, while for larger $M_{200}$ preventative feedback is also important. The effects of recycling and the efficiency of star formation are small. We show that if, instead of $M_{200}$, we use the cumulative mass of dark matter that fell in for the first time, the evolution of the SHM relation nearly vanishes. This suggests that the evolution is due to the definition of halo mass rather than to an evolving physical efficiency of galaxy formation. Finally, we demonstrate that the mass in the circum-galactic medium is much more sensitive to gas flows, especially recycling, than is the case for stars and the interstellar medium.

preprint2021arXiv

Morphological and Rotation Structures of Circumgalactic Mg II Gas in the EAGLE Simulation and the Dependence on Galaxy Properties

Low-ionization-state Mg II gas has been extensively studied in quasar sightline observations to understand the cool, $\sim$$10^4$ K gas in the circumgalactic medium. Motivated by recent observations showing that the Mg II gas around low-redshift galaxies has significant angular momentum, we use the high-resolution EAGLE cosmological simulation to analyze the morphological and rotation structures of the $z\approx0.3$ circumgalactic Mg II gas and examine how they change with the host galaxy properties. Around star-forming galaxies, we find that the Mg II gas has an axisymmetric instead of a spherical distribution, and the axis of symmetry aligns with that of the Mg II gas rotation. A similar rotating structure is less commonly found in the small sample of simulated quiescent galaxies. We also examine how often Mg II gas around galaxies selected using a line-of-sight velocity cut includes gas physically outside of the virial radius ($r_\mathrm{vir}$). For example, we show that at an impact parameter of 100 pkpc, a $\pm500$ km s$^{-1}$ velocity cut around galaxies with stellar masses of $10^9$-$10^{9.5}\mathrm{M_\odot}$ ($10^{10}$-$10^{10.5}\mathrm{M_\odot}$) selects Mg II gas beyond the virial radius 80% (6%) of the time. Because observers typically select Mg II gas around target galaxies using such a velocity cut, we discuss how this issue affects the study of circumgalactic Mg II gas properties, including the detection of corotation. While the corotating Mg II gas generally extends beyond $0.5r_\mathrm{vir}$, the Mg II gas outside of the virial radius contaminates the corotation signal and makes observers less likely to conclude that gas at large impact parameters (e.g., $\gtrsim0.25 r_\mathrm{vir}$) is corotating.

preprint2021arXiv

SEAGLE--II: Constraints on feedback models in galaxy formation from massive early type strong lens galaxies

We use nine different galaxy formation scenarios in ten cosmological simulation boxes from the EAGLE suite of ΛCDM hydrodynamical simulations to assess the impact of feedback mechanisms in galaxy formation and compare these to observed strong gravitational lenses. To compare observations with simulations, we create strong lenses with $M_\star$ > $10^{11}$ $M_\odot$ with the appropriate resolution and noise level, and model them with an elliptical power-law mass model to constrain their total mass density slope. We also obtain the mass-size relation of the simulated lens-galaxy sample. We find significant variation in the total mass density slope at the Einstein radius and in the projected stellar mass-size relation, mainly due to different implementations of stellar and AGN feedback. We find that for lens selected galaxies, models with either too weak or too strong stellar and/or AGN feedback fail to explain the distribution of observed mass-density slopes, with the counter-intuitive trend that increasing the feedback steepens the mass density slope around the Einstein radius ($\approx$ 3-10 kpc). Models in which stellar feedback becomes inefficient at high gas densities, or weaker AGN feedback with a higher duty cycle, produce strong lenses with total mass density slopes close to isothermal (i.e. -d log(ρ)/d log(r) $\approx$ 2.0) and slope distributions statistically agreeing with observed strong lens galaxies in SLACS and BELLS. Agreement is only slightly worse with the more heterogeneous SL2S lens galaxy sample. Observations of strong-lens selected galaxies thus appear to favor models with relatively weak feedback in massive galaxies.

preprint2021arXiv

The MUSE Extremely Deep Field: Evidence for SFR-induced cores in dark-matter dominated galaxies at z=1

Disc-halo decomposition on rotationally supported star-forming galaxies (SFGs) at $z>1$ are often limited to massive galaxies ($M_\star>10^{10}~M_\odot$) and rely on either deep Integral Field Spectroscopy data or stacking analyses. We present a study of the dark matter (DM) content of nine $z\approx1$ SFGs selected Using the brightest [OII] emitters in the deepest Multi-Unit Spectrograph Explorer (MUSE) field to date, namely the 140hr MUSE Extremely Deep Field, we perform disk-halo decompositions on 9 low-mass SFGs (with $10^{8.5}<M_\star<10^{10.5}~M_\odot$) using a novel 3D modeling approach, which together with the exquisite S/N allows us to measure individual rotation curves to $3\times R_e$. The disk-halo decomposition includes a stellar, DM, gas, and occasionally a bulge component. The DM component primarily uses the generalized $α,β,γ$ profile or a Navarro-Frenk-White (NFW) profile. The disk stellar masses $M_\star$ obtained from the [OII] disk-halo decomposition agree with the values inferred from the spectral energy distributions. While the rotation curves show diverse shapes, ranging from rising to declining at large radii, the DM fractions within the half-light radius $f_{\rm DM}(<R_e)$ are found to be 60\% to 95\%, extending to lower masses (densities) recent results on massive SFGs with $M_\star>10^{10}~M_\odot$. The DM halos show constant surface densities of $\sim100~M_\odot$ pc$^{-2}$. Half of the sample shows a strong preference for cored over cuspy DM profiles. The presence of DM cores appears to be related to galaxies with stellar-to-halo mass $\log M_\star/M_{\rm vir}\approx-2.5$. In addition, the cuspiness of the DM profiles is found to be a strong function of the recent star-formation activity. Both of these results are interpreted as evidence for feedback-induced core formation in the Cold Dark Matter context.

preprint2021arXiv

The surprising accuracy of isothermal Jeans modelling of self-interacting dark matter density profiles

Recent claims of observational evidence for self-interacting dark matter (SIDM) have relied on a semi-analytic method for predicting the density profiles of galaxies and galaxy clusters containing SIDM. We present a thorough description of this method, known as isothermal Jeans modelling, and then test it with a large ensemble of haloes taken from cosmological simulations. Our simulations were run with cold and collisionless dark matter (CDM) as well as two different SIDM models, all with dark matter only variants as well as versions including baryons and relevant galaxy formation physics. Using a mix of different box sizes and resolutions, we study haloes with masses ranging from 3e10 to 3e15 Msun. Overall, we find that the isothermal Jeans model provides as accurate a description of simulated SIDM density profiles as the Navarro-Frenk-White profile does of CDM halos. We can use the model predictions, compared with the simulated density profiles, to determine the input DM-DM scattering cross-sections used to run the simulations. This works especially well for large cross-sections, while with CDM our results tend to favour non-zero (albeit fairly small) cross-sections, driven by a bias against small cross-sections inherent to our adopted method of sampling the model parameter space. The model works across the whole halo mass range we study, although including baryons leads to DM profiles of intermediate-mass (10^12 - 10^13 Msun) haloes that do not depend strongly on the SIDM cross-section. The tightest constraints will therefore come from lower and higher mass haloes: dwarf galaxies and galaxy clusters.

preprint2020arXiv

An EAGLE view of the missing baryons

Context. A significant fraction of the predicted baryons remains undetected in the local universe. We adopted the common assumption that a large fraction of the missing baryons corresponds to the hot (log T(K) = 5.5-7) phase of the Warm Hot Intergalactic Medium (WHIM). We base our missing baryons search on the scenario whereby the WHIM has been heated up via accretion shocks and galactic outflows, and is concentrated towards the filaments of the Cosmic Web. Aims. Our aim is to improve the observational search of the poorly detected hot WHIM. Methods. We detect the filamentary structure within the EAGLE simulation by applying the Bisous formalism to the galaxy distribution. In addition, we use the MMF/NEXUS+ classification of the large scale environment of the dark matter component in EAGLE. We then study the spatio-thermal distribution of the hot baryons within the extracted filaments. Results. While the filaments occupy only 5% of the full simulation volume, the diffuse hot intergalactic medium in filaments amounts to 23% $-$ 25% of the total baryon budget, or 79% $-$ 87% of all the hot WHIM. The most optimal filament sample, with a missing baryon mass fraction of 82%, is obtained by selecting Bisous filaments with a high galaxy luminosity density. For these filaments we derived analytic formulae for the radial gas density and temperature profiles, consistent with recent Planck SZ and CMB lensing observations within the central $r$~ 1 Mpc. Conclusions. Results from EAGLE suggest that the missing baryons are strongly concentrated towards the filament axes. Since the filament finding methods used here are applicable to galaxy surveys, a large fraction of the missing baryons can be localised by focusing the observational efforts on the central 1 Mpc regions of the filaments. Moreover, focusing on high galaxy luminosity density regions will optimise the observational signal.

preprint2020arXiv

EAGLE and Illustris-TNG predictions for resolved eROSITA X-ray observations of the circumgalactic medium around normal galaxies

We simulate stacked observations of nearby hot X-ray coronae associated with galaxies in the EAGLE and Illustris-TNG hydrodynamic simulations. A forward modeling pipeline is developed to predict 4-year eROSITA observations and stacked image analysis, including the effects of instrumental and astrophysical backgrounds. We propose an experiment to stack z~0.01 galaxies separated by specific star-formation rate (sSFR) to examine how the hot (T>=10^6 K) circumgalactic medium (CGM) differs for high- and low-sSFR galaxies. The simulations indicate that the hot CGM of low-mass (M_*~10^{10.5} Msol), high-sSFR (defined as the top one-third ranked by sSFR) central galaxies will be detectable to a galactocentric radius r~30-50 kpc. Both simulations predict lower luminosities at fixed stellar mass for the low-sSFR galaxies (the lower third of sSFR) with Illustris-TNG predicting 3x brighter coronae around high-sSFR galaxies than EAGLE. Both simulations predict detectable emission out to r~150-200 kpc for stacks centered on high-mass (M_*~10^{11.0} Msol) galaxies, with EAGLE predicting brighter X-ray halos. The extended soft X-ray luminosity correlates strongly and positively with the mass of circumgalactic gas within the virial radius (f_{CGM}). Prior analyses of both simulations have established that f_{CGM} is reduced by expulsive feedback driven mainly by black hole growth, which quenches galaxy growth by inhibiting replenishment of the ISM. Both simulations predict that eROSITA stacks should not only conclusively detect and resolve the hot CGM around L^* galaxies for the first time, but provide a powerful probe of how the baryon cycle operates, for which there remains an absence of consensus between state-of-the-art simulations.

preprint2020arXiv

Elevated ionizing photon production efficiency in faint high-equivalent-width Lyman-alpha emitters

While low-luminosity galaxies dominate number counts at all redshifts, their contribution to cosmic Reionization is poorly understood due to a lack of knowledge of their physical properties. We isolate a sample of 35 z~4-5 continuum-faint Lyman-alpha emitters from deep VLT/MUSE spectroscopy and directly measure their Halpha emission using stacked Spitzer/IRAC Ch. 1 photometry. Based on Hubble Space Telescope imaging, we determine that the average UV continuum magnitude is fainter than -16 (~0.01 L_star), implying a median Lyman-alpha equivalent width of 249 Angstroms. By combining the Halpha measurement with the UV magnitude we determine the ionizing photon production efficiency, xi_ion, a first for such faint galaxies. The measurement of log (xi_ion [Hz/erg]) = 26.28 (+0.28; -0.40) is in excess of literature measurements of both continuum- and emission line-selected samples, implying a more efficient production of ionizing photons in these lower-luminosity, Lyman-alpha-selected systems. We conclude that this elevated efficiency can be explained by stellar populations with metallicities between 4e-4 and 0.008, with light-weighted ages less than 3 Myr.

preprint2020arXiv

Galactic inflow and wind recycling rates in the EAGLE simulations

The role of galactic wind recycling represents one of the largest unknowns in galaxy evolution, as any contribution of recycling to galaxy growth is largely degenerate with the inflow rates of first-time infalling material, and the rates with which outflowing gas and metals are driven from galaxies. We present measurements of the efficiency of wind recycling from the EAGLE cosmological simulation project, leveraging the statistical power of large-volume simulations that reproduce a realistic galaxy population. We study wind recycling at the halo scale, i.e. gas that has been ejected beyond the halo virial radius, and at the galaxy scale, i.e. gas that has been ejected from the ISM to at least $\approx 10 \, \%$ of the virial radius (thus excluding smaller-scale galactic fountains). Galaxy-scale wind recycling is generally inefficient, with a characteristic return timescale that is comparable or longer than a Hubble time, and with an efficiency that clearly peaks at the characteristic halo mass of $M_{200} = 10^{12} \, \mathrm{M_\odot}$. Correspondingly, the majority of gas being accreted onto galaxies in EAGLE is infalling for the first time. At the halo scale, the efficiency of recycling onto haloes differs by orders of magnitude from values assumed by semi-analytic galaxy formation models. Differences in the efficiency of wind recycling with other hydrodynamical simulations are currently difficult to assess, but are likely smaller. We are able to show that the fractional contribution of wind recycling to galaxy growth is smaller in EAGLE than in some other simulations. We find that cumulative first-time gas accretion rates at the virial radius are reduced relative to the expectation from dark matter accretion for haloes with mass, $M_{200} < 10^{12} \, \mathrm{M_\odot}$, indicating efficient preventative feedback on halo scales.

preprint2020arXiv

Galactic outflow rates in the EAGLE simulations

We present measurements of galactic outflow rates from the EAGLE suite of cosmological simulations. We find that gas is removed from the interstellar medium (ISM) of central galaxies with a dimensionless mass loading factor that scales approximately with circular velocity as $V_{\mathrm{c}}^{-3/2}$ in the low-mass regime where stellar feedback dominates. Feedback from active galactic nuclei (AGN) causes an upturn in the mass loading for halo masses $> 10^{12} \, \mathrm{M_\odot}$. We find that more gas outflows through the halo virial radius than is removed from the ISM of galaxies, particularly at low redshift, implying substantial mass loading within the circum-galactic medium (CGM). Outflow velocities span a wide range at a given halo mass/redshift, and on average increase positively with redshift and halo mass up to $M_{200} \sim 10^{12} \, \mathrm{M_\odot}$. Outflows exhibit a bimodal flow pattern on circum-galactic scales, aligned with the galactic minor axis. We present a number of like-for-like comparisons to outflow rates from other recent cosmological hydrodynamical simulations, and show that comparing the propagation of galactic winds as a function of radius reveals substantial discrepancies between different models. Relative to some other simulations, EAGLE favours a scenario for stellar feedback where agreement with the galaxy stellar mass function is achieved by removing smaller amounts of gas from the ISM, but with galactic winds that then propagate and entrain ambient gas out to larger radii.

preprint2020arXiv

Numerical convergence of hydrodynamical simulations of galaxy formation: the abundance and internal structure of galaxies and their cold dark matter haloes

We address the issue of numerical convergence in cosmological smoothed particle hydrodynamics simulations using a suite of runs drawn from the EAGLE project. Our simulations adopt subgrid models that produce realistic galaxy populations at a fiducial mass and force resolution, but systematically vary the latter in order to study their impact on galaxy properties. We provide several analytic criteria that help guide the selection of gravitational softening for hydrodynamical simulations, and present results from runs that both adhere to and deviate from them. Unlike dark matter-only simulations, hydrodynamical simulations exhibit a strong sensitivity to gravitational softening, and care must be taken when selecting numerical parameters. Our results--which focus mainly on star formation histories, galaxy stellar mass functions and sizes--illuminate three main considerations. First, softening imposes a minimum resolved escape speed, $v_ε$, due to the binding energy between gas particles. Runs that adopt such small softening lengths that $v_ε\gt 10\,{\rm km s^{-1}}$ (the sound speed in ionised $\sim 10^4\,{\rm K}$ gas) suffer from reduced effects of photo-heating. Second, feedback from stars or active galactic nuclei may suffer from numerical over-cooling if the gravitational softening length is chosen below a critical value, $ε_{\rm eFB}$. Third, we note that small softening lengths exacerbate the segregation of stars and dark matter particles in halo centres, often leading to the counter-intuitive result that galaxy sizes {\em increase} as softening is reduced. The structure of dark matter haloes in hydrodynamical runs respond to softening in a way that reflects the sensitivity of their galaxy populations to numerical parameters.

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

Radiative cooling rates, ion fractions, molecule abundances and line emissivities including self-shielding and both local and metagalactic radiation fields

We use the spectral synthesis code Cloudy to tabulate the properties of gas for an extensive range in redshift (z=0 to 9), temperature (log T [K] = 1 to 9.5), metallicity (log Z/$\mathrm{Z}_{\odot}$ = -4 to +0.5, Z=0), and density (log $n_{\mathrm{H}}$ [$\mathrm{cm}^{-3}$] = -8 to +6). This therefore includes gas with properties characteristic of the interstellar, circumgalactic and intergalactic media. The gas is exposed to a redshift-dependent UV/X-ray background, while for the self-shielded lower-temperature gas (i.e. ISM gas) an interstellar radiation field and cosmic rays are added. The radiation field is attenuated by a density- and temperature-dependent column of gas and dust. Motivated by the observed star formation law, this gas column density also determines the intensity of the interstellar radiation field and the cosmic ray density. The ionization balance, molecule fractions, cooling rates, line emissivities, and equilibrium temperatures are calculated self-consistently. We include dust, cosmic rays, and the interstellar radiation field step-by-step to study their relative impact. These publicly available tables are ideal for hydrodynamical simulations. They can be used stand alone or coupled to a non-equilibrium network for a subset of elements. The release includes a C routine to read in and interpolate the tables, as well as an easy to use python graphical user interface to explore the tables.

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 effect of gas accretion on the radial gas metallicity profile of simulated galaxies

We study the effect of the gas accretion rate ($\dot M_{\rm accr}$) on the radial gas metallicity profile (RMP) of galaxies using the EAGLE cosmological hydrodynamic simulations, focusing on central galaxies of stellar mass $M_\star \gtrsim 10^9 \, {\rm M_\odot}$ at $z \le 1$. We find clear relations between $\dot M_{\rm accr}$ and the slope of the RMP (measured within an effective radius), where higher $\dot M_{\rm accr}$ are associated with more negative slopes. The slope of the RMPs depends more strongly on $\dot M_{\rm accr}$ than on stellar mass, star formation rate or gas fraction, suggesting $\dot M_{\rm accr}$ to be a more fundamental driver of the RMP slope of galaxies. We find that eliminating the dependence on stellar mass is essential for pinning down the properties that shape the slope of the RMP. Although $\dot M_{\rm accr}$ is the main property modulating the slope of the RMP, we find that it causes other correlations that are more easily testable observationally: at fixed stellar mass, galaxies with more negative RMP slopes tend to have higher gas fractions and SFRs, while galaxies with lower gas fractions and SFRs tend to have flatter metallicity profiles within an effective radius.

preprint2020arXiv

The MUSE Hubble Ultra Deep Field Survey XIV. The evolution of the Lya emitter fraction from z=3 to z=6

The Lya emitter (LAE) fraction, X_LAE, is a potentially powerful probe of the evolution of the intergalactic neutral hydrogen gas fraction. However, uncertainties in the measurement of X_LAE are still debated. Thanks to deep data obtained with MUSE, we can measure the evolution of X_LAE homogeneously over a wide redshift range of z~3-6 for UV-faint galaxies (down to M_1500~-17.75). This is significantly fainter than in former studies, and allows us to probe the bulk of the population of high-z star-forming galaxies. We construct a UV-complete photo-redshift sample following UV luminosity functions and measure the Lya emission with MUSE using the second data release from the MUSE HUDF Survey. We derive the redshift evolution of X_LAE for M_1500 in [-21.75;-17.75] for the first time with a equivalent width range EW(Lya)>=65 A and find low values of X_ LAE<~30% at z<~6. For M_1500 in [-20.25;-18.75] and EW(Lya)<~25 A, our X_LAE values are consistent with those in the literature within 1sigma at z<~5, but our median values are systematically lower than reported values over the whole redshift range. In addition, we do not find a significant dependence of X_LAE on M_1500 for EW(Lya)>~50 A at z~3-4, in contrast with previous work. The differences in X_LAE mainly arise from selection biases for Lyman Break Galaxies (LBGs) in the literature: UV-faint LBGs are more easily selected if they have strong Lya emission, hence X_LAE is biased towards higher values. Our results suggest either a lower increase of X_LAE towards z~6 than previously suggested, or even a turnover of X_LAE at z~5.5, which may be the signature of a late or patchy reionization process. We compared our results with predictions from a cosmological galaxy evolution model. We find that a model with a bursty star formation (SF) can reproduce our observed X_LAE much better than models where SF is a smooth function of time.

preprint2020arXiv

The MUSE Hubble Ultra Deep Field Survey. XV. The mean rest-UV spectra of Ly-alpha emitters at z>3

We investigate the ultraviolet (UV) spectral properties of faint Lyman-$α$ emitters (LAEs) in the redshift range 2.9<z<4.6 and provide material to prepare future observations of the faint Universe. We use data from the MUSE Hubble Ultra Deep Survey to construct mean rest-frame spectra of continuum-faint (median M$_{UV}$ of -18 and down to M$_{UV}$ of -16), low stellar mass (median value of $10^{8.4}$ and down to $10^{7}M_{\odot}$) LAEs at redshift z>3. We compute various averaged spectra of LAEs sub-sampled on the basis of their observational (e.g., Ly$α$ strength, UV magnitude and spectral slope) and physical (e.g., stellar mass and star-formation rate) properties. We search for UV spectral features other than Ly$α$, such as higher-ionization nebular emission lines and absorption features. We successfully observe the OIII]1666 and CIII]909 collisionally excited emission lines and the HeII1640 recombination feature, as well as the resonant CIV1550 doublet either in emission or P-Cygni. We compare the observed spectral properties of the different mean spectra and find the emission lines to vary with the observational and physical properties of the LAEs. In particular, the mean spectra of LAEs with larger Ly$α$ equivalent widths, fainter UV magnitudes, bluer UV spectral slopes and lower stellar masses show the strongest nebular emission. The line ratios of these lines are similar to those measured in the spectra of local metal-poor galaxies, while their equivalent widths are weaker compared to the handful of extreme values detected in individual spectra of z>2 galaxies. This suggests that weak UV features are likely ubiquitous in high z, low-mass and faint LAEs. We publicly release the stacked spectra as they can serve as empirical templates for the design of future observations, such as those with the James Webb Space Telescope and the Extremely Large Telescope.

preprint2020arXiv

The warm-hot circumgalactic medium around EAGLE-simulation galaxies and its detection prospects with X-ray and UV line absorption

We use the EAGLE (Evolution and Assembly of GaLaxies and their Environments) cosmological simulation to study the distribution of baryons, and far-ultraviolet (O VI), extreme-ultraviolet (Ne VIII) and X-ray (O VII, O VIII, Ne IX, and Fe XVII) line absorbers, around galaxies and haloes of mass $\mathrm{M}_{200c}=10^{11}$-$10^{14.5}\,\mathrm{M}_{\odot}$ at redshift 0.1. EAGLE predicts that the circumgalactic medium (CGM) contains more metals than the interstellar medium across halo masses. The ions we study here trace the warm-hot, volume-filling phase of the CGM, but are biased towards temperatures corresponding to the collisional ionization peak for each ion, and towards high metallicities. Gas well within the virial radius is mostly collisionally ionized, but around and beyond this radius, and for O VI, photoionization becomes significant. When presenting observables we work with column densities, but quantify their relation with equivalent widths by analysing virtual spectra. Virial-temperature collisional ionization equilibrium ion fractions are good predictors of column density trends with halo mass, but underestimate the diversity of ions in haloes. Halo gas dominates the highest column density absorption for X-ray lines, but lower density gas contributes to strong UV absorption lines from O VI and Ne VIII. Of the O VII (O VIII) absorbers detectable in an Athena X-IFU blind survey, we find that 41 (56) per cent arise from haloes with $\mathrm{M}_{200c}=10^{12.0}$-$10^{13.5}\,\mathrm{M}_{\odot}$. We predict that the X-IFU will detect O VII (O VIII) in 77 (46) per cent of the sightlines passing $\mathrm{M}_{\star}=10^{10.5}$-$10^{11.0}\,\mathrm{M}_{\odot}$ galaxies within 100 pkpc (59 (82) per cent for $\mathrm{M}_{\star}>10^{11.0}\,\mathrm{M}_{\odot}$). Hence, the X-IFU will probe covering fractions comparable to those detected with the Cosmic Origins Spectrograph for O VI.

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

Feedback from supermassive black holes transforms centrals into passive galaxies by ejecting circumgalactic gas

Davies et al. (2019) established that for L^* galaxies the fraction of baryons in the circumgalactic medium (CGM) is inversely correlated with the mass of their central supermassive black holes (BHs) in the EAGLE hydrodynamic simulation. The interpretation is that, over time, a more massive BH has provided more energy to transport baryons beyond the virial radius, which additionally reduces gas accretion and star formation. We continue this research by focusing on the relationship between the 1) BH masses, 2) physical and observational properties of the CGM, and 3) galaxy colours for Milky Way-mass systems. The ratio of the cumulative BH feedback energy over the gaseous halo binding energy is a strong predictor of the CGM gas content, with BHs injecting >~10x the binding energy resulting in gas-poor haloes. Observable tracers of the CGM, including CIV, OVI, and HI absorption line measurements, are found to be effective tracers of the total z~0 CGM halo mass. We use high-cadence simulation outputs to demonstrate that BH feedback pushes baryons beyond the virial radius within 100 Myr timescales, but that CGM metal tracers take longer (0.5-2.5 Gyr) to respond. Secular evolution of galaxies results in blue, star-forming or red, passive populations depending on the cumulative feedback from BHs. The reddest quartile of galaxies with M_*=10^{10.2-10.7} M_solar (median u-r = 2.28) has a CGM mass that is 2.5x lower than the bluest quartile (u-r=1.59). We propose strategies for observing the predicted lower CGM column densities and covering fractions around galaxies hosting more massive BHs using the Cosmic Origins Spectrograph on Hubble.

preprint2019arXiv

Hot WHIM counterparts of FUV OVI absorbers: Evidence in the line-of-sight towards quasar 3C 273

We explore the high spectral resolution X-ray data towards the quasar 3C273 to search for signals of hot ($\sim10^{6-7}$ K) X-ray-absorbing gas co-located with two established intergalactic FUV OVI absorbers. We analyze the soft X-ray band grating data of all XMM-Newton and Chandra instruments to search for the hot phase absorption lines at the FUV predicted redshifts. The viability of potential line detections is examined by adopting the constraints of a physically justified absorption model. The WHIM hypothesis is investigated with a complementary 3D galaxy distribution analysis, and by comparison of the measurement results to the WHIM properties in the EAGLE cosmological, hydrodynamical simulation. At FUV redshift z=0.09017, we measured signals of two hot ion species, OVIII and NeIX, with a $3.9σ$ combined significance level. Considering the line features in all instruments collectively and assuming collisional equilibrium for absorbing gas, we were able to constrain the temperature ($kT=0.26\pm0.03$ keV) and the column density ($N_H\times{Z_\odot/Z}=1.3_{-0.5}^{+0.6}\times10^{19}$ cm$^{-2}$) of the absorber. Thermal analysis indicates that FUV and X-ray absorption relate to different phases, with estimated temperatures $T_{FUV}\approx3\times10^5$ and $T_{X-ray}\approx3\times10^6$ K, which match the EAGLE predictions for WHIM at the FUV/X-ray measured $N_{ion}$-ranges. We detected a large scale galactic filament crossing the sightline at the redshift of the absorption, linking the absorption to this structure. This study provides insights into co-existing warm and hot gas within a WHIM filament and estimates the ratio of the hot and warm phases. Because the hot phase is thermally distinct from the OVI gas, the estimated baryon content of the absorber is increased, conveying the promise of X-ray follow-up studies of FUV detected WHIM in refining the picture of the missing baryons.

preprint2019arXiv

Measuring the temperature and profiles of Lyman-$α$ absorbers

The distribution of the absorption line broadening observed in the Ly$α$ forest carries information about the temperature and widths of the filaments in the intergalactic medium (IGM). In this work, we present and test a new method for inferring the temperature of the IGM, the size of the absorbing filaments and the photo-ionization rate of hydrogen from the properties of absorption line broadening in the Ly$α$ forest. We apply our method to mock spectra from the reference model of the EAGLE cosmological simulation, and we demonstrate that we are able to reconstruct the IGM properties. Our method explicitly takes into account the degeneracy between IGM temperature, the photo-ionization rate of hydrogen and the widths of the intergalactic filaments.

preprint2019arXiv

MusE GAs FLOw and Wind (MEGAFLOW) IV: A two sightline tomography of a galactic wind

Galactic outflows are thought to eject baryons back out to the circum-galactic medium (CGM). Studies based on metal absorption lines (MgII in particular) in the spectra of background quasars indicate that the gas is ejected anisotropically, with galactic winds likely leaving the host in a bi-conical flow perpendicular to the galaxy disk. In this paper, we present a detailed analysis of an outflow from a z = 0.7 "green-valley" galaxy (log($M_*$/$\mathrm{M}_\odot$) = 9.9; SFR = 0.5 $\mathrm{M}_\odot\,\mathrm{yr}^{-1}$) probed by two background sources part of the MUSE Gas Flow and Wind (MEGAFLOW) survey. Thanks to a fortuitous configuration with a background quasar (SDSSJ1358+1145) and a bright background galaxy at $z = 1.4$, both at impact parameters of $\approx 15\,\mathrm{kpc}$, we can - for the first time - probe both the receding and approaching components of a putative galactic outflow around a distant galaxy. We measure a significant velocity shift between the MgII absorption from the two sightlines ($84\pm17\,\mathrm{km}\,\mathrm{s}^{-1}$), which is consistent with the expectation from our simple fiducial wind model, possibly combined with an extended disk contribution.

preprint2019arXiv

The impact of the observed baryon distribution in haloes on the total matter power spectrum

The interpretation of upcoming weak gravitational lensing surveys depends critically on our understanding of the matter power spectrum on scales $k < 10 h/\mathrm{Mpc}$, where baryonic processes are important. We study the impact of galaxy formation processes on the matter power spectrum using a halo model that treats the stars and gas separately from the dark matter distribution. We use empirical constraints from X-ray observations (hot gas) and halo occupation distribution modelling (stars) for the baryons. Since X-ray observations cannot generally measure the hot gas content outside $r_\mathrm{500c}$, we vary the gas density profiles beyond this radius. Compared with dark matter only models, we find a total power suppression of $1 \%$ ($5 \%$) on scales $0.2-1 h/\mathrm{Mpc}$ ($0.5-2h/\mathrm{Mpc}$), where lower baryon fractions result in stronger suppression. We show that groups of galaxies ($10^{13} < m_\mathrm{500c} / (\mathrm{M_\odot}/h) < 10^{14}$) dominate the total power at all scales $k \lesssim 10 h/\mathrm{Mpc}$. We find that a halo mass bias of $30 \%$ (similar to what is expected from the hydrostatic equilibrium assumption) results in an underestimation of the power suppression of up to $4 \%$ at $k = 1 h/\mathrm{Mpc}$, illustrating the importance of measuring accurate halo masses. Contrary to work based on hydrodynamical simulations, our conclusion that baryonic effects can no longer be neglected is not subject to uncertainties associated with our poor understanding of feedback processes. Observationally, probing the outskirts of groups and clusters will provide the tightest constraints on the power suppression for $k \lesssim 1 h/\mathrm{Mpc}$.

preprint2016arXiv

A chronicle of galaxy mass assembly in the EAGLE simulation

We analyse the mass assembly of central galaxies in the EAGLE hydrodynamical simulations. We build merger trees to connect galaxies to their progenitors at different redshifts and characterize their assembly histories by focusing on the time when half of the galaxy stellar mass was assembled into the main progenitor. We show that galaxies with stellar mass $M_*<10^{10.5}M_{\odot}$ assemble most of their stellar mass through star formation in the main progenitor (`in-situ' star formation). This can be understood as a consequence of the steep rise in star formation efficiency with halo mass for these galaxies. For more massive galaxies, however, an increasing fraction of their stellar mass is formed outside the main progenitor and subsequently accreted. Consequently, while for low-mass galaxies the assembly time is close to the stellar formation time, the stars in high-mass galaxies typically formed long before half of the present-day stellar mass was assembled into a single object, giving rise to the observed anti-hierarchical downsizing trend. In a typical present-day $M_*\geq10^{11}M_{\odot}$ galaxy, around $20\%$ of the stellar mass has an external origin. This fraction decreases with increasing redshift. Bearing in mind that mergers only make an important contribution to the stellar mass growth of massive galaxies, we find that the dominant contribution comes from mergers with galaxies of mass greater than one tenth of the main progenitor's mass. The galaxy merger fraction derived from our simulations agrees with recent observational estimates.

preprint2016arXiv

A large difference in the progenitor masses of active and passive galaxies in the EAGLE simulation

Cumulative number density matching of galaxies is a method to observationally connect descendent galaxies to their typical main progenitors at higher redshifts and thereby to assess the evolution of galaxy properties. The accuracy of this method is limited due to galaxy merging and scatter in the stellar mass growth history of individual galaxies. Behroozi et al. (2013) have introduced a refinement of the method, based on abundance matching of observed galaxies to the Bolshoi dark-matter-only simulation. The EAGLE cosmological hydro-simulation is well suited to test this method, because it reproduces the observed evolution of the galaxy stellar mass function and the passive fraction. We find agreement with the Behroozi et al. (2013) method for the complete sample of main progenitors of z = 0 galaxies, but we also find a strong dependence on the current star formation rate. Passive galaxies with a stellar mass up to 10^10.75 Msun have a completely different median mass history than active galaxies of the same mass. This difference persists if we only select central galaxies. This means that the cumulative number density method should be applied separately to active and passive galaxies. Even then, the typical main progenitor of a z = 0 galaxy already spans two orders of magnitude in stellar mass at z = 2.

preprint2016arXiv

Alignments between galaxies, satellite systems and haloes

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

preprint2016arXiv

Being WISE II: Reducing the Influence of Star Formation History on the Mass-to-Light Ratio of Quiescent Galaxies

Stellar population synthesis models can now reproduce the photometry of old stellar systems (age $>$ 2 Gyr) in the near-infrared (NIR) bands at 3.4 and 4.6$μ$m (WISE W1 $\&$ W2 or IRAC 1 $\&$ 2). In this paper we derive stellar mass-to-light ratios for these and optical bands, and confirm that the NIR M/L shows dramatically reduced sensitivity to both age and metallicity compared to optical bands, and further, that this behavior leads to significantly more robust stellar masses for quiescent galaxies with [Fe/H] > -0.5 regardless of star formation history (SFH). We then use realistic early-type galaxy SFHs and metallicity distributions from the EAGLE simulations of galaxy formation to investigate two methods to determine the appropriate M/L for a galaxy: 1) We show that the uncertainties introduced by an unknown SFH can be largely removed using a spectroscopically inferred luminosity-weighted age and metallicity for the population to select the appropriate single stellar population (SSP) equivalent M/L. Using this method, the maximum systematic error due to SFH on the M/L of an early-type galaxy is $<$ 4$\%$ at 3.4 $μ$m and typical uncertainties due to errors in the age and metallicity create scatter of $\lesssim$13$\%$. The equivalent values for optical bands are more than 2-3 times greater, even before considering uncertainties associated with internal dust extinction. 2) We demonstrate that if the EAGLE SFHs and metallicities accurately reproduce the true properties of early-type galaxies, the use of an iterative approach to select a mass dependent M/L can provide even more accurate stellar masses for early-type galaxies, with typical uncertainties $<$ 9$\%$.

preprint2016arXiv

Bimodality of low-redshift circumgalactic O VI in non-equilibrium EAGLE zoom simulations

We introduce a series of 20 cosmological hydrodynamical simulations of Lstar (M_200 =10^11.7 - 10^12.3 Msol) and group-sized (M_200 = 10^12.7 - 10^13.3 Msol) haloes run with the model used for the EAGLE project, which additionally includes a non-equilibrium ionization and cooling module that follows 136 ions. The simulations reproduce the observed correlation, revealed by COS-Halos at z~0.2, between O VI column density at impact parameters b < 150 kpc and the specific star formation rate (sSFR=SFR/Mstar) of the central galaxy at z~0.2. We find that the column density of circumgalactic O VI is maximal in the haloes associated with Lstar galaxies, because their virial temperatures are close to the temperature at which the ionization fraction of O VI peaks (T~10^5.5 K). The higher virial temperature of group haloes (> 10^6 K) promotes oxygen to higher ionization states, suppressing the O VI column density. The observed NO VI-sSFR correlation therefore does not imply a causal link, but reflects the changing characteristic ionization state of oxygen as halo mass is increased. In spite of the mass-dependence of the oxygen ionization state, the most abundant circumgalactic oxygen ion in both Lstar and group haloes is O VII; O VI accounts for only 0.1% of the oxygen in group haloes and 0.9-1.3% with Lstar haloes. Nonetheless, the metals traced by O VI absorbers represent a fossil record of the feedback history of galaxies over a Hubble time; their characteristic epoch of ejection corresponds to z > 1 and much of the ejected metal mass resides beyond the virial radius of galaxies. For both Lstar and group galaxies, more of the oxygen produced and released by stars resides in the circumgalactic medium (within twice the virial radius) than in the stars and ISM of the galaxy.

preprint2016arXiv

Chemical evolution of giant molecular clouds in simulations of galaxies

We present an analysis of Giant Molecular Clouds (GMCs) within hydrodynamic simulations of isolated, low-mass (M* ~ 10^9 M_sol) disc galaxies. We study the evolution of molecular abundances and the implications for CO emission and the X_CO conversion factor in individual clouds. We define clouds either as regions above a density threshold n_H,min = 10 cm^-3, or using an observationally motivated CO intensity threshold of 0.25 K km s^-1. Our simulations include a non-equilibrium chemical model with 157 species, including 20 molecules. We also investigate the effects of resolution and pressure floors (i.e. Jeans limiters). We find cloud lifetimes up to ~40 Myr, with a median of 13 Myr, in agreement with observations. At one tenth solar metallicity, young clouds (<10-15 Myr) are underabundant in H2 and CO compared to chemical equilibrium, by factors of ~3 and 1-2 orders of magnitude, respectively. At solar metallicity, GMCs reach chemical equilibrium faster (within ~1 Myr). We also compute CO emission from individual clouds. The mean CO intensity, I_CO, is strongly suppressed at low dust extinction, A_v, and possibly saturates towards high A_v, in agreement with observations. The I_CO - A_v relation shifts towards higher A_v for higher metallicities and, to a lesser extent, for stronger UV radiation. At one tenth solar metallicity, CO emission is weaker in young clouds (<10-15 Myr), consistent with the underabundance of CO. Consequently, X_CO decreases by an order of magnitude from 0 to 15 Myr, albeit with a large scatter.

preprint2016arXiv

Cosmic distribution of highly ionized metals and their physical conditions in the EAGLE simulations

We study the distribution and evolution of highly ionised intergalactic metals in the Evolution and Assembly of Galaxies and their Environment (EAGLE) cosmological, hydrodynamical simulations. EAGLE has been shown to reproduce a wide range of galaxy properties while its subgrid feedback was calibrated without considering gas properties. We compare the predictions for the column density distribution functions (CDDFs) and cosmic densities of SiIV, CIV, NV, OVI and NeVIII absorbers with observations at redshift z = 0 to ~ 6 and find reasonable agreement, although there are some differences. We show that the typical physical densities of the absorbing gas increase with column density and redshift, but decrease with the ionization energy of the absorbing ion. The typical metallicity increases with both column density and time. The fraction of collisionally ionized metal absorbers increases with time and ionization energy. While our results show little sensitivity to the presence or absence of AGN feedback, increasing/decreasing the efficiency of stellar feedback by a factor of two substantially decreases/increases the CDDFs and the cosmic densities of the metal ions. We show that the impact of the efficiency of stellar feedback on the CDDFs and cosmic densities is largely due to its effect on the metal production rate. However, the temperatures of the metal absorbers, particularly those of strong OVI, are directly sensitive to the strength of the feedback.

preprint2016arXiv

Far-infrared and dust properties of present-day galaxies in the EAGLE simulations

The EAGLE cosmological simulations reproduce the observed galaxy stellar mass function and many galaxy properties. In this work, we study the dust-related properties of present-day EAGLE galaxies through mock observations in the far-infrared and submm wavelength ranges obtained with the 3D dust radiative transfer code SKIRT. To prepare an EAGLE galaxy for radiative transfer processing, we derive a diffuse dust distribution from the gas particles and we re-sample the star-forming gas particles and the youngest star particles into star-forming regions that are assigned dedicated emission templates. We select a set of redshift-zero EAGLE galaxies that matches the K-band luminosity distribution of the galaxies in the Herschel Reference Survey (HRS), a volume-limited sample of about 300 normal galaxies in the Local Universe. We find overall agreement of the EAGLE dust scaling relations with those observed in the HRS, such as the dust-to-stellar mass ratio versus stellar mass and versus NUV-r colour relations. A discrepancy in the f_250/f_350 versus f_350/f_500 submm colour-colour relation implies that part of the simulated dust is insufficiently heated, likely because of limitations in our sub-grid model for star-forming regions. We also investigate the effect of adjusting the metal-to-dust ratio and the covering factor of the photodissociation regions surrounding the star-forming cores. We are able to constrain the important dust-related parameters in our method, informing the calculation of dust attenuation for EAGLE galaxies in the UV and optical domain.

preprint2016arXiv

Galaxies in the EAGLE hydrodynamical simulation and in the Durham and Munich semi-analytical models

We compare global predictions from the EAGLE hydrodynamical simulation, and two semi-analytic (SA) models of galaxy formation, L-GALAXIES and GALFORM. All three models include the key physical processes for the formation and evolution of galaxies and their parameters are calibrated against a small number of observables at $z\approx 0$. The two SA models have been applied to merger trees constructed from the EAGLE dark matter only simulation. We find that at $z\leq 2$, both the galaxy stellar mass functions for stellar masses $M_*<10^{10.5} M_{\odot}$ and the median specific star formation rates (sSFRs) in the three models agree to better than $0.4$~dex. The evolution of the sSFR predicted by the three models closely follows the mass assembly history of dark matter haloes. In both EAGLE and L-GALAXIES there are more central passive galaxies with $M_*<10^{9.5} M_{\odot}$ than in L-GALAXIES. This difference is related to galaxies that have entered and then left a larger halo and which are treated as satellites in GALFORM. In the range $0<z<1$, the slope of the evolution of the star formation rate density in EAGLE is a factor of $\approx 1.5$ steeper than for the two SA models. The median sizes for galaxies with $M_*>10^{9.5} M_{\odot}$ differ in some instances by an order of magnitude, while the stellar mass-size relation in EAGLE is a factor of $\approx 2$ tighter than for the two SA models. Our results suggest the need for a revision of how SA models treat the effect of baryonic self-gravity on the underlying dark matter. The treatment of gas flows in the models needs to be revised based on detailed comparison with observations to understand in particular the evolution of the stellar mass-metallicity relation.

preprint2016arXiv

Implications of a variable IMF for the interpretation of observations of galaxy populations

We investigate the effect of a metallicity-dependent stellar initial mass function (IMF), as deduced observationally by Martin-Navarro et al. (2015c), on the inferred stellar masses and star formation rates (SFRs) of a representative sample of 186,886 SDSS galaxies. Relative to a Chabrier IMF, for which we show the implied masses to be close to minimal, the inferred masses increase in both the low- and high-metallicity regimes due to the addition of stellar remnants and dwarf stars, respectively. The resulting galaxy stellar mass function (GSMF) shifts toward higher masses by 0.5 dex, without affecting the high-mass slope (and thus the need for effective quenching). The implied low-redshift SFR density increases by an order of magnitude. However, these results depend strongly on the assumed IMF parametrisation, which is not directly constrained by the observations. Varying the low-end IMF slope instead of the high-end IMF slope, while maintaining the same dwarf-to-giant ratio, results in a much more modest GSMF shift of 0.2 dex and a 10 per cent increase in the SFR density relative to the Chabrier IMF. A bottom-heavy IMF during the late, metal-rich evolutionary stage of a galaxy would help explain the rapid quenching and the bimodality in the galaxy population by on the one hand making galaxies less quenched (due to the continued formation of dwarf stars) and on the other hand reducing the gas consumption timescale. We conclude that the implications of the observational evidence for a variable IMF could vary from absolutely dramatic to mild but significant.

preprint2016arXiv

It's not easy being green: The evolution of galaxy colour in the EAGLE simulation

We examine the evolution of intrinsic u-r colours of galaxies in the EAGLE cosmological hydrodynamical simulations, which has been shown to reproduce the observed redshift z=0.1 colour-magnitude distribution well. The median u-r of star-forming ('blue cloud') galaxies reddens by 1 mag from z=2 to 0 at fixed stellar mass, as their specific star formation rates decrease with time. A red sequence starts to build-up around z=1, due to the quenching of low-mass satellite galaxies at the faint end, and due to the quenching of more massive central galaxies by their active galactic nuclei (AGN) at the bright end. This leaves a dearth of intermediate-mass red sequence galaxies at z=1, which is mostly filled in by z=0. We quantify the time-scales of colour transition due to satellite and AGN quenching, finding that most galaxies spend less than 2 Gyr in the 'green valley'. On examining the trajectories of galaxies in a colour-stellar mass diagram, we identify three characteristic tracks that galaxies follow (quiescently star-forming, quenching and rejuvenating galaxies) and quantify the fraction of galaxies that follow each track.

preprint2016arXiv

MusE GAs FLOw and Wind (MEGAFLOW) I: First MUSE results on background quasars

The physical properties of galactic winds are one of the keys to understand galaxy formation and evolution. These properties can be constrained thanks to background quasar lines of sight (LOS) passing near star-forming galaxies (SFGs). We present the first results of the MusE GAs FLOw and Wind (MEGAFLOW) survey obtained of 2 quasar fields which have 8 MgII absorbers of which 3 have rest-equivalent width greater than 0.8 Å. With the new Multi Unit Spectroscopic Explorer (MUSE) spectrograph on the Very Large Telescope (VLT), we detect 6 (75$\%$) MgII host galaxy candidates withing a radius of 30 arcsec from the quasar LOS. Out of these 6 galaxy--quasar pairs, from geometrical arguments, one is likely probing galactic outflows, two are classified as "ambiguous", two are likely probing extended gaseous disks and one pair seems to be a merger. We focus on the wind$-$pair and constrain the outflow using a high resolution quasar spectra from Ultraviolet and Visual Echelle Spectrograph (UVES). Assuming the metal absorption to be due to gas flowing out of the detected galaxy through a cone along the minor axis, we find outflow velocities of the order of $\approx$ 150 km/s (i.e. smaller than the escape velocity) with a loading factor, $η=\dot M_{\rm out}/$SFR, of $\approx$ 0.7. We see evidence for an open conical flow, with a low-density inner core. In the future, MUSE will provide us with about 80 multiple galaxy$-$quasar pairs in two dozen fields.

preprint2016arXiv

Music from the heavens - Gravitational waves from supermassive black hole mergers in the EAGLE simulations

We estimate the expected event rate of gravitational wave signals from mergers of supermassive black holes that could be resolved by a space-based interferometer, such as the Evolved Laser Interferometer Space Antenna (eLISA), utilising the reference cosmological hydrodynamical simulation from the EAGLE suite. These simulations assume a $Λ$CDM cosmogony with state-of-the-art subgrid models for radiative cooling, star formation, stellar mass loss, and feedback from stars and accreting black holes. They have been shown to reproduce the observed galaxy population with unprecedented fidelity. We combine the merger rates of supermassive black holes in EAGLE with the latest phenomenological waveform models to calculate the gravitational waves signals from the intrinsic parameters of the merging black holes. The EAGLE models predict $\sim2$ detections per year by a gravitational wave detector such as eLISA. We find that these signals are largely dominated by mergers between seed mass black holes merging at redshifts between $z\sim2$ and $z\sim1$. In order to investigate the dependence on the assumed black hole seed mass, we introduce an additional model with a black hole seed mass an order of magnitude smaller than in our reference model. We also consider a variation of the reference model where a prescription for the expected delays in the black hole merger timescale has been included after their host galaxies merge. We find that the merger rate is similar in all models, but that the initial black hole seed mass could be distinguished through their detected gravitational waveforms. Hence, the characteristic gravitational wave signals detected by eLISA will provide profound insight into the origin of supermassive black holes and the initial mass distribution of black hole seeds.

preprint2016arXiv

Observations of metals in the $z\approx3.5$ intergalactic medium and comparison to the EAGLE simulations

We study the $z\approx3.5$ intergalactic medium (IGM) by comparing new, high-quality absorption spectra of eight QSOs with $\langle z_{\rm QSO} \rangle=3.75$, to virtual observations of the EAGLE cosmological hydrodynamical simulations. We employ the pixel optical depth method and uncover strong correlations between various combinations of HI, CIII, CIV, SiIII, SiIV, and OVI. We find good agreement between many of the simulated and observed correlations, including OVI(HI). However, the observed median optical depths for the CIV(HI) and SiIV(HI) relations are higher than those measured from the mock spectra. The discrepancy increases from up to $\approx0.1$ dex at $τ_{\rm HI}=1$ to $\approx1$ dex at $τ_{\rm HI}=10^2$, where we are likely probing dense regions at small galactocentric distances. As possible solutions, we invoke (a) models of ionizing radiation softened above 4 Ryd to account for delayed completion of HeII reionization; (b) simulations run at a higher resolution; (c) the inclusion of additional line broadening due to unresolved turbulence; and (d) increased elemental abundancess; however, none of these factors can fully explain the observed differences. Enhanced photoionization of HI by local sources, which was not modelled, could offer a solution. However, the much better agreement with the observed OVI(HI) relation, which we find probes a hot and likely collisionally-ionized gas phase, indicates that the simulations are not in tension with the hot phase of the IGM, and suggests that the simulated outflows may entrain insufficient cool gas.

preprint2016arXiv

On the connection between the metal-enriched intergalactic medium and galaxies: an OVI-galaxy cross-correlation study at $z < 1$

We present new results on the auto- and cross-correlation functions of galaxies and OVI absorbers in a $\sim 18~\textrm{Gpc}^3$ comoving volume at $z < 1$. We use a sample of 51,296 galaxies and 140 OVI absorbers in the column density range $13 \lesssim \log N \lesssim 15$ to measure two-point correlation functions in the two dimensions transverse and orthogonal to the line-of-sight $ξ(r_{\perp}, r_{\parallel})$. We furthermore infer the corresponding 'real-space' correlation functions, $ξ(r)$, by projecting $ξ(r_{\perp}, r_{\parallel})$ along $r_{\parallel}$, and assuming a power-law form, $ξ(r) = (r / r_0)^{-γ}$. Comparing the results from the absorber-galaxy cross-correlation function, $ξ_{\textrm{ag}}$, the galaxy auto-correlation function, $ξ_{\textrm{gg}}$, and the absorber auto-correlation function, $ξ_{\textrm{aa}}$, we constrain the statistical connection between galaxies and the metal-enriched intergalactic medium as a function of star-formation activity. We also compare these results to predictions from the EAGLE cosmological hydrodynamical simulation and find a reasonable agreement. We find that: (i) OVI absorbers show very little velocity dispersion with respect to galaxies on $\sim$ Mpc scales, likely $\lesssim$ 100 \kms; (ii) OVI absorbers and galaxies may not linearly trace the same underlying distribution of matter in general. In particular, our results demonstrate that OVI absorbers are less clustered, and potentially more extended around galaxies than galaxies are around themselves; (iii) On $\gtrsim 100$ kpc scales, the likelihood of finding OVI absorbers around star-forming galaxies is similar to the likelihood of finding OVI absorbers around non star-forming galaxies (abridged)

preprint2016arXiv

Simulated Milky Way analogues: implications for dark matter direct searches

We study the implications of galaxy formation on dark matter direct detection using high resolution hydrodynamic simulations of Milky Way-like galaxies simulated within the EAGLE and APOSTLE projects. We identify Milky Way analogues that satisfy observational constraints on the Milky Way rotation curve and total stellar mass. We then extract the dark matter density and velocity distribution in the Solar neighbourhood for this set of Milky Way analogues, and use them to analyse the results of current direct detection experiments. For most Milky Way analogues, the event rates in direct detection experiments obtained from the best fit Maxwellian distribution (with peak speed of 223 - 289 km/s) are similar to those obtained directly from the simulations. As a consequence, the allowed regions and exclusion limits set by direct detection experiments in the dark matter mass and spin-independent cross section plane shift by a few GeV compared to the Standard Halo Model, at low dark matter masses. For each dark matter mass, the halo-to-halo variation of the local dark matter density results in an overall shift of the allowed regions and exclusion limits for the cross section. However, the compatibility of the possible hints for a dark matter signal from DAMA and CDMS-Si and null results from LUX and SuperCDMS is not improved.

preprint2016arXiv

Size matters: abundance matching, galaxy sizes, and the Tully-Fisher relation in EAGLE

The Tully-Fisher relation (TFR) links the stellar mass of a disk galaxy, $M_{\rm str}$, to its rotation speed: it is well approximated by a power law, shows little scatter, and evolves weakly with redshift. The relation has been interpreted as reflecting the mass-velocity scaling ($M\propto V^3$) of dark matter halos, but this interpretation has been called into question by abundance-matching (AM) models, which predict the galaxy-halo mass relation to be non-monotonic and rapidy evolving. We study the TFR of luminous spirals and its relation to AM using the EAGLE set of $Λ$CDM cosmological simulations. Matching both relations requires disk sizes to satisfy constraints given by the concentration of halos and their response to galaxy assembly. EAGLE galaxies approximately match these constraints and show a tight mass-velocity scaling that compares favourably with the observed TFR. The TFR is degenerate to changes in galaxy formation efficiency and the mass-size relation; simulations that fail to match the galaxy stellar mass function may fit the observed TFR if galaxies follow a different mass-size relation. The small scatter in the simulated TFR results because, at fixed halo mass, galaxy mass and rotation speed correlate strongly, scattering galaxies along the main relation. EAGLE galaxies evolve with lookback time following approximately the prescriptions of AM models and the observed mass-size relation of bright spirals, leading to a weak TFR evolution consistent with observation out to $z=1$. $Λ$CDM models that match both the abundance and size of galaxies as a function of stellar mass have no difficulty reproducing the observed TFR and its evolution.

preprint2016arXiv

Subhalo abundance matching and assembly bias in the EAGLE simulation

Subhalo abundance matching (SHAM) is a widely-used method to connect galaxies with dark matter structures in numerical simulations. SHAM predictions agree remarkably well with observations, yet they still lack strong theoretical support. We examine the performance, implementation, and assumptions of SHAM using the EAGLE project simulations. We find that $V_{\rm relax}$, the highest value of the circular velocity attained by a subhalo while it satisfies a relaxation criterion, is the subhalo property that correlates most strongly with galaxy stellar mass ($M_{\rm star}$). Using this parameter in SHAM, we retrieve the real-space clustering of EAGLE to within our statistical uncertainties on scales greater than $2$ Mpc for galaxies with $8.77<\log_{10}(M_{\rm star}[M_\odot])<10.77$. Conversely, clustering is overestimated by $30\%$ on scales below $2$ Mpc for galaxies with $8.77<\log_{10}(M_{\rm star}[M_\odot])<9.77$ because SHAM slightly overpredicts the fraction of satellites in massive haloes compared to EAGLE. The agreement is even better in redshift-space, where the clustering is recovered to within our statistical uncertainties for all masses and separations. Additionally, we analyse the dependence of galaxy clustering on properties other than halo mass, i.e. the assembly bias. We demonstrate assembly bias alters the clustering in EAGLE by $20\%$ and that $V_{\rm relax}$ captures its effect to within $15\%$. We trace small differences in the clustering to the failure of SHAM as typically implemented, i.e. the $M_{\rm star}$ assigned to a subhalo does not depend on i) its host halo mass, ii) whether it is a central or a satellite. In EAGLE we find that these assumptions are not completely satisfied.

preprint2016arXiv

Supermassive black holes in the EAGLE Universe. Revealing the observables of their growth

We investigate the evolution of supermassive black holes in the `Evolution and Assembly of GaLaxies and their Environments' (EAGLE) cosmological hydrodynamic simulations. The largest of the EAGLE volumes covers a $(100 \,\rm cMpc)^3$ and includes state-of-the-art physical models for star formation and black hole growth that depend only on local gas properties. We focus on the black hole mass function, Eddington ratio distribution and the implied duty cycle of nuclear activity. The simulation is broadly consistent with observational constraints on these quantities. In order to make a more direct comparison with observational data, we calculate the soft and hard X-ray luminosity functions of the active galactic nuclei (AGN). Between redshifts $0$ and $1$, the simulation is in agreement with data. At higher redshifts, the simulation tends to underpredict the luminosities of the brightest observed AGN. This may be due to the limited volume of the simulation, or a fundamental deficiency of the underlying model. It seems unlikely that additional unresolved variability can account for this difference. The simulation shows a similar `downsizing' of the AGN population as seen in observational surveys.

preprint2016arXiv

The Athena X-ray Integral Field Unit (X-IFU)

The X-ray Integral Field Unit (X-IFU) on board the Advanced Telescope for High-ENergy Astrophysics (Athena) will provide spatially resolved high-resolution X-ray spectroscopy from 0.2 to 12 keV, with 5 arc second pixels over a field of view of 5 arc minute equivalent diameter and a spectral resolution of 2.5 eV up to 7 keV. In this paper, we first review the core scientific objectives of Athena, driving the main performance parameters of the X-IFU, namely the spectral resolution, the field of view, the effective area, the count rate capabilities, the instrumental background. We also illustrate the breakthrough potential of the X-IFU for some observatory science goals. Then we briefly describe the X-IFU design as defined at the time of the mission consolidation review concluded in May 2016, and report on its predicted performance. Finally, we discuss some options to improve the instrument performance while not increasing its complexity and resource demands (e.g. count rate capability, spectral resolution). The X-IFU will be provided by an international consortium led by France, The Netherlands and Italy, with further ESA member state contributions from Belgium, Finland, Germany, Poland, Spain, Switzerland and two international partners from the United States and Japan.

preprint2016arXiv

The brighter galaxies reionised the Universe

Hydrogen in the Universe was (re)ionised between redshifts $z \approx 10$ and $z \approx 6$. The nature of the sources of the ionising radiation is hotly debated, with faint galaxies below current detection limits regarded as prime candidates. Here we consider a scenario in which ionising photons escape through channels punctured in the interstellar medium by outflows powered by starbursts. We take account of the observation that strong outflows occur only when the star formation density is sufficiently high, and estimate the galaxy-averaged escape fraction as a function of redshift and luminosity from the resolved star formation surface densities in the EAGLE cosmological hydrodynamical simulation. We find that the fraction of ionising photons that escape from galaxies increases rapidly with redshift, reaching values of 5-20 percent at $z > 6$, with the brighter galaxies having higher escape fractions. Combining the dependence of escape fraction on luminosity and redshift with the observed luminosity function, we demonstrate that galaxies emit enough ionising photons to match the existing constraints on reionisation while also matching the observed UV-background post-reionisation. Our findings suggest that galaxies above the current Hubble Space Telescope detection limit emit half of the ionising radiation required to reionise the Universe.

preprint2016arXiv

The dark nemesis of galaxy formation: why hot haloes trigger black hole growth and bring star formation to an end

Galaxies fall into two clearly distinct types: `blue-sequence' galaxies that are rapidly forming young stars, and `red-sequence' galaxies in which star formation has almost completely ceased. Most galaxies more massive than $3\times10^{10} M_\odot$ follow the red-sequence while less massive central galaxies lie on the blue sequence. We show that these sequences are created by a competition between star formation-driven outflows and gas accretion on to the supermassive black hole at the galaxy's center. We develop a simple analytic model for this interaction. In galaxies less massive than $3\times10^{10} M_\odot$, young stars and supernovae drive a high entropy outflow that is more buoyant that any diffuse corona. The outflow balances the rate of gas inflow, preventing high gas densities building up in the central regions. More massive galaxies, however, are surrounded by a hot corona. We argue that above a halo mass of $\sim 10^{12} M_\odot$, the supernova-driven outflow is no longer buoyant and star formation is unable to prevent the build up of gas in the central regions. This triggers a strongly non-linear response from the black hole. Its accretion rate rises rapidly, heating the galaxy's corona, disrupting the incoming supply of cool gas and starving the galaxy of the fuel for star formation. The host galaxy makes a transition to the red sequence, and further growth predominantly occurs through galaxy mergers. We show that the analytic model provides a good description of galaxy evolution in the EAGLE hydrodynamic simulations, and demonstrate that, so long as star formation-driven outflows are present, the transition mass scale is almost independent of subgrid parameter choice. The transition mass disappears entirely, however, if star formation driven outflows are absent.

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 effect of baryons on redshift space distortions and cosmic density and velocity fields in the EAGLE simulation

We use the EAGLE galaxy formation simulation to study the effects of baryons on the power spectrum of the total matter and dark matter distributions and on the velocity fields of dark matter and galaxies. On scales $k{\stackrel{>}{{}_\sim}} 4{h\,{\rm Mpc}^{-1}}$ the effect of baryons on the amplitude of the total matter power spectrum is greater than $1\%$. The back-reaction of baryons affects the density field of the dark matter at the level of $\sim3\%$ on scales of $1\leq k/({h\,{\rm Mpc}^{-1}})\leq 5$. The dark matter velocity divergence power spectrum at $k{\stackrel{<}{{}_\sim}}0.5{h\,{\rm Mpc}^{-1}}$ is changed by less than $1\%$. The 2D redshift-space power spectrum is affected at the level of $\sim6\%$ at $|\vec{k}|{\stackrel{>}{{}_\sim}} 1{h\,{\rm Mpc}^{-1}}$ (for $μ>0.5$), but for $|\vec{k}|\leq 0.4{h\,{\rm Mpc}^{-1}}$ it differs by less than $1\%$. We report vanishingly small baryonic velocity bias for haloes: the peculiar velocities of haloes with $M_{200}>3\times10^{11}{{\rm M}_{\odot}}$ (hosting galaxies with $M_{*}>10^9{{\rm M}_{\odot}}$) are affected at the level of at most $1~$km/s, which is negligible for $1\%$-precision cosmology. We caution that since EAGLE overestimates cluster gas fractions it may also underestimate the impact of baryons, particularly for the total matter power spectrum. Nevertheless, our findings suggest that for theoretical modelling of redshift space distortions and galaxy velocity-based statistics, baryons and their back-reaction can be safely ignored at the current level of observational accuracy. However, we confirm that the modelling of the total matter power spectrum in weak lensing studies needs to include realistic galaxy formation physics in order to achieve the accuracy required in the precision cosmology era.

preprint2016arXiv

The effects of metallicity, UV radiation and non-equilibrium chemistry in high-resolution simulations of galaxies

We present a series of hydrodynamic simulations of isolated galaxies with stellar mass of $10^{9} \, \rm{M}_{\odot}$. The models use a resolution of $750 \, \rm{M}_{\odot}$ per particle and include a treatment for the full non-equilibrium chemical evolution of ions and molecules (157 species in total), along with gas cooling rates computed self-consistently using the non-equilibrium abundances. We compare these to simulations evolved using cooling rates calculated assuming chemical (including ionisation) equilibrium, and we consider a wide range of metallicities and UV radiation fields, including a local prescription for self-shielding by gas and dust. We find higher star formation rates and stronger outflows at higher metallicity and for weaker radiation fields, as gas can more easily cool to a cold (few hundred Kelvin) star forming phase under such conditions. Contrary to variations in the metallicity and the radiation field, non-equilibrium chemistry generally has no strong effect on the total star formation rates or outflow properties. However, it is important for modelling molecular outflows. For example, the mass of H$_{2}$ outflowing with velocities $> 50 \, \rm{km} \, \rm{s}^{-1}$ is enhanced by a factor $\sim 20$ in non-equilibrium. We also compute the observable line emission from CII and CO. Both are stronger at higher metallicity, while CII and CO emission are higher for stronger and weaker radiation fields respectively. We find that CII is generally unaffected by non-equilibrium chemistry. However, emission from CO varies by a factor of $\sim 2 - 4$. This has implications for the mean $X_{\rm{CO}}$ conversion factor between CO emission and H$_{2}$ column density, which we find is lowered by up to a factor $\sim 2.3$ in non-equilibrium, and for the fraction of CO-dark molecular gas.

preprint2016arXiv

The environmental dependence of HI in galaxies in the EAGLE simulations

We use the EAGLE suite of cosmological hydrodynamical simulations to study how the HI content of present-day galaxies depends on their environment. We show that EAGLE reproduces observed HI mass-environment trends very well, while semi-analytic models typically overpredict the average HI masses in dense environments. The environmental processes act primarily as an on/off switch for the HI content of satellites with stellar mass Mstar>10^9 Msun. At a fixed Mstar, the fraction of HI-depleted satellites increases with increasing host halo mass M200 in response to stronger environmental effects, while at a fixed M200 it decreases with increasing satellite Mstar as the gas is confined by deeper gravitational potentials. HI-depleted satellites reside mostly, but not exclusively, within the virial radius r200 of their host halo. We investigate the origin of these trends by focussing on three environmental mechanisms: ram pressure stripping by the intra-group medium, tidal stripping by the host halo, and satellite-satellite encounters. By tracking back in time the evolution of the HI-depleted satellites, we find that the most common cause of HI removal is satellite encounters. The timescale for HI removal is typically less than 0.5 Gyr. Tidal stripping occurs in halos of M200<10^{14} Msun within half r200, while the other processes act also in more massive halos, generally within r200. Conversely, we find that ram pressure stripping is the most common mechanism that disturbs the HI morphology of galaxies at redshift z=0. This implies that HI removal due to satellite-satellite interactions occurs on shorter timescales than the other processes.

preprint2016arXiv

The fundamental plane of star formation in galaxies revealed by the EAGLE hydrodynamical simulations

We investigate correlations between different physical properties of star-forming galaxies in the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) cosmological hydrodynamical simulation suite over the redshift range $0\le z\le 4.5$. A principal component analysis reveals that neutral gas fraction ($f_{\rm gas, neutral}$), stellar mass ($M_{\rm stellar}$) and star formation rate (SFR) account for most of the variance seen in the population, with galaxies tracing a two-dimensional, nearly flat, surface in the three-dimensional space of $f_{\rm gas, neutral}-M_{\rm stellar}-\rm SFR$ with little scatter. The location of this plane varies little with redshift, whereas galaxies themselves move along the plane as their $f_{\rm gas, neutral}$ and SFR drop with redshift. The positions of galaxies along the plane are highly correlated with gas metallicity. The metallicity can therefore be robustly predicted from $f_{\rm gas, neutral}$, or from the $M_{\rm stellar}$ and SFR. We argue that the appearance of this "fundamental plane of star formation" is a consequence of self-regulation, with the plane's curvature set by the dependence of the SFR on gas density and metallicity. We analyse a large compilation of observations spanning the redshift range $0\lesssim \rm z\lesssim 2.5$, and find that such a plane is also present in the data. The properties of the observed fundamental plane of star formation are in good agreement with EAGLE's predictions.

preprint2016arXiv

The low-mass end of the baryonic Tully-Fisher relation

The scaling of disk galaxy rotation velocity with baryonic mass (the "Baryonic Tully-Fisher" relation, BTF) has long confounded galaxy formation models. It is steeper than the M ~ V^3 scaling relating halo virial masses and circular velocities and its zero point implies that galaxies comprise a very small fraction of available baryons. Such low galaxy formation efficiencies may in principle be explained by winds driven by evolving stars, but the tightness of the BTF relation argues against the substantial scatter expected from such vigorous feedback mechanism. We use the APOSTLE/EAGLE simulations to show that the BTF relation is well reproduced in LCDM simulations that match the size and number of galaxies as a function of stellar mass. In such models, galaxy rotation velocities are proportional to halo virial velocity and the steep velocity-mass dependence results from the decline in galaxy formation efficiency with decreasing halo mass needed to reconcile the CDM halo mass function with the galaxy luminosity function. Despite the strong feedback, the scatter in the simulated BTF is smaller than observed, even when considering all simulated galaxies and not just rotationally-supported ones. The simulations predict that the BTF should become increasingly steep at the faint end, although the velocity scatter at fixed mass should remain small. Observed galaxies with rotation speeds below ~40 km/s seem to deviate from this prediction. We discuss observational biases and modeling uncertainties that may help to explain this disagreement in the context of LCDM models of dwarf galaxy formation.

preprint2016arXiv

The oldest and most metal poor stars in the APOSTLE Local Group simulations

We examine the spatial distribution of the oldest and most metal poor stellar populations of Milky Way-sized galaxies using the APOSTLE cosmological hydrodynamical simulations of the Local Group. In agreement with earlier work, we find strong radial gradients in the fraction of the oldest (tform < 0.8 Gyr) and most metal poor ([Fe/H]< -2.5) stars, both of which increase outwards. The most metal poor stars form over an extended period of time; half of them form after z = 5.3, and the last 10% after z = 2.8. The age of the metal poor stellar population also shows significant variation with environment; a high fraction of them are old in the galaxy's central regions and an even higher fraction in some individual dwarf galaxies, with substantial scatter from dwarf to dwarf. Overall, over half of the stars that belong to both the oldest and most metal-poor population are found outside the solar circle. Somewhat counter-intuitively, we find that dwarf galaxies with a large fraction of metal poor stars that are very old are systems where metal poor stars are relatively rare, but where a substantial old population is present. Our results provide guidance for interpreting the results of surveys designed to hunt for the earliest and most pristine stellar component of our Milky Way.

preprint2016arXiv

The origin of compact galaxies with anomalously high black hole masses

Observations of local galaxies harbouring supermassive black holes (BHs) of anomalously high mass, M_BH, relative to their stellar mass, M_star, appear to be at odds with simple models of the co-evolution between galaxies and their central BHs. We study the origin of such outliers in a Lambda cold dark matter context using the EAGLE cosmological, hydrodynamical simulation. We find 15 'M_BH(M_star)-outlier' galaxies, defined as having M_BH more than 1.5 dex above the median M_BH(M_star) relation in the simulation, M_{BH,med}. All M_BH(M_star)-outliers are satellite galaxies, typically with M_star ~ 10^10 M_sun and M_BH ~ 10^8 M_sun. They have all become outliers due to a combination of tidal stripping of their outer stellar component acting over several Gyr and early formation times leading to rapid BH growth at high redshift, with the former mechanism being most important for 67 per cent of these outliers. The same mechanisms also cause the M_BH(M_star)-outlier satellites to be amongst the most compact galaxies in the simulation, making them ideal candidates for ultracompact dwarf galaxy progenitors. The 10 most extreme central galaxies found at z=0 (with log_{10}(M_BH/M_{BH,med}) in [1.2, 1.5]) grow rapidly in M_BH to lie well above the present-day M_BH-M_star relation at early times (z > 2), and either continue to evolve parallel to the z=0 relation or remain unchanged until the present day, making them 'relics' of the high-redshift universe. This high-z formation mechanism may help to explain the origin of observed M_BH(M_star)-outliers with extended dark matter haloes and undisturbed morphologies.

preprint2016arXiv

The origin of the $α$-enhancement of massive galaxies

We study the origin of the stellar $α$-element-to-iron abundance ratio, $[α/\mathrm{Fe}]_{\ast}$, of present-day central galaxies, using cosmological, hydrodynamical simulations from the Evolution and Assembly of GaLaxies and their Environments (EAGLE) project. For galaxies with stellar masses of $M_{\ast} > 10^{10.5}$ M$_{\odot}$, $[α/\mathrm{Fe}]_{\ast}$ increases with increasing galaxy stellar mass and age. These trends are in good agreement with observations of early-type galaxies, and are consistent with a `downsizing' galaxy formation scenario: more massive galaxies have formed the bulk of their stars earlier and more rapidly, hence from an interstellar medium that was mostly $α$-enriched by massive stars. In the absence of feedback from active galactic nuclei (AGN), however, $[α/\mathrm{Fe}]_{\ast}$ in $M_{\ast} > 10^{10.5}$ M$_{\odot}$ galaxies is roughly constant with stellar mass and decreases with mean stellar age, extending the trends found for lower-mass galaxies in both simulations with and without AGN. We conclude that AGN feedback can account for the $α$-enhancement of massive galaxies, as it suppresses their star formation, quenching more massive galaxies at earlier times, thereby preventing the iron from longer-lived intermediate-mass stars (supernova Type Ia) from being incorporated into younger stars.

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

Ubiquitous giant Ly $α$ nebulae around the brightest quasars at $z\sim3.5$ revealed with MUSE

Direct Ly $α$ imaging of intergalactic gas at $z\sim2$ has recently revealed giant cosmological structures around quasars, e.g. the Slug Nebula (Cantalupo et al. 2014). Despite their high luminosity, the detection rate of such systems in narrow-band and spectroscopic surveys is less than 10%, possibly encoding crucial information on the distribution of gas around quasars and the quasar emission properties. In this study, we use the MUSE integral-field instrument to perform a blind survey for giant Ly $α$ nebulae around 17 bright radio-quiet quasars at $3<z<4$ that does not suffer from most of the limitations of previous surveys. After data reduction and analysis performed with specifically developed tools, we found that each quasar is surrounded by giant Ly $α$ nebulae with projected sizes larger than 100 physical kpc and, in some cases, extending up to 320 kpc. The circularly averaged surface brightness profiles of the nebulae appear very similar to each other despite their different morphologies and are consistent with power laws with slopes $\approx-1.8$. The similarity between the properties of all these nebulae and the Slug Nebula suggests a similar origin for all systems and that a large fraction of gas around bright quasars could be in a relatively "cold" (T$\sim$10$^4$K) and dense phase. In addition, our results imply that such gas is ubiquitous within at least 50 kpc from bright quasars at $3<z<4$ independently of the quasar emission opening angle, or extending up to 200 kpc for quasar isotropic emission.

preprint2015arXiv

An assessment of the evidence from ATLAS3D for a variable initial mass function

The ATLAS3D Survey has reported evidence for a non-universal stellar initial mass function (IMF) for early type galaxies (ETGs) (Cappellari et al. 2012, 2013b,a). The IMF was constrained by comparing stellar mass measurements from kinematic data with those from spectral energy distribution (SED) fitting. Here we investigate possible effects of scatter in the reported stellar mass measurements and their potential impact on the IMF determination. We find that a trend of the IMF mismatch parameter with the kinematic mass to light ratio, comparable to the trend observed by Cappellari et al. (2012), could arise if the Gaussian errors of the kinematic mass determination are typically 30%. Without additional data, it is hard to separate between the option that the IMF has a true large intrinsic variation or the option that the errors in the determination are larger than anticipated. A correlation of the IMF with other properties would help to make this distinction, but no strong correlation has been found yet. The strongest correlation is with velocity dispersion. However, it has a large scatter and the correlation depends on sample selection and distance measurements. The correlation with velocity dispersion could be partly caused by the colour-dependent calibration of the surface brightness fluctuation distances of Tonry et al. (2001). We find that the K-band luminosity limited ATLAS3D Survey is incomplete for the highest M/L galaxies below 10^10.3 M_sun. There is a significant IMF - velocity dispersion trend for galaxies with SED masses above this limit, but no trend for galaxies with kinematic masses above this limit. We also find an IMF trend with distance, but no correlation between nearest neighbour ETGs, which excludes a large environmental dependence. Our findings do not rule out the reported IMF variations, but they suggest that further study is needed.

preprint2015arXiv

Baryon effects on the internal structure of LCDM halos in the EAGLE simulations

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

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

Dark matter annihilation radiation in hydrodynamic simulations of Milky Way haloes

We obtain predictions for the properties of cold dark matter annihilation radiation using high resolution hydrodynamic zoom-in cosmological simulations of Milky Way-like galaxies (APOSTLE project) carried out as part of the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) programme. Galactic halos in the simulation have significantly different properties from those assumed in the "standard halo model" often used in dark matter detection studies. The formation of the galaxy causes a contraction of the dark matter halo, whose density profile develops a steeper slope than the Navarro-Frenk-White (NFW) profile between $r\approx1.5$ kpc and $r\approx10$ kpc. At smaller radii, $r\lesssim1.5$ kpc, the halos develop a flatter than NFW slope. This unexpected feature may be specific to our particular choice of subgrid physics model but nevertheless the dark matter density profiles agree within 30% as the mass resolution is increased by a factor 150. The inner regions of the halos are almost perfectly spherical (axis ratios $b/a > 0.97$ within $r=1$ kpc) and there is no offset larger than 45 pc between the centre of the stellar distribution and the centre of the dark halo. The morphology of the predicted dark matter annihilation radiation signal is in broad agreement with $γ$-ray observations at large Galactic latitudes ($b\gtrsim3^\circ$). At smaller angles, the inferred signal in one of our four galaxies is similar to that which is observed but it is significantly weaker in the other three.

preprint2015arXiv

Detection of hot, metal-enriched outflowing gas around $z\approx\,$2.3 star-forming galaxies in the Keck Baryonic Structure Survey

We use quasar absorption lines to study the physical conditions in the circumgalactic medium of redshift $z\approx 2.3$ star-forming galaxies taken from the Keck Baryonic Structure Survey (KBSS). In Turner et al. 2014 we used the pixel optical depth technique to show that absorption by HI and the metal ions OVI, NV, CIV, CIII and SiIV is strongly enhanced within $|Δv|\lesssim170$ km/s and projected distances $|d|\lesssim180$ proper kpc from sightlines to the background quasars. Here we demonstrate that the OVI absorption is also strongly enhanced at fixed HI, CIV, and SiIV optical depths, and that this enhancement extends out to $\sim350$ km/s. At fixed HI the increase in the median OVI optical depth near galaxies is 0.3-0.7 dex and is detected at 2--3-$σ$ confidence for all seven HI bins that have $\log_{10}τ_{\rm HI}\ge-1.5$. We use ionization models to show that the observed strength of OVI as a function of HI is consistent with enriched, photoionized gas for pixels with $τ_{\rm HI}\gtrsim10$. However, for pixels with $τ_{\rm HI} \lesssim 1$ this would lead to implausibly high metallicities at low densities if the gas were photoionized by the background radiation. This indicates that the galaxies are surrounded by gas that is sufficiently hot to be collisionally ionized ($T > 10^5\,$K) and that a substantial fraction of the hot gas has a metallicity $\gtrsim 10^{-1}$ of solar. Given the high metallicity and large velocity extent (out to $\sim1.5\times v_{\rm circ}$) of this gas, we conclude that we have detected hot, metal enriched outflows arising from star-forming galaxies.

preprint2015arXiv

Galaxies that Shine: radiation-hydrodynamical simulations of disk galaxies

Radiation feedback is typically implemented using subgrid recipes in hydrodynamical simulations of galaxies. Very little work has so far been performed using radiation-hydrodynamics (RHD), and there is no consensus on the importance of radiation feedback in galaxy evolution. We present RHD simulations of isolated galaxy disks of different masses with a resolution of 18 pc. Besides accounting for supernova feedback, our simulations are the first galaxy-scale simulations to include RHD treatments of photo-ionisation heating and radiation pressure, from both direct optical/UV radiation and multi-scattered, re-processed infrared (IR) radiation. Photo-heating smooths and thickens the disks and suppresses star formation about as much as the inclusion of ("thermal dump") supernova feedback does. These effects decrease with galaxy mass and are mainly due to the prevention of the formation of dense clouds, as opposed to their destruction. Radiation pressure, whether from direct or IR radiation, has little effect, but for the IR radiation we show that its impact is limited by our inability to resolve the high optical depths for which multi-scattering becomes important. While artificially boosting the IR optical depths does reduce the star formation, it does so by smoothing the gas rather than by generating stronger outflows. We conclude that although higher-resolution simulations, and potentially also different supernova implementations, are needed for confirmation, our findings suggest that radiation feedback is more gentle and less effective than is often assumed in subgrid prescriptions.

preprint2015arXiv

Intrinsic alignments of galaxies in the EAGLE and cosmo-OWLS simulations

We report results for the alignments of galaxies in the EAGLE and cosmo-OWLS simulations as a function of galaxy separation and halo mass. The combination of these hydro-cosmological simulations enables us to span four orders of magnitude in halo mass ($10.7<log_{10}(M_{200}/[h^{-1}M_\odot])<15$) and a large range of separations ($-1<log_{10}(r/[h^{-1}Mpc])< 2$). We focus on two classes of alignments: the orientations of galaxies with respect to either the directions to, or the orientations of, surrounding galaxies. We find that the strength of the alignment is a strongly decreasing function of the distance between galaxies. The orientation-direction alignment can remain significant up to ~100 Mpc, for galaxies hosted by the most massive haloes in our simulations. Galaxies hosted by more massive subhaloes show stronger alignment. At a fixed halo mass, more aspherical or prolate galaxies exhibit stronger alignments. The spatial distribution of satellites is anisotropic and significantly aligned with the major axis of the main host halo. The major axis of satellite galaxies, when all stars are considered, are preferentially aligned towards the centre of the main host halo. The predicted projected direction-orientation alignment, $ε_{g+}(r_{p})$, is in broad agreement with recent observations when only stars within the typical observable extent of a galaxy are used to define galaxy orientations. We find that the orientation-orientation alignment is weaker than the orientation-direction alignment on all scales. Overall, the strength of galaxy alignments depends strongly on the subset of stars that are used to measure the orientations of galaxies and it is always weaker than the alignment of the dark matter haloes. Thus, alignment models that use halo orientation as a direct proxy for galaxy orientation will overestimate the impact of intrinsic alignments on weak lensing analyses.

preprint2015arXiv

Metal-line absorption around $z\approx$2.4 star-forming galaxies in the Keck Baryonic Structure Survey

(Abridged) We study metal absorption around 854 $z\approx$2.4 star-forming galaxies taken from the Keck Baryonic Structure Survey (KBSS). The galaxies examined in this work lie in the fields of 15 hyper-luminous background QSOs, with galaxy impact parameters ranging from 35 proper kpc (pkpc) to 2 proper Mpc (pMpc). Using the pixel optical depth technique, we present the first galaxy-centred 2-D maps of the median absorption by OVI, NV, CIV, CIII, and SIV, as well as updated results for HI. At small galactocentric radii we detect a strong enhancement of the absorption relative to randomly located regions that extend out to at least 180 pkpc in the transverse direction, and $\pm$240 km/s along the line-of-sight (LOS, ~1 pMpc in the case of pure Hubble flow) for all ions except NV. For CIV (and HI) we detect a significant enhancement of the absorption signal out to 2 pMpc in the transverse direction, corresponding to the maximum impact parameter in our sample. After normalising the median absorption profiles to account for variations in line strengths and detection limits, in the transverse direction we find no evidence for a sharp drop-off in metals distinct from that of HI. We argue instead that non-detection of some metal line species in the extended circumgalactic medium is consistent with differences in the detection sensitivity. We also present measurements of covering fractions and equivalent widths as a function of projected galaxy distance. Limiting the sample to the 340 galaxies with redshifts measured from nebular emission lines does not decrease the extent of the enhancement along the LOS compared to that in the transverse direction. This rules out redshift errors as the source of the observed redshift-space anisotropy and thus implies that we have detected the signature of gas peculiar velocities from infall, outflows, or virial motions for HI, OVI, CIV, CIII, and CIV.

preprint2015arXiv

Molecular hydrogen abundances of galaxies in the EAGLE simulations

We investigate the abundance of galactic molecular hydrogen (H$_2$) in the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) cosmological hydrodynamic simulations. We assign H$_2$ masses to gas particles in the simulations in post-processing using two different prescriptions that depend on the local dust-to-gas ratio and the interstellar radiation field. Both result in H$_2$ galaxy mass functions that agree well with observations in the local and high-redshift Universe. The simulations reproduce the observed scaling relations between the mass of H$_2$ and the stellar mass, star formation rate and stellar surface density. Towards high edshifts, galaxies in the simulations display larger H$_2$ mass fractions, and correspondingly lower H$_2$ depletion timescales, also in good agreement with observations. The comoving mass density of H$_2$ in units of the critical density, $Ω_{\rm H_2}$, peaks at $z\approx 1.2-1.5$, later than the predicted peak of the cosmic star formation rate activity, at $z\approx 2$. This difference stems from the decrease in gas metallicity and increase in interstellar radiation field with redshift, both of which hamper H$_2$ formation. We find that the cosmic H$_2$ budget is dominated by galaxies with $M_{\rm H_2}>10^9\,\rm M_{\odot}$, star formation rates $>10\,\rm M_{\odot}\,\rm yr^{-1}$ and stellar masses $M_{\rm stellar}>10^{10}\,\rm M_{\odot}$, which are readily observable in the optical and near-IR. The match between the H$_2$ properties of galaxies that emerge in the simulations and observations is remarkable, particularly since H$_2$ observations were not used to adjust parameters in EAGLE.

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

Recycled stellar ejecta as fuel for star formation and implications for the origin of the galaxy mass-metallicity relation

We use cosmological, hydrodynamical simulations from the EAGLE and OWLS projects to assess the significance of recycled stellar ejecta as fuel for star formation. The fractional contributions of stellar mass loss to the cosmic star formation rate (SFR) and stellar mass densities increase with time, reaching $35 \%$ and $19 \%$, respectively, at $z=0$. The importance of recycling increases steeply with galaxy stellar mass for $M_{\ast} < 10^{10.5}$ M$_{\odot}$, and decreases mildly at higher mass. This trend arises from the mass dependence of feedback associated with star formation and AGN, which preferentially suppresses star formation fuelled by recycling. Recycling is more important for satellites than centrals and its contribution decreases with galactocentric radius. The relative contribution of AGB stars increases with time and towards galaxy centers. This is a consequence of the more gradual release of AGB ejecta compared to that of massive stars, and the preferential removal of the latter by star formation-driven outflows and by lock up in stellar remnants. Recycling-fuelled star formation exhibits a tight, positive correlation with galaxy metallicity, with a secondary dependence on the relative abundance of alpha elements (which are predominantly synthesized in massive stars), that is insensitive to the subgrid models for feedback. Hence, our conclusions are directly relevant for the origin of the mass-metallicity relation and metallicity gradients. Applying the relation between recycling and metallicity to the observed mass-metallicity relation yields our best estimate of the mass-dependent contribution of recycling. For centrals with a mass similar to that of the Milky Way, we infer the contributions of recycled stellar ejecta to the SFR and stellar mass to be $35 \%$ and $20 \%$, respectively.

preprint2015arXiv

Spatially adaptive radiation-hydrodynamical simulations of galaxy formation during cosmological reionization

We present a suite of cosmological radiation-hydrodynamical simulations of the assembly of galaxies driving the reionization of the intergalactic medium (IGM) at z >~ 6. The simulations account for the hydrodynamical feedback from photoionization heating and the explosion of massive stars as supernovae (SNe). Our reference simulation, which was carried out in a box of size 25 comoving Mpc/h using 2 x 512^3 particles, produces a reasonable reionization history and matches the observed UV luminosity function of galaxies. Simulations with different box sizes and resolutions are used to investigate numerical convergence, and simulations in which either SNe or photoionization heating or both are turned off, are used to investigate the role of feedback from star formation. Ionizing radiation is treated using accurate radiative transfer at the high spatially adaptive resolution at which the hydrodynamics is carried out. SN feedback strongly reduces the star formation rates (SFRs) over nearly the full mass range of simulated galaxies and is required to yield SFRs in agreement with observations. Photoheating helps to suppress star formation in low-mass galaxies, but its impact on the cosmic SFR is small. Because the effect of photoheating is masked by the strong SN feedback, it does not imprint a signature on the UV galaxy luminosity function, although we note that our resolution is insufficient to model star-forming minihaloes cooling through molecular hydrogen transitions. Photoheating does provide a strong positive feedback on reionization because it smooths density fluctuations in the IGM, which lowers the IGM recombination rate substantially. Our simulations demonstrate a tight non-linear coupling of galaxy formation and reionization, motivating the need for the accurate and simultaneous inclusion of photoheating and SN feedback in models of the early Universe.

preprint2015arXiv

The accretion history of dark matter halos I: The physical origin of the universal function

Understanding the universal accretion history of dark matter halos is the first step towards determining the origin of their structure. We use the extended Press-Schechter formalism to derive the halo mass accretion history from the growth rate of initial density perturbations. We show that the halo mass history is well described by an exponential function of redshift in the high-redshift regime. However, in the low-redshift regime the mass history follows a power law because the growth of density perturbations is halted in the dark energy dominated era due to the accelerated expansion of the Universe. We provide an analytic model that follows the expression $M(z) = M_{0}(1+z)^{af(M_{0})}e^{-f(M_{0})z}$, where $M_{0} = M(z=0)$, $a$ depends on cosmology and $f(M_{0})$ depends only on the linear matter power spectrum. The analytic model does not rely on calibration against numerical simulations and is suitable for any cosmology. We compare our model with the latest empirical models for the mass accretion history in the literature and find very good agreement. We provide numerical routines for the model online.

preprint2015arXiv

The accretion history of dark matter halos II: The connections with the mass power spectrum and the density profile

We explore the relation between the structure and mass accretion histories of dark matter halos using a suite of cosmological simulations. We confirm that the formation time, defined as the time when the virial mass of the main progenitor equals the mass enclosed within the scale radius, correlates strongly with concentration. We provide a semi-analytic model for halo mass history that combines analytic relations with fits to simulations. This model has the functional form, $M(z) = M_{0}(1+z)^αe^{βz}$, where the parameters $α$ and $β$ are directly correlated with concentration. We then combine this model for the halo mass history with the analytic relations between $α$, $β$ and the linear power spectrum derived by Correa et al. (2014) to establish the physical link between halo concentration and the initial density perturbation field. Finally, we provide fitting formulas for the halo mass history as well as numerical routines, we derive the accretion rate as a function of halo mass, and we demonstrate how the halo mass history depends on cosmology and the adopted definition of halo mass.

preprint2015arXiv

The accretion history of dark matter halos III: A physical model for the concentration-mass relation

We present a semi-analytic, physically motivated model for dark matter halo concentration as a function of halo mass and redshift. The semi-analytic model combines an analytic model for the halo mass accretion history (MAH), based on extended Press Schechter (EPS) theory, with an empirical relation between concentration and formation time obtained through fits to the results of numerical simulations. Because the semi-analytic model is based on EPS theory, it can be applied to wide ranges in mass, redshift and cosmology. The resulting concentration-mass (c-M) relations are found to agree with the simulations, and because the model applies only to relaxed halos, they do not exhibit the upturn at high masses or high redshifts found by some recent works. We predict a change of slope in the z=0 c-M relation at a mass scale of $10^{11}\rm{M}_{\odot}$. We find that this is due to the change in the functional form of the halo MAH, which goes from being dominated by an exponential (for high-mass halos) to a power-law (for low-mass halos). During the latter phase, the core radius remains approximately constant, and the concentration grows due to the drop of the background density. We also analyse how the c-M relation predicted by this work affects the power produced by dark matter annihilation, finding that at z = 0 the power is two orders of magnitude lower than that obtained from extrapolating best-fitting c-M relations. We provide fits to the c-M relations as well as numerical routines to compute concentrations and MAHs.

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 APOSTLE simulations: solutions to the Local Group's cosmic puzzles

The Local Group of galaxies offer some of the most discriminating tests of models of cosmic structure formation. For example, observations of the Milky Way (MW) and Andromeda satellite populations appear to be in disagreement with N-body simulations of the "Lambda Cold Dark Matter" (ΛCDM) model: there are far fewer satellite galaxies than substructures in cold dark matter halos (the "missing satellites" problem); dwarf galaxies seem to avoid the most massive substructures (the "too-big-to-fail" problem); and the brightest satellites appear to orbit their host galaxies on a thin plane (the "planes of satellites" problem). Here we present results from APOSTLE (A Project Of Simulating The Local Environment), a suite of cosmological hydrodynamic simulations of twelve volumes selected to match the kinematics of the Local Group (LG) members. Applying the Eagle code to the LG environment, we find that our simulations match the observed abundance of LG galaxies, including the satellite galaxies of the MW and Andromeda. Due to changes to the structure of halos and the evolution in the LG environment, the simulations reproduce the observed relation between stellar mass and velocity dispersion of individual dwarf spheroidal galaxies without necessitating the formation of cores in their dark matter profiles. Satellite systems form with a range of spatial anisotropies, including one similar to that of the MW, confirming that such a configuration is not unexpected in ΛCDM. Finally, based on the observed velocity dispersion, size, and stellar mass, we provide new estimates of the maximum circular velocity for the halos of nine MW dwarf spheroidals.

preprint2015arXiv

The broadening of Lyman-alpha forest absorption lines

We provide an analytical description of the line broadening of HI absorbers in the Lyman-alpha forest resulting from Doppler broadening and Jeans smoothing. We demonstrate that our relation captures the dependence of the line-width on column density for narrow lines in z~3 mock spectra remarkably well. Broad lines at a given column density arise when the underlying density structure is more complex, and such clustering is not captured by our model. Our understanding of the line broadening opens the way to a new method to characterise the thermal state of the intergalactic medium and to determine the sizes of the absorbing structures.

preprint2015arXiv

The contributions of matter inside and outside of haloes to the matter power spectrum

Halo-based models have been successful in predicting the clustering of matter. However, the validity of the postulate that the clustering is fully determined by matter inside haloes remains largely untested, and it is not clear a priori whether non-virialised matter might contribute significantly to the non-linear clustering signal. Here, we investigate the contribution of haloes to the matter power spectrum as a function of both scale and halo mass by combining a set of cosmological N-body simulations to calculate the contributions of different spherical overdensity regions, Friends-of-Friends (FoF) groups and matter outside haloes to the power spectrum. We find that matter inside spherical overdensity regions of size R200,mean cannot account for all power for 1<k<100 h/Mpc, regardless of the minimum halo mass. At most, it accounts for 95% of the power (k>20 h/Mpc). For 2<k<10 h/Mpc, haloes with mass M200,mean<10^11 Msun/h contribute negligibly to the power spectrum, and our results appear to be converged with decreasing halo mass. When haloes are taken to be regions of size R200,crit, the amount of power unaccounted for is larger on all scales. Accounting also for matter inside FoF groups but outside R200,mean increases the contribution of halo matter on most scales probed here by 5-15%. Matter inside FoF groups with M200,mean>10^9 Msun/h accounts for essentially all power for 3<k<100 h/Mpc. We therefore expect halo models that ignore the contribution of matter outside R200,mean to overestimate the contribution of haloes of any mass to the power on small scales (k>1 h/Mpc).

preprint2015arXiv

The distribution of atomic hydrogen in EAGLE galaxies: morphologies, profiles, and HI holes

We compare the mass and internal distribution of atomic hydrogen (HI) in 2200 present-day central galaxies with M_star > 10^10 M_Sun from the 100 Mpc EAGLE Reference simulation to observational data. Atomic hydrogen fractions are corrected for self-shielding using a fitting formula from radiative transfer simulations and for the presence of molecular hydrogen using an empirical or a theoretical prescription from the literature. The resulting neutral hydrogen fractions, M_(HI+H2)/M_star, agree with observations to better than 0.1 dex for galaxies with M_star between 10^10 and 10^11 M_Sun. Our fiducial, empirical H2 model based on gas pressure results in galactic HI mass fractions, M_HI/M_star, that agree with observations from the GASS survey to better than 0.3 dex, but the alternative theoretical H2 formula leads to a negative offset in M_HI/M_star of up to 0.5 dex. Visual inspection reveals that most HI disks in simulated HI-rich galaxies are vertically disturbed, plausibly due to recent accretion events. Many galaxies (up to 80 per cent) contain spuriously large HI holes, which are likely formed as a consequence of the feedback implementation in EAGLE. The HI mass-size relation of all simulated galaxies is close to (but 16 per cent steeper than) observed, and when only galaxies without large holes in the HI disc are considered, the agreement becomes excellent (better than 0.1 dex). The presence of large HI holes also makes the radial HI surface density profiles somewhat too low in the centre, at Σ_HI > 1 M_Sun pc^-2 (by a factor of <~ 2 compared to data from the Bluedisk survey). In the outer region (Σ_HI < 1 M_Sun pc^-2), the simulated profiles agree quantitatively with observations. Scaled by HI size, the simulated profiles of HI-rich (M_HI > 10^9.8 M_Sun) and control galaxies (10^9.1 M_Sun > M_HI > 10^9.8 M_Sun) follow each other closely, as observed. (Abridged)

preprint2015arXiv

The distribution of neutral hydrogen around high-redshift galaxies and quasars in the EAGLE simulation

The observed high covering fractions of neutral hydrogen (HI) with column densities above $\sim 10^{17} \rm{cm}^{-2}$ around Lyman-Break Galaxies (LBGs) and bright quasars at redshifts z ~ 2-3 has been identified as a challenge for simulations of galaxy formation. We use the EAGLE cosmological, hydrodynamical simulation, which has been shown to reproduce a wide range of galaxy properties and for which the subgrid feedback was calibrated without considering gas properties, to study the distribution of HI around high-redshift galaxies. We predict the covering fractions of strong HI absorbers ($N_{\rm{HI}} \gtrsim 10^{17} \rm{cm}^{-2}$) inside haloes to increase rapidly with redshift but to depend only weakly on halo mass. For massive ($M_{200} \gtrsim 10^{12} {\rm M_{\odot}}$) halos the covering fraction profiles are nearly scale-invariant and we provide fitting functions that reproduce the simulation results. While efficient feedback is required to increase the HI covering fractions to the high observed values, the distribution of strong absorbers in and around halos of a fixed mass is insensitive to factor of two variations in the strength of the stellar feedback. In contrast, at fixed stellar mass the predicted HI distribution is highly sensitive to the feedback efficiency. The fiducial EAGLE simulation reproduces both the observed global column density distribution function of HI and the observed radial covering fraction profiles of strong HI absorbers around LBGs and bright quasars.

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 effect of baryons on the inner density profiles of rich clusters

We use the "Evolution and assembly of galaxies and their environments" (EAGLE) cosmological simulation to investigate the effect of baryons on the density profiles of rich galaxy clusters. We focus on EAGLE clusters with $M_{200}>10^{14}~M_\odot$ of which we have six examples. The central brightest cluster galaxies (BCGs) in the simulation have steep stellar density profiles, $ρ_*(r) \propto r^{-3}$. Stars dominate the mass density for $r<10~\rm{kpc}$, and, as a result, the total mass density profiles are steeper than the Navarro-Frenk-White (NFW) profile, in remarkable agreement with observations. The dark matter halo itself closely follows the NFW form at all resolved radii ($r\gtrsim3.0~\rm{kpc}$). The EAGLE BCGs have similar surface brightness and line-of-sight velocity dispersion profiles as the BCGs in the sample of Newman et al., which have the most detailed measurements currently available. After subtracting the contribution of the stars to the central density, Newman et al. infer significantly shallower slopes than the NFW value, in contradiction with the EAGLE results. We discuss possible reasons for this discrepancy, and conclude that an inconsistency between the kinematical model adopted by Newman et al. for their BCGs, which assumes isotropic stellar orbits, and the kinematical structure of the EAGLE BCGs, in which the orbital stellar anisotropy varies with radius and tends to be radially biased, could explain at least part of the discrepancy.

preprint2015arXiv

The link between the assembly of the inner dark matter halo and the angular momentum evolution of galaxies in the EAGLE simulation

We explore the co-evolution of the specific angular momentum of dark matter haloes and the cold baryons that comprise the galaxies within. We study over two thousand central galaxies within the reference cosmological hydrodynamical simulation of the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) project. We employ a methodology within which the evolutionary history of a system is specified by the time-evolving properties of the Lagrangian particles that define it at z=0. We find a strong correlation between the evolution of the specific angular momentum of today's stars (cold gas) and that of the inner (whole) dark matter halo they are associated with. This link is particularly strong for the stars formed before the epoch of maximum expansion and subsequent collapse of the central dark matter halo (turnaround). Spheroids are typically assembled primarily from stars formed prior to turnaround, and are therefore destined to suffer a net loss of angular momentum associated with the strong merging activity during the assembly of the inner dark matter halo. Stellar discs retain their specific angular momentum since they are comprised of stars formed mainly after turnaround, from gas that mostly preserves the high specific angular momentum it acquired by tidal torques during the linear growth of the halo. Since the specific angular momentum loss of the stars is tied to the galaxy's morphology today, it may be possible to use our results to predict, statistically, the assembly history of a halo given the morphology of the galaxy it hosts.

preprint2015arXiv

The unexpected diversity of dwarf galaxy rotation curves

We examine the circular velocity profiles of galaxies in ΛCDM cosmological hydrodynamical simulations from the EAGLE and LOCAL GROUPS projects and compare them with a compilation of observed rotation curves of galaxies spanning a wide range in mass. The shape of the circular velocity profiles of simulated galaxies varies systematically as a function of galaxy mass, but shows remarkably little variation at fixed maximum circular velocity. This is especially true for low-mass dark matter-dominated systems, reflecting the expected similarity of the underlying cold dark matter haloes. This is at odds with observed dwarf galaxies, which show a large diversity of rotation curve shapes, even at fixed maximum rotation speed. Some dwarfs have rotation curves that agree well with simulations, others do not. The latter are systems where the inferred mass enclosed in the inner regions is much lower than expected for cold dark matter haloes and include many galaxies where previous work claims the presence of a constant density "core". The "cusp vs core" issue is thus better characterized as an "inner mass deficit" problem than as a density slope mismatch. For several galaxies the magnitude of this inner mass deficit is well in excess of that reported in recent simulations where cores result from baryon-induced fluctuations in the gravitational potential. We conclude that one or more of the following statements must be true: (i) the dark matter is more complex than envisaged by any current model; (ii) current simulations fail to reproduce the effects of baryons on the inner regions of dwarf galaxies; and/or (iii) the mass profiles of "inner mass deficit" galaxies inferred from kinematic data are incorrect.

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

Constraining the epoch of reionization with the variance statistic: simulations of the LOFAR case

Several experiments are underway to detect the cosmic redshifted 21-cm signal from neutral hydrogen from the Epoch of Reionization (EoR). Due to their very low signal-to-noise ratio, these observations aim for a statistical detection of the signal by measuring its power spectrum. We investigate the extraction of the variance of the signal as a first step towards detecting and constraining the global history of the EoR. Signal variance is the integral of the signal's power spectrum, and it is expected to be measured with a high significance. We demonstrate this through results from a simulation and parameter estimation pipeline developed for the Low Frequency Array (LOFAR)-EoR experiment. We show that LOFAR should be able to detect the EoR in 600 hours of integration using the variance statistic. Additionally, the redshift ($z_r$) and duration ($Δz$) of reionization can be constrained assuming a parametrization. We use an EoR simulation of $z_r = 7.68$ and $Δz = 0.43$ to test the pipeline. We are able to detect the simulated signal with a significance of 4 standard deviations and extract the EoR parameters as $z_r = 7.72^{+0.37}_{-0.18}$ and $Δz = 0.53^{+0.12}_{-0.23}$ in 600 hours, assuming that systematic errors can be adequately controlled. We further show that the significance of detection and constraints on EoR parameters can be improved by measuring the cross-variance of the signal by cross-correlating consecutive redshift bins.

preprint2014arXiv

Local Group galaxies emerge from the dark

The "Lambda Cold Dark Matter" (LCDM) model of cosmic structure formation is eminently falsifiable: once its parameters are fixed on large scales, it becomes testable in the nearby Universe. Observations within our Local Group of galaxies, including the satellite populations of the Milky Way and Andromeda, appear to contradict LCDM predictions: there are far fewer satellite galaxies than dark matter halos (the "missing satellites" problem), galaxies seem to avoid the largest substructures (the "too big to fail" problem), and the brightest satellites appear to orbit their host galaxies on a thin plane (the "planes of satellites" problem). We present results from the first hydrodynamic simulations of the Local Group that match the observed abundance of galaxies. We find that when baryonic and dark matter are followed simultaneously in the context of a realistic galaxy formation model, all three "problems" are resolved within the LCDM paradigm.

preprint2014arXiv

Non-equilibrium chemistry and cooling in the diffuse interstellar medium - I. Optically thin regime

An accurate treatment of the multiphase interstellar medium (ISM) in hydrodynamic galaxy simulations requires that we follow not only the thermal evolution of the gas, but also the evolution of its chemical state, including its molecular chemistry, without assuming chemical (including ionisation) equilibrium. We present a reaction network that can be used to solve for this thermo-chemical evolution. Our model follows the evolution of all ionisation states of the 11 elements that dominate the cooling rate, along with important molecules such as H2 and CO, and the intermediate molecular species that are involved in their formation (20 molecules in total). We include chemical reactions on dust grains, thermal processes involving dust, cosmic ray ionisation and heating and photochemical reactions. We focus on conditions typical for the diffuse ISM, with densities of 10^-2 cm^-3 < nH < 10^4 cm^-3 and temperatures of 10^2 K < T < 10^4 K, and we consider a range of radiation fields, including no UV radiation. In this paper we consider only gas that is optically thin, while paper II considers gas that becomes shielded from the radiation field. We verify the accuracy of our model by comparing chemical abundances and cooling functions in chemical equilibrium with the photoionisation code Cloudy. We identify the major coolants in diffuse interstellar gas to be CII, SiII and FeII, along with OI and H2 at densities nH > 10^2 cm^-3. Finally, we investigate the impact of non-equilibrium chemistry on the cooling functions of isochorically or isobarically cooling gas. We find that, at T < 10^4 K, recombination lags increase the electron abundance above its equilibrium value at a given temperature, which can enhance the cooling rate by up to two orders of magnitude. The cooling gas also shows lower H2 abundances than in equilibrium, by up to an order of magnitude.

preprint2014arXiv

Non-equilibrium chemistry and cooling in the diffuse interstellar medium - II. Shielded gas

We extend the non-equilibrium model for the chemical and thermal evolution of diffuse interstellar gas presented in Richings et al. (2014) to account for shielding from the UV radiation field. We attenuate the photochemical rates by dust and by gas, including absorption by HI, H2, HeI, HeII and CO where appropriate. We then use this model to investigate the dominant cooling and heating processes in interstellar gas as it becomes shielded from the UV radiation. We consider a one-dimensional plane-parallel slab of gas irradiated by the interstellar radiation field, either at constant density and temperature or in thermal and pressure equilibrium. The dominant thermal processes tend to form three distinct regions in the clouds. At low column densities cooling is dominated by ionised metals such as SiII, FeII, FeIII and CII, which are balanced by photoheating, primarily from HI. Once the hydrogen-ionising radiation becomes attenuated by neutral hydrogen, photoelectric dust heating dominates, while CII becomes dominant for cooling. Finally, dust shielding triggers the formation of CO and suppresses photoelectric heating. The dominant coolants in this fully shielded region are H2 and CO. The column density of the HI-H2 transition predicted by our model is lower at higher density (or at higher pressure for gas clouds in pressure equilibrium) and at higher metallicity, in agreement with previous PDR models. We also compare the HI-H2 transition in our model to two prescriptions for molecular hydrogen formation that have been implemented in hydrodynamic simulations.

preprint2014arXiv

The chosen few: the low mass halos that host faint galaxies

Since reionization prevents star formation in most halos below 3 x 10^9 solar masses, dwarf galaxies only populate a fraction of existing dark matter halos. We use hydrodynamic cosmological simulations of the Local Group to study the discriminating factors for galaxy formation in the early Universe and connect them to the present-day properties of galaxies and halos. A combination of selection effects related to reionization, and the subsequent evolution of halos in different environments, introduces strong biases between the population of halos that host dwarf galaxies, and the total halo population. Halos that host galaxies formed earlier and are more concentrated. In addition, halos more affected by tidal stripping are more likely to host a galaxy for a given mass or maximum circular velocity, vmax, today. Consequently, satellite halos are populated more frequently than field halos, and satellite halos of 10^8 - 10^9 solar masses or vmax of 12 - 20 km/s, similar to the Local Group dwarf spheroidals, have experienced a greater than average reduction in both mass and vmax after infall. They are on closer, more radial orbits with higher infall velocities and earlier infall times. Together, these effects make dwarf galaxies highly biased tracers of the underlying dark matter distribution.

preprint2014arXiv

The EAGLE project: Simulating the evolution and assembly of galaxies and their environments

We introduce the Virgo Consortium's EAGLE project, a suite of hydrodynamical simulations that follow the formation of galaxies and black holes in representative volumes. We discuss the limitations of such simulations in light of their finite resolution and poorly constrained subgrid physics, and how these affect their predictive power. One major improvement is our treatment of feedback from massive stars and AGN in which thermal energy is injected into the gas without the need to turn off cooling or hydrodynamical forces, allowing winds to develop without predetermined speed or mass loading factors. Because the feedback efficiencies cannot be predicted from first principles, we calibrate them to the z~0 galaxy stellar mass function and the amplitude of the galaxy-central black hole mass relation, also taking galaxy sizes into account. The observed galaxy mass function is reproduced to $\lesssim 0.2$ dex over the full mass range, $10^8 < M_*/M_\odot \lesssim 10^{11}$, a level of agreement close to that attained by semi-analytic models, and unprecedented for hydrodynamical simulations. We compare our results to a representative set of low-redshift observables not considered in the calibration, and find good agreement with the observed galaxy specific star formation rates, passive fractions, Tully-Fisher relation, total stellar luminosities of galaxy clusters, and column density distributions of intergalactic CIV and OVI. While the mass-metallicity relations for gas and stars are consistent with observations for $M_* \gtrsim 10^9 M_\odot$, they are insufficiently steep at lower masses. The gas fractions and temperatures are too high for clusters of galaxies, but for groups these discrepancies can be resolved by adopting a higher heating temperature in the subgrid prescription for AGN feedback. EAGLE constitutes a valuable new resource for studies of galaxy formation.

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

A high molecular fraction in a sub-damped absorber at z=0.56

Measuring rest-frame ultraviolet rotational transitions from the Lyman and Werner bands in absorption against a bright background continuum is one of the few ways to directly measure molecular hydrogen (H2). Here we report the detection of Lyman-Werner absorption from H2 at z=0.56 in a sub-damped Ly-alpha system with neutral hydrogen column density N(HI) = 10^(19.5 +/- 0.2) cm^-2. This is the first H2 system analysed at a redshift < 1.5 beyond the Milky Way halo. It has a surprisingly high molecular fraction: log f(H2) > -1.93 +/- 0.36 based on modelling the line profiles, with a robust model-independent lower limit of f(H2) > 10^-3. This is higher than f(H2) values seen along sightlines with similar N(HI) through the Milky Way disk and the Magellanic clouds. The metallicity of the absorber is 0.19 +0.21 -0.10 solar, with a dust-to-gas ratio < 0.36 times the value in the solar neighbourhood. Absorption from associated low-ionisation metal transitions such as OI and FeII is observed in addition to OVI. Using Cloudy models we show that there are three phases present; a ~100 K phase giving rise to H2, a ~10^4 K phase where most of the low-ionisation metal absorption is produced; and a hotter phase associated with OVI. Based on similarities to high velocity clouds in the Milky Way halo showing H2 and the presence of two nearby galaxy candidates with impact parameters of ~10 kpc, we suggest that the absorber may be produced by a tidally-stripped structure similar to the Magellanic Stream.

preprint2013arXiv

A measurement of galaxy halo mass from the surrounding HI Lyα absorption

We measure the dark matter halo masses of z~2.36 UV color-selected star-forming galaxies by matching the observed median HI Lyα absorption around them, as observed in the spectra of background QSOs, to the absorption around haloes above a given mass in cosmological simulations. Focusing on transverse separations 0-2 pMpc and line of sight separations 154-616 km/s, we find a minimum halo mass of log10(Mmin/Msun)=11.6(+0.2)(-0.2), which is in good agreement with published halo mass estimates from clustering analyses. We verified that the measured halo mass is insensitive to a change in the cosmological parameters (WMAP1 vs. WMAP3) and to the inclusion of strong AGN feedback. One unique strength of this method is that it can be used in narrow field galaxy-QSO surveys, i.e. ~30 x 30 arcseconds. In addition, we find that the observed anisotropy in the 2-D HI Lyα absorption distribution on scales of 1.5-2 pMpc is consistent with being a consequence of large-scale gas infall into the potential wells occupied by galaxies.

preprint2013arXiv

Absorption signatures of warm-hot gas at low redshift: NeVIII

At z<1 a large fraction of the baryons is thought to reside in diffuse gas that has been shock-heated to high temperatures (log T=5-6). Absorption by the 770.41, 780.32 Å doublet of NeVIII in quasar spectra represents a unique tool to study this elusive warm-hot phase. We have developed an analytic model for the properties of NeVIII absorbers that allows for an inhomogeneous metal distribution. Our model agrees with the predictions of a simulation from the OWLS project indicating that the average line-of-sight metal-filling fraction within the absorbing gas is low (c_L~0.1). Most of the NeVIII in our model is produced in low-density, collisionally ionized gas (log n=-6 to -4 cm^-3, log T=5-6). Strong NeVIII absorbers (log N(NeVIII)>14), like those recently detected by HST/COS, are found to arise in higher density gas (log n>-4, log T=5.75). NeVIII cloudlets harbour only 1 per cent of the cosmic baryon budget. The baryon content of the surrounding gas (which has similar densities and temperatures as the NeVIII cloudlets) is a factor 1/c_L higher. We conclude that NeVIII absorbers are robust probes of shock-heated diffuse gas, but that spectra with signal-to-noise ratios S/N>100 would be required to detect the bulk of the baryons in warm-hot gas.

preprint2013arXiv

AGN proximity zone fossils and the delayed recombination of metal lines

We model the time-dependent evolution of metal-enriched intergalactic and circumgalactic gas exposed to the fluctuating radiation field from an active galactic nucleus (AGN). We consider diffuse gas densities (n_H=10^-5-10^-3 cm^-3) exposed to the extra-galactic background (EGB) and initially in thermal equilibrium (T \sim 10^4-10^4.5 K). Once the proximate AGN field turns on, additional photo-ionisation rapidly ionises the HI and metals. The enhanced AGN radiation field turns off after a typical AGN lifetime (tau_AGN=1-20 Myr) and the field returns to the EGB intensity, but the metals remain out of ionisation equilibrium for timescales that can significantly exceed tau_AGN. We define this phase as the AGN proximity zone "fossil" phase and show that high ionisation stages (e.g. OVI, NeVIII, MgX) are in general enhanced, while the abundances of low ions are reduced. In contrast, HI re-equilibrates rapidly (<<tau_AGN) owing to its low neutral fraction at diffuse densities. We demonstrate that metal column densities of intervening gas observed in absorption in quasar sight lines are significantly affected by delayed recombination for a wide range of densities, metallicities, and AGN strengths, lifetimes, and duty cycles. We model the exceptionally strong z=0.9 NeVIII absorbers observed by Tripp et al. (2011) as arising in a possible fossil zone or near a recently turned-on AGN and we demonstrate that at low redshift even moderate strength AGN could significantly enhance the high-ion metal columns in the circumgalactic media of galaxies observed without active AGN. Fossil proximity zones may be very important during the quasar era, z \sim 2-5. AGN proximity zone fossils allow a whole new class of non-equilibrium solutions that may be applicable to a large fraction of observed metal absorbers and which could potentially change the inferred physical conditions and masses of diffuse gases.

preprint2013arXiv

Effect of baryonic feedback on two- and three-point shear statistics: prospects for detection and improved modelling

Accurate knowledge of the effect of feedback from galaxy formation on the matter distribution is a key requirement for future weak lensing experiments. Recent studies using hydrodynamic simulations have shown that different baryonic feedback scenarios lead to significantly different two-point shear statistics. In this paper we extend earlier work to three-point shear statistics. We show that, relative to the predictions of dark matter only models, the amplitude of the signal can be reduced by as much as 30-40% on scales of a few arcminutes. We find that baryonic feedback may affect two- and three-point shear statistics differently and demonstrate that this can be used to assess the fidelity of various feedback models. In particular, upcoming surveys such as Euclid might be able to discriminate between different feedback models by measuring both second- and third-order statistics. Because it will likely remain impossible to predict baryonic feedback with high accuracy from first principles, we argue in favour of phenomenological models that can capture the relevant effects of baryonic feedback processes in addition to changes in cosmology. We construct such a model by modifying the dark matter-only halo model to characterise the generic effects of energetic feedback using a small number of parameters. We use this model to perform a likelihood analysis in a simplified case in which two- and three-point shear statistics are measured between 0.5 and 20 arcmin and in which the amplitude of fluctuations, sigma8, the matter density parameter, Om, and the dark energy parameter, w0, are the only unknown free parameters. We demonstrate that for weak lensing surveys such as Euclid, marginalising over the feedbac parameters describing the effects of baryonic processes, such as outflows driven by feedback from star formation and AGN, may be able to mitigate the bias affecting Om, sigma8 and w0.

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

How the diffuse Universe cools

In this work we investigate the cooling channels of diffuse gas (i.e. n_H<0.1 cm^-3) in cosmology. We aim to identify the wavelengths where most of the energy is radiated in the form of emission lines or continuum radiation, and the main elements and ions responsible for the emission. We use a subset of cosmological, hydrodynamical runs from the OWLS project to calculate the emission of diffuse gas and its evolution with time. We find that at z=0 (z=2) about 70 (80) per cent of the energy emitted by diffuse gas is carried by emission lines, with the continuum radiation contributing the remainder. Hydrogen lines in the Lyman series are the primary contributors to the line emission, with a share of 16 (20) per cent. Oxygen lines are the main metal contributors at high redshift, while silicon, carbon and iron lines are strongest at low redshift, when the contributions of AGB stars and supernova Ia explosions to the metal budget become important and when there is more hot gas. The ionic species carrying the most energy are OIII, CII, CIII, SiII, SiIII, FeII and SIII. The great majority of energy is emitted in the UV band (lambda=100-4000 A), both as continuum radiation and line emission. With almost no exception, all the strongest lines fall in this band. At high energies, continuum radiation is dominant (e.g., 80 per cent in the X-ray band), while lines contribute progressively more at lower energies. While the results do depend on the details of the numerical implementation of the physical processes modelled in the simulations, the comparison of results from different simulations demonstrates that the variations are overall small, and that the conclusions are fairly robust. Given the overwhelming importance of UV emission for the cooling of diffuse gas, it is desirable to build instruments dedicated to the detection and characterisation of diffuse UV emission.

preprint2013arXiv

Non-Equilibirum ionisation and cooling of metal-enriched gas in the presence of a photo-ionisation background

Simulations of the formation of galaxies, as well as ionisation models used to interpret observations of quasar absorption lines, generally either assume ionisation equilibrium or ignore the presence of the extra-galactic background (EGB) radiation. We introduce a method to compute the non-equilibrium ionisation and cooling of diffuse gas exposed to the EGB. Our method iterates the ionisation states of the 11 elements that dominate the cooling and uses tabulated ion-by-ion cooling and photo-heating efficiencies to update the temperature of the gas. Our reaction network includes radiative and di-electric recombination, collisional ionisation, photo-ionisation, Auger ionisation, and charge transfer. We verify that our method reproduces published results for collisional equilibrium, collisional non-equilibrium, and photo-ionised equilibrium. Non-equilibrium effects can become very important in cooling gas, particularly below 10^6 K. Photo-ionisation and non-equilibrium effects both tend to boost the degree of ionisation and to reduce cooling efficiencies. The effect of the EGB is larger for lower densities. Hence, photo-ionisation affects cooling more under isochoric than under isobaric conditions. Non-equilibrium effects are smaller in the presence of the EGB and are thus overestimated when using collisional-only processes. The inclusion of the EGB alters the observational diagnostics of diffuse, metal-enriched gas (e.g. metal absorption lines) even more significantly than the cooling efficiencies. We argue that the cooling efficiency should be considered if ionisation models are used to infer physical conditions from observed line ratios, as the a priori probability of observing gas is lower if its cooling time is shorter. We provide on-line tables of ionisation fractions and cooling efficiencies for equilibrium and non-equilibrium scenarios. (abridged)

preprint2013arXiv

On the connection between the intergalactic medium and galaxies: The HI-galaxy cross-correlation at z < 1

We present a new optical spectroscopic survey of 1777 'star-forming' ('SF') and 366 'non-star-forming' ('non-SF') galaxies at redshifts z < 1 (2143 in total), 22 AGN and 423 stars, observed by instruments such as DEIMOS, VIMOS and GMOS, in 3 fields containing 5 quasi-stellar objects (QSOs) with HST UV spectroscopy. We also present a new spectroscopic survey of 165 'strong' (10^14 < NHI < 10^17 cm^-2), and 489 'weak' (10^13 < NHI < 10^14 cm^-2) intervening HI absorption line systems at z < 1 (654 in total), observed in the spectra of 8 QSOs by COS and FOS on the HST. Combining these new data with previously published galaxy catalogs such as VVDS and GDDS, we have gathered a sample of 654 HI absorption systems and 17509 galaxies at transverse scales < 50 Mpc. We present observational results on the HI-galaxy and galaxy-galaxy correlations at transverse scales r < 10 Mpc, and the HI-HI auto-correlation at transverse scales r < 2 Mpc. The two-point correlation functions are measured both along and transverse to the line-of-sight. We constrain the HI-galaxy statistical connection, as a function of both HI column density and galaxy star-forming activity. Our results are consistent with the following conclusions: (1) the bulk of HI systems on Mpc scales have little velocity dispersion (<120 km/s) with respect to the bulk of galaxies; (2) the vast majority of strong HI systems and SF galaxies are distributed in the same locations, together with 75+-15% of non-SF galaxies, all of which typically reside in dark matter haloes of similar masses; (3) 25+-15% of non-SF galaxies reside in galaxy clusters and are not correlated with strong HI systems at scales < 2 Mpc; and (4) 50% of weak HI systems reside within galaxy voids (hence not correlated with galaxies), and are confined in dark matter haloes of masses smaller than those hosting... [abridged]

preprint2013arXiv

On the evolution of the HI column density distribution in cosmological simulations

We use a set of cosmological simulations combined with radiative transfer calculations to investigate the distribution of neutral hydrogen in the post-reionization Universe. We assess the contributions from the metagalactic ionizing background, collisional ionization and diffuse recombination radiation to the total ionization rate at redshifts z=0-5. We find that the densities above which hydrogen self-shielding becomes important are consistent with analytic calculations and previous work. However, because of diffuse recombination radiation, whose intensity peaks at the same density, the transition between highly ionized and self-shielded regions is smoother than what is usually assumed. We provide fitting functions to the simulated photoionization rate as a function of density and show that post-processing simulations with the fitted rates yields results that are in excellent agreement with the original radiative transfer calculations. The predicted neutral hydrogen column density distributions agree very well with the observations. In particular, the simulations reproduce the remarkable lack of evolution in the column density distribution of Lyman limit and weak damped Lyα systems below z = 3. The evolution of the low column density end is affected by the increasing importance of collisional ionization with decreasing redshift. On the other hand, the simulations predict the abundance of strong damped Lyα systems to broadly track the cosmic star formation rate density.

preprint2013arXiv

Physical properties of simulated galaxy populations at z=2 - I. Effect of metal-line cooling and feedback from star formation and AGN

We use hydrodynamical simulations from the OWLS project to investigate the dependence of the physical properties of galaxy populations at redshift 2 on metal-line cooling and feedback from star formation and active galactic nuclei (AGN). We find that if the sub-grid feedback from star formation is implemented kinetically, the feedback is only efficient if the initial wind velocity exceeds a critical value. This critical velocity increases with galaxy mass and also if metal-line cooling is included. This suggests that radiative losses quench the winds if their initial velocity is too low. If the feedback is efficient, then the star formation rate is inversely proportional to the amount of energy injected per unit stellar mass formed (which is proportional to the initial mass loading for a fixed wind velocity). This can be understood if the star formation is self-regulating, i.e. if the star formation rate (and thus the gas fraction) increase until the outflow rate balances the inflow rate. Feedback from AGN is efficient at high masses, while increasing the initial wind velocity with gas pressure or halo mass allows one to generate galaxy-wide outflows at all masses. Matching the observed galaxy mass function requires efficient feedback. In particular, the predicted faint-end slope is too steep unless we resort to highly mass loaded winds for low-mass objects. Such efficient feedback from low-mass galaxies (M_* << 10^10 Msun) also reduces the discrepancy with the observed specific star formation rates, which are higher than predicted unless the feedback transitions from highly efficient to inefficient just below the observed stellar mass range.

preprint2013arXiv

Physical properties of simulated galaxy populations at z=2 - II. Effects of cosmology, reionization and ISM physics

We use hydrodynamical simulations from the OWLS project to investigate the dependence of the physical properties of galaxy populations at redshift 2 on the assumed star formation law, the equation of state imposed on the unresolved interstellar medium, the stellar initial mass function, the reionization history, and the assumed cosmology. This work complements that of Paper I, where we studied the effects of varying models for galactic winds driven by star formation and AGN. The normalisation of the matter power spectrum strongly affects the galaxy mass function, but has a relatively small effect on the physical properties of galaxies residing in haloes of a fixed mass. Reionization suppresses the stellar masses and gas fractions of low-mass galaxies, but by z = 2 the results are insensitive to the timing of reionization. The stellar initial mass function mainly determines the physical properties of galaxies through its effect on the efficiency of the feedback, while changes in the recycled mass and metal fractions play a smaller role. If we use a recipe for star formation that reproduces the observed star formation law independently of the assumed equation of state of the unresolved ISM, then the latter is unimportant. The star formation law, i.e. the gas consumption time scale as a function of surface density, determines the mass of dense, star-forming gas in galaxies, but affects neither the star formation rate nor the stellar mass. This can be understood in terms of self-regulation: the gas fraction adjusts until the outflow rate balances the inflow rate.

preprint2013arXiv

Soft X-ray and ultra-violet metal-line emission from the gas around galaxies

(Abridged) A large fraction of the gas in galactic haloes has temperatures between 10^4.5 and 10^7 K. At these temperatures, cooling is dominated by metal-line emission if the metallicity Z>~0.1 Zsun. We explore the detectability of several lines using large cosmological, hydrodynamical simulations. We stack surface brightness maps centred on galaxies to calculate the expected mean surface brightness profiles for different halo masses. Assuming a detection limit of 10^-1 photon s^-1 cm^-2 sr^-1, proposed X-ray telescopes can detect O VIII emission from z=0.125 out to 80% of the virial radius (Rvir) of groups and clusters and out to 0.4Rvir for haloes with masses Mhalo=10^12-13 Msun. Emission lines from C VI, N VII, O VII, and Ne X can be detected out to smaller radii, 0.1-0.5Rvir. With a detection limit of 10^-20 erg s^-1 cm^-2 arcsec^-2, future UV telescopes can detect C III emission out to 0.2-0.6Rvir at z=0.25. C IV, O VI, Si III, and Si IV can be seen out to 0.1-0.2Rvir for Mhalo>10^12 Msun. Optical HI H-alpha emission is comparable in strength to C III emission. At z=3 it may be possible to observe C III out to 0.2-0.3Rvir and other rest-frame UV lines out to ~0.1Rvir for Mhalo>10^11 Msun with upcoming optical instruments. Metal-line emission is typically biased towards high density and metallicity and towards the temperature at which the emissivity curve of the corresponding metal line peaks. The bias is similar for the different soft X-ray lines considered, whereas it varies strongly between different UV lines. Active galactic nucleus (AGN) feedback can change the inner surface brightness profiles significantly, but it generally does not change the radius out to which the emission can be observed. Metal-line emission is a promising probe of the warm and hot, enriched gas around galaxies and provides a unique window into the interactions between galaxies and their gaseous haloes.

preprint2013arXiv

The abundance of (not just) dark matter haloes

We study the effect of baryons on the abundance of structures and substructures in a Lambda-CDM cosmology, using a pair of high resolution cosmological simulations from the GIMIC project. Both simulations use identical initial conditions, but while one contains only dark matter, the other also includes baryons. We find that gas pressure, reionisation, supernova feedback, stripping, and truncated accretion systematically reduce the total mass and the abundance of structures below ~10^12 solar masses compared to the pure dark matter simulation. Taking this into account and adopting an appropriate detection threshold lowers the abundance of observed galaxies with maximum circular velocities below 100 km/s, significantly reducing the reported discrepancy between Lambda-CDM and the measured HI velocity function of the ALFALFA survey. We also show that the stellar-to-total mass ratios of galaxies with stellar masses of ~10^5 - 10^7 solar masses inferred from abundance matching of the (sub)halo mass function to the observed galaxy mass function increase by a factor of ~2. In addition, we find that an important fraction of low-mass subhaloes are completely devoid of stars. Accounting for the presence of dark subhaloes below 10^10 solar masses further reduces the abundance of observable objects, and leads to an additional increase in the inferred stellar-to-total mass ratio by factors of 2 - 10 for galaxies in haloes of 10^9 - 10^10 solar masses. This largely reconciles the abundance matching results with the kinematics of individual dwarf galaxies in Lambda-CDM. We propose approximate corrections to the masses of objects derived from pure dark matter calculations to account for baryonic effects.

preprint2013arXiv

The effect of recombination radiation on the temperature and ionization state of partially ionized gas

A substantial fraction of all ionizing photons originate from radiative recombinations. However, in radiative transfer calculations this recombination radiation is often assumed to be absorbed 'on-the-spot' because for most methods the computational cost associated with the inclusion of gas elements as sources is prohibitive. We present a new, CPU and memory efficient implementation for the transport of ionizing recombination radiation in the TRAPHIC radiative transfer scheme. TRAPHIC solves the radiative transfer equation by tracing photon packets at the speed of light and in a photon-conserving manner in spatially adaptive smoothed particle hydrodynamics simulations. Our new implementation uses existing features of the TRAPHIC scheme to add recombination radiation at no additional cost in the limit in which the fraction of the simulation box filled with radiation approaches 1. We test the implementation by simulating an HII region in photoionization equilibrium and comparing to reference solutions presented in the literature, finding excellent agreement. We apply our implementation to discuss the evolution of the HII region to equilibrium. We show that the widely used case A and B approximations yield accurate ionization profiles only near the source and near the ionization front, respectively. We also discuss the impact of recombination radiation on the geometry of shadows behind optically thick absorbers. We demonstrate that the shadow region may be completely ionized by the diffuse recombination radiation field and discuss the important role of heating by recombination radiation in the shadow region.

preprint2013arXiv

The host-galaxy response to the afterglow of GRB 100901A

For Gamma-Ray Burst 100901A, we have obtained Gemini-North and Very Large Telescope optical afterglow spectra at four epochs: one hour, one day, three days and one week after the burst, thanks to the afterglow remaining unusually bright at late times. Apart from a wealth of metal resonance lines, we also detect lines arising from fine-structure levels of the ground state of Fe II, and from metastable levels of Fe II and Ni II at the host redshift (z = 1.4084). These lines are found to vary significantly in time. The combination of the data and modelling results shows that we detect the fall of the Ni II 4 F9/2 metastable level population, which to date has not been observed. Assuming that the population of the excited states is due to the UV-radiation of the afterglow, we estimate an absorber distance of a few hundred pc. This appears to be a typical value when compared to similar studies. We detect two intervening absorbers (z = 1.3147, 1.3179). Despite the wide temporal range of the data, we do not see significant variation in the absorption lines of these two intervening systems.

preprint2013arXiv

The Hot and Energetic Universe: A White Paper presenting the science theme motivating the Athena+ mission

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

preprint2013arXiv

The Hot and Energetic Universe: The missing baryons and the warm-hot intergalactic medium

The backbone of the large-scale structure of the Universe is determined by processes on a cosmological scale and by the gravitational interaction of the dominant dark matter. However, the mobile baryon population shapes the appearance of these structures. Theory predicts that most of the baryons reside in vast unvirialized filamentary structures that connect galaxy groups and clusters, but the observational evidence is currently lacking. Because the majority of the baryons are supposed to exist in a large-scale, hot and dilute gaseous phase, X-rays provide the ideal tool to progress our understanding. Observations with the Athena+ X-ray Integral Field Unit will reveal the location, chemical composition, physical state and dynamics of the active population of baryons.

preprint2013arXiv

The Impact of Different Physical Processes on the Statistics of Lyman-limit and Damped Lyman-alpha Absorbers

We compute the z = 3 neutral hydrogen column density distribution function f(NHI) for 19 simulations drawn from the OWLS project using a post-processing correction for self-shielding calculated with full radiative transfer of the ionising background radiation. We investigate how different physical processes and parameters affect the abundance of Lyman-limit systems (LLSs) and damped Lyman-alpha absorbers (DLAs) including: i) metal-line cooling; ii) the efficiency of feedback from SNe and AGN; iii) the effective equation of state for the ISM; iv) cosmological parameters; v) the assumed star formation law and; vi) the timing of hydrogen reionization . We find that the normalisation and slope, D = d log10 f /d log10 NHI, of f(NHI) in the LLS regime are robust to changes in these physical processes. Among physically plausible models, f(NHI) varies by less than 0.2 dex and D varies by less than 0.18 for LLSs. This is primarily due to the fact that these uncertain physical processes mostly affect star-forming gas which contributes less than 10% to f(NHI) in the the LLS column density range. At higher column densities, variations in f(NHI) become larger (approximately 0.5 dex at NHI = 10^22 cm^-2 and 1.0 dex at NHI = 10^23 cm^-2) and molecular hydrogen formation also becomes important. Many of these changes can be explained in the context of self-regulated star formation in which the amount of star forming gas in a galaxy will adjust such that outflows driven by feedback balance inflows due to accretion. Data and code to reproduce all figures can be found at the following url: https://bitbucket.org/galtay/hi-cddf-owls-1

preprint2013arXiv

The impact of local stellar radiation on the HI column density distribution

It is often assumed that local sources of ionizing radiation have little impact on the distribution of HI in the post-reionization Universe. While this is a good assumption for the IGM, analytic arguments suggest that local sources may typically be more important than the background radiation for high column density absorbers (N_HI > 10^17 /cm^2). We post-process cosmological simulations with accurate radiation transport to investigate the impact of local stellar sources on the HI distribution. We demonstrate that the limited numerical resolution and the simplified treatment of the ISM that are typical of current cosmological simulations provide significant challenges, but that many of the problems can be overcome by taking two steps. First, using ISM particles rather than stellar particles as sources results in a much better sampling of the source distribution. Second, by rescaling the source luminosities so that the amount of radiation escaping into the IGM agrees with that required to produce the observed background, many of the results become insensitive to errors in the predicted fraction of the radiation that escapes the immediate vicinity of the sources. By adopting this strategy and by varying the assumptions about the structure of the unresolved ISM, we conclude that we can robustly estimate the effect of local sources for column densities N_HI << 10^21 /cm^2. However, neither the escape fraction of ionizing radiation nor the effect of local sources on the abundance of N_HI >~ 10^21 systems can be predicted with confidence. We find that local stellar radiation is unimportant for N_HI << 10^17, but that it can affect Lyman Limit and DLA systems. For 10^18 < N_HI < 10^21 the impact of local sources increases with redshift. At z = 5 the abundance of absorbers with N_HI >> 10^17 is substantially reduced, but at z = 0 the reduction only becomes significant for N_HI >~ 10^21 /cm^2.

preprint2013arXiv

The interaction between feedback from active galactic nuclei and supernovae

Energetic feedback from supernovae (SNe) and from active galactic nuclei (AGN) are both important processes that are thought to control how much gas is able to condense into galaxies and form stars. We show that although both AGN and SNe suppress star formation, they mutually weaken one another's effect by up to an order of magnitude in haloes in the mass range for which both feedback processes are efficient (10^11.25 M_sun < m_200 < 10^12.5 M_sun). These results demonstrate the importance of the simultaneous, non-independent inclusion of these two processes in models of galaxy formation to estimate the total feedback strength. These results are of particular relevance to semi-analytic models, which implicitly assume the effects of the two feedback processes to be independent, and also to hydrodynamical simulations that model only one of the feedback processes.

preprint2012arXiv

Absorption signatures of warm-hot gas at low redshift: Broad HI Lyman-Alpha Absorbers

We investigate the physical state of HI absorbing gas at low redshift (z = 0.25) using a subset of cosmological, hydrodynamic simulations from the OWLS project, focusing in particular on broad (b_HI > 40 km/s) HI Lyman-Alpha absorbers (BLAs), which are believed to originate in shock-heated gas in the warm-hot intergalactic medium (WHIM). Our fiducial model, which includes radiative cooling by heavy elements and feedback by supernovae and active galactic nuclei, predicts that by z = 0.25 nearly 60 per cent of the gas mass ends up at densities and temperatures characteristic of the WHIM and we find that half of this fraction is due to outflows. The standard HI observables (distribution of HI column densities N_HI, distribution of Doppler parameters b_HI, b_HI - N_HI correlation) and the BLA line number density predicted by our simulations are in remarkably good agreement with observations. BLAs arise in gas that is hotter, more highly ionised and more enriched than the gas giving rise to typical Lyman-Alpha forest absorbers. The majority of the BLAs arise in warm-hot (log (T/K) ~ 5) gas at low (log Delta < 1.5) overdensities. On average, thermal broadening accounts for at least 60 per cent of the BLA line width, which in turn can be used as a rough indicator of the thermal state of the gas. Detectable BLAs account for only a small fraction of the true baryon content of the WHIM at low redshift. In order to detect the bulk of the mass in this gas phase, a sensitivity at least one order of magnitude better than achieved by current ultraviolet spectrographs is required. We argue that BLAs mostly trace gas that has been shock-heated and enriched by outflows and that they therefore provide an important window on a poorly understood feedback process.

preprint2012arXiv

Cold accretion flows and the nature of high column density HI absorption at redshift 3

Simulations predict that galaxies grow primarily through the accretion of gas that has not gone through an accretion shock near the virial radius and that this cold gas flows towards the central galaxy along dense filaments and streams. There is, however, little observational evidence for the existence of these cold flows. We use a large, cosmological, hydrodynamical simulation that has been post-processed with radiative transfer to study the contribution of cold flows to the observed z=3 column density distribution of neutral hydrogen, which our simulation reproduces. We find that nearly all of the HI absorption arises in gas that has remained colder than 10^5.5 K, at least while it was extragalactic. In addition, the majority of the HI is rapidly falling towards a nearby galaxy, with non-negligible contributions from outflowing and static gas. Above a column density of N_HI = 10^17 cm^-2, most of the absorbers reside inside haloes, but the interstellar medium only dominates for N_HI > 10^21 cm^-2. Haloes with total mass below 10^10 Msun dominate the absorption for 10^17<N_HI < 10^21 cm^-2, but the average halo mass increases sharply for higher column densities. Although very little of the HI in absorbers with N_HI <~ 10^20 cm^-2 resides inside galaxies, systems with N_HI > 10^17 cm^-2 are closely related to star formation: most of their HI either will become part of the interstellar medium before z=2 or has been ejected from a galaxy at z>3. Cold accretion flows are critical for the success of our simulation in reproducing the observed rate of incidence of damped Lyman-alpha and particularly that of Lyman limit systems. We therefore conclude that cold accretion flows exist and have already been detected in the form of high column density HI absorbers.

preprint2012arXiv

Hydrodynamical simulations and semi-analytic models of galaxy formation: two sides of the same coin

In this work we develop a new method to turn a state-of-the-art hydrodynamical cosmological simulation of galaxy formation (HYD) into a simple semi-analytic model (SAM). This is achieved by summarizing the efficiencies of accretion, cooling, star formation, and feedback given by the HYD, as functions of the halo mass and redshift. Surprisingly, by turning the HYD into a SAM, we conserve the mass of individual galaxies, with deviations at the level of 0.1 dex, on an object-by-object basis. This is true for all redshifts, and for the mass of stars and gas components, although the agreement reaches 0.2 dex for satellite galaxies at low redshift. We show that the same level of accuracy is obtained even in case the SAM uses only one phase of gas within each galaxy. Moreover, we demonstrate that the formation history of one massive galaxy provides sufficient information for the SAM to reproduce the population of galaxies within the entire cosmological box. The reasons for the small scatter between the HYD and SAM galaxies are: a) The efficiencies are matched as functions of the halo mass and redshift, meaning that the evolution within merger-trees agrees on average. b) For a given galaxy, efficiencies fluctuate around the mean value on time scales of 0.2-2 Gyr. c) The various mass components of galaxies are obtained by integrating the efficiencies over time, averaging out these fluctuations. We compare the efficiencies found here to standard SAM recipes and find that they often deviate significantly. For example, here the HYD shows smooth accretion that is less effective for low mass haloes, and is always composed of hot or dilute gas; cooling is less effective at high redshift; and star formation changes only mildly with cosmic time. The method developed here can be applied in general to any HYD, and can thus serve as a common language for both HYDs and SAM (Abridged).

preprint2012arXiv

Neutral hydrogen optical depth near star-forming galaxies at z~2.4 in the Keck Baryonic Structure Survey

[Abridged] We study the interface between galaxies and the intergalactic medium by measuring the absorption by neutral hydrogen in the vicinity of star-forming galaxies at z~2.4. Our sample consists of 679 rest-frame-UV selected galaxies with spectroscopic redshifts that have impact parameters < 2 (proper) Mpc to the line of sight of one of 15 bright, background QSOs. We present the first 2-D maps of the absorption around galaxies, plotting the median Ly-a pixel optical depth as a function of transverse and line of sight, LOS, separation from galaxies. The median optical depth, and hence the median density of atomic hydrogen, drops by more than an order of magnitude around 100 kpc, which is similar to the virial radius of the halos thought to host the galaxies. The median remains enhanced (>3 sigma level) out to at least 2.8 Mpc (i.e. > 9 comoving Mpc). Within 100 (200) kpc, and over \pm 165 km/s, the covering fraction of gas with Ly-a optical depth greater than unity is 100(+0)(-32)% (86(+14)(-18)%). Absorbers with tau(Ly-a)> 0.1 are typically closer to galaxies than random. The mean galaxy overdensity around absorbers increases with the optical depth and also as the length scale over which the galaxy overdensity is evaluated is decreased. Absorbers with tau(Ly-a)~1 reside in regions where the galaxy number density is close to the cosmic mean on scales > 0.25 Mpc. We detect two types of redshift space anisotropies. On scales < 200 km/s, or < 1 Mpc, the absorption is stronger along the LOS than in the transverse direction. This "finger of God" effect may be partly due to redshift errors, but is probably dominated by gas motions within or very close to the halos. On the other hand, on scales of 1.4 - 2.0 Mpc the absorption is compressed along the LOS (>3 sigma significance), which we attribute to large-scale infall (i.e. the Kaiser effect).

preprint2012arXiv

Properties of gas in and around galaxy haloes

We study the properties of gas inside and around galaxy haloes as a function of radius and halo mass, focussing mostly on z=2, but also showing some results for z=0. For this purpose, we use a suite of large cosmological, hydrodynamical simulations from the OverWhelmingly Large Simulations project. The properties of cold- and hot-mode gas, which we separate depending on whether the temperature has been higher than 10^5.5 K while it was extragalactic, are clearly distinguishable in the outer parts of massive haloes (virial temperatures >> 10^5 K. The differences between cold- and hot-mode gas resemble those between inflowing and outflowing gas. The cold-mode gas is mostly confined to clumpy filaments that are approximately in pressure equilibrium with the diffuse, hot-mode gas. Besides being colder and denser, cold-mode gas typically has a much lower metallicity and is much more likely to be infalling. However, the spread in the properties of the gas is large, even for a given mode and a fixed radius and halo mass, which makes it impossible to make strong statements about individual gas clouds. Metal-line cooling causes a strong cooling flow near the central galaxy, which makes it hard to distinguish gas accreted through the cold and hot modes in the inner halo. Stronger feedback results in larger outflow velocities and pushes hot-mode gas to larger radii. The gas properties evolve as expected from virial arguments, which can also account for the dependence of many gas properties on halo mass. We argue that cold streams penetrating hot haloes are observable as high-column density HI Lyman-alpha absorption systems in sightlines near massive foreground galaxies.

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

Sample variance and Lyman-alpha forest transmission statistics

We compare the observed probability distribution function of the transmission in the \HI\ Lyman-alpha forest, measured from the UVES 'Large Programme' sample at redshifts z=[2,2.5,3], to results from the GIMIC cosmological simulations. Our measured values for the mean transmission and its PDF are in good agreement with published results. Errors on statistics measured from high-resolution data are typically estimated using bootstrap or jack-knife resampling techniques after splitting the spectra into chunks. We demonstrate that these methods tend to underestimate the sample variance unless the chunk size is much larger than is commonly the case. We therefore estimate the sample variance from the simulations. We conclude that observed and simulated transmission statistics are in good agreement, in particular, we do not require the temperature-density relation to be 'inverted'.

preprint2012arXiv

Simulating galactic outflows with thermal supernova feedback

Cosmological simulations make use of sub-grid recipes for the implementation of galactic winds driven by massive stars because direct injection of supernova energy in thermal form leads to strong radiative losses, rendering the feedback inefficient. We argue that the main cause of the catastrophic cooling is a mismatch between the mass of the gas in which the energy is injected and the mass of the parent stellar population. Because too much mass is heated, the temperatures are too low and the cooling times too short. We use analytic arguments to estimate, as a function of the gas density and the numerical resolution, the minimum heating temperature that is required for the injected thermal energy to be efficiently converted into kinetic energy. We then propose and test a stochastic implementation of thermal feedback that uses this minimum temperature increase as an input parameter and that can be employed in both particle- and grid-based codes. We use smoothed particle hydrodynamics simulations to test the method on models of isolated disc galaxies in dark matter haloes with total mass 10^10 and 10^12 h^-1 solar masses. The thermal feedback strongly suppresses the star formation rate and can drive massive, large-scale outflows without the need to turn off radiative cooling temporarily. In accord with expectations derived from analytic arguments, for sufficiently high resolution the results become insensitive to the imposed temperature jump and also agree with high-resolution simulations employing kinetic feedback.

preprint2012arXiv

The impact of feedback from galaxy formation on the Lyman-alpha transmitted flux

The forest of Lyman-alpha absorption lines seen in the spectra of distant quasars has become an important probe of the distribution of matter in the Universe. We use large, hydrodynamical simulations from the OWLS project to investigate the effect of feedback from galaxy formation on the probability distribution function and the power spectrum of the Lyman-alpha transmitted flux. While metal-line cooling is unimportant, both galactic outflows from massive galaxies driven by active galactic nuclei and winds from low-mass galaxies driven by supernovae have a substantial impact on the flux statistics. At redshift z=2.25, the effects on the flux statistics are of a similar magnitude as the statistical uncertainties of published data sets. The changes in the flux statistics are not due to differences in the temperature-density relation of the photo-ionised gas. Instead, they are caused by changes in the density distribution and in the fraction of hot, collisionally ionised gas. It may be possible to disentangle astrophysical and cosmological effects by taking advantage of the fact that they induce different redshift dependencies. In particular, the magnitude of the feedback effects appears to decrease rapidly with increasing redshift. Analyses of Lyman-alpha forest data from surveys that are currently in process, such as BOSS/SDSS-III and X-Shooter/VLT, must take galactic winds into account.

preprint2012arXiv

The Origin of Disks and Spheroids in Simulated Galaxies

In the simplest scenario, disk galaxies form predominantly in halos with high angular momentum and quiet recent assembly history, whereas spheroids are the slowly-rotating remnants of repeated merging events. We explore these assumptions using one hundred systems with halo masses similar to that of the Milky Way, identified in a series of cosmological gasdynamical simulations GIMIC. At z=0, the simulated galaxies exhibit a wide variety of morphologies, from dispersion-dominated spheroids to pure disk galaxies. Surprisingly, these morphological features are very poorly correlated with their halo properties: disks form in halos with high and low net spin, and mergers play a negligible role in the formation of spheroid stars, most of which form in-situ. More important to morphology is the coherent alignment of the angular momentum of baryons that accrete over time to form a galaxy. Spheroids tend to form when the spin of newly-accreted gas is misaligned with that of the extant galaxy, leading to the episodic formation of stars with different kinematics that cancel out the net rotation of the system. Disks, on the other hand, form out of gas that flows in with similar angular momentum to that of earlier-accreted material. Gas accretion from a hot corona thus favours disk formation, whereas gas that flows "cold", often along separate, misaligned filaments, favours the formation of spheroids. In this scenario, most spheroids consist of superpositions of stellar components with distinct kinematics, age, and metallicity, an arrangement that might survive to the present day given the paucity of major mergers. Since angular momentum is acquired largely at turnaround, morphology is imprinted early by the interplay of the tidal field and the shape of the material destined to form the galaxy.

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

Calibrating Galaxy Redshifts Using Absorption by the Surrounding Intergalactic Medium

Rest-frame UV spectral lines of star-forming galaxies are systematically offset from the galaxies' systemic redshifts, probably because of large-scale outflows. We calibrate galaxy redshifts measured from rest-frame UV lines by utilizing the fact that the mean HI Ly-alpha absorption profiles around the galaxies, as seen in spectra of background objects, must be symmetric with respect to the true galaxy redshifts if the galaxies are oriented randomly with respect to the lines of sight to the background objects. We use 15 QSOs at z~2.5-3 and more than 600 foreground galaxies with spectroscopic redshifts at z~1.9-2.5. All galaxies are within 2 Mpc proper from the lines of sight to the background QSOs. We find that LyA emission and ISM absorption redshifts require systematic shifts of v_LyA=-295(+35)(-35) km/s and v_ISM=145(+70)(-35) km/s. Assuming a Gaussian distribution, we put 1-sigma upper limits on possible random redshift offsets of <220 km/s for LyA and <420 km/s for ISM redshifts. For the small subset (<10%) of galaxies for which near-IR spectra have been obtained, we can compare our results to direct measurements based on nebular emission lines which we confirm to mark the systemic redshifts. While our v_ISM agrees with the direct measurements, our v_LyA is significantly smaller. However, when we apply our method to the near-IR subsample which is characterized by slightly different selection effects, the best-fit velocity offset comes into agreement with the direct measurement. This confirms the validity of our approach, and implies that no single number appropriately describes the whole population of galaxies, in line with the observation that the line offset depends on galaxy spectral morphology. This method provides accurate redshift calibrations and will enable studies of circumgalactic matter around galaxies for which rest-frame optical observations are not available.

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

Dark matter halo concentrations in the Wilkinson Microwave Anisotropy Probe year 5 cosmology

We use a combination of three large N-body simulations to investigate the dependence of dark matter halo concentrations on halo mass and redshift in the WMAP year 5 cosmology. The median relation between concentration and mass is adequately described by a power-law for halo masses in the range 10^11 - 10^15 Msol/h and redshifts z < 2, regardless of whether the halo density profiles are fit using NFW or Einasto profiles. Compared with recent analyses of the Millennium Simulation, which uses a value of sigma_8 that is higher than allowed by WMAP5, z = 0 halo concentrations are reduced by factors ranging from 23 per cent at 10^11 Msol/h to 16 per cent at 10^14 Msol/h. The predicted concentrations are much lower than inferred from X-ray observations of groups and clusters.

preprint2011arXiv

Disentangling galaxy environment and host halo mass

[Abridged] The properties of observed galaxies and dark matter haloes in simulations depend on their environment. The term environment has been used to describe a wide variety of measures that may or may not correlate with each other. Popular measures of environment include the distance to the N'th nearest neighbour, the number density of objects within some distance, or the mass of the host dark matter halo. We use results from the Millennium simulation and a semi-analytic model for galaxy formation to quantify the relations between environment and halo mass. We show that the environmental parameters used in the observational literature are in effect measures of halo mass, even if they are measured for a fixed stellar mass. The strongest correlation between environment and halo mass arises when the number of objects is counted out to a distance of 1.5-2 times the virial radius of the host halo and when the galaxies/haloes are required to be relatively bright/massive. For observational studies the virial radius is not easily determined, but the number of neighbours out to 1-2 Mpc/h gives a similarly strong correlation. For the distance to the N'th nearest neighbour the correlation with halo mass is nearly as strong provided N>2. We demonstrate that this environmental parameter becomes insensitive to halo mass if it is constructed from dimensionless quantities. This can be achieved by scaling the minimum luminosity/mass of neighbours to that of the object in question and by dividing the distance to a length scale associated with either the neighbour or the galaxy under consideration. We show how such a halo mass independent environmental parameter can be defined for observational and numerical studies. The results presented here will help future studies to disentangle the effects of halo mass and external environment on the properties of galaxies and dark matter haloes.

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

Modelling neutral hydrogen in galaxies using cosmological hydrodynamical simulations

The characterisation of the atomic and molecular hydrogen content of high-redshift galaxies is a major observational challenge that will be addressed over the coming years with a new generation of radio telescopes. We investigate this important issue by considering the states of hydrogen across a range of structures within high-resolution cosmological hydrodynamical simulations. Additionally, our simulations allow us to investigate the sensitivity of our results to numerical resolution and to sub-grid baryonic physics (especially feedback from supernovae and active galactic nuclei). We find that the most significant uncertainty in modelling the neutral hydrogen distribution arises from our need to model a self-shielding correction in moderate density regions. Future simulations incorporating radiative transfer schemes will be vital to improve on our empirical self-shielding threshold. Irrespective of the exact nature of the threshold we find that while the atomic hydrogen mass function evolves only mildly from redshift two to zero, the molecular hydrogen mass function increases with increasing redshift, especially at the high-mass end. Interestingly, the weak evolution of the neutral hydrogen mass function is insensitive to the feedback scheme utilised, but the opposite is true for the molecular gas, which is more closely associated with the star formation in the simulations.

preprint2011arXiv

Rest-frame UV line emission from the intergalactic medium at 2<z<5

Rest-frame UV emission lines offer the possibility to directly image the gas around high-redshift galaxies with upcoming optical instruments. We use a suite of large, hydrodynamical simulations to predict the nature and detectability of emission lines from the intergalactic medium at 2<z<5. The brightest emission comes from HI Ly-alpha and the strongest metal line, CIII, is about an order of magnitude fainter, although HI Ly-alpha may be fainter if the gas is self-shielded to the UV background or if dust is important. The highest surface brightness regions for CIV, SiIII, SiIV and OVI are fainter than CIII by factors of a few. The NV and NeVIII lines, as well as HeII H-alpha, are substantially weaker but their maximum surface brightnesses still exceed 100 photon/cm^2/s/sr at z=2 (for 2" pixels). Lower ionisation lines arise in denser and colder gas that produces clumpier emission. The brightest HI Ly-alpha emission arises in highly overdense gas, but the highest surface brightness emission from high-ionisation metal lines traces a wider range of overdensities. Bright metal-line emission traces gas with temperatures close to the peak of the corresponding emissivity curve. While HI Ly-alpha, HeII H-alpha, CIII, SiIII, and SiIV are excellent probes of cold accretion flows and the colder parts of outflows, CIV, NV, OVI, and NeVIII are powerful tracers of the diffuse WHIM and galactic winds. A comparison of results from simulations with varying physical prescriptions demonstrates that the predictions for the brighter metal-line emission are robust to within factors of a few. Several emission lines from the high-redshift IGM will become detectable in the near future, possibly starting with the Cosmic Web Imager on Palomar. MUSE and the Keck Cosmic Web Imager have the potential to revolutionise studies of the interactions between high-redshift galaxies and their environment. (Abridged)

preprint2011arXiv

The correlation structure of dark matter halo properties

We investigate the correlation between nine different dark matter halo properties using a rank correlation analysis and a Principal Component Analysis for a sample of haloes spanning five orders of magnitude in mass. We consider mass and dimensionless measures of concentration, age, relaxedness, sphericity, triaxiality, substructure, spin, and environment, where the latter is defined in a way that makes it insensitive to mass. We find that concentration is the most fundamental property. Except for environment, all parameters are strongly correlated with concentration. Concentration, age, substructure, mass, sphericity and relaxedness can be considered a single family of parameters, albeit with substantial scatter. In contrast, spin, environment, and triaxiality are more independent, although spin does correlate strongly with substructure and both spin and triaxiality correlate substantially with concentration. Although mass sets the scale of a halo, all other properties are more sensitive to concentration.

preprint2011arXiv

The drop in the cosmic star formation rate below redshift 2 is caused by a change in the mode of gas accretion and by AGN feedback

The cosmic star formation rate is observed to drop sharply after redshift z=2. We use a large, cosmological, smoothed particle hydrodynamics simulation to investigate how this decline is related to the evolution of gas accretion and to outflows driven by active galactic nuclei (AGN). We find that the drop in the star formation rate follows a corresponding decline in the global cold-mode accretion rate density onto haloes, but with a delay of order the gas consumption time scale in the interstellar medium. Here we define cold-mode (hot-mode) accretion as gas that is accreted and whose temperature has never exceeded (did exceed) 10^5.5 K. In contrast to cold-mode accretion, which peaks at z~3, the hot mode continues to increase to z~1 and remains roughly constant thereafter. By the present time, the hot mode strongly dominates the global accretion rate onto haloes. Star formation does not track hot-mode halo accretion because most of the hot halo gas never accretes onto galaxies. AGN feedback plays a crucial role by preferentially preventing gas that entered haloes in the hot mode from accreting onto their central galaxies. Consequently, in the absence of AGN feedback, gas accreted in the hot mode would become the dominant source of fuel for star formation and the drop off in the cosmic star formation rate would be much less steep.

preprint2011arXiv

The effect of variations in the input physics on the cosmic distribution of metals predicted by simulations

[Abridged] We investigate how a range of physical processes affect the cosmic metal distribution using a suite of cosmological, hydrodynamical simulations. Focusing on z = 0 and 2, we study the metallicities and metal mass fractions for stars as well as for the ISM, and several more diffuse gas phases. We vary the cooling rates, star formation law, structure of the ISM, galactic winds, feedback from AGN, reionization history, stellar IMF, and cosmology. In all models stars and the warm-hot IGM (WHIM) constitute the dominant repository of metals, while for z > 2 the ISM is also important. In models with galactic winds, predictions for the metallicities of the various phases vary at the factor of two level and are broadly consistent with observations. The exception is the cold-warm IGM, whose metallicity varies at the order of magnitude level if the prescription for galactic winds is varied, even for a fixed wind energy per unit stellar mass formed, and falls far below the observed values if winds are not included. At the other extreme, the metallicity of the intracluster medium (ICM) is largely insensitive to the presence of galactic winds, indicating that its enrichment is regulated by other processes. The mean metallicities of stars (~ Z_sun), the ICM (~ 0.1 Z_sun), and the WHIM (~ 0.1 Z_sun) evolve only slowly, while those of the cold halo gas and the IGM increase by more than an order of magnitude from z = 5 to 0. Higher velocity outflows are more efficient at transporting metals to low densities, but actually predict lower metallicities for the cold-warm IGM since the winds shock-heat the gas to high temperatures, thereby increasing the fraction of the metals residing in, but not the metallicity of, the WHIM. Besides galactic winds driven by feedback from star formation, the metal distribution is most sensitive to the inclusion of metal-line cooling and feedback from AGN.

preprint2011arXiv

The effects of galaxy formation on the matter power spectrum: A challenge for precision cosmology

Upcoming weak lensing surveys, such as LSST, EUCLID, and WFIRST, aim to measure the matter power spectrum with unprecedented accuracy. In order to fully exploit these observations, models are needed that, given a set of cosmological parameters, can predict the non-linear matter power spectrum at the level of 1% or better for scales corresponding to comoving wave numbers 0.1<k<10 h/Mpc. We have employed the large suite of simulations from the OWLS project to investigate the effects of various baryonic processes on the matter power spectrum. In addition, we have examined the distribution of power over different mass components, the back-reaction of the baryons on the CDM, and the evolution of the dominant effects on the matter power spectrum. We find that single baryonic processes are capable of changing the power spectrum by up to several tens of per cent. Our simulation that includes AGN feedback, which we consider to be our most realistic simulation as, unlike those used in previous studies, it has been shown to solve the overcooling problem and to reproduce optical and X-ray observations of groups of galaxies, predicts a decrease in power relative to a dark matter only simulation ranging, at z=0, from 1% at k~0.3 h/Mpc to 10% at k~1 h/Mpc and to 30% at k~10 h/Mpc. This contradicts the naive view that baryons raise the power through cooling, which is the dominant effect only for k>70 h/Mpc. Therefore, baryons, and particularly AGN feedback, cannot be ignored in theoretical power spectra for k>0.3 h/Mpc. It will thus be necessary to improve our understanding of feedback processes in galaxy formation, or at least to constrain them through auxiliary observations, before we can fulfil the goals of upcoming weak lensing surveys.

preprint2011arXiv

The filling factor of intergalactic metals at redshift z=3

Observations of quasar absorption line systems reveal that the z=3 intergalactic medium (IGM) is polluted by heavy elements down to HI optical depths tau_HI<<10. What is not yet clear, however, is what fraction of the volume needs to be enriched by metals and whether it suffices to enrich only regions close to galaxies in order to reproduce the observations. We use gas density fields derived from large cosmological simulations, together with synthetic quasar spectra and imposed, model metal distributions to investigate what enrichment patterns can reproduce the observed median optical depth of CIV as a function of tau_HI. Our models can only satisfy the observational constraints if the z=3 IGM was primarily enriched by galaxies that reside in low-mass (m_tot<10^10 M_sun) haloes that can eject metals out to distances >10^2 kpc. Galaxies in more massive haloes cannot possibly account for the observations as they are too rare for their outflows to cover a sufficiently large fraction of the volume. Galaxies need to enrich gas out to distances that are much greater than the virial radii of their host haloes. Assuming the metals to be well mixed on small scales, our modeling requires that the fractions of the simulated volume and baryonic mass that are polluted with metals are, respectively, >10% and >50% in order to match observations.

preprint2011arXiv

The rates and modes of gas accretion on to galaxies and their gaseous haloes

(Abridged) We study the rate at which gas accretes onto galaxies and haloes and investigate whether the accreted gas was shocked to high temperatures before reaching a galaxy. For this purpose we use a suite of large cosmological, hydrodynamical simulations from the OWLS project. We improve on previous work by considering a wider range of halo masses and redshifts, by distinguishing accretion onto haloes and galaxies, by including important feedback processes, and by comparing simulations with different physics. The specific rate of gas accretion onto haloes is, like that for dark matter, only weakly dependent on halo mass. For halo masses Mhalo>>10^11 Msun it is relatively insensitive to feedback processes. In contrast, accretion rates onto galaxies are determined by radiative cooling and by outflows driven by supernovae and active galactic nuclei. Galactic winds increase the halo mass at which the central galaxies grow the fastest by about two orders of magnitude to Mhalo~10^12 Msun. Gas accretion is bimodal, with maximum past temperatures either of order the virial temperature or <~10^5 K. The fraction of gas accreted on to haloes in the hot mode is insensitive to feedback and metal-line cooling. It increases with decreasing redshift, but is mostly determined by halo mass, increasing gradually from less than 10% for ~10^11 Msun to greater than 90% at 10^13 Msun. In contrast, for accretion onto galaxies the cold mode is always significant and the relative contributions of the two accretion modes are more sensitive to feedback and metal-line cooling. The majority of stars present in any mass halo at any redshift were formed from gas accreted in the cold mode, although the hot mode contributes typically over 10% for Mhalo>~10^11 Msun. Galaxies, but not necessarily their gaseous haloes, are predominantly fed by gas that did not experience an accretion shock when it entered the host halo.

preprint2011arXiv

Through Thick and Thin - HI Absorption in Cosmological Simulations

We investigate the column density distribution function of neutral hydrogen at redshift z = 3 using a cosmological simulation of galaxy formation from the OverWhelmingly Large Simulations (OWLS) project. The base simulation includes gravity, hydrodynamics, star formation, supernovae feedback, stellar winds, chemodynamics, and element-by-element cooling in the presence of a uniform UV background. Self-shielding and formation of molecular hydrogen are treated in post-processing, without introducing any free parameters, using an accurate reverse ray-tracing algorithm and an empirical relation between gas pressure and molecular mass fraction. The simulation reproduces the observed z = 3 abundance of Ly-A forest, Lyman Limit and Damped Ly-A HI absorption systems probed by quasar sight lines over ten orders of magnitude in column density. Self-shielding flattens the column density distribution for NHI > 10^18 cm-2, while the conversion to fully neutral gas and conversion of HI to H2 steepen it around column densities of NHI = 10^20.3 cm-2 and NHI = 10^21.5 cm-2, respectively.

preprint2011arXiv

Towards an understanding of the evolution of the scaling relations for supermassive black holes

The growth of the supermassive black holes (BHs) that reside at the centres of most galaxies is intertwined with the physical processes that drive the formation of the galaxies themselves. The evolution of the relations between the mass of the BH, m_BH, and the properties of its host therefore represent crucial aspects of the galaxy formation process. We use a cosmological simulation, as well as an analytical model, to investigate how and why the scaling relations for BHs evolve with cosmic time. We find that a simulation that reproduces the observed redshift zero relations between m_BH and the properties of its host galaxy, as well as the thermodynamic profiles of the intragroup medium, also reproduces the observed evolution in the ratio m_BH/m_s for massive galaxies, although the evolution of the m_BH/sigma relation is in apparent conflict with observations. The simulation predicts that the relations between m_BH and the binding energies of both the galaxy and its dark matter halo do not evolve, while the ratio m_BH/m_halo increases with redshift. The simple, analytic model of Booth & Schaye (2010), in which the mass of the BH is controlled by the gravitational binding energy of its host halo, quantitatively reproduces the latter two results. Finally, we can explain the evolution in the relations between m_BH and the mass and binding energy of the stellar component of its host galaxy for massive galaxies (m_s~10^11 M_sun) at low redshift (z<1) if these galaxies grow primarily through dry mergers.

preprint2010arXiv

Absorption signatures of warm-hot gas at low redshift: OVI

We investigate the origin and physical properties of OVI absorbers at low redshift (z = 0.25) using a subset of cosmological, hydrodynamical simulations from the OverWhelmingly Large Simulations (OWLS) project. Intervening OVI absorbers are believed to trace shock-heated gas in the Warm-Hot Intergalactic Medium (WHIM) and may thus play a key role in the search for the missing baryons in the present-day Universe. When compared to observations, the predicted distributions of the different OVI line parameters (column density, Doppler parameter, rest equivalent width) from our simulations exhibit a lack of strong OVI absorbers. This suggests that physical processes on sub-grid scales (e.g. turbulence) may strongly influence the observed properties of OVI systems. We find that the intervening OVI absorption arises mainly in highly metal-enriched (0.1 << Z/Z_sun < 1) gas at typical overdensities of 1 << rho/<rho> < 100. One third of the OVI absorbers in our simulation are found to trace gas at temperatures T < 10^5 K, while the rest arises in gas at higher temperatures around T =10^5.3 K. The OVI resides in a similar region of (rho,T)-space as much of the shock-heated baryonic matter, but the vast majority of this gas has a lower metal content and does not give rise to detectable OVI absorption As a consequence of the patchy metal distribution, OVI absorbers in our simulations trace only a very small fraction of the cosmic baryons (<2 percent) and the cosmic metals. Instead, these systems presumably trace previously shock-heated, metal-rich material from galactic winds that is now cooling. The common approach of comparing OVI and HI column densities to estimate the physical conditions in intervening absorbers from QSO observations may be misleading, as most of the HI (and most of the gas mass) is not physically connected with the high-metallicity patches that give rise to the OVI absorption.

preprint2010arXiv

Dark matter haloes determine the masses of supermassive black holes

The energy and momentum deposited by the radiation from accretion onto the supermassive black holes (BHs) that reside at the centres of virtually all galaxies can halt or even reverse gas inflow, providing a natural mechanism for supermassive BHs to regulate their growth and to couple their properties to those of their host galaxies. However, it remains unclear whether this self-regulation occurs on the scale at which the BH is gravitationally dominant, on that of the stellar bulge, the galaxy, or that of the entire dark matter halo. To answer this question, we use self-consistent simulations of the co-evolution of the BH and galaxy populations that reproduce the observed correlations between the masses of the BHs and the properties of their host galaxies. We first confirm unambiguously that the BHs regulate their growth: the amount of energy that the BHs inject into their surroundings remains unchanged when the fraction of the accreted rest mass energy that is injected, is varied by four orders of magnitude. The BHs simply adjust their masses so as to inject the same amount of energy. We then use simulations with artificially reduced star formation rates to demonstrate explicitly that BH mass is not set by the stellar mass. Instead, we find that it is determined by the mass of the dark matter halo with a secondary dependence on the halo concentration, of the form that would be expected if the halo binding energy were the fundamental property that controls the mass of the BH. We predict that the logarithmic slope of the relation between dark matter halo mass and black hole mass is 1.55+/-0.05 and that the scatter around the mean relation in part reflects the scatter in the halo concentration-mass relation.

preprint2010arXiv

Feedback and the Structure of Simulated Galaxies at redshift z=2

We study the properties of simulated high-redshift galaxies using cosmological N-body/gasdynamical runs from the OverWhelmingly Large Simulations (OWLS) project. The runs contrast several feedback implementations of varying effectiveness: from no-feedback, to supernova-driven winds to powerful AGN-driven outflows. These different feedback models result in large variations in the abundance and structural properties of bright galaxies at z=2. We find that feedback affects the baryonic mass of a galaxy much more severely than its spin, which is on average roughly half that of its surrounding dark matter halo in our runs. Feedback induces strong correlations between angular momentum content and galaxy mass that leave their imprint on galaxy scaling relations and morphologies. Encouragingly, we find that galaxy disks are common in moderate-feedback runs, making up typically ~50% of all galaxies at the centers of haloes with virial mass exceeding 1e11 M_sun. The size, stellar masses, and circular speeds of simulated galaxies formed in such runs have properties that straddle those of large star-forming disks and of compact early-type galaxies at z=2. Once the detailed abundance and structural properties of these rare objects are well established it may be possible to use them to gauge the overall efficacy of feedback in the formation of high redshift galaxies.

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

Impact of baryon physics on dark matter structures: a detailed simulation study of halo density profiles

The back-reaction of baryons on the dark matter halo density profile is of great interest, not least because it is an important systematic uncertainty when attempting to detect the dark matter. Here, we draw on a large suite of high resolution cosmological hydrodynamical simulations, to systematically investigate this process and its dependence on the baryonic physics associated with galaxy formation. The inclusion of baryons results in significantly more concentrated density profiles if radiative cooling is efficient and feedback is weak. The dark matter halo concentration can in that case increase by as much as 30 (10) per cent on galaxy (cluster) scales. The most significant effects occur in galaxies at high redshift, where there is a strong anti-correlation between the baryon fraction in the halo centre and the inner slope of both the total and the dark matter density profiles. If feedback is weak, isothermal inner profiles form, in agreement with observations of massive, early-type galaxies. However, we find that AGN feedback, or extremely efficient feedback from massive stars, is necessary to match observed stellar fractions in groups and clusters, as well as to keep the maximum circular velocity similar to the virial velocity as observed for disk galaxies. These strong feedback models reduce the baryon fraction in galaxies by a factor of 3 relative to the case with no feedback. The AGN is even capable of reducing the baryon fraction by a factor of 2 in the inner region of group and cluster haloes. This in turn results in inner density profiles which are typically shallower than isothermal and the halo concentrations tend to be lower than in the absence of baryons.

preprint2010arXiv

Metal-line emission from the warm-hot intergalactic medium: I. Soft X-rays

Emission lines from metals offer one of the most promising ways to detect the elusive warm-hot intergalactic medium (WHIM; 10^5 K<T<10^7 K), which is thought to contain a substantial fraction of the baryons in the low-redshift Universe. We present predictions for the soft X-ray line emission from the WHIM using a subset of cosmological simulations from the OverWhelmingly Large Simulations (OWLS) project. We use the OWLS models to test the dependence of the predicted emission on a range of physical prescriptions, such as cosmology, gas cooling and feedback from star formation and accreting black holes. Provided that metal-line cooling is taken into account, the models give surprisingly similar results, indicating that the predictions are robust. Soft X-ray lines trace the hotter part of the WHIM (T>10^6 K). We find that the OVIII 18.97A is the strongest emission line, with a predicted maximum surface brightness of ~10^2 photon/s/cm^2/sr, but a number of other lines are only slightly weaker. All lines show a strong correlation between the intensity of the observed flux and the density and metallicity of the gas responsible for the emission. On the other hand, the potentially detectable emission consistently corresponds to the temperature at which the emissivity of the electronic transition peaks. The emission traces neither the baryonic nor the metal mass. In particular, the emission that is potentially detectable with proposed missions, traces overdense (rho>10^2rho_mean) and metal-rich (Z>0.1Z_sun) gas in and around galaxies and groups. While soft X-ray line emission is therefore not a promising route to close the baryon budget, it does offer the exciting possibility to image the gas accreting onto and flowing out of galaxies.

preprint2010arXiv

Metal-line emission from the warm-hot intergalactic medium: II. Ultraviolet

Approximately half the baryons in the local Universe are thought to reside in the warm-hot intergalactic medium (WHIM). Emission lines from metals in the UV band are excellent tracers of the cooler fraction of this gas. We present predictions for the surface brightness of a sample of UV lines that could potentially be observed by the next generation of UV telescopes at z<1. We use a subset of simulations from the OWLS project to create emission maps and to investigate the effect of varying the physical prescriptions for star formation, supernova and AGN feedback, chemodynamics and radiative cooling. Most models produce results in agreement within a factor of a few, indicating that the predictions are robust. Of the lines we consider, C III is the strongest line, but it typically traces gas colder than 10^5 K. The same is true for Si IV. The second strongest line, C IV, traces circum-galactic gas with T~10^5 K. O VI and Ne VIII probe the warmer (T~10^5.5 K and T~10^6 K, respectively) and more diffuse gas that may be a better tracer of the large scale structure. N V emission is intermediate between C IV and O VI. The intensity of all emission lines increases strongly with gas density and metallicity, and for the bright emission it is tightly correlated with the temperature for which the line emissivity is highest. In particular, the C III, C IV, Si IV and O VI emission that is sufficiently bright to be potentially detectable in the near future (>10^3 photon/s/cm^2/sr), comes from relatively dense (rho>10^2 rho_mean) and metal rich (Z>0.1 Z_sun) gas. As such, emission lines are highly biased tracers of the missing baryons and are not an optimal tool to close the baryon budget. However, they do provide a powerful means to detect the gas cooling onto or flowing out of galaxies and groups. (Abridged)

preprint2010arXiv

Multi-frequency, thermally coupled radiative transfer with TRAPHIC: Method and tests

We present an extension of TRAPHIC, the method for radiative transfer of ionising radiation in smoothed particle hydrodynamics simulations that we introduced in Pawlik & Schaye (2008). The new version keeps all advantages of the original implementation: photons are transported at the speed of light, in a photon-conserving manner, directly on the spatially adaptive, unstructured grid traced out by the particles, in a computation time that is independent of the number of radiation sources, and in parallel on distributed memory machines. We extend the method to include multiple frequencies, both hydrogen and helium, and to model the coupled evolution of the temperature and ionisation balance. We test our methods by performing a set of simulations of increasing complexity and including a small cosmological reionisation run. The results are in excellent agreement with exact solutions, where available, and also with results obtained with other codes if we make similar assumptions and account for differences in the atomic rates used. We use the new implementation to illustrate the differences between simulations that compute photoheating in the grey approximation and those that use multiple frequency bins. We show that close to ionising sources the grey approximation asymptotes to the multi-frequency result if photoheating rates are computed in the optically thin limit, but that the grey approximation breaks down everywhere if, as is often done, the optically thick limit is assumed.

preprint2010arXiv

Power spectrum extraction for redshifted 21-cm epoch of reionization experiments: the LOFAR case

One of the aims of the Low Frequency Array (LOFAR) Epoch of Reionization (EoR) project is to measure the power spectrum of variations in the intensity of redshifted 21-cm radiation from the EoR. The sensitivity with which this power spectrum can be estimated depends on the level of thermal noise and sample variance, and also on the systematic errors arising from the extraction process, in particular from the subtraction of foreground contamination. We model the extraction process using realistic simulations of the cosmological signal, the foregrounds and noise, and so estimate the sensitivity of the LOFAR EoR experiment to the redshifted 21-cm power spectrum. Detection of emission from the EoR should be possible within 360 hours of observation with a single station beam. Integrating for longer, and synthesizing multiple station beams within the primary (tile) beam, then enables us to extract progressively more accurate estimates of the power at a greater range of scales and redshifts. We discuss different observational strategies which compromise between depth of observation, sky coverage and frequency coverage. A plan in which lower frequencies receive a larger fraction of the time appears to be promising. We also study the nature of the bias which foreground fitting errors induce on the inferred power spectrum, and discuss how to reduce and correct for this bias. The angular and line-of-sight power spectra have different merits in this respect, and we suggest considering them separately in the analysis of LOFAR data.

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 enrichment history of cosmic metals

We use a suite of cosmological, hydrodynamical simulations to investigate the chemical enrichment history of the Universe. Specifically, we trace the origin of the metals back in time to investigate when various gas phases were enriched and by what halo masses. We find that the age of the metals decreases strongly with the density of the gas in which they end up. At least half of the metals that reside in the diffuse intergalactic medium (IGM) at redshift zero (two) were ejected from galaxies above redshift two (three). The mass of the haloes that last contained the metals increases rapidly with the gas density. More than half of the mass in intergalactic metals was ejected by haloes with total masses less than 1e11 solar masses and stellar masses less than 1e9 solar masses. The range of halo masses that contributes to the enrichment is wider for the hotter part of the IGM. By combining the `when' and `by what' aspects of the enrichment history, we show that metals residing in lower density gas were typically ejected earlier and by lower mass haloes.

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

Cosmological simulations of the growth of supermassive black holes and feedback from active galactic nuclei: method and tests

(Abridged) We present a method that tracks the growth of supermassive black holes (BHs) and the feedback from AGN in cosmological simulations. Our model is a substantially modified version of the one by Springel et al. (2005). Because cosmological simulations lack both the resolution and the physics to model the multiphase interstellar medium, they tend to strongly underestimate the Bondi-Hoyle accretion rate. To allow low-mass BHs to grow, it is therefore necessary to increase the predicted Bondi-Hoyle rates in dense gas by large, ad-hoc factors. We explore the physical regimes where the use of such factors is reasonable, and through this introduce a new prescription for gas accretion. Feedback from AGN is modeled by coupling a fraction of the rest-mass energy of the accreted gas thermally into the surrounding medium. We describe the implementation as well as the limitations of the model and motivate all the changes relative to previous work. We investigate the robustness of the predictions for the cosmic star formation history, the redshift zero cosmic BH density, BH scaling relations, and galaxy specific star formation rates. We find that the freedom introduced by the need to increase the predicted accretion rates, the standard procedure in the literature, is the most significant source of uncertainty. Our simulations demonstrate that supermassive BHs are able to regulate their growth by releasing a fixed amount of energy for a given halo mass, independent of the assumed efficiency of AGN feedback, which sets the normalization of the BH scaling relations. Regardless of whether BH seeds are initially placed above or below the BH scaling relations they grow onto the same relations. AGN feedback efficiently suppresses star formation in high-mass galaxies.

preprint2009arXiv

Galaxies-Intergalactic Medium Interaction Calculation --I. Galaxy formation as a function of large-scale environment

[Abridged] We present the first results of hydrodynamical simulations that follow the formation of galaxies to z=0 in spherical regions of radius ~20 Mpc/h drawn from the Millennium Simulation. The regions have overdensities that deviate by (-2, -1, 0, +1, +2)sigma from the cosmic mean, where sigma is the rms mass fluctuation on a scale of ~20Mpc/h at z=1.5. The simulations have mass resolution of up to 10^6 Msun/h, cover the entire range of large-scale environments and allow extrapolation of statistics to the entire 500 (Mpc/h)^3 Millennium volume. They include gas cooling, photoheating from an ionising background, SNe feedback and winds, but no AGN. We find that the specific SFR density at z <~ 10 varies systematically from region to region by up to an order of magnitude, but the global value, averaged over all volumes, reproduces observational data. Massive, compact galaxies, similar to those observed in the GOODS fields, form in the overdense regions as early as z=6, but do not appear in the underdense regions until z~3. These environmental variations are not caused by a dependence of the star formation properties on environment, but rather by a strong variation of the halo mass function from one environment to another, with more massive haloes forming preferentially in the denser regions. At all epochs, stars form most efficiently in haloes of circular velocity ~ 250 km/s. However, the star formation history exhibits a form of "downsizing" (even in the absence of AGN): the stars comprising massive galaxies at z=0 have mostly formed by z=1-2, whilst those comprising smaller galaxies typically form at later times. However, additional feedback is required to limit star formation in massive galaxies at late times.

preprint2009arXiv

The physics driving the cosmic star formation history

We investigate the physics driving the cosmic star formation (SF) history using the more than fifty large, cosmological, hydrodynamical simulations that together comprise the OverWhelmingly Large Simulations (OWLS) project. We systematically vary the parameters of the model to determine which physical processes are dominant and which aspects of the model are robust. Generically, we find that SF is limited by the build-up of dark matter haloes at high redshift, reaches a broad maximum at intermediate redshift, then decreases as it is quenched by lower cooling rates in hotter and lower density gas, gas exhaustion, and self-regulated feedback from stars and black holes. The higher redshift SF is therefore mostly determined by the cosmological parameters and to a lesser extent by photo-heating from reionization. The location and height of the peak in the SF history, and the steepness of the decline towards the present, depend on the physics and implementation of stellar and black hole feedback. Mass loss from intermediate-mass stars and metal-line cooling both boost the SF rate at late times. Galaxies form stars in a self-regulated fashion at a rate controlled by the balance between, on the one hand, feedback from massive stars and black holes and, on the other hand, gas cooling and accretion. Paradoxically, the SF rate is highly insensitive to the assumed SF law. This can be understood in terms of self-regulation: if the SF efficiency is changed, then galaxies adjust their gas fractions so as to achieve the same rate of production of massive stars. Self-regulated feedback from accreting black holes is required to match the steep decline in the observed SF rate below redshift two, although more extreme feedback from SF, for example in the form of a top-heavy IMF at high gas pressures, can help.

preprint2008arXiv

Simulating galactic outflows with kinetic supernova feedback

Feedback from star formation is thought to play a key role in the formation and evolution of galaxies, but its implementation in cosmological simulations is currently hampered by a lack of numerical resolution. We present and test a sub-grid recipe to model feedback from massive stars in cosmological smoothed particle hydrodynamics simulations. The energy is distributed in kinetic form among the gas particles surrounding recently formed stars. The impact of the feedback is studied using a suite of high-resolution simulations of isolated disc galaxies embedded in dark halos with total mass 10^{10} and 10^{12} Msol/h. We focus in particular on the effect of pressure forces on wind particles within the disc, which we turn off temporarily in some of our runs to mimic a recipe that has been widely used in the literature. We find that this popular recipe gives dramatically different results because (ram) pressure forces on expanding superbubbles determine both the structure of the disc and the development of large-scale outflows. Pressure forces exerted by expanding superbubbles puff up the disc, giving the dwarf galaxy an irregular morphology and creating a galactic fountain in the massive galaxy. Hydrodynamic drag within the disc results in a strong increase of the effective mass loading of the wind for the dwarf galaxy, but quenches much of the outflow in the case of the high-mass galaxy.

preprint2007arXiv

On the relation between the Schmidt and Kennicutt-Schmidt star formation laws and its implications for numerical simulations

When averaged over large scales, star formation in galaxies is observed to follow the empirical Kennicutt-Schmidt (KS) law for surface densities above a constant threshold. While the observed law involves surface densities, theoretical models and simulations generally work with volume density laws (i.e. Schmidt laws). We derive analytic relations between star formation laws expressed in terms of surface densities, volume densities, and pressures and we show how these relations depend on parameters such as the effective equation of state of the multiphase interstellar medium. Our analytic relations enable us to implement observed surface density laws into simulations. Because the parameters of our prescription for star formation are observables, we are not free to tune them to match the observations. We test our theoretical framework using high-resolution simulations of isolated disc galaxies that assume an effective equation of state for the multiphase interstellar medium. We are able to reproduce the star formation threshold and both the slope and the normalisation of arbitrary input KS laws without tuning any parameters and with very little scatter, even for unstable galaxies and even if we use poor numerical resolution. Moreover, we can do so for arbitrary effective equations of state. Our prescription therefore enables simulations of galaxies to bypass our current inability to simulate the formation of stars. On the other hand, the fact that we can reproduce arbitrary input thresholds and KS laws, rather than just the particular ones picked out by nature, indicates that simulations that lack the physics and/or resolution to simulate the multiphase interstellar medium can only provide limited insight into the origin of the observed star formation laws.

preprint2002arXiv

Metallicity of the Intergalactic Medium Using Pixel Statistics: I. Method

Studies of absorption spectra of high-z QSOs have revealed that the intergalactic medium at z ~ 2-3 is enriched to ~ 0.1%-1% solar for gas densities more than a few times the mean cosmic density, but have not yet produced an accurate metallicity estimate, nor constrained variations in the metallicity with density, redshift, or spatial location. This paper discusses the ``pixel optical depth'' (POD) method of QSO spectrum analysis, using realistic simulated spectra from cosmological simulations. In this method, absorption in Ly-alpha is compared to corresponding metal absorption on a pixel-by-pixel basis, yielding for each analyzed spectrum a single statistical correlation encoding metal enrichment information. Our simulations allow testing and optimization of each step of the technique's implementation. Tests show that previous studies have probably been limited by C IV self-contamination and O VI contamination by HI lines; we have developed and tested an effective method of correcting for both contaminants. We summarize these and other findings, and provide a useful recipe for the POD technique's application to observed spectra. Our tests reveal that the POD technique applied to spectra of presently available quality is effective in recovering useful metallicity information even in underdense gas. We present an extension of the POD technique to directly recover the intergalactic metallicity as a function of gas density. For a given ionizing background, both the oxygen and carbon abundance can be measured with errors of at most a factor of a few over at least an order of magnitude in density, using a single high-quality spectrum.

preprint2001arXiv

Observational signatures of feedback in QSO absorption spectra

Models for the formation of galaxies and clusters of galaxies require strong feedback in order to explain the observed properties of these systems. We investigate whether such feedback has observational consequences for the intergalactic medium, as probed in absorption towards background quasars. A typical quasar sight-line intersects one proto-cluster per unit redshift, and significant feedback from forming galaxies or AGN, heating the proto-cluster gas, will result in a large clearing of reduced absorption in the Lyman alpha-forest. Such a gap could be detected at redshift > 3 when the mean opacity is high. Feedback from Lyman-break galaxies in proto-clusters can be probed by the absorption lines produced in their winds. Strong feedback from galaxies has a major impact on the number and properties of absorption lines with neutral hydrogen column densities of 10^16/ cm^2. This feedback can be probed with multiple sight-lines and by studying the unsaturated higher-order lines of the Lyman series. Galactic winds from dwarf galaxies should break-up into clouds, in order not to over produce the number of absorption lines. These clouds can then coast to large distances.

Joop Schaye

What is connected

Connect this record

See the researcher in context

Building this map preview

183 published item(s)

Cosmological back-reaction of baryons on dark matter in the CAMELS simulations

Luminosity-Dependent Assembly Bias of Central Galaxies from Weak Lensing and Clustering

Modelling the evolution and influence of dust in cosmological simulations that include the cold phase of the interstellar medium

MUSEQuBES: Investigating the Physical and Chemical Properties of Circumgalactic Gas Around Ly$α$ Emitters at $z \approx 3.3$

Galaxy cluster aperture masses are more robust to baryonic effects than 3D halo masses

Galaxy cluster photons alter the ionisation state of the nearby warm-hot intergalactic medium

How baryons can significantly bias cluster count cosmology

Spin-driven jet feedback in idealised simulations of galaxy groups and clusters

The effects of local stellar radiation and dust depletion on non-equilibrium interstellar chemistry

The importance of black hole repositioning for galaxy formation simulations

The importance of the way in which supernova energy is distributed around young stellar populations in simulations of galaxies

The interplay between AGN feedback and precipitation of the intracluster medium in simulations of galaxy groups and clusters

The MUSE Extremely Deep Field: a first panoramic view of an Mg II emitting intragroup medium

The MUSE eXtremely Deep Field: Individual detections of Lyα haloes around rest-frame UV-selected galaxies at z~2.9-4.4

The warm-hot circumgalactic medium around EAGLE-simulation galaxies and its detection prospects with X-ray line emission

Why are we still using 3D masses for cluster cosmology?

Baryonic mass budgets for haloes in the EAGLE simulation, including ejected and prevented gas

How gas flows shape the stellar-halo mass relation in the EAGLE simulation

Morphological and Rotation Structures of Circumgalactic Mg II Gas in the EAGLE Simulation and the Dependence on Galaxy Properties

SEAGLE--II: Constraints on feedback models in galaxy formation from massive early type strong lens galaxies

The MUSE Extremely Deep Field: Evidence for SFR-induced cores in dark-matter dominated galaxies at z=1

The surprising accuracy of isothermal Jeans modelling of self-interacting dark matter density profiles

An EAGLE view of the missing baryons

EAGLE and Illustris-TNG predictions for resolved eROSITA X-ray observations of the circumgalactic medium around normal galaxies

Elevated ionizing photon production efficiency in faint high-equivalent-width Lyman-alpha emitters

Galactic inflow and wind recycling rates in the EAGLE simulations

Galactic outflow rates in the EAGLE simulations

Numerical convergence of hydrodynamical simulations of galaxy formation: the abundance and internal structure of galaxies and their cold dark matter haloes

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

Radiative cooling rates, ion fractions, molecule abundances and line emissivities including self-shielding and both local and metagalactic radiation fields

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 effect of gas accretion on the radial gas metallicity profile of simulated galaxies

The MUSE Hubble Ultra Deep Field Survey XIV. The evolution of the Lya emitter fraction from z=3 to z=6

The MUSE Hubble Ultra Deep Field Survey. XV. The mean rest-UV spectra of Ly-alpha emitters at z>3

The warm-hot circumgalactic medium around EAGLE-simulation galaxies and its detection prospects with X-ray and UV line absorption

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

Feedback from supermassive black holes transforms centrals into passive galaxies by ejecting circumgalactic gas

Hot WHIM counterparts of FUV OVI absorbers: Evidence in the line-of-sight towards quasar 3C 273

Measuring the temperature and profiles of Lyman-$α$ absorbers

MusE GAs FLOw and Wind (MEGAFLOW) IV: A two sightline tomography of a galactic wind

The impact of the observed baryon distribution in haloes on the total matter power spectrum

A chronicle of galaxy mass assembly in the EAGLE simulation

A large difference in the progenitor masses of active and passive galaxies in the EAGLE simulation

Alignments between galaxies, satellite systems and haloes

Being WISE II: Reducing the Influence of Star Formation History on the Mass-to-Light Ratio of Quiescent Galaxies

Bimodality of low-redshift circumgalactic O VI in non-equilibrium EAGLE zoom simulations

Chemical evolution of giant molecular clouds in simulations of galaxies

Cosmic distribution of highly ionized metals and their physical conditions in the EAGLE simulations

Far-infrared and dust properties of present-day galaxies in the EAGLE simulations

Galaxies in the EAGLE hydrodynamical simulation and in the Durham and Munich semi-analytical models

Implications of a variable IMF for the interpretation of observations of galaxy populations

It's not easy being green: The evolution of galaxy colour in the EAGLE simulation

MusE GAs FLOw and Wind (MEGAFLOW) I: First MUSE results on background quasars

Music from the heavens - Gravitational waves from supermassive black hole mergers in the EAGLE simulations

Observations of metals in the $z\approx3.5$ intergalactic medium and comparison to the EAGLE simulations

On the connection between the metal-enriched intergalactic medium and galaxies: an OVI-galaxy cross-correlation study at $z < 1$

Simulated Milky Way analogues: implications for dark matter direct searches

Size matters: abundance matching, galaxy sizes, and the Tully-Fisher relation in EAGLE

Subhalo abundance matching and assembly bias in the EAGLE simulation

Supermassive black holes in the EAGLE Universe. Revealing the observables of their growth

The Athena X-ray Integral Field Unit (X-IFU)

The brighter galaxies reionised the Universe

The dark nemesis of galaxy formation: why hot haloes trigger black hole growth and bring star formation to an end

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

The EAGLE simulations: atomic hydrogen associated with galaxies

The effect of baryons on redshift space distortions and cosmic density and velocity fields in the EAGLE simulation

The effects of metallicity, UV radiation and non-equilibrium chemistry in high-resolution simulations of galaxies

The environmental dependence of HI in galaxies in the EAGLE simulations

The fundamental plane of star formation in galaxies revealed by the EAGLE hydrodynamical simulations

The low-mass end of the baryonic Tully-Fisher relation

The oldest and most metal poor stars in the APOSTLE Local Group simulations

The origin of compact galaxies with anomalously high black hole masses

The origin of the $α$-enhancement of massive galaxies