Deep observations have detected extended Ly$\alpha$ 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$\alpha$, while CIII] shows a shallower profile similar to Ly$\alpha$ 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$\alpha$ haloes, which has much weaker metal emission lines relative to Ly$\alpha$. 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.

We compare a sample of five high-resolution, high S/N Ly$\alpha$ 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.

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.

Theories of gravity that incorporate new scalar degrees of freedom typically require screening mechanisms to ensure consistency with Solar System tests. One widely-studied mechanism -- the chameleon -- can lead to violations of the equivalence principle (EP), as screened and unscreened objects fall differently. If the stars are screened but the surrounding dark matter is not, this leads to asymmetry between leading and trailing streams. We provide analytic estimates of the magnitude of this effect for realistic Galactic mass distributions. Using a restricted N-body code, we simulate 4 satellites with a range of masses and orbits, together with a variety of strengths of the fifth force and screening levels of the Milky Way and satellite. The ratio of the cumulative number function of stars in the leading and trailing stream as a function of longitude from the satellite is computable from simulations, measurable from the stellar data and can provide a direct test. We forecast constraints for streams at large Galactocentric distances, using the specific example case of Hu-Sawicki gravity. Streams with apocentres between 100 and 200 kpc provide attainable constraints at the level of $|f_{R0}| = 10^{-7}$. Still more stringent constraints at the level of $10^{-7.5}$ or even $10^{-8}$ are plausible provided the environmental screening of the satellite is accounted for. These would be among the tightest astrophysical constraints to date. We note further signatures of chameleon gravity: (i) the trailing stellar stream may become detached from the dark matter progenitor if all the stars are lost, (ii) in the extreme fifth force regime, striations in the stellar trailing tail may develop, (iii) if the satellite is fully screened, its orbital frequency is lower than that of the associated dark matter, which is preferentially liberated into the leading tidal tail.

Pair creation on the cosmic infrared background and subsequent inverse-Compton scattering on the CMB potentially reprocesses the TeV emission of blazars into faint GeV halos with structures sensitive to intergalactic magnetic fields (IGMF). Previous work has shown that these halos are then highly-anisotropic and extended. If the coherence length of the IGMF is greater than the inverse-Compton cooling length of the pairs, then the orientation of the gamma-ray halo will be correlated with the direction of the magnetic field which is unknown and expected to change for each source. In order to constructively add each source we then use angular power spectra which are insensitive to the jet orientation. By looking at known GeV blazars detected by Fermi, we exclude the existence of an IGMF with coherence lengths $>100$ Mpc at greater than $3.9\sigma$ with current-day strengths in the range $10^{-16}$ to $10^{-15}$ G, and at 2$\sigma$ from $10^{-17}$ to $10^{-14}$ G. This provides a direct measurement of the non-existence of gamma-ray halos, providing an important check on previous results.

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.

Over the last decades, cosmological simulations of galaxy formation have been instrumental for advancing our understanding of structure and galaxy formation in the Universe. These simulations follow the non-linear evolution of galaxies modeling a variety of physical processes over an enormous range of scales. A better understanding of the physics relevant for shaping galaxies, improved numerical methods, and increased computing power have led to simulations that can reproduce a large number of observed galaxy properties. Modern simulations model dark matter, dark energy, and ordinary matter in an expanding space-time starting from well-defined initial conditions. The modeling of ordinary matter is most challenging due to the large array of physical processes affecting this matter component. Cosmological simulations have also proven useful to study alternative cosmological models and their impact on the galaxy population. This review presents a concise overview of the methodology of cosmological simulations of galaxy formation and their different applications.

We study the gas and stellar mass content of galaxy groups and clusters in the FABLE suite of cosmological hydrodynamical simulations, including the evolution of their central brightest cluster galaxies (BCGs), satellite galaxies and intracluster light (ICL). The total gas and stellar mass of FABLE clusters are in very good agreement with observations and show negligible redshift evolution at fixed halo mass for $M_{500} \gtrsim 3 \times 10^{14} M_{\odot}$ at $z \lesssim 1$, in line with recent findings from Sunyaev-Zel'dovich (SZ)-selected cluster samples. Importantly, the simulations predict significant redshift evolution in these quantities in the low mass ($M_{500} \sim 10^{14} M_{\odot}$) regime, which will be testable with upcoming SZ surveys such as SPT-3G. While the stellar masses of FABLE BCGs are in reasonable agreement with observations, the total stellar mass in satellite galaxies is lower than observed and the total mass in ICL is somewhat higher. This may be caused by enhanced tidal stripping of satellite galaxies due to their large sizes. BCGs are characterised by moderate stellar mass growth at $z < 1$ coincident with a late-time development of the ICL. The level of BCG mass growth is in good agreement with recent observations, however, we caution that the inferred growth depends sensitively on the mass definition. We further show that in-situ star formation contributes more than half the mass of a BCG over its lifetime, the bulk of which is gained at $z > 1$ where star formation rates are highest. The stellar mass profiles of the BCG+ICL component are similar to observed profiles out to $\sim 100$ kpc at $z \approx 0$ and follow a close to power law shape out to several hundred kpc. We further demonstrate that the inferred size growth of BCGs can be severely biased by the choice of parametric model and the outer radius of the fit.

Galaxy cluster outskirts mark the transition region from the mildly non-linear cosmic web to the highly non-linear, virialised, cluster interior. It is in this transition region that the intra-cluster medium (ICM) begins to influence the properties of accreting galaxies and groups, as ram pressure impacts a galaxy's cold gas content and subsequent star formation rate. Conversely, the thermodynamical properties of the ICM in this transition region should also feel the influence of accreting substructure (i.e. galaxies and groups), whose passage can drive shocks. In this paper, we use a suite of cosmological hydrodynamical zoom simulations of a single galaxy cluster, drawn from the nIFTy comparison project, to study how the dynamics of substructure accreted from the cosmic web influences the thermodynamical properties of the ICM in the cluster's outskirts. We demonstrate how features evident in radial profiles of the ICM (e.g. gas density and temperature) can be linked to strong shocks, transient and short-lived in nature, driven by the passage of substructure. The range of astrophysical codes and galaxy formation models in our comparison are broadly consistent in their predictions (e.g. agreeing when and where shocks occur, but differing in how strong shocks will be); this is as we would expect of a process driven by large-scale gravitational dynamics and strong, inefficiently radiating, shocks. This suggests that mapping such shock structures in the ICM in a cluster's outskirts (via e.g. radio synchrotron emission) could provide a complementary measure of its recent merger and accretion history.

Jet feedback from active galactic nuclei (AGN) harboured by brightest cluster galaxies is expected to play a fundamental role in regulating cooling in the intracluster medium (ICM). While observations and theory suggest energy within jet lobes balances ICM radiative losses, the modus operandi of energy communication with the ICM remains unclear. We present simulations of very high resolution AGN-driven jets launching in a live, cosmological galaxy cluster, within the moving mesh code AREPO. As the jet propagates through the ICM the majority of its energy, which is initially in the kinetic form, thermalizes quickly through internal shocks and inflates lobes of very hot gas. The jets effectively heat the cluster core, with $PdV$ work and weather-aided mixing being the main channels of energy transfer from the lobes to the ICM, while strong shocks and turbulence are subdominant. We additionally present detailed mock X-ray maps at different stages of evolution, revealing clear cavities surrounded by X-ray bright rims, with lobes being detectable for up to $\sim 10^8$ yr even when magnetic draping is ineffective. We find bulk motions in the cluster can significantly affect lobe propagation, offsetting them from the jet direction and imparting bulk velocities that can dominate over the buoyantly rising motion.

We study the redshift evolution of the X-ray and Sunyaev-Zel'dovich (SZ) scaling relations for galaxy groups and clusters in the FABLE suite of cosmological hydrodynamical simulations. Using an expanded sample of $27$ high-resolution zoom-in simulations, together with a uniformly-sampled cosmological volume to sample low-mass systems, we find very good agreement with the majority of observational constraints up to $z \sim 1$. We predict significant deviations of all examined scaling relations from the simple self-similar expectations. While the slopes are approximately independent of redshift, the normalisations evolve positively with respect to self-similarity, even for commonly-used mass proxies such as the $Y_{\mathrm{X}}$ parameter. These deviations are due to a combination of factors, including more effective AGN feedback in lower mass haloes, larger binding energy of gas at a given halo mass at higher redshifts and larger non-thermal pressure support from kinetic motions at higher redshifts. Our results have important implications for cluster cosmology from upcoming SZ surveys such as SPT-3G, ACTpol and CMB-S4, as relatively small changes in the observable--mass scaling relations (within theoretical uncertainties) have a large impact on the predicted number of high-redshift clusters and hence on our ability to constrain cosmology using cluster abundances. In addition, we find that the intrinsic scatter of the relations, which agrees well with most observational constraints, increases at lower redshifts and for lower mass systems. This calls for a more complex parametrization than adopted in current observational studies to be able to accurately account for selection biases.

In chameleon $f(R)$-gravity, the fifth force will lead to `upturns' in galaxy rotation curves near the screening radius. The location of the upturn depends on the cosmic background value of the scalar field $f_{R0}$, as well as the mass, size and environment of the galaxy. We search for this signature of modified gravity in the SPARC sample of measured rotation curves, using an MCMC technique to derive constraints on $f_{R0}$. Assuming NFW dark matter haloes and with $f_{R0}$ freely varying for each galaxy, most galaxies prefer $f(R)$ gravity to $\Lambda$CDM, but there is a large spread of inferred $f_{R0}$ values, inconsistent with a single global value. Requiring instead a consistent $f_{R0}$ value for the whole sample, models with $\log_{10}|f_{R0}| > -6.1$ are excluded. On the other hand, models in the range $-7.5<\log_{10}|f_{R0}|<-6.5$ seem to be favoured with respect to $\Lambda$CDM, with a significant peak at -7. However, this signal is largely a result of galaxies for which the $f(R)$ signal is degenerate with the core/cusp problem, and when the NFW profile is replaced with a cored halo profile, $\Lambda$CDM gives better fits than any given $f(R)$ model. Thus, we find no convincing evidence of $f(R)$ gravity down to the level of $|f_{R0}| \sim 6 \times 10^{-8}$, with the caveat that if cored halo density profiles cannot ultimately be explained within $\Lambda$CDM, a screened modified gravity theory could possibly provide an alternative solution for the core/cusp problem. However, the $f(R)$ models studied here fall short of achieving this.

Mapping the intergalactic medium (IGM) in Lyman-$\alpha$ emission would yield unprecedented tomographic information on the large-scale distribution of baryons and potentially provide new constraints on the UV background and various feedback processes relevant for galaxy formation. Here, we examine the Lyman-$\alpha$ emission of the moderate-density IGM due to collisional excitations and recombinations in the presence of a UV background in the Sherwood simulation suite. We focus on large-scale-structure filaments in which Lyman-$\alpha$ radiative transfer effects are expected to be moderate. At low density the emission is primarily due to fluorescent re-emission of the ionising UV background due to recombinations, while collisional excitations dominate at higher densities. We discuss prospects of current and future observational facilities to detect this emission and focus on VLT/MUSE for a more detailed sensitivity analysis. We construct mock MUSE observations resembling the MUSE Hubble Deep Field South from our simulations and show that our predictions are consistent with recent analyses of diffuse Lyman-$\alpha$ emission using MUSE at redshifts $3 < z < 6$. We find that it should be most feasible to detect the Lyman-$\alpha$ emission from filaments in the IGM in overdense regions, somewhat surprisingly towards the high-redshift end ($z \gtrsim 4.5$) accessible by MUSE, and if narrowband widths close to the MUSE spectral resolution limit are used. This is due to the higher densities and lower temperatures in the IGM at higher redshift. High-redshift protoclusters therefore appear to be the ideal environment to observe filamentary structures in the IGM in Lyman-$\alpha$ emission.

High-redshift QSO spectra show large spatial fluctuations in the Ly-alpha opacity of the intergalactic medium on surprisingly large scales at z>~5.5. We present a radiative transfer simulation of cosmic reionization driven by galaxies that reproduces this large scatter and the rapid evolution of the Ly-alpha opacity distribution at 5<z<6. The simulation also reproduces the low Thomson scattering optical depth reported by the latest CMB measurement and is consistent with the observed short near-zones and strong red damping wings in the highest-redshift QSOs. It also matches the rapid disappearance of observed Ly-alpha emission by galaxies at z>~6. Reionization is complete at z=5.3 in our model, and 50% of the volume of the Universe is ionized at z=7. Agreement with the Ly-alpha forest data in such a late reionization model requires a rapid evolution of the ionizing emissivity of galaxies that peaks at z~6.8. The late end of reionization results in a large scatter in the photoionisation rate and the neutral hydrogen fraction at redshifts as low as z<~5.5 with large residual neutral 'islands' that can produce very long Gunn-Peterson troughs resembling those seen in the data.

An intergalactic magnetic field stronger than $3\times10^{-13}$~G would explain the lack of a bright, extended degree-scale, GeV-energy inverse Compton component in the gamma-ray spectra of TeV-blazars. A robustly predicted consequence of the presence of such a field is the existence of degree-scale GeV-energy gamma-ray halos -- gamma-ray bow ties -- about TeV-bright active galactic nuclei, corresponding to more than half of all radio galaxies. However, the emitting regions of these halos are confined to and aligned with the direction of the relativistic jets associated with gamma-ray sources. Based on the orientation of radio jets, we align and stack corresponding degree-scale gamma-ray images of isolated Fanaroff-Riley class I and II objects and exclude the existence of these halos at overwhelming confidence, limiting the intergalactic field strength to $<10^{-15}$~G for large-scale fields and progressively larger in the diffusive regime when the correlation length of the field becomes small in comparison to 1 Mpc. When combined with prior limits on the strength of the intergalactic magnetic field, this excludes a purely magnetic explanation for the absence of halos. Thus, it requires the existence of novel physical processes that preempt the creation of halos, e.g., the presence of beam-plasma instabilities in the intergalactic medium or a drastic cutoff of the very high energy spectrum of these sources.

Current constraints on gravity are relatively weak on galactic and intergalactic scales. Screened modified gravity models can exhibit complex behaviour there without violating stringent tests of gravity within our Solar System. They might hence provide viable extensions of the theory of gravity. Here, we use galaxy kinematics to constrain screened modified gravity models. We focus on chameleon $f(R)$ gravity and predict its impact on galaxy rotation curves and radial acceleration relations. This is achieved by post-processing state-of-the-art galaxy formation simulations from the \textsc{auriga project}, using the \textsc{mg-gadget} code. For a given galaxy, the surface dividing screened and un-screened regions adopts an oblate shape, reflecting the disc morphology of the galaxy's mass distribution. At the `screening radius'---the radius at which screening is triggered in the disc plane---characteristic `upturns' are present in both rotation curves and radial acceleration relations. The locations of these features depend on various factors, such as the galaxy mass, the concentration of the density profile and the value of the background field amplitude $f_{R0}$. Self-screening of stars and environmental screening also play a role. For Milky Way-size galaxies, we find that a model with $|f_{R0}|=10^{-7}$ results in rotation curves that are indistinguishable from $\Lambda$CDM, while for $|f_{R0}| \geq 2 \times 10^{-6}$ the simulated galaxies are entirely unscreened, violating Solar System constraints. For intermediate values, distinct upturns are present. With a careful statistical analysis of existing samples of observed rotation curves, including lower mass objects, constraints on $f(R)$ gravity with a sensitivity down to $|f_{R0}|\sim10^{-7}$ should be possible.

We present the Feedback Acting on Baryons in Large-scale Environments (FABLE) suite of cosmological hydrodynamical simulations of galaxies, groups and clusters. The simulations use the AREPO moving-mesh code with a set of physical models for galaxy formation based on the successful Illustris simulation, but with updated AGN and supernovae feedback models. This allows us to simultaneously reproduce the observed redshift evolution of the galaxy stellar mass function together with the stellar and gas mass fractions of local groups and clusters across a wide range of halo masses. Focusing on the properties of groups and clusters, we find very good agreement with a range of observed scaling relations, including the X-ray luminosity--total mass and gas mass relations as well as the total mass--temperature and Sunyaev-Zel'dovich flux--mass relations. Careful comparison of our results with scaling relations based on X-ray hydrostatic masses as opposed to weak lensing-derived masses reveals some discrepancies, which hint towards a non-negligible X-ray mass bias in observed samples. We further show that radial profiles of density, pressure and temperature of the simulated intracluster medium are in very good agreement with observations, in particular for $r > 0.3\ r_{500}$. In the innermost regions however we find too large entropy cores, which indicates that a more sophisticated modelling of the physics of AGN feedback may be required to accurately reproduce the observed populations of cool-core and non-cool-core clusters.

We model the transmission of the Lyman-alpha line through the circum- and intergalactic media around dark matter haloes expected to host Lyman-alpha emitters (LAEs) at z > 5.7, using the high-dynamic-range Sherwood simulations. We find very different CGM environments around more massive haloes (~10^11 M_sun) compared to less massive haloes (~10^9 M_sun) at these redshifts, which can contribute to a different evolution of the Lyman-alpha transmission from LAEs within these haloes. Additionally we confirm that part of the differential evolution could result from bright LAEs being more likely to reside in larger ionized regions. We conclude that a combination of the CGM environment and the IGM ionization structure is likely to be responsible for the differential evolution of the bright and faint ends of the LAE luminosity function at z > 6. More generally, we confirm the suggestion that the self-shielded neutral gas in the outskirts of the host halo can strongly attenuate the Lyman-alpha emission from high redshift galaxies. We find that this has a stronger effect on the more massive haloes hosting brighter LAEs. The faint-end of the LAE luminosity function is thus a more reliable probe of the average ionization state of the IGM. Comparing our model for LAEs with a range of observational data we find that the favoured reionization histories are our previously advocated `Late' and `Very Late' reionization histories, in which reionization finishes rather rapidly at around z ~ 6.

We introduce a set of four very high resolution cosmological simulations for exploring $f(R)$-gravity, with $2048^3$ particles in $768\,h^{-1}\textrm{Mpc}$ and $1536\,h^{-1}\textrm{Mpc}$ simulation boxes, both for a $|\overline{f_{R0}}| = 10^{-5}$ model and a $\Lambda$CDM comparison universe, making the set the largest simulations of $f(R)$-gravity to date. In order to mimic real observations, the simulations include a continuous 2D and 3D lightcone output which is dedicated to study lensing and clustering statistics in modified gravity. In this work, we present a detailed analysis and resolution study for the matter power spectrum in $f(R)$-gravity over a wide range of scales. We also analyse the angular matter power spectrum and lensing convergence on the lightcone. In addition, we investigate the impact of modified gravity on the halo mass function, matter and halo auto-correlation functions, linear halo bias and the concentration-mass relation. We find that the impact of $f(R)$-gravity is generally larger on smaller scales and increases with decreasing redshift. Comparing our simulations to state-of-the-art hydrodynamical simulations we confirm a degeneracy between $f(R)$-gravity and baryonic feedback in the matter power spectrum on small scales, but also find that scales around $k = 1\, h\, {\rm Mpc}^{-1}$ are promising to distinguish both effects. The lensing convergence power spectrum is increased in $f(R)$-gravity. Interestingly available numerical fits are in good agreement overall with our simulations for both standard and modified gravity, but tend to overestimate their relative difference on non-linear scales by few percent. We also find that the halo bias is lower in $f(R)$-gravity compared to general relativity, whereas halo concentrations are increased for unscreened halos.

We calibrate here cosmological radiative transfer simulation with ATON/RAMSES with a range of measurements of the Lyman alpha opacity from QSO absorption spectra. We find the Lyman alpha opacity to be very sensitive to the exact timing of hydrogen reionisation. Models reproducing the measured evolution of the mean photoionisation rate and average mean free path reach overlap at z ~ 7 and predict an accelerated evolution of the Lyman alpha opacity at z > 6 consistent with the rapidly evolving luminosity function of Lyman alpha emitters in this redshift range. Similar to "optically thin" simulations our full radiative transfer simulations fail, however, to reproduce the high-opacity tail of the Lyman alpha opacity PDF at z > 5. We argue that this is due to spatial UV fluctuations in the post-overlap phase of reionisation on substantially larger scales than predicted by our source model, where the ionising emissivity is dominated by large numbers of sub-L* galaxies. We further argue that this suggests a significant contribution to the ionising UV background by much rarer bright sources at high redshift.

The thermal state of the post-reionization IGM is sensitive to the timing of reionization and the nature of the ionizing sources. We have modelled here the thermal state of the IGM in cosmological radiative transfer simulations of a realistic, extended, spatially inhomogeneous hydrogen reionization process, carefully calibrated with Ly-alpha forest data. We compare these with cosmological simulations run using a spatially homogeneous ionizing background. The simulations with a realistic growth of ionized regions and a realistic spread in reionization redshifts show, as expected, significant spatial fluctuations in the temperature-density relation (TDR) of the post-reionization IGM. The most recently ionized regions are hottest and exhibit a flatter TDR. In simulations consistent with the average TDR inferred from Ly-alpha forest data, these spatial fluctuations have a moderate but noticeable effect on the statistical properties of the Ly-alpha opacity of the IGM at z ~ 4-6. This should be taken into account in accurate measurements of the thermal properties of the IGM and the free-streaming of dark matter from Ly-alpha forest data in this redshift range. The spatial variations of the TDR predicted by our simulations are, however, smaller by about a factor two than would be necessary to explain the observed large spatial opacity fluctuations on large (> 50 comoving Mpc/h) scales at z > 5.5.

Understanding how baryonic processes shape the intracluster medium (ICM) is of critical importance to the next generation of galaxy cluster surveys. However, most models of structure formation neglect potentially important physical processes, like anisotropic thermal conduction (ATC). In this letter, we explore the impact of ATC on the prevalence of cool-cores (CCs) using 12 pairs of magnetohydrodynamical galaxy cluster simulations, simulated using the IllustrisTNG model with and without ATC. Although the impact of ATC varies from cluster to cluster and with CC criterion, its inclusion produces a systematic shift to larger CC fractions at z = 0 for all CC criteria considered. Additionally, the inclusion of ATC yields a flatter CC fraction redshift evolution, easing the tension with the observed evolution. With ATC included, the energy required for the central black hole to achieve self-regulation is reduced and the gas fraction in the cluster core increases, resulting in larger CC fractions. ATC makes the ICM unstable to perturbations and the increased efficiency of AGN feedback suggests that its inclusion results in a greater level of mixing in the ICM. Therefore, ATC is potentially an important physical process in reproducing the thermal structure of the ICM.

High-signal to noise observations of the Ly$\alpha$ forest transmissivity in the z = 7.085 QSO ULAS J1120+0641 show seven narrow transmission spikes followed by a long 240 cMpc/h trough. Here we use radiative transfer simulations of cosmic reionization previously calibrated to match a wider range of Ly$\alpha$ forest data to show that the occurrence of seven transmission spikes in the narrow redshift range z = 5.85 - 6.1 is very sensitive to the exact timing of reionization. Occurrence of the spikes requires the most under dense regions of the IGM to be already fully ionised. The rapid onset of a long trough at z = 6.12 requires a strong decrease of the photo-ionisation rate at z$\sim$6.1 in this line-of-sight, consistent with the end of percolation at this redshift. The narrow range of reionisation histories that we previously found to be consistent with a wider range of Ly$\alpha$ forest data have a reasonable probability of showing seven spikes and the mock absorption spectra provide an excellent match to the spikes and the trough in the observed spectrum of ULAS J1120+0641. Despite the large overall opacity of Ly$\alpha$ at z > 5.8, larger samples of high signal-to-noise observations of rare transmission spikes should therefore provide important further insights into the exact timing of the percolation of HII bubbles at the tail-end of reionization

We explore how inhomogeneity in the background plasma number density alters the growth of electrostatic unstable wavemodes of beam plasma systems. This is particularly interesting for blazar-driven beam-plasma instabilities, which may be suppressed by inhomogeneities in the intergalactic medium as was recently claimed in the literature. Using high resolution Particle-In-Cell simulations with the SHARP code, we show that the growth of the instability is local, i.e., regions with almost homogeneous background density will support the growth of the Langmuir waves, even when they are separated by strongly inhomogeneous regions, resulting in an overall slower growth of the instability. We also show that if the background density is continuously varying, the growth rate of the instability is lower; though in all cases, the system remains within the linear regime longer and the instability is not extinguished. In all cases, the beam loses approximately the same fraction of its initial kinetic energy in comparison to the uniform case at non-linear saturation. Thus, inhomogeneities in the intergalactic medium are unlikely to suppress the growth of blazar-driven beam-plasma instabilities.

In principle, the angular anisotropy in the extragalactic gamma-ray background (EGRB) places severe constraints upon putative populations of unresolved gamma-ray point sources. Existing estimates of the EGRB anisotropy have been constructed by excising known point sources, e.g., taken from the First or 2 Year Fermi-LAT Source Catalog (1FGL or 2FGL, respectively) and statistically analyzing the residual gamma-ray sky maps. We perform an independent check of the EGRB anisotropy limits by comparing the values obtained from the 1FGL-masked sky maps to the signal implied by sources that lie below the 1FGL detection threshold in the more sensitive 2FGL and 1FHL (First Fermi-LAT catalog of >10 GeV sources). As such, our analysis provides an internal consistency check of implications for source counts and spectral index distributions of gamma-ray bright active galactic nuclei obtained from Fermi-LAT data. Based on this, we find evidence for substantially larger anisotropies than those previously reported at energies above 5 GeV, where BL Lac objects are likely to provide the bulk of their contribution to the EGRB. This uncertainty in the EGRB anisotropy cautions against using it as an independent constraint for the high-redshift gamma-ray universe. Moreover, this would suggest that contrary to previous claims, smooth extensions of the resolved point-source population may be able to simultaneously explain both the isotropic and anisotropic components of the EGRB.

Recent observations suggest that hydrogen reionization ends late ($z \simeq 6$) and proceeds quickly. We present here a new model of the meta-galactic UV/X-ray background (UVB) that is consistent with this. It adopts the most recent determinations of the ionizing emissivity due to stars and AGN, as well as of the HI absorber column density distribution. Another major improvement is a new treatment of the intergalactic medium (IGM) opacity for ionizing photons that is able to consistently capture the transition from a neutral to an ionized IGM. Previous synthesis models of the UVB, when used in simulations, yield reionization and thermal histories that are inconsistent with the assumed ionizing emissivities. With our new treatment, this discrepancy is fully resolved. In our fiducial model, galaxies leaking $\lesssim 18\%$ of their Lyman continuum emission drive HI reionization, while AGN drive HeII reionization (completing at $z \simeq 6.2$ and $2.8$ respectively). Due to the limited time available for cooling between HI and HeII reionization, higher IGM temperatures are predicted for late reionization scenarios. In our fiducial model, the predicted temperatures agree well with observational constraints at $z \lesssim 4$, while being slightly high compared to (somewhat uncertain) data above that. Models with a larger contribution of AGN are instead disfavoured by the temperature data, as well as by measurements of the HI and HeII Lyman-$\alpha$ forest opacities. We also present "equivalent-equilibrium" ionization/heating rates that mimic our fiducial UVB model for use in simulation codes that assume ionization equilibrium.

Many astrophysical plasmas are prone to beam-plasma instabilities. For relativistic and dilute beams, the {\it spectral} support of the beam-plasma instabilities is narrow, i.e., the linearly unstable modes that grow with rates comparable to the maximum growth rate occupy a narrow range of wave numbers. This places stringent requirements on the box-sizes when simulating the evolution of the instabilities. We identify the implied lower limits on the box size imposed by the longitudinal beam plasma instability, i.e., typically the most stringent condition required to correctly capture the linear evolution of the instabilities in multidimensional simulations. We find that sizes many orders of magnitude larger than the resonant wavelength are typically required. Using one-dimensional particle-in-cell simulations, we show that the failure to sufficiently resolve the spectral support of the longitudinal instability yields slower growth and lower levels of saturation, potentially leading to erroneous physical conclusion.

We study the effect of different feedback prescriptions on the properties of the low redshift ($z\leq1.6$) Ly$\alpha$ forest using a selection of hydrodynamical simulations drawn from the Sherwood simulation suite. The simulations incorporate stellar feedback, AGN feedback and a simplified scheme for efficiently modelling the low column density Ly$\alpha$ forest. We confirm a discrepancy remains between Cosmic Origins Spectrograph (COS) observations of the Ly$\alpha$ forest column density distribution function (CDDF) at $z \simeq 0.1$ for high column density systems ($N_{\rm HI}>10^{14}\rm\,cm^{-2}$), as well as Ly$\alpha$ velocity widths that are too narrow compared to the COS data. Stellar or AGN feedback -- as currently implemented in our simulations -- have only a small effect on the CDDF and velocity width distribution. We conclude that resolving the discrepancy between the COS data and simulations requires an increase in the temperature of overdense gas with $\Delta=4$--$40$, either through additional He$\,\rm \scriptstyle II\ $ photo-heating at $z>2$ or fine-tuned feedback that ejects overdense gas into the IGM at just the right temperature for it to still contribute significantly to the Ly$\alpha$ forest. Alternatively a larger, currently unresolved turbulent component to the line width could resolve the discrepancy.

Numerical heating in particle-in-cell (PIC) codes currently precludes the accurate simulation of cold, relativistic plasma over long periods, severely limiting their applications in astrophysical environments. We present a spatially higher-order accurate relativistic PIC algorithm in one spatial dimension, which conserves charge and momentum exactly. We utilize the smoothness implied by the usage of higher-order interpolation functions to achieve a spatially higher-order accurate algorithm (up to fifth order). We validate our algorithm against several test problems -- thermal stability of stationary plasma, stability of linear plasma waves, and two-stream instability in the relativistic and non-relativistic regimes. Comparing our simulations to exact solutions of the dispersion relations, we demonstrate that SHARP can quantitatively reproduce important kinetic features of the linear regime. Our simulations have a superior ability to control energy non-conservation and avoid numerical heating in comparison to common second-order schemes. We provide a natural definition for convergence of a general PIC algorithm: the complement of physical modes captured by the simulation, i.e., those that lie above the Poisson noise, must grow commensurately with the resolution. This implies that it is necessary to simultaneously increase the number of particles per cell and decrease the cell size. We demonstrate that traditional ways for testing for convergence fail, leading to plateauing of the energy error. This new PIC code enables us to faithfully study the long-term evolution of plasma problems that require absolute control of the energy and momentum conservation.

We present here predictions for the spatial distribution of 21 cm brightness temperature fluctuations from high-dynamic-range simulations for AGN-dominated reionization histories that have been tested against available Lyman-alpha and CMB data. We model AGN by extrapolating the observed M-sigma relation to high redshifts and assign them ionizing emissivities consistent with recent UV luminosity function measurements. 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 21 cm spin temperature is significantly larger than the CMB temperature. Our AGN-dominated reionization histories increase the variance of the 21 cm emission by a factor of up to ten compared to similar reionization histories dominated by faint galaxies, to values close to 100 mK^2 at scales accessible to experiments (k < 1 h/cMpc). This is lower than the sensitivity claimed to have been already reached by ongoing experiments by only a factor of about two or less. When reionization is dominated by AGN, the 21 cm power spectrum is enhanced on all scales due to the enhanced bias of the clustering of the more massive haloes and the peak in the large scale 21 cm power is strongly enhanced and moved to larger scales due to bigger characteristic bubble sizes. AGN dominated reionization should be easily detectable by LOFAR (and later HERA and SKA1) at their design sensitivity, assuming successful foreground subtraction and instrument calibration. Conversely, these could become the first non-trivial reionization scenarios to be ruled out by 21 cm experiments, thereby constraining the contribution of AGN to reionization.

We derive constraints on feedback by active galactic nuclei (AGN) by setting limits on their thermal Sunyaev-Zel'dovich (SZ) imprint on the cosmic microwave background (CMB). The amplitude of any SZ signature is small and degenerate with the poorly known sub-mm spectral energy distribution of the AGN host galaxy and other unresolved dusty sources along the line of sight. Here we break this degeneracy by combining microwave and sub-mm data from Planck with all-sky far-infrared maps from the AKARI satellite. We first test our measurement pipeline using the Sloan Digital Sky Survey (SDSS) redMaPPer catalogue of galaxy clusters, finding a highly significant detection ($>$$20\sigma$) of the SZ effect together with correlated dust emission. We then constrain the SZ signal associated with spectroscopically confirmed quasi-stellar objects (QSOs) from SDSS data release 7 (DR7) and the Baryon Oscillation Spectroscopic Survey (BOSS) DR12. We obtain a low-significance ($1.6\sigma$) hint of an SZ signal, pointing towards a mean thermal energy of $\simeq 5 \times 10^{60}$ erg, lower than reported in some previous studies. A comparison of our results with high-resolution hydrodynamical simulations including AGN feedback suggests QSO host masses of $M_{200c} \sim 4 \times 10^{12}~h^{-1}M_\odot$, but with a large uncertainty. Our analysis provides no conclusive evidence for an SZ signal specifically associated with AGN feedback.

We compare the low redshift (z ~ 0.1) Lyman-alpha forest from hydrodynamical simulations with data from the Cosmic Origin Spectrograph (COS). We find tension between the observed number of lines with b-parameters in the range 25-45 km/s and the predictions from simulations that incorporate either vigorous feedback from active galactic nuclei or that exclude feedback altogether. The gas in these simulations is, respectively, either too hot to contribute to the Lyman-alpha absorption or too cold to produce the required line widths. Matching the observed b-parameter distribution therefore requires feedback processes that thermally or turbulently broaden the absorption features without collisionally (over-)ionising hydrogen. This suggests the Lyman-alpha forest b-parameter distribution is a valulable diagnostic of galactic feedback in the low redshift Universe. We furthermore confirm the low redshift Lyman-alpha forest column density distribution is better reproduced by an ultraviolet background with an HI photo-ionisation rate a factor 1.5-3 higher than predicted by Haardt & Madau (2012).

Lyman-alpha forest data probing the post-reionization Universe shows surprisingly large opacity fluctuations over rather large ($\ge$50 comoving Mpc/h) spatial scales. We model these fluctuations using a hybrid approach utilizing the large volume Millennium simulation to predict the spatial distribution of QSOs combined with smaller scale full hydrodynamical simulation performed with RAMSES and post-processed with the radiative transfer code ATON. We produce realictic mock absorption spectra that account for the contribution of galaxies and QSOs to the ionising UV background. This improved models confirm our earlier findings that a significant ($\ge$50%) contribution of ionising photons from QSOs can explain the large reported opacity fluctuations on large scales. The inferred QSO luminosity function is thereby consistent with recent estimates of the space density of QSOs at this redshift. Our simulations still somewhat struggle, however, to reproduce the very long (110 comoving Mpc/h) high opacity absorption through observed in ULAS J0148+0600, perhaps suggesting an even later end of reionization than assumed in our previously favoured model. Medium-deep/medium area QSO surveys as well as targeted searches for the predicted strong transverse QSO proximity effect whould illuminate the origin of the observed large scale opacity fluctuations. They would allow to substantiate whether UV fluctuations due to QSO are indeed primarily responsible, or whether significant contributions from other recently proposed mechansims such as large scale fluctuations in temperature and mean free path (even in the absence of rare bright sources) are required.

Many-degree-scale gamma-ray halos are expected to surround extragalactic high-energy gamma ray sources. These arise from the inverse Compton emission of an intergalactic population of relativistic electron/positron pairs generated by the annihilation of >100 GeV gamma rays on the extragalactic background light. These are typically anisotropic due to the jetted structure from which they originate or the presence of intergalactic magnetic fields. Here we propose a novel method for detecting these inverse-Compton gamma-ray halos based upon this anisotropic structure. Specifically, we show that by stacking suitably defined angular power spectra instead of images it is possible to robustly detect gamma-ray halos with existing Fermi Large Area Telescope (LAT) observations for a broad class of intergalactic magnetic fields. Importantly, these are largely insensitive to systematic uncertainties within the LAT instrumental response or associated with contaminating astronomical sources.

Modern theories of galaxy formation predict that galaxies impact on their gaseous surroundings, playing the fundamental role of regulating the amount of gas converted into stars. While star-forming galaxies are believed to provide feedback through galactic winds, Quasi-Stellar Objects (QSOs) are believed instead to provide feedback through the heat generated by accretion onto a central supermassive black hole. A quantitative difference in the impact of feedback on the gaseous environments of star-forming galaxies and QSOs has not been established through direct observations. Using the Sherwood cosmological simulations, we demonstrate that measurements of neutral hydrogen in the vicinity of star-forming galaxies and QSOs during the era of peak galaxy formation show excess LyA absorption extending up to comoving radii of about 150 kpc for star-forming galaxies and 300 - 700 kpc for QSOs. Simulations including supernovae-driven winds with the wind velocity scaling like the escape velocity of the halo account for the absorption around star-forming galaxies but not QSOs.

Feedback from central supermassive blackholes is often invoked to explain the low star formation rates in massive galaxies at the centers of galaxy clusters. However, the detailed physics of the coupling of the injected feedback energy with the intracluster medium is still unclear. Using high-resolution magnetohydrodynamic cosmological simulations of galaxy cluster formation, we investigate the role of anisotropic thermal conduction in shaping the thermodynamic structure of clusters, and, in particular, in modifying the impact of black hole feedback. Stratified anisotropically conducting plasmas are formally always unstable, and thus more prone to mixing, an expectation borne out by our results. The increased mixing efficiently isotropizes the injected feedback energy which in turn significantly improves the coupling between the feedback energy and the intracluster medium. This facilitates an earlier disruption of the cool core, reduces the star formation rate by more than an order of magnitude, and results in earlier quenching despite an overall lower amount of feedback energy injected into the cluster core. With conduction, the metallicity gradients and dispersions are lowered, aligning them better with observational constraints. These results highlight the important role of thermal conduction in establishing and maintaining quiescence of massive galaxies.

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{\"u}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.

We examine the properties of the galaxies and dark matter haloes residing in the cluster infall region surrounding the simulated $\Lambda$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.

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.

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%.

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.

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

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.

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.

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.

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.

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.

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}.

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 $\Lambda$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.

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.

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.

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.

We have simulated the formation of a galaxy cluster in a $\Lambda$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.

In this work, we analyze the Lyman-$\alpha$ 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-$\alpha$ 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-$\alpha$ forest is rather small. Even models with strong modifications of gravity, like $|\bar{f}_{R0}| = 10^{-4}$, do not change the statistical Lyman-$\alpha$ 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.

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.

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.

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.

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 $\Lambda$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 $\Lambda$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 $\Sigma _{i}m_{\nu _{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].

The extragalactic gamma-ray sky at TeV energies is dominated by blazars, a subclass of accreting super-massive black holes with powerful relativistic outflows directed at us. Only constituting a small fraction of the total power output of black holes, blazars were thought to have a minor impact on the universe at best. As we argue here, the opposite is true and the gamma-ray emission from TeV blazars can be thermalized via beam-plasma instabilities on cosmological scales with order unity efficiency, resulting in a potentially dramatic heating of the low-density intergalactic medium. Here, we review this novel heating mechanism and explore the consequences for the formation of structure in the universe. In particular, we show how it produces an inverted temperature-density relation of the intergalactic medium that is in agreement with observations of the Lyman-alpha forest. This suggests that blazar heating can potentially explain the paucity of dwarf galaxies in galactic halos and voids, and the bimodality of galaxy clusters. This also transforms our understanding of the evolution of blazars, their contribution to the extra-galactic gamma-ray background, and how their individual spectra can be used in constraining intergalactic magnetic fields.

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.

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.

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.

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.

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]

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.

We present relations between X-ray luminosity and velocity dispersion (L-sigma), X-ray luminosity and gas mass (L-Mgas), and cluster radius and velocity dispersion (r500-sigma) for 62 galaxy clusters in the HIFLUGCS, an X-ray flux-limited sample minimizing bias toward any cluster morphology. Our analysis in total is based on ~1.3Ms of clean X-ray XMM-Newton data and 13439 cluster member galaxies with redshifts. Cool cores are among the major contributors to the scatter in the L-sigma relation. When the cool-core-corrected X-ray luminosity is used the intrinsic scatter decreases to 0.27 dex. Even after the X-ray luminosity is corrected for the cool core, the scatter caused by the presence of cool cores dominates for the low-mass systems. The scatter caused by the non-cool-core clusters does not strongly depend on the mass range, and becomes dominant in the high-mass regime. The observed L-sigma relation agrees with the self-similar prediction, matches that of a simulated sample with AGN feedback disregarding six clusters with <45 cluster members with spectroscopic redshifts, and shows a common trend of increasing scatter toward the low-mass end, i.e., systems with sigma<500km/s. A comparison of observations with simulations indicates an AGN-feedback-driven impact in the low-mass regime. The best fits to the $L-M_{\rm gas}$ relations for the disturbed clusters and undisturbed clusters in the observational sample closely match those of the simulated samples with and without AGN feedback, respectively. This suggests that one main cause of the scatter is AGN activity providing feedback in different phases, e.g., during a feedback cycle. The slope and scatter in the observed r500-sigma relation is similar to that of the simulated sample with AGN feedback except for a small offset but still within the scatter.

We investigate the impact of AGN feedback on the gravitational lensing properties of a sample of galaxy clusters with masses in the range 10^14 - 10^15 Msol, using state-of-the-art simulations. Adopting a ray-tracing algorithm, we compute the cross-section of giant arcs from clusters simulated with dark-matter only physics (DM), dark matter plus gas with cooling and star formation (CSF), and dark matter plus gas with cooling, star formation and AGN feedback (CSFBH). Once AGN feedback is included, baryonic physics boosts the strong lensing cross-section by much less than previously estimated using clusters simulated with only cooling and star formation. For a cluster with a virial mass of 7.4 X 10^14 Msol, inclusion of baryonic physics without feedback can boost the cross-section by as much as a factor of 3, in agreement with previous studies, whereas once AGN feedback is included this maximal figure falls to a factor of 2 at most. Typically, clusters simulated with DM and CSFBH physics have similar cross-sections for the production of giant arcs. We also investigate how baryonic physics affects the weak lensing properties of the simulated clusters by fitting NFW profiles to synthetic weak lensing data sets using a Markov Chain Monte Carlo approach, and by performing non-parametric mass reconstructions. Without the inclusion of AGN feedback, measured concentration parameters can be much larger than those obtained with AGN feedback, which are similar to the dark-matter only case.

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.

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.

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.

Recently, rapid observational and theoretical progress has established that black holes (BHs) play a decisive role in the formation and evolution of individual galaxies as well as galaxy groups and clusters. In particular, there is compelling evidence that BHs vigorously interact with their surroundings in the central regions of galaxy clusters, indicating that any realistic model of cluster formation needs to account for these processes. This is also suggested by the failure of previous generations of hydrodynamical simulations without BH physics to simultaneously account for the paucity of strong cooling flows in clusters, the slope and amplitude of the observed cluster scaling relations, and the high-luminosity cut-off of central cluster galaxies. 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. By comparing a large sample of simulated AGN-heated clusters with observations, our new simulation technique should make it possible to reliably calibrate observational biases in cluster surveys, thereby enabling various high-precision cosmological studies of the dark matter and dark energy content of the universe.

We show how hydrostatic equilibrium in galaxy clusters can be quantitatively probed combining X-ray, SZ, and gravitational-lensing data. Our previously published method for recovering three-dimensional cluster gas distributions avoids the assumption of hydrostatic equilibrium. Independent reconstructions of cumulative total-mass profiles can then be obtained from the gas distribution, assuming hydrostatic equilibrium, and from gravitational lensing, neglecting it. Hydrostatic equilibrium can then be quantified comparing the two. We describe this procedure in detail and show that it performs well on progressively realistic synthetic data. An application to a cluster merger demonstrates how hydrostatic equilibrium is violated and restored as the merger proceeds.

We propose and test a new method based on Richardson-Lucy deconvolution to reconstruct three-dimensional gas density and temperature distributions in galaxy clusters from combined X-ray and thermal Sunyaev-Zel'dovich observations. Clusters are assumed to be axially symmetric and arbitrarily inclined with respect to the line-of-sight. No equilibrium assumption other than local thermal equilibrium is needed. We test the algorithm with synthetic observations of analytically modeled and numerically simulated galaxy clusters and discuss the quality of the density and temperature reconstructions in idealised situations and in presence of observational noise, deviations from axial symmetry and cluster substructure. We find that analytic and numerical gas density and temperature distributions can be accurately reconstructed in three dimensions, even if observational noise is present. We also discuss methods for determining the inclination angle from data and show that it can be constrained using X-ray temperature maps. For a realistic cluster and including observational noise the three-dimensional reconstructions reach a level of accuracy of about 15%.

Previous studies of strong gravitational lensing by galaxy clusters neglected the potential impact of the intracluster gas. Here, we compare simulations of strong cluster lensing including gas physics at increasing levels of complexity, i.e. with adiabatic, cooling, star-forming, feedback-receiving, and thermally conducting gas, and with different implementations of the artificial viscosity in the SPH simulations. Each cluster was simulated starting from the same initial conditions such as to allow directly comparing the simulated clusters. We compare the clusters' shapes, dynamics and density profiles and study their strong-lensing cross sections computed by means of ray-tracing simulations. We find that the impact of adiabatic gas depends on the amount of turbulence that builds up, which means that the artificial viscosity plays an important role. With the common viscosity implementation, adiabatic gas has little effect on strong cluster lensing, while lower viscosity allows stronger turbulence, thus higher non-thermal pressure and a generally broader gas distribution which tends to lower lensing cross sections. Conversely, cooling and star formation steepen the core density profiles and can thus increase the strong-lensing efficiency considerably.

We calculate the one-loop vacuum polarization tensor of scalar electrodynamics in a locally de Sitter space-time, endowed with a nearly minimally coupled, light scalar field. We show that the photon dynamics is well approximated by a (local) Proca Lagrangean. Since the photon mass can be much larger than the Hubble parameter, the photons may propagate slowly during inflation. Finally, we briefly discuss magnetic field generation on cosmological scales, and point out that, while the spectrum of the magnetic field is identical to that obtained from the massless scalar, B_l ~ B_0/l, the amplitude B_0 may be significantly enhanced, implying that the seed field bound for the galactic dynamo can be easily met.

We propose a new mechanism for magnetic field generation from inflation, by which strong magnetic fields can be generated on cosmological scales. These fields may be observable by cosmic microwave background radition measurements, and may have a dynamical impact on structure formation. The mechanism is based on the observation that a light nearly minimally coupled charged scalar may be responsible for the creation of a negative photon mass-squared (provided the scalar field coupling to the curvature scalar is small but negative), which in turn results in abundant photon production - and thus in growing magnetic fields - during inflation.