Supermassive black holes (SMBHs) which reside at the centres of galaxies can inject vast amounts of energy into the surrounding gas and are thought to be a viable mechanism to quench star-formation in massive galaxies. Here we study the $10^{9-12.5}\,\mathrm{M_\odot}$ stellar mass central galaxy population of the IllustrisTNG simulation, specifically the TNG100 and TNG300 volumes at $z=0$, and show how the three components - SMBH, galaxy, and circumgalactic medium (CGM) - are interconnected in their evolution. We find that gas entropy is a sensitive diagnostic of feedback injection. In particular, we demonstrate how the onset of the low-accretion BH feedback mode, realised in the IllustrisTNG model as a kinetic, BH-driven wind, leads not only to star-formation quenching at stellar masses $\gtrsim10^{10.5}\,\mathrm{M_\odot}$ but also to a change in thermodynamic properties of the (non-star-forming) gas, both within the galaxy and beyond. The IllustrisTNG kinetic feedback from SMBHs increases the average gas entropy, within the galaxy and in the CGM, lengthening typical gas cooling times from $10-100\, \mathrm{Myr}$ to $1-10\,\mathrm{Gyr}$, effectively ceasing ongoing star-formation and inhibiting radiative cooling and future gas accretion. In practice, the same AGN feedback channel is simultaneously `ejective' and `preventative' and leaves an imprint on the temperature, density, entropy, and cooling times also in the outer reaches of the gas halo, up to distances of several hundred kiloparsecs. In the IllustrisTNG model, a long-lasting quenching state can occur for a heterogeneous CGM, whereby the hot and dilute CGM gas of quiescent galaxies contains regions of low-entropy gas with short cooling times.

We study the stellar-to-halo mass relation (SHMR) for central and satellite galaxies with total dynamical masses above 10^10.5 Msun using the suite of cosmological magneto-hydrodynamical simulations IllustrisTNG. In particular, we quantify environmental effects on satellite populations from TNG50, TNG100, and TNG300 located within the virial radius of group- and cluster-like hosts with total masses of 10^12-15.2 Msun. At fixed stellar mass, the satellite SHMR exhibits a distinct shift towards lower dynamical mass compared to the SHMR of centrals. Conversely, at fixed dynamical mass, satellite galaxies appear to have larger stellar-to-total mass fractions than centrals by up to a factor of a few. The systematic deviation from the central SHMR is larger for satellites in more massive hosts, at smaller cluster-centric distances, with earlier infall times, and that inhabit higher local density environments; moreover, it is in place already at early times (z < 2). Systematic environmental effects contribute to the perceived galaxy-to-galaxy variation in the measured SHMR when galaxies cannot be separated into satellites and centrals. The SHMR of satellites exhibits a larger scatter than centrals, over the whole range of dynamical mass (up to 0.8 dex). The shift of the satellite SHMR results mostly from tidal stripping of their dark matter, which affects satellites in an outside-in fashion: the departure of the satellite SHMR from the centrals' relation diminishes for measurements of dynamical mass in progressively smaller apertures. Finally, we provide a family of fitting functions for the SHMR predicted by IllustrisTNG.

In recent years, the outskirts of galaxy clusters have emerged as one of the new frontiers and unique laboratories for studying the growth of large scale structure in the universe. Modern cosmological hydrodynamical simulations make firm and testable predictions of the thermodynamic and chemical evolution of the X-ray emitting intracluster medium. However, recent X-ray and Sunyaev-Zeldovich effect observations have revealed enigmatic disagreements with theoretical predictions, which have motivated deeper investigations of a plethora of astrophysical processes operating in the virialization region in the cluster outskirts. Much of the physics of cluster outskirts is fundamentally different from that of cluster cores, which has been the main focus of X-ray cluster science over the past several decades. A next-generation X-ray telescope, equipped with sub-arcsecond spatial resolution over a large field of view along with a low and stable instrumental background, is required in order to reveal the full story of the growth of galaxy clusters and the cosmic web and their applications for cosmology.

We use IllustrisTNG, a suite of gravity and MHD simulations, to study the demographics and properties of jellyfish galaxies in the full cosmological context. By jellyfish galaxies, we mean satellites orbiting in massive groups and clusters that exhibit highly asymmetric distributions of gas and gas tails. We use the TNG100 run and select galaxies at redshifts $z\le0.6$ with stellar mass exceeding $10^{9.5}{\rm M_\odot}$ and with host halo masses of $10^{13}-10^{14.6}\,{\rm M_\odot}$. Among more than about 6000 (2600) galaxies with stars (and some gas), we identify 800 jellyfish galaxies by visually inspecting their gas and stellar mass maps in random projections. About $31\%$ of cluster satellites are found with signatures of ram-pressure stripping and gaseous tails stemming from the main luminous bodies. This is a lower limit, since the random orientation entails a loss of about $30\%$ of galaxies that in an optimal projection would otherwise be identified as jellyfish. The connection with ram-pressure stripping is further confirmed by a series of findings: jellyfish galaxies are more frequent at intermediate and large cluster-centric distances ($r/R_{\rm 200c}\gtrsim 0.25$); they move through the ICM with larger bulk velocities and Mach numbers than the general cluster population, typically orbiting supersonically and experiencing larger ram pressures. Furthermore, the gaseous tails usually extend in opposite directions to the galaxy trajectory, with no relation between tail orientation and the host's center. The frequency of jellyfish galaxies shows a very weak dependence on redshift $(0\le z\le0.6)$ but larger fractions of disturbed gaseous morphologies occur in more massive hosts and at smaller satellite masses. Finally, jellyfish galaxies are late infallers ($< 2.5-3$ Gyrs ago, at $z=0$) and the emergence of gaseous tails correlates well with the presence of bow shocks in the ICM.

The mass of galaxy clusters is not a direct observable, nonetheless it is commonly used to probe cosmological models. Based on the combination of all main cluster observables, that is, the X-ray emission, the thermal Sunyaev-Zel'dovich (SZ) signal, the velocity dispersion of the cluster galaxies, and gravitational lensing, the gravitational potential of galaxy clusters can be jointly reconstructed. We derive the two main ingredients required for this joint reconstruction: the potentials individually reconstructed from the observables and their covariance matrices, which act as a weight in the joint reconstruction. We show here the method to derive these quantities. The result of the joint reconstruction applied to a real cluster will be discussed in a forthcoming paper. We apply the Richardson-Lucy deprojection algorithm to data on a two-dimensional (2D) grid. We first test the 2D deprojection algorithm on a $\beta$-profile. Assuming hydrostatic equilibrium, we further reconstruct the gravitational potential of a simulated galaxy cluster based on synthetic SZ and X-ray data. We then reconstruct the projected gravitational potential of the massive and dynamically active cluster Abell 2142, based on the X-ray observations collected with XMM-Newton and the SZ observations from the Planck satellite. Finally, we compute the covariance matrix of the projected reconstructed potential of the cluster Abell 2142 based on the X-ray measurements collected with XMM-Newton. The gravitational potentials of the simulated cluster recovered from synthetic X-ray and SZ data are consistent, even though the potential reconstructed from X-rays shows larger deviations from the true potential. Regarding Abell 2142, the projected gravitational cluster potentials recovered from SZ and X-ray data reproduce well the projected potential inferred from gravitational-lensing observations. (abridged)

Cold Fronts and shocks are hallmarks of the complex intra-cluster medium (ICM) in galaxy clusters. They are thought to occur due to gas motions within the ICM and are often attributed to galaxy mergers within the cluster. Using hydro-cosmological simulations of clusters of galaxies, we show that collisions of inflowing gas streams, seen to penetrate to the very centre of about half the clusters, offer an additional mechanism for the formation of shocks and cold fronts in cluster cores. Unlike episodic merger events, a gas stream inflow persists over a period of several Gyrs and it could generate a particular pattern of multiple cold fronts and shocks.

We find, using cosmological simulations of galaxy clusters, that the hot X-ray emitting intra-cluster medium (ICM) enclosed within the outer accretion shock extends out to $R_{\rm shock}\sim(2 - 3) R_{\rm vir}$, where $R_{\rm vir}$ is the standard virial radius of the halo. Using a simple analytic model for satellite galaxies in the cluster, we evaluate the effect of ram pressure stripping on the gas in the inner discs and in the haloes at different distances from the cluster centre. We find that significant removal of star-forming disc gas occurs only at $r \lesssim 0.5 R_{\rm vir}$, while gas removal from the satellite halo is more effective and can occur when the satellite is found between $ R_{\rm vir}$ and $R_{\rm shock}$. Removal of halo gas sets the stage for quenching of the star formation by starvation over $2\textrm{--}3\,\mathrm{Gyr}$, prior to the satellite entry to the inner cluster halo. This scenario explains the presence of quenched galaxies, preferentially discs, at the outskirts of galaxy clusters, and the delayed quenching of satellites compared to central galaxies.

Virial shocks at edges of cosmic-web structures are a clear prediction of standard structure formation theories. We derive a criterion for the stability of the post-shock gas and of the virial shock itself in spherical, filamentary and planar infall geometries. When gas cooling is important, we find that shocks become unstable, and gas flows uninterrupted towards the center of the respective halo, filament or sheet. For filaments, we impose this criterion on self-similar infall solutions. We find that instability is expected for filament masses between $10^{11}-10^{13}M_\odot Mpc^{-1}.$ Using a simplified toy model, we then show that these filaments will likely feed halos with $10^{10}M_{\odot}\lesssim M_{halo}\lesssim 10^{13}M_{\odot}$ at redshift $z=3$, as well as $10^{12}M_{\odot}\lesssim M_{halo}\lesssim 10^{15}M_{\odot}$ at $z=0$. The instability will affect the survivability of the filaments as they penetrate gaseous halos in a non-trivial way. Additionally, smaller halos accreting onto non-stable filaments will not be subject to ram-pressure inside the filaments. The instreaming gas will continue towards the center, and stop either once its angular momentum balances the gravitational attraction, or when its density becomes so high that it becomes self-shielded to radiation.

We utilize cosmological simulations of 16 galaxy clusters at redshifts \zeq{0} and \zeq{0.6} to study the effect of inflowing streams on the properties of the X-ray emitting intracluster medium. We find that the mass accretion occurs predominantly along streams that originate from the cosmic web and consist of heated gas. Clusters that are unrelaxed in terms of their X-ray morphology are characterized by higher mass inflow rates and deeper penetration of the streams, typically into the inner third of the virial radius. The penetrating streams generate elevated random motions, bulk flows and cold fronts. The degree of penetration of the streams may change over time such that clusters can switch from being unrelaxed to relaxed over a time-scale of several giga years.

The massive galaxies in the young universe, ten billion years ago, formed stars at surprising intensities. Although this is commonly attributed to violent mergers, the properties of many of these galaxies are incompatible with such events, showing gas-rich, clumpy, extended rotating disks not dominated by spheroids (Genzel et al. 2006, 2008). Cosmological simulations and clustering theory are used to explore how these galaxies acquired their gas. Here we report that they are stream-fed galaxies, formed from steady, narrow, cold gas streams that penetrate the shock-heated media of massive dark matter haloes (Dekel & Birnboim 2006; Keres et al. 2005). A comparison with the observed abundance of star-forming galaxies implies that most of the input gas must rapidly convert to stars. One-third of the stream mass is in gas clumps leading to mergers of mass ratio greater than 1:10, and the rest is in smoother flows. With a merger duy cycle of 0.1, three-quarters of the galaxies forming stars at a given rate are fed by smooth streams. The rarer, submillimetre galaxies that form stars even more intensely are largely merger-induced starbursts. Unlike destructive mergers, the streams are likely to keep the rotating disk configuration intact, although turbulent and broken into giant star-forming clumps that merge into a central spheroid (Noguchi 1999; Genzel et al. 2008, Elmegreen, Bournaud & Elmegreen 2008, Dekel, Sari & Ceverino 2009). This stream-driven scenario for the formation of disks and spheroids is an alternative to the merger picture.