BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Ewald Puchwein

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

ResearcherUnclaimed source record
astro-ph.COastro-ph.GAastro-ph.HEgr-qcastro-ph.IMastro-ph.SRhep-phphysics.comp-phphysics.plasm-ph

Catalog footprint

What is connected

39works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

39 published item(s)

preprint2026arXiv

Characterizing turbulence in galaxy clusters: defining turbulent energies and assessing multi-scale versus fixed-scale filters

Disentangling turbulence and bulk motions in the intracluster medium (ICM) of galaxy clusters is inherently ambiguous, as the plasma is continuously stirred by different processes on disparate scales. This poses a serious problem in the interpretation of both observations and numerical simulations. In this paper, we use filtering operators in real space to separate bulk motion from turbulence at different scales. We show how filters can be used to define consistent kinetic and magnetic energies for the bulk and turbulent component. We apply our GPU-accelerated filtering pipeline to a simulation of a major galaxy cluster merger, which is part of the PICO-Clusters suite of zoom-in cosmological simulations of massive clusters using the moving mesh code Arepo and the IllustrisTNG galaxy formation model. We find that during the merger the turbulent pressure fraction on physical scales $\lesssim$160 kpc reaches a maximum of 5%, before decreasing to 2% after $\sim$1.3 Gyr from the core passage. These low values are consistent with recent observations of clusters with XRISM, and suggest that unless a cluster was recently perturbed by a major merger, turbulence levels are low. We then re-examine the popular multiscale iterative filter method. In our tests, we find that its use can introduce artifacts, and that it does not reliably disentangle fluctuations living on widely separated length scales. Rather, we believe it is more fruitful to use fixed-scale filters and turbulent energies to compare between simulations and observations. This work significantly improves our understanding of turbulence generation by major mergers in galaxy clusters, which can be probed by XRISM and next-generation X-ray telescopes, allowing us to connect high-resolution cosmological simulations to observations.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

The THESAN-ZOOM project: Star-formation efficiencies in high-redshift galaxies

Recent JWST observations hint at unexpectedly intense cosmic star-formation in the early Universe, often attributed to enhanced star-formation efficiencies (SFEs). Here, we analyze the SFE in THESAN-ZOOM, a novel zoom-in radiation-hydrodynamic simulation campaign of high-redshift ($z \gtrsim 3$) galaxies employing a state-of-the-art galaxy formation model resolving the multiphase interstellar medium (ISM). The halo-scale SFE ($ε^{\ast}_{\rm halo}$) - the fraction of baryons accreted by a halo that are converted to stars - follows a double power-law dependence on halo mass, with a mild redshift evolution above $M_{\rm halo} \gtrsim 10^{9.5}\,{\rm M}_{\odot}$. The power-law slope is roughly $1/3$ at large halo masses, consistent with expectations when gas outflows are momentum-driven. At lower masses, the slope is roughly $2/3$ and is more aligned with the energy-driven outflow scenario. $ε^{\ast}_{\rm halo}$ is a factor of $2-3$ larger than commonly assumed in empirical galaxy-formation models at $M_{\rm halo} \lesssim 10^{11}\,{\rm M}_{\odot}$. On galactic (kpc) scales, the Kennicutt-Schmidt (KS) relation of neutral gas is universal in THESAN-ZOOM, following $Σ_{\rm SFR} \propto Σ_{\rm gas}^2$, indicative of a turbulent energy balance in the ISM maintained by stellar feedback. The rise of $ε^{\ast}_{\rm halo}$ with halo mass can be traced primarily to increasing gas surface densities in massive galaxies, while the underlying KS relation and neutral, star-forming gas fraction remain unchanged. Although the increase in $ε^{\ast}_{\rm halo}$ with redshift is relatively modest, it is sufficient to explain the large observed number density of UV-bright galaxies at $z \gtrsim 12$. However, reproducing the brightest sources at $M_{\rm UV} \lesssim -21$ may require extrapolating the SFE beyond the halo mass range directly covered by THESAN-ZOOM.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Limits on non-canonical heating and turbulence in the intergalactic medium from the low redshift Lyman-alpha forest

We examine the column density distribution function (CDDF) and Doppler parameter distribution from hydrodynamical simulations and Cosmic Origins Spectrograph (COS) observations of the Lyman-alpha forest at redshift $0\leq z\leq 0.2$. Allowing for a factor of two uncertainty in the metagalactic HI photoionisation rate, our hydrodynamical simulations are in good agreement ($1$-$1.5σ$) with the shape and amplitude of the observed CDDF at HI column densities $10^{13.3}\rm\,cm^{-2}\leq N_{\rm HI}\leq 10^{14.5}\rm\,cm^{-2}$. However, the Doppler widths of the simulated lines remain too narrow with respect to the COS data. We argue that invoking AGN feedback does not resolve this discrepancy. We also disfavour enhanced photoheating rates as a potential solution, as this requires an unphysically hard UV background spectrum. If instead appealing to a non-canonical source of heating, an additional specific heat injection of $u \lesssim 6.9\rm\,eV\,m_{\rm p}^{-1}$ is required at $z\lesssim 2.5$ for gas that has $N_{\rm HI}\simeq 10^{13.5}\rm\,cm^{-2}$ by $z=0.1$. Alternatively, there may be an unresolved line of sight turbulent velocity component of $v_{\rm turb}\lesssim 8.5\rm\,km\,s^{-1}(N_{\rm HI}/10^{13.5}\rm\,cm^{-2})^{0.21}$ for the coldest gas in the diffuse IGM.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The effect of inhomogeneous reionisation on the Lyman-$α$ forest power spectrum at redshift $z>4$: implications for thermal parameter recovery

We use the Sherwood-Relics suite of hybrid hydrodynamical and radiative transfer simulations to model the effect of inhomogeneous reionisation on the 1D power spectrum of the \Lya forest transmitted flux at redshifts $4.2\leq z \leq 5$. Relative to models that assume a homogeneous UV background, reionisation suppresses the power spectrum at small scales, $k \sim 0.1\rm\,km^{-1}\,s$ by $\sim 10$ per cent because of spatial variations in the thermal broadening kernel and the divergent peculiar velocity field associated with over-pressurised intergalactic gas. On larger scales, $k<0.03\rm\,km^{-1}\,s$, the power spectrum is instead enhanced by $10$--$50$ per cent by large scale spatial variations in the neutral hydrogen fraction. The effect of inhomogeneous reionisation must therefore be accounted for in analyses of forthcoming high precision measurements. We provide a correction for the Ly$α$ forest power spectrum at $4.1\leq z \leq 5.4$ that can be easily applied within other parameter inference frameworks using similar reionisation models. We perform a Bayesian analysis of mock data to assess the extent of systematic biases that may arise in measurements of the intergalactic medium if ignoring this correction. At the scales probed by current high resolution Ly$α$ forest data at $z>4$, $0.006 \rm \,km^{-1}\,s\leq k \leq 0.2 \rm\, km^{-1}\,s$, we find inhomogeneous reionisation does not introduce any significant bias in thermal parameter recovery for the current measurement uncertainties of $\sim 10$ per cent. However, for $5$ per cent uncertainties, $\sim 1σ$ shifts between the estimated and true parameters occur.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Concerns about Modelling of the EDGES Data

It is predicted that the spectrum of radio emission from the whole sky should show a dip arising from the action of the light from the first stars on the hydrogen atoms in the surrounding gas, which causes the 21-cm line to appear in absorption against the cosmic microwave background. Bowman et al. 2018 identified a broad flat-bottomed absorption profile centred at 78 MHz, which could be this feature, although the depth of the profile is much larger than expected. We have examined the modelling process they used and find that their data implies unphysical parameters for the foreground emission and also that their solution is not unique in the sense that we found other simple formulations for the signal that are different in shape but that also fit their data. We argue that this calls into question the interpretation of these data as an unambiguous detection of the cosmological 21-cm absorption signature.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Metal enrichment in the circumgalactic medium and Lyα haloes around quasars at z $\sim$3

Deep observations have detected extended Ly$α$ emission nebulae surrounding tens of quasars at redshift 2 to 6. However, the metallicity of such extended haloes is still poorly understood. We perform a detailed analysis on a large sample of 80 quasars at $z\sim3$ based on MUSE-VLT data. We find clear evidence of extended emission of the UV nebular lines such as CIV or HeII for about 20$\%$ of the sample, while CIII] is only marginally detected in a few objects. By stacking the cubes we detect emission of CIV, HeII and CIII] out to a radius of about 45 kpc. CIV and HeII show a radial decline much steeper than Ly$α$, while CIII] shows a shallower profile similar to Ly$α$ in the inner 45 kpc. We infer that the average metallicity of the circumgalactic gas within the central 30-50~kpc is $\sim$0.5 solar, or even higher. However, we also find evidence of a component of the Ly$α$ haloes, which has much weaker metal emission lines relative to Ly$α$. We suggest that the high metallicity of the circumgalactic medium within the central 30-50 kpc is associated with chemical pre-enrichment by past quasar-driven outflows and that there is a more extended component of the CGM that has much lower metallicity and likely associated with near-pristine gas accreted from the intergalactic medium. We show that our observational results are in good agreement with the expectations of the FABLE zoom-in cosmological simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Probing the thermal state of the intergalactic medium at $z>5$ with the transmission spikes in high-resolution Ly$α$ forest spectra

We compare a sample of five high-resolution, high S/N Ly$α$ forest spectra of bright $6<z \lesssim 6.5$ QSOs aimed at spectrally resolving the last remaining transmission spikes at $z>5$ with those obtained from mock absorption spectra from the Sherwood and Sherwood-Relics suites of hydrodynamical simulations of the intergalactic medium (IGM). We use a profile fitting procedure for the inverted transmitted flux, $1-F$, similar to the widely used Voigt profile fitting of the transmitted flux $F$ at lower redshifts, to characterise the transmission spikes that probe predominately underdense regions of the IGM. We are able to reproduce the width and height distributions of the transmission spikes, both with optically thin simulations of the post-reionization Universe using a homogeneous UV background and full radiative transfer simulations of a late reionization model. We find that the width of the fitted components of the simulated transmission spikes is very sensitive to the instantaneous temperature of the reionized IGM. The internal structures of the spikes are more prominant in low temeperature models of the IGM. The width distribution of the observed transmission spikes, which require high spectral resolution ($\leq $ 8 km/s) to be resolved, is reproduced for optically thin simulations with a temperature at mean density of $T_0= (11000 \pm 1600,10500\pm 2100,12000 \pm 2200)$ K at $z= (5.4,5.6,5.8)$. This is weakly dependent on the slope of the temperature-density relation, which is favoured to be moderately steeper than isothermal. In the inhomogeneous, late reionization, full radiative transfer simulations where islands of neutral hydrogen persist to $z\sim5.3$, the width distribution of the observed transmission spikes is consistent with the range of $T_0$ caused by spatial fluctuations in the temperature-density relation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The growth of the longitudinal beam-plasma instability in the presence of an inhomogeneous background

We study the longitudinal stability of beam-plasma systems in the presence of a density inhomogeneity in the background plasma. Previous works have focused on the non-relativistic regime where hydrodynamical models are used to evolve pre-existing Langmuir waves within inhomogeneous background plasmas. Here, for the first time we study the problem with kinetic equations in a fully-relativistic way. We do not assume the existence of Langmuir waves, and we focus on the rate and the mechanism by which waves are excited in such systems from an initial perturbation. We derive the structure of the unstable modes and compute an analytical approximation for their growth rates. Our computation is limited to dilute and cold beams, and shows an excellent agreement with particle-in-cell simulations performed using the SHARP code. We show that, due to such an inhomogeneity, the virulent beam-plasma instabilities in the intergalactic medium are not suppressed but their counterparts in the solar wind can be suppressed as evidenced by propagating type-III solar radio bursts.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Bow Ties in the Sky I: The Angular Structure of Inverse Compton Gamma-ray Halos in the Fermi Sky

Extended inverse Compton halos are generally anticipated around extragalactic sources of gamma rays with energies above 100 GeV. These result from inverse Compton scattered cosmic microwave background photons by a population of high-energy electron/positron pairs produced by the annihilation of the high-energy gamma rays on the infrared background. Despite the observed attenuation of the high-energy gamma rays, the halo emission has yet to be directly detected. Here, we demonstrate that in most cases these halos are expected to be highly anisotropic, distributing the up-scattered gamma rays along axes defined either by the radio jets of the sources or oriented perpendicular to a global magnetic field. We present a pedagogical derivation of the angular structure in the inverse Compton halo and provide an analytic formalism that facilitates the generation of mock images. We discuss exploiting this fact for the purpose of detecting gamma-ray halos in a set of companion papers.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Models of the Cosmological 21 cm Signal from the Epoch of Reionization Calibrated with Lyman-alpha and CMB Data

We present here 21 cm predictions from high dynamic range simulations for a range of reionization histories that have been tested against available Lyman-alpha and CMB data. We assess the observability of the predicted spatial 21 cm fluctuations by ongoing and upcoming experiments in the late stages of reionization in the limit in which the hydrogen spin temperature is significantly larger than the CMB temperature. Models consistent with the available Lyman-alpha data and CMB measurement of the Thomson optical depth predict typical values of 10--20 mK^2 for the variance of the 21 cm brightness temperature at redshifts z=7--10 at scales accessible to ongoing and upcoming experiments (k < 1 h/cMpc). This is only a factor of a few below the sensitivity claimed to have been already reached by ongoing experiments in the signal rms value. Our different models for the reionization history make markedly different predictions for the redshift evolution and thus frequency dependence of the 21 cm power spectrum and should be easily discernible by LOFAR (and later HERA and SKA1) at their design sensitivity. Our simulations have sufficient resolution to assess the effect of high-density Lyman limit systems that can self-shield against ionizing radiation and stay 21 cm bright even if the hydrogen in their surroundings is highly ionized. Our simulations predict that including the effect of the self-shielded gas in highly ionized regions reduces the large scale 21 cm power by about 30%.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Probing the end of reionization with the near-zones of z > 6 QSOs

QSO near-zones are an important probe of the the ionization state of the IGM at z ~ 6-7, at the end of reionization. We present here high-resolution cosmological 3D radiative transfer simulations of QSO environments for a wide range of host halo masses, 10$^{10-12.5}$ M_sun. Our simulated near-zones reproduce both the overall decrease of observed near-zone sizes at 6 < z < 7 and their scatter. The observable near-zone properties in our simulations depend only very weakly on the mass of the host halo. The size of the H II region expanding into the IGM is generally limited by (super-)Lyman Limit systems loosely associated with (low-mass) dark matter haloes. This leads to a strong dependence of near-zone size on direction and drives the large observed scatter. In the simulation centred on our most massive host halo, many sightlines show strong red damping wings even for initial volume averaged neutral hydrogen fractions as low as ~ 10$^{-3}$. For QSO lifetimes long enough to allow growth of the central supermassive black hole while optically bright, we can reproduce the observed near-zone of ULAS J1120+0641 only with an IGM that is initially neutral. Our results suggest that larger samples of z > 7 QSOs will provide important constraints on the evolution of the neutral hydrogen fraction and thus on how late reionization ends.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Testing the effect of galactic feedback on the IGM at z ~ 6 with metal-line absorbers

We present models of low- and high-ionization metal-line absorbers (O I, C II, C IV and Mg II) during the end of the reionization epoch, at z ~ 6. Using four cosmological hydrodynamical simulations with different feedback schemes (including the Illustris and Sherwood simulations) and two different choices of hydro-solver, we investigate how the overall incidence rate and equivalent width distribution of metal-line absorbers varies with the galactic wind prescription. We find that the O I and C II absorbers are reasonably insensitive to the feedback scheme. All models, however, struggle to reproduce the observations of C IV and Mg II, which are probing down to lower overdensities than O I and C II at z ~ 6, suggesting that the metals in the simulations are not being transported out into the IGM efficiently enough. The situation is improved but not resolved if we choose a harder (but still reasonable) and/or (locally) increased UV background at z ~ 6.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The effect of AGN feedback on the X-ray morphologies of clusters -- simulations vs. observations

We study the effect of Active Nuclei Galaxy (AGN) feedback as one of the major mechanisms modifying the cluster morphology influencing scaling relations, which are the most uncertain factor in constraining cosmology with clusters of galaxies. Using cosmological hydrodynamical simulations we investigate how the AGN feedback changes the X-ray morphology of the simulated systems, and compare to the observed REXCESS (Representative XMM-Newton Cluster Structure Survey) clusters. We apply centre shifts and power ratios to characterise the cluster morphology, and find that our simulated clusters are more substructured than the observed ones. We show that the degree of this discrepancy is affected by the inclusion of AGN feedback. While the clusters simulated with the AGN feedback are in much better agreement with the REXCESS L_X-T relation, they are also more substructured, which increases the tension with observations. This suggests that not only global cluster properties such as L_X and T and radial profiles should be used to compare and to calibrate simulations with observations, but also substructure measures such as centre shifts and power ratios. We discuss what changes in the simulations might ease the tension with observational constraints on these quantities.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Linear Instability of Dilute Ultrarelativistic $e^{\pm}$ Pair Beams

The annihilation of TeV photons from extragalactic TeV sources and the extragalactic background light produces ultrarelativistic $e^{\pm}$ beams, which are subject to powerful plasma instabilities that sap their kinetic energy. Here we study the linear phase of the plasma instabilities that these pair beams drive. To this end, we calculate the linear growth rate of the beam plasma and oblique instability in the electrostatic approximation in both the reactive and kinetic regimes, assuming a Maxwell-J{ü}ttner distribution for the pair beam. We reproduce the well-known reactive and kinetic growth rates for both the beam plasma and oblique mode. We demonstrate for the oblique instability that there is a broad spectrum of unstable modes that grow at the maximum rate for a wide range of beam temperatures and wave vector orientations relative to the beam. We also delineate the conditions for applicability for the reactive and kinetic regimes and find that the beam plasma mode transitions to the reactive regime at a lower Lorentz factor than the oblique mode due to a combination of their different scalings and the anisotropy of the velocity dispersions. Applying these results to the ultrarelativistic $e^{\pm}$ beams from TeV blazars, we confirm that these beams are unstable to both, the kinetic oblique mode and the reactive beam-plasma mode. These results are important in understanding how powerful plasma instabilities may sap the energy of the ultrarelativistic $e^{\pm}$ beams as they propagate through intergalactic space.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Sherwood simulation suite: overview and data comparisons with the Lyman-alpha forest at redshifts 2 < z < 5

We introduce a new set of large scale, high resolution hydrodynamical simulations of the intergalactic medium: the Sherwood simulation suite. These are performed in volumes 10^3-160^3 h^-3 comoving Mpc^{3}, span almost four orders of magnitude in mass resolution with up to 17.2 billion particles, and employ a variety of physics variations including warm dark matter and galactic outflows. We undertake a detailed comparison of the simulations to high resolution, high signal-to-noise observations of the Lyman-alpha forest over the redshift range 2 < z < 5. The simulations are in very good agreement with the observational data, lending further support to the paradigm that the Lyman-alpha forest is a natural consequence of the web-like distribution of matter arising in LCDM cosmological models. Only a small number of minor discrepancies remain with respect to the observational data. Saturated Lyman-alpha absorption lines with column densities N_HI > 10^14.5 cm^-2 at 2 < z < 2.5 are underpredicted in the models. An uncertain correction for continuum placement bias is required to match the distribution and power spectrum of the transmitted flux, particularly at z > 4. Finally, the temperature of intergalactic gas in the simulations may be slightly too low at z=2.7 and a flatter temperature-density relation is required at z=2.4, consistent with the expected effects of non-equilibrium ionisation during He-II reionisation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Zoomed cosmological simulations of Milky Way sized halos in f(R)-gravity

We investigate the impact of f(R) modified gravity on the internal properties of Milky Way sized dark matter halos in a set of cosmological zoom simulations of seven halos from the Aquarius suite, carried out with our code MG-GADGET in the Hu & Sawicki f(R) model. Also, we calculate the fifth forces in ideal NFW-halos as well as in our cosmological simulations and compare them against analytic model predictions for the fifth force inside spherical objects. We find that these theoretical predictions match the forces in the ideal halos very well, whereas their applicability is somewhat limited for realistic cosmological halos. Our simulations show that f(R) gravity significantly affects the dark matter density profile of Milky Way sized objects as well as their circular velocities. In unscreened regions, the velocity dispersions are increased by up to 40% with respect to LCDM for viable f(R) models. This difference is larger than reported in previous works. The Solar circle is fully screened in $f_{R0} = -10^{-6}$ models for Milky Way sized halos, while this location is unscreened for slightly less massive objects. Within the scope of our limited halo sample size, we do not find a clear dependence of the concentration parameter of dark matter halos on $f_{R0}$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Lyman-alpha emitters gone missing: evidence for late reionization?

We combine high resolution hydrodynamical simulations with an intermediate resolution, dark matter only simulation and an analytical model for the growth of ionized regions to estimate the large scale distribution and redshift evolution of the visibility of Lyman-alpha emission in 6<=z<=8 galaxies. The inhomogeneous distribution of neutral hydrogen during the reionization process results in significant fluctuations in the Lyman-alpha transmissivity on large scales. The transmissivity depends not only on the ionized fraction of the intergalactic medium by volume and the amplitude of the local ionizing background, but is also rather sensitive to the evolution of the relative velocity shift of the Lyman-alpha emission line due to resonant scattering. We reproduce a decline in the space density of Lyman-alpha emitting galaxies as rapid as observed with a rather rapidly evolving neutral fraction between z=6-8, and a typical Lyman-alpha line velocity offset of 100 km/s redward of systemic at z=6 which decreases toward higher redshift. The new (02/2015) Planck results indicate such a recent end to reionization is no longer disfavoured by constraints from the cosmic microwave background.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Modified Gravity N-body Code Comparison Project

Self-consistent ${\it N}$-body simulations of modified gravity models are a key ingredient to obtain rigorous constraints on deviations from General Relativity using large-scale structure observations. This paper provides the first detailed comparison of the results of different ${\it N}$-body codes for the $f(R)$, DGP, and Symmetron models, starting from the same initial conditions. We find that the fractional deviation of the matter power spectrum from $Λ$CDM agrees to better than $1\%$ up to $k \sim 5-10~h/{\rm Mpc}$ between the different codes. These codes are thus able to meet the stringent accuracy requirements of upcoming observational surveys. All codes are also in good agreement in their results for the velocity divergence power spectrum, halo abundances and halo profiles. We also test the quasi-static limit, which is employed in most modified gravity ${\it N}$-body codes, for the Symmetron model for which the most significant non-static effects among the models considered are expected. We conclude that this limit is a very good approximation for all of the observables considered here.

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Patchy blazar heating: diversifying the thermal history of the intergalactic medium

TeV-blazars potentially heat the intergalactic medium (IGM) as their gamma rays interact with photons of the extragalactic background light to produce electron-positron pairs, which lose their kinetic energy to the surrounding medium through plasma instabilities. This results in a heating mechanism that is only weakly sensitive to the local density, and therefore approximately spatially uniform, naturally producing an inverted temperature-density relation in underdense regions. In this paper we go beyond the approximation of uniform heating and quantify the heating rate fluctuations due to the clustered distribution of blazars and how this impacts on the thermal history of the IGM. We analytically compute a filtering function that relates the heating rate fluctuations to the underlying dark matter density field. We implement it in the cosmological code GADGET-3 and perform large scale simulations to determine the impact of inhomogeneous heating. We show that, because of blazar clustering, blazar heating is inhomogeneous for z>= 2. At high redshift, the temperature-density relation shows an important scatter and presents a low temperature envelope of unheated regions, in particular at low densities and within voids. However, the median temperature of the IGM is close to that in the uniform case, albeit slightly lower at low redshift. We find that blazar heating is more complex than initially assumed and that the temperature-density relation is not unique. Our analytic model for the heating rate fluctuations couples well with large scale simulations and provides a cost-effective alternative to subgrid models.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Imprint of f(R) Gravity on Non-Linear Structure Formation

We test the imprint of f(R) modified gravity on the halo mass function, using N-body simulations and a theoretical model developed in (Kopp et al. 2013). We find a very good agreement between theory and simulations. We extend the theoretical model to the conditional mass function and apply it to the prediction of the linear halo bias in f(R) gravity. Using the halo model we obtain a prediction for the non-linear matter power spectrum accurate to ~10% at z=0 and up to k=2h/Mpc. We also study halo profiles for the f(R) models and find a deviation from the standard general relativity result up to 40%, depending on the halo masses and redshift. This has not been pointed out in previous analysis. Finally we study the number density and profiles of voids identified in these f(R) N-body simulations. We underline the effect of the bias and the sampling to identify voids. We find significant deviation from GR when measuring the f(R) void profiles with fR0<-10^{-6}.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The photoheating of the intergalactic medium in synthesis models of the UV background

We compare cosmological hydrodynamical simulations combined with the homogeneous metagalactic UV background (UVB) of Haardt & Madau (2012) (HM2012) to observations of the Lyman-alpha forest that are sensitive to the thermal and ionization state of the intergalactic medium (IGM). The transition from optically thick to thin photoheating predicted by the simple one-zone, radiative transfer model implemented by HM2012 predicts a thermal history that is in remarkably good agreement with the observed rise of the IGM temperature at z~3 if we account for the expected evolution of the volume filling factor of HeIII. Our simulations indicate that there may be, however, some tension between the observed peak in the temperature evolution and the rather slow evolution of the HeII opacities suggested by recent Hubble Space Telescope/COS measurements. The HM2012 UVB also underpredicts the metagalactic hydrogen photoionization rate required by our simulations to match the observed opacity of the forest at z>4 and z<2.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Cosmic Degeneracies I: Joint N-body Simulations of Modified Gravity and Massive Neutrinos

We present the first suite of cosmological N-body simulations that simultaneously include the effects of two different and theoretically independent extensions of the standard $Λ$CDM cosmological scenario - namely an $f(R)$ theory of Modified Gravity (MG) and a cosmological background of massive neutrinos - with the aim to investigate their possible observational degeneracies. We focus on three basic statistics of the large-scale matter distribution, more specifically the nonlinear matter power spectrum, the halo mass function, and the halo bias, for which we determine the deviation with respect to the fiducial $Λ$CDM cosmology in the context of both separate and combined simulations of $f(R)$ MG and massive neutrinos scenarios. Our results show that while these two extended models separately determine very prominent and potentially detectable features in all the three statistics, when we allow them to be simultaneously at work these features are strongly suppressed, resulting in much weaker deviations from the standard model's predictions. In particular, when an $f(R)$ gravity model with $f_{R0}=-1\times 10^{-4}$ is combined with a total neutrino mass of $Σ_{i}m_{ν_{i}}=0.4$ eV, the resulting matter power spectrum, halo mass function, and bias at z=0 are found to be consistent with the standard model's predictions at the 10%, 20%, and 5% accuracy levels, respectively. Therefore, our results imply an intrinsic theoretical limit to the effective discriminating power of present and future observational data sets with respect to these widely considered extensions of the standard cosmological scenario in the absence of independent measurements of the neutrino masses from laboratory experiments, even though the high-redshift evolution might still allow to partially break the degeneracy [Abridged].

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nonlinear growing neutrino cosmology

The energy scale of Dark Energy, $\sim 2 \times 10^{-3}$ eV, is a long way off compared to all known fundamental scales - except for the neutrino masses. If Dark Energy is dynamical and couples to neutrinos, this is no longer a coincidence. The time at which Dark Energy starts to behave as an effective cosmological constant can be linked to the time at which the cosmic neutrinos become nonrelativistic. This naturally places the onset of the Universe's accelerated expansion in recent cosmic history, addressing the why-now problem of Dark Energy. We show that these mechanisms indeed work in the Growing Neutrino Quintessence model - even if the fully nonlinear structure formation and backreaction are taken into account, which were previously suspected of spoiling the cosmological evolution. The attractive force between neutrinos arising from their coupling to Dark Energy grows as large as $10^6$ times the gravitational strength. This induces very rapid dynamics of neutrino fluctuations which are nonlinear at redshift $z \approx 2$. Nevertheless, a nonlinear stabilization phenomenon ensures only mildly nonlinear oscillating neutrino overdensities with a large-scale gravitational potential substantially smaller than that of cold dark matter perturbations. Depending on model parameters, the signals of large-scale neutrino lumps may render the cosmic neutrino background observable.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Scaling relations and mass bias in hydrodynamical f(R) gravity simulations of galaxy clusters

We investigate the impact of chameleon-type f(R) gravity models on the properties of galaxy clusters and groups. Our f(R) simulations follow for the first time also the hydrodynamics of the intracluster and intragroup medium. This allows us to assess how f(R) gravity alters the X-ray scaling relations of clusters and how hydrostatic and dynamical mass estimates are biased when modifications of gravity are ignored in their determination. We find that velocity dispersions and intracluster medium temperatures are both increased by up to 1/3 in f(R) gravity in low-mass halos, while the difference disappears in massive objects. The mass scale of the transition depends on the background value f_R0 of the scalar degree of freedom. These changes in temperature and velocity dispersion alter the mass-temperature and X-ray luminosity-temperature scaling relations and bias dynamical and hydrostatic mass estimates that do not explicitly account for modified gravity towards higher values. Recently, a relative enhancement of X-ray compared to weak lensing masses was found by the Planck Collaboration (2013). We demonstrate that an explanation for this offset may be provided by modified gravity and the associated bias effects, which interestingly are of the required size. Finally, we find that the abundance of subhalos at fixed cluster mass is only weakly affected by f(R) gravity.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Effect of Nonlinear Landau Damping on Ultrarelativistic Beam Plasma Instabilities

Very-high energy gamma-rays from extragalactic sources pair-produce off of the extragalactic background light, yielding an electron-positron pair beam. This pair beam is unstable to various plasma instabilities, especially the "oblique" instability, which can be the dominant cooling mechanism for the beam. However, recently, it has been claimed that nonlinear Landau damping renders it physically irrelevant by reducing the effective damping rate to a low level. Here, we show with numerical calculations that the effective damping rate is $8\times 10^{-4}$ of the growth rate of the linear instability, which is sufficient for the "oblique" instability to be the dominant cooling mechanism of these pair beams. In particular, we show that previous estimates of this rate ignored the exponential cutoff in the scattering amplitude at large wavenumber and assumed that the damping of scattered waves entirely depends on collisions, ignoring collisionless processes. We find that the total wave energy eventually grows to approximate equipartition with the beam by increasingly depositing energy into long wavelength modes. As we have not included the effect of nonlinear wave-wave interactions on these long wavelength modes, this scenario represents the "worst-case" scenario for the oblique instability. As it continues to drain energy from the beam at a faster rate than other processes, we conclude that the "oblique" instability is sufficiently strong to make it the physically dominant cooling mechanism for high-energy pair beams in the intergalactic medium.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Lyman-alpha forest in f(R) modified gravity

In this work, we analyze the Lyman-$α$ forest in cosmological hydrodynamical simulations of chameleon-type f(R) gravity with the goal to assess whether the impact of such models is detectable in absorption line statistics. We carry out a set of hydrodynamical simulations with the cosmological simulation code MG-GADGET, including star formation and cooling effects, and create synthetic Lyman-$α$ absorption spectra from the simulation outputs. We statistically compare simulations with f(R) and ordinary general relativity, focusing on flux probability distribution functions (PDFs) and flux power-spectra, an analysis of the column density and line width distributions, as well as the matter power spectrum. We find that the influence of f(R) gravity on the Lyman-$α$ forest is rather small. Even models with strong modifications of gravity, like $|\bar{f}_{R0}| = 10^{-4}$, do not change the statistical Lyman-$α$ properties by more than 10%. The column density and line width distributions are hardly affected at all. It is therefore not possible to get competitive constraints on the background field $f_R$ using current observational data. An improved understanding of systematics in the observations and a more accurate modeling of the baryonic/radiative physics would be required to allow this in the future. The impact of f(R) on the matter power spectrum in our results is consistent with previous works.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Implications of Plasma Beam Instabilities for the Statistics of the Fermi Hard Gamma-ray Blazars and the Origin of the Extragalactic Gamma-Ray Background

Fermi has been instrumental in constraining the luminosity function and redshift evolution of gamma-ray bright blazars. This includes limits upon the spectrum and anisotropy of the extragalactic gamma-ray background (EGRB), redshift distribution of nearby Fermi active galactic nuclei (AGN), and the construction of a log(N)-log(S) relation. Based upon these, it has been argued that the evolution of the gamma-ray bright blazar population must be much less dramatic than that of other AGN. However, critical to such claims is the assumption that inverse Compton cascades reprocess emission above a TeV into the Fermi energy range, substantially enhancing the strength of the observed limits. Here we demonstrate that in the absence of such a process, due, e.g., to the presence of virulent plasma beam instabilities that preempt the cascade, a population of TeV-bright blazars that evolve similarly to quasars is consistent with the population of hard gamma-ray blazars observed by Fermi. Specifically, we show that a simple model for the properties and luminosity function is simultaneously able to reproduce their log(N)-log(S) relation, local redshift distribution, and contribution to the EGRB and its anisotropy without any free parameters. Insofar the naturalness of a picture in which the hard gamma-ray blazar population exhibits the strong redshift evolution observed in other tracers of the cosmological history of accretion onto halos is desirable, this lends support for the absence of the inverse Compton cascades and the existence of the beam plasma instabilities.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Modified Gravity-GADGET: A new code for cosmological hydrodynamical simulations of modified gravity models

We present a new massively parallel code for N-body and cosmological hydrodynamical simulations of modified gravity models. The code employs a multigrid-accelerated Newton-Gauss-Seidel relaxation solver on an adaptive mesh to efficiently solve for perturbations in the scalar degree of freedom of the modified gravity model. As this new algorithm is implemented as a module for the P-Gadget3 code, it can at the same time follow the baryonic physics included in P-Gadget3, such as hydrodynamics, radiative cooling and star formation. We demonstrate that the code works reliably by applying it to simple test problems that can be solved analytically, as well as by comparing cosmological simulations to results from the literature. Using the new code, we perform the first non-radiative and radiative cosmological hydrodynamical simulations of an f(R)-gravity model. We also discuss the impact of AGN feedback on the matter power spectrum, as well as degeneracies between the influence of baryonic processes and modifications of gravity.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Semi-analytic galaxy formation in f(R)-gravity cosmologies

Modifications of the equations of general relativity at large distances offer one possibility to explain the observed properties of our Universe without invoking a cosmological constant. Numerous proposals for such modified gravity cosmologies exist, but often their consequences for structure formation in the non-linear sector are not yet accurately known. In this work, we employ high-resolution numerical simulations of f(R)-gravity models coupled with a semi-analytic model (SAM) for galaxy formation to obtain detailed predictions for the evolution of galaxy properties. The f(R)-gravity models imply the existence of a `fifth-force', which is however locally suppressed, preserving the successes of general relativity on solar system scales. We show that dark matter haloes in f(R)-gravity models are characterized by a modified virial scaling with respect to the LCDM scenario, reflecting a higher dark matter velocity dispersion at a given mass. This effect is taken into account in the SAM by an appropriate modification of the mass--temperature relation. We find that the statistical properties predicted for galaxies (such as the stellar mass function and the cosmic star formation rate) in f(R)-gravity show generally only very small differences relative to LCDM, smaller than the dispersion between the results of different SAM models, which can be viewed as a measure of their systematic uncertainty. We also demonstrate that galaxy bias is not able to disentangle between f(R)-gravity and the standard cosmological scenario. However, f(R)-gravity imprints modifications in the linear growth rate of cosmic structures at large scale, which can be recovered from the statistical properties of large galaxy samples.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Shaping the galaxy stellar mass function with supernova- and AGN-driven winds

Cosmological hydrodynamical simulations of galaxy formation in representative regions of the Universe typically need to resort to subresolution models to follow some of the feedback processes crucial for galaxy formation. Here, we show that an energy-driven outflow model in which the wind velocity decreases and the wind mass loading increases in low-mass galaxies, as suggested by observations, can produce a good match to the low-mass end of the observed galaxy stellar mass function. The high-mass end can be recovered simultaneously if feedback from active galactic nuclei (AGN) and a correction for diffuse stellar light plausibly missed in observations are included. At the same time, our model is in good agreement with the stellar mass functions at redshifts z=1 and z=2, and with the observed redshift evolution of the cosmic star formation rate density. In addition, it accurately reproduces the observed gas to stellar mass ratios and specific star formation rates of galaxies as a function of their stellar mass. This agreement with a diverse set of data marks significant progress in hydrodynamically modelling the formation of a representative galaxy population. It also suggests that the mass flux in real galactic winds should strongly increase towards low-mass galaxies. Without this assumption, an overproduction of galaxies at the faint-end of the galaxy luminosity function seems inevitable in our models.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Lyman-alpha forest in a blazar-heated Universe

It has been realised only recently that TeV emission from blazars can significantly heat the intergalactic medium (IGM) by pair-producing high-energy electrons and positrons, which in turn excite vigorous plasma instabilities, leading to a local dissipation of the pairs' kinetic energy. In this work, we use cosmological hydrodynamical simulations to model the impact of this blazar heating on the Lyman-alpha forest at redshifts z~2-3. We find that blazar heating produces an inverted temperature-density relation in the IGM and naturally resolves many of the problems present in previous simulations of the forest that included photoheating alone. In particular, our simulations with blazar heating simultaneously reproduce the observed effective optical depth and temperature as a function of redshift, the observed probability distribution functions of the transmitted flux, and the observed flux power spectra, over the full redshift range 2<z<3 analysed here. Additionally, by deblending the Lyman-alpha forest into a sum of thermally broadened individual lines, we find superb agreement with the observed lower cutoff of the line-width distribution and abundances of neutral hydrogen column densities. Using the most recent constraints on the cosmic ultraviolet (UV) background, this excellent agreement with observations does not require rescaling the amplitude of the UV background; a procedure that was routinely used in the past to match the observed level of transmitted flux. We also show that our blazar-heated model matches the data better than standard simulations even when such a rescaling is allowed. This concordance between Lyman-alpha data and simulations, which are based on the most recent cosmological parameters, suggests that the inclusion of blazar heating alleviates previous tensions on constraints for sigma_8 derived from Lyman-alpha measurements and other cosmological data. [abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Star-formation efficiency and metal enrichment of the intracluster medium in local massive clusters of galaxies

We have investigated the baryon-mass content in a subsample of 19 clusters of galaxies extracted from the X-ray flux-limited sample HIFLUGCS according to their positions in the sky. For these clusters, we measured total masses and characteristic radii on the basis of a rich optical spectroscopic data set, the physical properties of the intracluster medium (ICM) using XMM-Newton and ROSAT X-ray data, and total (galaxy) stellar masses utilizing the SDSS DR7 multi-band imaging. The observed (hot) gas-mass fractions are almost constant in this mass range. We confirm that the stellar mass fraction decreases as the total mass increases and shows (20+/-4)% scatter; in addition, we show that it decreases as the central entropy increases. The latter behavior supports a twofold interpretation where heating from merging quenches the star-formation activity of galaxies in massive systems, and feedback from supernovae and/or radio galaxies drives a significant amount of gas to the regions beyond r_{500} or, alternatively, a substantially large amount of intracluster light (ICL) is associated with galaxies in nonrelaxed systems. Furthermore, less massive clusters are confirmed to host less gas per unit total mass; however, they exhibit higher mass fractions in metals, so that their ICM is more metal-rich. This again supports the interpretation that in the potential wells of low-mass systems the star-formation efficiency of galaxies was high or, alternatively, some gas is missing from the hot phase of the ICM. The former hypothesis is preferred as the main driver of the mass-dependent metal enrichment since the total mass-to-optical luminosity ratio increases as the total mass increases.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Intracluster stars in simulations with AGN feedback

We use a set of high-resolution hydrodynamical simulations of clusters of galaxies to study the build-up of the intracluster light (ICL), an interesting and likely significant component of their total stellar mass. Our sample of groups and clusters includes AGN feedback and is of high enough resolution to accurately resolve galaxy populations down to the smallest galaxies that are expected to significantly contribute to the stellar mass budget. We describe and test four different methods to identify the ICL in simulations, thereby allowing us to assess the reliability of the measurements. For all of the methods, we consistently find a very significant ICL stellar fraction (~45%) which exceeds the values typically inferred from observations. However, we show that this result is robust with respect to numerical resolution and integration accuracy, remarkably insensitive to changes in the star formation model, and almost independent of halo mass. It is also almost invariant when black hole growth is included, even though AGN feedback successfully prevents excessive overcooling in clusters and leads to a drastically improved agreement of the simulated cluster galaxy population with observations. In particular, the luminosities of central galaxies and the ages of their stellar populations are much more realistic when including AGN. In the light of these findings, it appears challenging to construct a simulation model that simultaneously matches the cluster galaxy population and at the same time produces a low ICL component. We find that intracluster stars are preferentially stripped in a cluster's densest region from massive galaxies that fall into the cluster at z>1. Surprisingly, some of the intracluster stars also form in the intracluster medium inside cold gas clouds that are stripped out of infalling galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Cluster strong lensing in the Millennium Simulation: The effect of galaxies and structures along the line-of-sight

We use ray-tracing through the Millennium simulation to study how secondary matter structures along the line-of-sight and the stellar mass in galaxies affect strong cluster lensing, in particular the cross-section for giant arcs. Furthermore, we investigate the distribution of the cluster Einstein radii and the radial distribution of giant arcs. We find that additional structures along the line-of-sight increase the strong-lensing optical depth by ~10-25%, while strong-lensing cross-sections of individual clusters are frequently boosted by as much as ~50%. The enhancement is mainly due to structures that are not correlated with the lens. Cluster galaxies increase the strong-lensing optical depth by up to a factor of 2, while interloping galaxies are not significant. We conclude that these effects need to be taken into account for predictions of the giant arc abundance, but they are not large enough to fully account for the reported discrepancy between predicted and observed abundances. Furthermore, we find that Einstein radii defined via the area enclosed by the critical curve are 10-30% larger than those defined via radial surface mass density profiles. The contributions of radial and tangential arcs to the radial distribution of arcs can be clearly distinguished. The radial distribution of tangential arcs is very broad and extends out to several Einstein radii. Thus, individual arcs are not well suited for constraining Einstein radii.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Simulations of AGN feedback in galaxy clusters and groups

There is compelling evidence that black holes (BHs) in cluster centers vigorously interact with their surroundings, indicating that any realistic model of cluster formation needs to account for these processes. Here we use high-resolution cosmological simulations of a large cluster and group sample to study how BHs affect their host systems. We focus on two specific properties, the halo gas fraction and the X-ray luminosity-temperature scaling relation, both of which are notoriously difficult to reproduce in self-consistent hydrodynamical simulations. We show that BH feedback can solve both of these issues, bringing them in excellent agreement with observations, without alluding to the `cooling only' solution that produces unphysically bright central galaxies.

0 citations0 reviewsNo code linkedOpen access