BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Joop Schaye

Joop Schaye contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.GAastro-ph.COastro-ph.HE

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
42works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

42 published item(s)

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 &#34;back-reaction&#34; 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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 &#34;repositioning&#34; 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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&#39;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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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 &#34;ejected&#34; and &#34;prevented&#34; 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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&#39;). 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.

0 citations0 reviewsNo code linkedOpen access
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&#39;, 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&#39; models, while enhancing it for `thawing&#39; 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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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 &#34;green-valley&#34; 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.

0 citations0 reviewsNo code linkedOpen access
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}$.

0 citations0 reviewsNo code linkedOpen access