Supernova (SN) feedback plays a vital role in the evolution of galaxies. While modern cosmological simulations capture the leading structures within galaxies, they struggle to provide sufficient resolution to study small-scale stellar feedback, such as the detailed evolution of SN remnants. It is thus common practice to assume subgrid models that are rarely extended to low metallicities, and which routinely use the standard solar abundance pattern. With the aid of 1-d hydrodynamical simulations, we extend these models to consider low metallicities and non-solar abundance patterns as derived from spectra of Milky Way stars. For that purpose, a simple, yet effective framework has been developed to generate non-solar abundance pattern cooling functions. We find that previous treatments markedly over-predict SN feedback at low metallicities and show that non-negligible changes in the evolution of SN remnants of up to $\approx 50\%$ in $cooling\; mass$ and $\approx 27\%$ in $momentum\; injection\; from\; SN\; remnants$ arise from non-solar abundance patterns. We use our simulations to quantify these results as a function of metallicity and abundance pattern variations and present analytic formulae to accurately describe the trends. These formulae have been designed to serve as subgrid models for SN feedback in cosmological hydrodynamical simulations.

Dark matter dominates the properties of large cosmological structures such as galaxy clusters, and the mass profiles of the dark matter have been measured for these equilibrated structures for years using X-rays, lensing or galaxy velocities. A new method has been proposed, which should allow us to estimate a dynamical property of the dark matter, namely the velocity anisotropy. For the gas a similar velocity anisotropy is zero due to frequent collisions, however, the collisionless nature of dark matter allows it to be non-trivial. Numerical simulations have for years found non-zero and radially varying dark matter velocity anisotropies. Here we employ the method proposed by Hansen and Pifaretti (2007), and developed by Host et al. (2009) to estimate the dark matter velocity anisotropy in the bright galaxy cluster Perseus, to near 5 times the radii previously obtained. We find the dark matter velocity anisotropy to be consistent with the results of numerical simulations, however, still with large error-bars. At half the virial radius we find the velocity anisotropy to be non-zero at 1.7 standard deviations, lending support to the collisionless nature of dark matter.

A suite of idealised, global, gravitationally unstable, star-forming galactic disc simulations with 2 pc spatial resolution, performed with the adaptive mesh refinement code {\sc ramses} is used in this paper to predict the emergent effects of supernova feedback. The simulations include a simplified prescriptions for formation of single stellar populations of mass $\sim 100 \, M_{\odot}$, radiative cooling, photoelectric heating, an external gravitational potential for a dark matter halo and an old stellar disc, self-gravity, and a novel implementation of supernova feedback. The results of these simulations show that gravitationally unstable discs can generate violent supersonic winds with mass loading factors $\eta \gtrsim 10$, followed by a galactic fountain phase. These violent winds are generated by highly clustered supernovae exploding in dense environments created by gravitational instability, and they are not produced in simulation without self-gravity. The violent winds significantly perturb the vertical structure of the disc, which is later re-established during the galactic fountain phase. Gas resettles into a quasi-steady, highly turbulent disc with volume-weighted velocity dispersion $\sigma > 50 \, {\rm km/s}$. The new configuration drives weaker galactic winds with mass loading factor $\eta \leq 0.1$. The whole cycle takes place in $\leq 10$ dynamical times. Such high time variability needs to be taken into account when interpreting observations of galactic winds from starburst and post-starburst galaxies.

The connections among galaxies, the dark matter halos where they form and the properties of the large-scale Cosmic Web still need to be completely disentangled. We use the cosmological hydrodynamical simulation TNG100 of the IllustrisTNG suite to quantify the effects played by the large-scale density field and the Cosmic Web morphology on the relation between halo mass and galaxy stellar mass. We select objects with total dynamical mass in the range $\geq 6.3\times 10^{10} h ^{-1}\, M_{\odot}$ up to a few $10^{14} h^{-1} \, M_{\odot}$ between redshift $z=4$ and redshift $z=0$. A Cosmic Web class (knot, filament, sheet, void) is assigned to each region of the volume using a density field deformation tensor-based method. We find that galaxy stellar mass strongly correlates with total dynamical mass and formation time, and more weakly with large-scale overdensity and Cosmic Web class. The latter two quantities correlate with each other, but are not entirely degenerate. Furthermore, we find that at fixed halo mass, galaxies with stellar mass lower than the median value are more likely to be found in voids and sheets, whereas galaxies with stellar mass higher than the median are more likely to be found in filaments and knots. Finally, we find that the dependence on environment is stronger for satellites than for centrals, and discuss the physical implications of these results.

The recently discovered discrepancy between galaxy mass measurements from weak lensing and predictions from abundance matching questions our understanding of cosmology, or of the galaxy-halo connection, or of both. We re-examined this tension by considering, as models, different cosmological simulations in the Illustris suite. We produced excess profiles $R\Delta\Sigma$ from subhalo snapshots at different redshifts in Illustris-1 and IllustrisTNG (TNG100 and TNG300) simulations, enabling a direct comparison with weak-lensing measurements. We separate the individual contributions of stars, dark matter and gas within $\approx1$ Mpc (comoving length), beyond which correlated two-halo terms dominate. The mismatch between measurements and predictions is more severe than in previous studies: $R\Delta\Sigma$ profiles from IllustrisTNG are $\approx2$ times higher than the measured ones. Contrary to abundance matching results, the mismatch is mostly unchanged with increasing redshifts. The contribution of gas to the $R\Delta\Sigma$ profiles is $5-10\%$ over the scales dominated by one-halo terms. Different procedures to link stellar and halo masses (abundance matching, cosmological simulations) are still significantly discrepant with weak lensing measurements, but their trends are different. Therefore, the change in cosmological parameters advocated through abundance-matching arguments may not resolve this tension. Also, current criteria to select isolated massive galaxies in simulations are susceptible to resolution issues and may not correspond to observational criteria. The (currently subdominant) contribution of gas is non-negligible, and even if the major discrepancy within stellar and halo masses is resolved, it will be an appreciable source of systematics in the LSST era, when uncertainties on the $R\Delta\Sigma$ profiles are expected to be $\approx10$ times smaller.

We analyze the IllustrisTNG simulations to study the mass, volume fraction and phase distribution of gaseous baryons embedded in the knots, filaments, sheets and voids of the Cosmic Web from redshift $z=8$ to redshift $z=0$. We find that filaments host more star-forming gas than knots, and that filaments also have a higher relative mass fraction of gas in this phase than knots. We also show that the cool, diffuse Intergalactic Medium (IGM; $T<10^5 \, {\rm K}$, $ n_{\rm H}<10^{-4}(1+z) \, {\rm cm^{-3}}$) and the Warm-Hot Intergalactic Medium (WHIM; $ 10^5 \, {\rm K} <T<10^7 \, {\rm K}$, $ n_{\rm H} <10^{-4}(1+z)\, {\rm cm^{-3}}$) constitute $\sim 39\%$ and $\sim 46\%$ of the baryons at redshift $z=0$, respectively. Our results indicate that the WHIM may constitute the largest reservoir of {\it missing} baryons at redshift $z=0$. Using our Cosmic Web classification, we predict the WHIM to be the dominant baryon mass contribution in filaments and knots at redshift $z=0$, but not in sheets and voids where the cool, diffuse IGM dominates. We also characterise the evolution of WHIM and IGM from redshift $z=4$ to redshift $z=0$, and find that the mass fraction of WHIM in filaments and knots evolves only by a factor $\sim 2$ from redshift $z=0$ to $z=1$, but declines faster at higher redshift. The WHIM only occupies $4-11\%$ of the volume at redshift $0\leq z \leq 1$. We predict the existence of a significant number of currently undetected OVII and NeIX absorption systems in cosmic filaments which could be detected by future X-ray telescopes like Athena.

We study how jets driven by active galactic nuclei influence the cooling flow in Perseus-like galaxy cluster cores with idealised, non-relativistic, hydrodynamical simulations performed with the Eulerian code ATHENA using high-resolution Godunov methods with low numerical diffusion. We use novel analysis methods to measure the cooling rate, the heating rate associated to multiple mechanisms, and the power associated with adiabatic compression/expansion. A significant reduction of the cooling rate and cooling flow within 20 kpc from the centre can be achieved with kinetic jets. However, at larger scales and away from the jet axis, the system relaxes to a cooling flow configuration. Jet feedback is anisotropic and is mostly distributed along the jet axis, where the cooling rate is reduced and a significant fraction of the jet power is converted into kinetic power of heated outflowing gas. Away from the jet axis weak shock heating represents the dominant heating source. Turbulent heating is significant only near the cluster centre, but it becomes inefficient at 50 kpc scales where it only represents a few percent of the total heating rate. Several details of the simulations depend on the choice made for the hydro solver, a consequence of the difficulty of achieving proper numerical convergence for this problem: current physics implementations and resolutions do not properly capture multi-phase gas that develops as a consequence of thermal instability. These processes happen at the grid scale and leave numerical solutions sensitive to the properties of the chosen hydro solver.

We use three-dimensional hydrodynamic simulations of vertically stratified patches of galactic discs to study how the spatio-temporal clustering of supernovae (SNe) enhances the power of galactic winds. SNe that are randomly distributed throughout a galactic disc drive inefficient galactic winds because most supernova remnants lose their energy radiatively before breaking out of the disc. Accounting for the fact that most star formation is clustered alleviates this problem. Super-bubbles driven by the combined effects of clustered SNe propagate rapidly enough to break out of galactic discs well before the clusters' SNe stop going off. The radiative losses post-breakout are reduced dramatically and a large fraction ($\gtrsim 0.2$) of the energy released by SNe vents into the halo powering a strong galactic wind. These energetic winds are capable of providing strong preventative feedback and eject substantial mass from the galaxy with outflow rates on the order of the star formation rate. The momentum flux in the wind is only of order that injected by the SNe, because the hot gas vents before doing significant work on the surroundings. We show that our conclusions hold for a range of galaxy properties, both in the local Universe (e.g., M82) and at high redshift (e.g., $z \sim 2$ star forming galaxies). We further show that if the efficiency of forming star clusters increases with increasing gas surface density, as suggested by theoretical arguments, the condition for star cluster-driven super-bubbles to break out of galactic discs corresponds to a threshold star formation rate surface density for the onset of galactic winds $\sim 0.03$ M$_\odot$ yr$^{-1}$ kpc$^{-2}$, of order that observed.

We use the IllustrisTNG simulations to investigate the evolution of the mass-metallicity relation (MZR) for star-forming cluster galaxies as a function of the formation history of their cluster host. The simulations predict an enhancement in the gas-phase metallicities of star-forming cluster galaxies (10^9< M_star<10^10 M_sun) at z<1.0 in comparisons to field galaxies. This is qualitatively consistent with observations. We find that the metallicity enhancement of cluster galaxies appears prior to their infall into the central cluster potential, indicating for the first time a systematic "chemical pre-processing" signature for {\it infalling} cluster galaxies. Namely, galaxies which will fall into a cluster by z=0 show a ~0.05 dex enhancement in the MZR compared to field galaxies at z<0.5. Based on the inflow rate of gas into cluster galaxies and its metallicity, we identify that the accretion of pre-enriched gas is the key driver of the chemical evolution of such galaxies, particularly in the stellar mass range (10^9< M_star<10^10 M_sun). We see signatures of an environmental dependence of the ambient/inflowing gas metallicity which extends well outside the nominal virial radius of clusters. Our results motivate future observations looking for pre-enrichment signatures in dense environments.

Sub-structures in the hot gas of galaxy clusters are related to their formation history and to the astrophysical processes at play in the intracluster medium (ICM). The thermal Sunyaev-Zel'dovich (tSZ) effect is directly sensitive to the line-of-sight integrated ICM pressure, and is thus particularly adapted to study ICM sub-structures. We apply structure-enhancement filtering algorithms to high resolution tSZ observations of distant clusters, in order to search for pressure discontinuities, compressions, as well as secondary peaks in the ICM. The same filters are applied to synthetic tSZ images extracted from RHAPSODY-G hydrodynamic simulations, in order to better interpret the extracted features. We also study the noise propagation trough the filters and quantify the impact of systematic effects, point source residuals being identified as the dominant potential contaminant. In 3 of our 6 NIKA clusters we identify features at high S/N that show clear evidence for merger events. In MACSJ0717 (z=0.55), three strong pressure gradients are observed on the E, SE and W sectors, and two main peaks in the pressure distribution are identified. We observe a lack of tSZ compact structure in the cool-core cluster PSZ1G045.85 (z=0.61), and a tSZ gradient ridge dominates in the SE. In the highest z cluster, CLJ1227 (z=0.89), we detect a ~45" (360 kpc) long ridge pressure gradient associated with a secondary pressure peak in the W region. Our results show that current tSZ facilities have now reached the angular resolution and sensitivity to allow an exploration of the details of pressure sub-structures in clusters, even at high z. This opens the possibility to quantify the impact of the dynamical state on the relation between the tSZ signal and the mass of clusters, which is important when using tSZ clusters to test cosmological models. This work also marks the first NIKA cluster sample data release.

We consider the dynamics in and near galaxy clusters. Gas, dark matter and galaxies are presently falling into the clusters between approximately 1 and 5 virial radii. At very large distances, beyond 10 virial radii, all matter is following the Hubble flow, and inside the virial radius the matter particles have on average zero radial velocity. The cosmological parameters are imprinted on the infall profile of the gas, however, no method exists, which allows a measurement of it. We consider the results of two cosmological simulations (using the numerical codes RAMSES and Gadget) and find that the gas and dark matter radial velocities are very similar. We derive the relevant dynamical equations, in particular the generalized hydrostatic equilibrium equation, including both the expansion of the Universe and the cosmological background. This generalized gas equation is the main new contribution of this paper. We combine these generalized equations with the results of the numerical simulations to estimate the contribution to the measured cluster masses from the radial velocity: inside the virial radius it is negligible, and inside two virial radii the effect is below 40%, in agreement the earlier analyses for DM. We point out how the infall velocity in principle may be observable, by measuring the gas properties to distance of about two virial radii, however, this is practically not possible today.

Recent observations of galaxies in a cluster at z=0.35 show that their integrated gas-phase metallicities increase with decreasing cluster-centric distance. To test if ram pressure stripping (RPS) is the underlying cause, we use a semi-analytic model to quantify the "observational bias" that RPS introduces into the aperture-based metallicity measurements. We take integral field spectroscopy of local galaxies, remove gas from their outer galactic disks via RPS, and then conduct mock slit observations of cluster galaxies at z=0.35. Our RPS model predicts a typical cluster-scale metallicity gradient of -0.03 dex/Mpc. By removing gas from the outer galactic disks, RPS introduces a mean metallicity enhancement of +0.02 dex at a fixed stellar mass. This gas removal and subsequent quenching of star formation preferentially removes low mass cluster galaxies from the observed star-forming population. As only the more massive star-forming galaxies survive to reach the cluster core, RPS produces a cluster-scale stellar mass gradient of -0.05 log(M_*/M_sun)/Mpc. This mass segregation drives the predicted cluster-scale metallicity gradient of -0.03 dex/Mpc. However, the effects of RPS alone can not explain the higher metallicities measured in cluster galaxies at z=0.35. We hypothesize that additional mechanisms including steep internal metallicity gradients and self-enrichment due to gas strangulation are needed to reproduce our observations at z=0.35.

We use idealized three-dimensional hydrodynamic simulations of global galactic discs to study the launching of galactic winds by supernovae (SNe). The simulations resolve the cooling radii of the majority of supernova remnants (SNRs) and thus self-consistently capture how SNe drive galactic winds. We find that SNe launch highly supersonic winds with properties that agree reasonably well with expectations from analytic models. The energy loading ($\eta_E = \dot{E}_{\rm wind} / \dot{E}_{\rm SN}$) of the winds in our simulations are well converged with spatial resolution while the wind mass loading ($\eta_M = \dot{M}_{\rm wind} / \dot{M}_\star$) decreases with resolution at the resolutions we achieve. We present a simple analytic model based on the concept that SNRs with cooling radii greater than the local scale height breakout of the disc and power the wind. This model successfully explains the dependence (or lack thereof) of $\eta_E$ (and by extension $\eta_M$) on the gas surface density, star formation efficiency, disc radius, and the clustering of SNe. The winds in the majority of our simulations are weaker than expected in reality, likely due to the fact that we seed SNe preferentially at density peaks. Clustering SNe in time and space substantially increases the wind power.

We compare galaxy scaling relations as a function of environment at $z\sim2$ with our ZFIRE survey where we have measured H$\alpha$ fluxes for 90 star-forming galaxies selected from a mass-limited [$\log(M_{\star}/M_{\odot})>9$] sample based on ZFOURGE. The cluster galaxies (37) are part of a confirmed system at z=2.095 and the field galaxies (53) are at $1.9<z<2.4$; all are in the COSMOS legacy field. There is no statistical difference between H$\alpha$-emitting cluster and field populations when comparing their star formation rate (SFR), stellar mass ($M_{\star}$), galaxy size ($r_{eff}$), SFR surface density [$\Sigma$(H$\alpha_{star}$)], and stellar age distributions. The only difference is that at fixed stellar mass, the H$\alpha$-emitting cluster galaxies are $\log(r_{eff})\sim0.1$ larger than in the field. Approximately 19% of the H$\alpha$-emitters in the cluster and 26% in the field are IR-luminous ($L_{IR}>2\times10^{11} L_{\odot}$). Because the LIRGs in our combined sample are $\sim5$ times more massive than the low-IR galaxies, their radii are $\sim70$% larger. To track stellar growth, we separate galaxies into those that lie above, on, and below the H$\alpha$ star-forming main sequence (SFMS) using $\Delta$SFR$(M_{\star})=\pm0.2$ dex. Galaxies above the SFMS (starbursts) tend to have higher H$\alpha$ SFR surface densities and younger light-weighted stellar ages compared to galaxies below the SFMS. Our results indicate that starbursts (+SFMS) in the cluster and field at $z\sim2$ are growing their stellar cores. Lastly, we compare to the (SFR-$M_{\star}$) relation from RHAPSODY cluster simulations and find the predicted slope is nominally consistent with the observations. However, the predicted cluster SFRs tend to be too low by a factor of $\sim2$ which seems to be a common problem for simulations across environment.

Galaxy cluster mass determinations achieved using X-ray and Sunyaev-Zel'dovich data combined with the assumption of hydrostatic equilibrium are generally biased. The bias exists for two main reasons: non-thermal pressure forces are expected to contribute to the overall pressure balance and deviations from spherical symmetry and hydrostatic equilibrium can be present. In this paper, we use a sample of zoom-in hydrodynamical simulations of galaxy clusters to measure the magnitude of hydrostatic bias and the non-thermal contribution to the total pressure. We propose a new empirical model for non-thermal pressure based on our simulations that can be applied to observations. We show that our model can be successfully applied to remove most of the bias related to neglection of non-thermal pressure, which is usually not included in hydrostatic cluster mass profile reconstructions. The use of this model may significantly improve the calibration of cluster scaling relations that are a key tool for cluster cosmology.

We present the first observation of cluster-scale radial metallicity gradients from star-forming galaxies. We use the DEIMOS spectrograph on the Keck II telescope to observe two CLASH clusters at z~0.35: MACS1115+0129 and RXJ1532+3021. Based on our measured interstellar medium (ISM) properties of star-forming galaxies out to a radius of 2.5 Mpc from the cluster centre, we find that the galaxy metallicity decreases as a function of projected cluster-centric distance (-0.15+/-0.08 dex/Mpc) in MACS1115+01. On the mass-metallicity relation (MZR), star-forming galaxies in MACS1115+01 are offset to higher metallicity (~0.2 dex) than the local SDSS galaxies at a fixed mass range. In contrast, the MZR of RXJ1532+30 is consistent with the local comparison sample. RXJ1532+30 exhibits a bimodal radial metallicity distribution, with one branch showing a similar negative gradient as MACS1115+01 (-0.14+/-0.05 dex/Mpc) and the other branch showing a positive radial gradient. The positive gradient branch in RXJ1532+30 is likely caused by either interloper galaxies or an in-plane merger, indicating that cluster-scale abundance gradients probe cluster substructures and thus the dynamical state of a cluster. Most strikingly, we discover that neither the radial metallicity gradient nor the offset from the MZR is driven by the stellar mass. We compare our observations with Rhapsody-G cosmological hydrodynamical zoom-in simulations of relaxed galaxy clusters and find that the simulated galaxy cluster also exhibits a negative abundance gradient, albeit with a shallower slope (-0.04+/-0.03 dex/Mpc). Our observations suggest that the negative radial gradient originates from ram-pressure stripping and/or strangulation processes in the cluster environments.

We use the first release of the VImos Public Extragalactic Redshift Survey of galaxies (VIPERS) of ~50,000 objects to measure the biasing relation between galaxies and mass in the redshift range z=[0.5,1.1]. We estimate the 1-point distribution function [PDF] of VIPERS galaxies from counts in cells and, assuming a model for the mass PDF, we infer their mean bias relation. The reconstruction of the bias relation from PDFs is performed through a novel method that accounts for Poisson noise, redshift distortions, inhomogeneous sky coverage and other selection effects. With this procedure we constrain galaxy bias and its deviations from linearity down to scales as small as 4 Mpc/h and out to z=1.1. We detect small (~3%) but significant deviations from linear bias. The mean biasing function is close to linear in regions above the mean density. The mean slope of the biasing relation is a proxy to the linear bias parameter. It increases both with luminosity, in agreement with results of previous analyses, and with redshift. However, we detect a strong bias evolution only for z>0.9 in agreement with some, but not all, previous studies. We also detected a significant increase of the bias with the scale, from 4 to 8 Mpc/h, now seen for the first time out to z=1. The amplitude of nonlinearity depends on redshift, luminosity and on scales but no clear trend is detected. Thanks to the large cosmic volume probed by VIPERS we find that the mismatch between the previous estimates of bias at z~1 from zCOSMOS and VVDS-Deep galaxy samples is fully accounted for by cosmic variance. The results of our work confirm the importance of going beyond the over-simplistic linear bias hypothesis showing that non-linearities can be accurately measured through the applications of the appropriate statistical tools to existing datasets like VIPERS.

We study the evolution of the stellar component and the metallicity of both the intracluster medium and of stars in massive ($M_{\rm vir}\approx 6\times 10^{14}$ M$_{\odot}/h$) simulated galaxy clusters from the Rhapsody-G suite in detail and compare them to observational results. The simulations were performed with the AMR code RAMSES and include the effect of AGN feedback at the sub-grid level. AGN feedback is required to produce realistic galaxy and cluster properties and plays a role in mixing material in the central regions and regulating star formation in the central galaxy. In both our low and high resolution runs with fiducial stellar yields, we find that stellar and ICM metallicities are a factor of two lower than in observations. We find that cool core clusters exhibit steeper metallicity gradients than non-cool core clusters, in qualitative agreement with observations. We verify that the ICM metallicities measured in the simulation can be explained by a simple "regulator" model in which the metallicity is set by a balance of stellar yield and gas accretion. It is plausible that a combination of higher resolution and higher metal yield in AMR simulation would allow the metallicity of simulated clusters to match observed values; however this hypothesis needs to be tested with future simulations. Comparison to recent literature highlights that results concerning the metallicity of clusters and cluster galaxies might depend sensitively on the scheme chosen to solve the hydrodynamics.

We use local Cartesian simulations with a vertical gravitational potential to study how supernova (SN) feedback in stratified galactic discs drives turbulence and launches galactic winds. Our analysis includes three disc models with gas surface densities ranging from Milky Way-like galaxies to gas-rich ultra-luminous infrared galaxies (ULIRGs), and two different SN driving schemes (random and correlated with local gas density). In order to isolate the physics of SN feedback, we do not include additional feedback processes. We find that, in these local box calculations, SN feedback excites relatively low mass-weighted gas turbulent velocity dispersions ~3-7 km/s and low wind mass loading factors < 1 in all the cases we study. The low turbulent velocities and wind mass loading factors predicted by our local box calculations are significantly below those suggested by observations of gas-rich and rapidly star-forming galaxies; they are also in tension with global simulations of disc galaxies regulated by stellar feedback. Using a combination of numerical tests and analytic arguments, we argue that local Cartesian boxes cannot predict the properties of galactic winds because they do not capture the correct global geometry and gravitational potential of galaxies. The wind mass loading factors are in fact not well-defined in local simulations because they decline significantly with increasing box height. More physically realistic calculations (e.g., including a global galactic potential and disc rotation) will likely be needed to fully understand disc turbulence and galactic outflows, even for the idealized case of feedback by SNe alone.

Supernova (SN) feedback is one of the key processes shaping the interstellar medium (ISM) of galaxies. SNe contribute to (and in some cases may dominate) driving turbulence in the ISM and accelerating galactic winds. Modern cosmological simulations have sufficient resolution to capture the main structures in the ISM of galaxies, but are typically still not capable of explicitly resolving all of the small-scale stellar feedback processes, including the expansion of supernova remnants (SNRs). We perform a series of controlled three-dimensional hydrodynamic (adaptive mesh refinement) simulations of single SNRs expanding in an inhomogeneous density field with statistics motivated by those of the turbulent ISM. We use these to quantify the momentum and thermal energy injection from SNe as a function of spatial scale and the density, metallicity, and structure of the ambient medium. We develop a series of analytic formulae that we fit to the simulations. These formulae can be used as a basis for a more predictive sub-resolution model for SN feedback for galaxy formation simulations. We then use simulations of multiple, stochastically driven SNe that resolve the key phases of SNRs to test the sub-resolution model, and show that it accurately captures the turbulent kinetic energy and thermal energy in the ISM. By contrast, proposed SN feedback models in the literature based on `delayed cooling' significantly overpredict the late-time thermal energy and momentum in SNRs.

We present the Rhapsody-G suite of cosmological hydrodynamic AMR zoom simulations of ten massive galaxy clusters at the $M_{\rm vir}\sim10^{15}\,{\rm M}_\odot$ scale. These simulations include cooling and sub-resolution models for star formation and stellar and supermassive black hole feedback. The sample is selected to capture the whole gamut of assembly histories that produce clusters of similar final mass. We present an overview of the successes and shortcomings of such simulations in reproducing both the stellar properties of galaxies as well as properties of the hot plasma in clusters. In our simulations, a long-lived cool-core/non-cool core dichotomy arises naturally, and the emergence of non-cool cores is related to low angular momentum major mergers. Nevertheless, the cool-core clusters exhibit a low central entropy compared to observations, which cannot be alleviated by thermal AGN feedback. For cluster scaling relations we find that the simulations match well the $M_{500}-Y_{500}$ scaling of Planck SZ clusters but deviate somewhat from the observed X-ray luminosity and temperature scaling relations in the sense of being slightly too bright and too cool at fixed mass, respectively. Stars are produced at an efficiency consistent with abundance matching constraints and central galaxies have star formation rates consistent with recent observations. While our simulations thus match various key properties remarkably well, we conclude that the shortcomings strongly suggest an important role for non-thermal processes (through feedback or otherwise) or thermal conduction in shaping the intra-cluster medium.

Within a sufficiently large cosmic volume, conservation of baryons implies a simple `closed box' view in which the sum of the baryonic components must equal a constant fraction of the total enclosed mass. We present evidence from Rhapsody-G hydrodynamic simulations of massive galaxy clusters that the closed-box expectation may hold to a surprising degree within the interior, non-linear regions of haloes. At a fixed halo mass, we find a significant anti-correlation between hot gas mass fraction and galaxy mass fraction (cold gas + stars), with a rank correlation coefficient of -0.69 within $R_{500c}$. Because of this anti-correlation, the total baryon mass serves as a low-scatter proxy for total cluster mass. The fractional scatter of total baryon fraction scales approximately as $0.02 (\Delta_c/100)^{0.6}$, while the scatter of either gas mass or stellar mass is larger in magnitude and declines more slowly with increasing radius. We discuss potential observational tests using cluster samples selected by optical and hot gas properties; the simulations suggest that joint selection on stellar and hot gas has potential to achieve 5% scatter in total halo mass.

We introduce PHEW (Parallel HiErarchical Watershed), a new segmentation algorithm to detect structures in astrophysical fluid simulations, and its implementation into the adaptive mesh refinement (AMR) code ramses. PHEW works on the density field defined on the adaptive mesh, and can thus be used on the gas density or the dark matter density after a projection of the particles onto the grid. The algorithm is based on a "watershed" segmentation of the computational volume into dense regions, followed by a merging of the segmented patches based on the saddle point topology of the density field. PHEW is capable of automatically detecting connected regions above the adopted density threshold, as well as the entire set of substructures within. Our algorithm is fully parallel and uses the MPI library. We describe in great detail the parallel algorithm and perform a scaling experiment which proves the capability of phew to run efficiently on massively parallel systems.

We develop an extension of \textit{the Halo Model} that describes analytically the corrections to the matter power spectrum due to the physics of baryons. We extend these corrections to the weak-lensing shear angular power spectrum. Within each halo, our baryonic model accounts for: 1) a central galaxy, the major stellar component whose properties are derived from abundance matching techniques; 2) a hot plasma in hydrostatic equilibrium and 3) an adiabatically-contracted dark matter component. This analytic approach allows us to compare our model to the dark-matter-only case. Our basic assumptions are tested against the hydrodynamical simulations of Martizzi et. al. (2014), with which a remarkable agreement is found. Our baryonic model has only one free parameter, $M_{\rm crit}$, the critical halo mass that marks the transition between feedback-dominated halos, mostly devoid of gas, and gas rich halos, in which AGN feedback effects become weaker. We explore the entire cosmological parameter space, using the angular power spectrum in three redshift bins as the observable, assuming a Euclid-like survey. We derive the corresponding constraints on the cosmological parameters, as well as the possible bias introduced by neglecting the effects of baryonic physics. We find that, up to $\ell_{max}$=4000, baryonic physics plays very little role in the cosmological parameters estimation. However, if one goes up to $\ell_{max}$=8000, the marginalized errors on the cosmological parameters can be significantly reduced, but neglecting baryonic physics can lead to bias in the recovered cosmological parameters up to 10$\sigma$. These biases are removed if one takes into account the main baryonic parameter, $M_{\rm crit}$, which can also be determined up to 1-2\%, along with the other cosmological parameters.

A large sample of cosmological hydrodynamical zoom-in simulations with Adaptive Mesh Refinement (AMR) is analysed to study the properties of simulated Brightest Cluster Galaxies (BCGs). Following the formation and evolution of BCGs requires modeling an entire galaxy cluster, because the BCG properties are largely influenced by the state of the gas in the cluster and by interactions and mergers with satellites. BCG evolution is also deeply influenced by the presence of gas heating sources such as Active Galactic Nuclei (AGNs) that prevent catastrophic cooling of large amounts of gas. We show that AGN feedback is one of the most important mechanisms in shaping the properties of BCGs at low redshift by analysing our statistical sample of simulations with and without AGN feedback. When AGN feedback is included BCG masses, sizes, star formation rates and kinematic properties are closer to those of the observed systems. Some small discrepancies are observed only for the most massive BCGs and in the fraction of star-forming BCGs, effects that might be due to physical processes that are not included in our model.

We report the discovery of an object near M87 in the Virgo Cluster with an extraordinary blueshift of -1025 km/s, offset from the systemic velocity by >2300 km/s. Evaluation of photometric and spectroscopic data provides strong evidence that this object is a distant massive globular cluster, which we call HVGC-1 in analogy to Galactic hypervelocity stars. We consider but disfavor more exotic interpretations, such as a system of stars bound to a recoiling black hole. The odds of observing an outlier as extreme as HVGC-1 in a virialized distribution of intracluster objects are small; it appears more likely that the cluster was (or is being) ejected from Virgo following a three-body interaction. The nature of the interaction is unclear, and could involve either a subhalo or a binary supermassive black hole at the center of M87.

We study the effect of baryonic processes on the halo mass function in the galaxy cluster mass range using a catalogue of 153 high resolution cosmological hydrodynamical simulations performed with the AMR code ramses. We use the results of our simulations within a simple analytical model to gauge the effects of baryon physics on the halo mass function. Neglect of AGN feedback leads to a significant boost in the cluster mass function similar to that reported by other authors. However, including AGN feedback not only gives rise to systems that are similar to observed galaxy clusters, but they also reverse the global baryonic effects on the clusters. The resulting mass function is closer to the unmodified dark matter halo mass function but still contains a mass dependent bias at the 5-10% level. These effects bias measurements of the cosmological parameters, such as $\sigma_8$ and $\Omega_m$. For current cluster surveys baryonic effects are within the noise for current survey volumes, but forthcoming and planned large SZ, X-ray and multi-wavelength surveys will be biased at the percent level by these processes. The predictions for the halo mass function including baryonic effects need to be carefully studied with larger and improved simulations. However, simulations of full cosmological boxes with the resolution we achieve and including AGN feedback are still computationally challenging.

In a recent study (Martizzi et al. 2012), we used cosmological simulations to show that active galactic nuclei (AGN) feedback on the gas distribution in clusters of galaxies can be important in determining the spatial distribution of stars and dark matter in the central regions of these systems. The hierarchical assembly of dark matter, baryons and black holes obscures the physical mechanism behind the restructuring process. Here we use idealized simulations to follow the response of a massive dark matter halo as we feed the central black hole with a controlled supply of cold gas. This removes most of the complexity taking place in the cosmological simulations that may have biased our previous study. We confirm our previous results: gas heated and expelled from the central regions of the halo by AGN feedback can return after cooling; repeated cycles generate gravitational potential fluctuations responsible for irreversible modifications of the dark matter mass profile. The main result is the expulsion of large amounts of baryons and dark matter from the central regions of the halo. According to the work presented here, outflow induced fluctuations represent the only mechanism able to efficiently create dark matter cores in clusters of galaxies.

The spatial distribution of matter in clusters of galaxies is mainly determined by the dominant dark matter component, however, physical processes involving baryonic matter are able to modify it significantly. We analyse a set of 500 pc resolution cosmological simulations of a cluster of galaxies with mass comparable to Virgo, performed with the AMR code RAMSES. We compare the mass density profiles of the dark, stellar and gaseous matter components of the cluster that result from different assumptions for the subgrid baryonic physics and galaxy formation processes. First, the prediction of a gravity only N-body simulation is compared to that of a hydrodynamical simulation with standard galaxy formation recipes, then all results are compared to a hydrodynamical simulation which includes thermal AGN feedback from Super Massive Black Holes (SMBH). We find the usual effects of overcooling and adiabatic contraction in the run with standard galaxy formation physics, but very different results are found when implementing SMBHs and AGN feedback. Star formation is strongly quenched, producing lower stellar densities throughout the cluster, and much less cold gas is available for star formation at low redshifts. At redshift z = 0 we find a flat density core of radius 10 kpc in both of the dark and stellar matter density profiles. We specu- late on the possible formation mechanisms able to produce such cores and we conclude that they can be produced through the coupling of different processes: (I) dynamical friction from the decay of black hole orbits during galaxy mergers; (II) AGN driven gas outflows producing fluctuations of the gravitational potential causing the removal of collisionless matter from the central region of the cluster; (III) adiabatic expansion in response to the slow expulsion of gas from the central region of the cluster during the quiescent mode of AGN activity.

We use 500 pc resolution cosmological simulations of a Virgo-like galaxy cluster to study the properties of the brightest cluster galaxy (BCG) that forms at the center of the halo. We compared two simulations; one incorporating only supernovae feedback and a second that also includes prescriptions for black hole growth and the resulting AGN feedback from gas accretion. As previous work has shown, with supernovae feedback alone we are unable to reproduce any of the observed properties of massive cluster ellipticals. The resulting BCG is rotating quickly, has a high Sersic index, a strong mass excess in the center and a total central density profile falling more steeply than isothermal. Furthermore, it is far too efficient at converting most of the available baryons into stars which is strongly constrained by abundance matching. With a treatment of black hole dynamics and AGN feedback the BCG properties are in good agreement with data: they rotate slowly, have a cored surface density profile, a flat or rising velocity dispersion profile and a low stellar mass fraction. The AGN provides a new mechanism to create cores in luminous elliptical galaxies; the core expands due to the combined effects of heating from dynamical friction of sinking massive black holes and AGN feedback that ejects gaseous material from the central regions.

We use 1 kpc resolution cosmological AMR simulations of a Virgo-like galaxy cluster to investigate the effect of feedback from supermassive black holes (SMBH) on the mass distribution of dark matter, gas and stars. We compared three different models: (i) a standard galaxy formation model featuring gas cooling, star formation and supernovae feedback, (ii) a "quenching" model for which star formation is artificially suppressed in massive halos and finally (iii) the recently proposed AGN feedback model of Booth & Schaye (2009). Without AGN feedback (even in the quenching case), our simulated cluster suffers from a strong overcooling problem, with a stellar mass fraction significantly above observed values in M87. The baryon distribution is highly concentrated, resulting in a strong adiabatic contraction (AC) of dark matter. With AGN feedback, on the contrary, the stellar mass in the bright central galaxy (BCG) lies below observational estimates and the overcooling problem disappears. The stellar mass of the BCG is seen to increase with increasing mass resolution, suggesting that our stellar masses converges to the correct value from below. The gas and total mass distributions are in striking agreement with observations. We also find a slight deficit (~10%) of baryons at the virial radius, due to the effect of AGN-driven shock waves pushing gas to Mpc scales and beyond. This baryon deficit results in a slight adiabatic expansion of the dark matter distribution, that can be explained quantitatively by AC theory.

We assess the possibility to detect and characterize the physical state of the missing baryons at low redshift by analyzing the X-ray absorption spectra of the Gamma Ray Burst [GRB] afterglows, measured by a micro calorimeters-based detector with 3 eV resolution and 1000 cm2 effective area and capable of fast re-pointing, similar to that on board of the recently proposed X-ray satellites EDGE and XENIA. For this purpose we have analyzed mock absorption spectra extracted from different hydrodynamical simulations used to model the properties of the Warm Hot Intergalactic Medium [WHIM]. These models predict the correct abundance of OVI absorption lines observed in UV and satisfy current X-ray constraints. According to these models space missions like EDGE and XENIA should be able to detect about 60 WHIM absorbers per year through the OVII line. About 45 % of these have at least two more detectable lines in addition to OVII that can be used to determine the density and the temperature of the gas. Systematic errors in the estimates of the gas density and temperature can be corrected for in a robust, largely model-independent fashion. The analysis of the GRB absorption spectra collected in three years would also allow to measure the cosmic mass density of the WHIM with about 15 % accuracy, although this estimate depends on the WHIM model. Our results suggest that GRBs represent a valid, if not preferable, alternative to Active Galactic Nuclei to study the WHIM in absorption. The analysis of the absorption spectra nicely complements the study of the WHIM in emission that the spectrometer proposed for EDGE and XENIA would be able to carry out thanks to its high sensitivity and large field of view.