X-ray and microwave cluster scaling relations are immensely valuable for cosmological analysis. However, their power is limited by astrophysical systematics that bias mass estimates and introduce additional scatter. Turbulence injected into the intracluster medium via mass assembly contributes substantially to cluster non-thermal pressure support, a significant source of such uncertainties. We use an analytical model to compute the assembly-driven non-thermal pressure profiles of haloes based on Monte Carlo-generated accretion histories. We introduce a fitting function for the average non-thermal pressure fraction profile, which exhibits minimal dependence on redshift at fixed peak height. Using the model, we predict deviations from self-similarity and the intrinsic scatter in the Sunyaev--Zel'dovich effect observable-mass scaling relation ($Y_\mathrm{SZ}-M$) due solely to inter-cluster variation in mass accretion histories. We study the dependence of $Y_\mathrm{SZ}-M$ on aperture radius, cosmology, redshift, and mass limit. The model predicts $5-9\%$ scatter in $Y_\mathrm{SZ}-M$ at $z=0$, increasing as the aperture used to compute $Y_\mathrm{SZ}$ increases from $R_\mathrm{500c}$ to $5R_\mathrm{500c}$. The predicted scatter lies slightly below that of studies based on non-radiative hydro-simulations, illustrating that assembly history variance is likely responsible for a substantial fraction of scatter in $Y_\mathrm{SZ}-M$. This should be regarded as a lower bound, which will likely increase with the use of an updated gas density model that incorporates a more realistic response to halo assembly. As redshift increases, $Y_\mathrm{SZ}-M$ deviates more from self-similarity and scatter increases. We show that the $Y_\mathrm{SZ}-M$ residuals correlate strongly with the recent halo mass accretion rate, potentially providing an opportunity to infer the latter.

The diffuse plasma inside clusters of galaxies has X-ray emitting temperatures of a few keV. The physical mechanisms that heat this intracluster medium (ICM) to such temperatures include the accretion shock at the periphery of a galaxy cluster, the shocks driven by merger events, as well as a somewhat overlooked mechanism -- the dissipation of intracluster turbulent motions. We study the relative role of these heating mechanisms using galaxy clusters in Lagrangian tracer particle re-simulations of the Omega500 cosmological simulation. We adopt a novel analysis method of decomposing the temperature increase at each time step into the contribution from dissipative heating and that from adiabatic heating. In the high-resolution spatial-temporal map of these heating rates, merger tracks are clearly visible, demonstrating the dominant role of merger events in heating the ICM. The dissipative heating contributed by each merger event is extended in time and also occurs in the rarefaction regions, suggesting the importance of heating by the dissipation of merger-induced turbulence. Quantitative analysis shows that turbulence heating, rather than direct heating at merger shocks, dominates the temperature increase of the ICM, especially at inner radii $r < r_{\rm 500c}$. In addition, we find that many merger shocks can propagate with almost constant velocity to very large radii $r \gg r_{\rm 500c}$, some even reach and join with the accretion shock and becoming the outer boundary of the ICM. Altogether, these results suggest that the ICM is heated more in an `inside-out' fashion rather than `outside-in' as depicted in the classical smooth accretion picture.

We study the distribution of line-of-sight velocities of galaxies in the vicinity of SDSS redMaPPer galaxy clusters. Based on their velocities, galaxies can be split into two categories: galaxies that are dynamically associated with the cluster, and random line-of-sight projections. Both the fraction of galaxies associated with the galaxy clusters, and the velocity dispersion of the same, exhibit a sharp feature as a function of radius. The feature occurs at a radial scale $R_{\rm edge} \approx 2.2R_{\rm{\lambda}}$, where $R_{\rm{\lambda}}$ is the cluster radius assigned by redMaPPer. We refer to $R_{\rm edge}$ as the "edge radius." These results are naturally explained by a model that further splits the galaxies dynamically associated with a galaxy cluster into a component of galaxies orbiting the halo and an infalling galaxy component. The edge radius $R_{\rm edge}$ constitutes a true "cluster edge", in the sense that no orbiting structures exist past this radius. A companion paper (Aung et al. 2020) tests whether the "halo edge" hypothesis holds when investigating the full three-dimensional phase space distribution of dark matter substructures in numerical simulations, and demonstrates that this radius coincides with a suitably defined splashback radius.

The phase space structure of dark matter halos can be used to measure the mass of the halo, infer mass accretion rates, and probe the effects of modified gravity. Previous studies showed that the splashback radius can be measured in position space using the slope of the density profile. Using N-body simulations, we show that the phase space structure of the dark matter halo does not end at this splashback radius. Instead, there exists a region where infalling, splashback, and virialized halos are mixed spatially. We model the distribution of the three kinematically distinct populations and show that there exists an "edge radius" beyond which a dark matter halo has no orbiting substructures. This radius is a fixed multiple of the splashback radius as defined in previous works, and can be interpreted as a radius which contains a fixed fraction of the apocenters of dark matter particles. Our results provide a firm theoretical foundation to the satellite galaxy model adopted in the companion paper by Tomooka et al., where we analyzed the phase space distribution of SDSS redMaPPer clusters.

We study the effects of Kelvin Helmholtz Instability (KHI) on the cold streams that feed massive halos at high redshift, generalizing our earlier results to include the effects of radiative cooling and heating from a UV background, using analytic models and high resolution idealized simulations. We currently do not consider self-shielding, thermal conduction or gravity. A key parameter in determining the fate of the streams is the ratio of the cooling time in the turbulent mixing layer which forms between the stream and the background following the onset of the instability, $t_{\rm cool,mix}$, to the time in which the mixing layer expands to the width of the stream in the non-radiative case, $t_{\rm shear}$. This can be converted into a critical stream radius, $R_{\rm s,crit}$, such that $R_{\rm s}/R_{\rm s,crit}=t_{\rm shear}/t_{\rm cool,mix}$. If $R_{\rm s}<R_{\rm s,crit}$, the non-linear evolution proceeds similarly to the non-radiative case studied by Mandelker et al. 2019a. If $R_{\rm s}>R_{\rm s,crit}$, which we find to almost always be the case for astrophysical cold streams, the stream is not disrupted by KHI. Rather, background mass cools and condenses onto the stream, and can increase the mass of cold gas by a factor of $\sim 3$ within 10 stream sound crossing times. The mass entrainment induces thermal energy losses from the background and kinetic energy losses from the stream, which we model analytically. Roughly half of the dissipated energy is radiated away from gas with $T<5\times 10^4{\rm K}$, suggesting much of it will be emitted in Ly$\alpha$.

We present an analytic toy model for the radiation produced by the interaction between the cold streams thought to feed massive halos at high redshift and their hot CGM. We begin by deriving cosmologically motivated parameters for the streams as they enter the halo virial radius, $R_{\rm v}$, as a function of halo mass and redshift. For $10^{12}M_{\odot}$ halos at $z=2$, we find the Hydrogen number density in streams to be $n_{\rm H,s}\sim (0.1-5)\times 10^{-2}{\rm cm}^{-3}$, a factor of $\delta \sim (30-300)$ times denser than the hot CGM density, while the stream radii are in the range $R_{\rm s}\sim (0.03-0.50)R_{\rm v}$. As the streams accelerate towards the halo centre, they become denser and narrower. The stream-CGM interaction induces Kelvin-Helmholtz Instability (KHI), which leads to entrainment of CGM mass by the stream and therefore to stream deceleration by momentum conservation. Assuming that the entrainment rates derived by Mandelker et al. 2019 in the absence of gravity can be applied locally at each halocentric radius, we derive equations of motion for the stream in the halo. Using these, we derive the net acceleration, mass growth, and energy dissipation induced by the stream-CGM interaction, as a function of halo mass and redshift, for different CGM density profiles. For the range of model parameters considered, we find that the interaction can induce dissipation luminosities $L_{\rm diss}>10^{42}~{\rm erg~s^{-1}}$ within $\le 0.6 R_{\rm v}$ of halos with $M_{\rm v}>10^{12}M_{\odot}$ at $z=2$, with the emission scaling with halo mass and redshift approximately as $\propto M_{\rm v}\,(1+z)^2$. The magnitude and spatial extent of the emission produced in massive halos at high redshift is consistent with observed Ly$\alpha$ blobs, though better treatment of the UV background and self-shielding is needed to solidify this conclusion.

Self-gravitating gaseous filaments exist on many astrophysical scales, from sub-pc filaments in the interstellar medium to Mpc scale streams feeding galaxies from the cosmic web. These filaments are often subject to Kelvin-Helmotz Instability (KHI) due to shearing against a confining background medium. We study the nonlinear evolution of KHI in pressure-confined self-gravitating gas streams initially in hydrostatic equilibrium, using analytic models and hydrodynamic simulations, not including radiative cooling. We derive a critical line-mass as a function of the stream Mach number and density contrast with respect to the background, $\mu_{cr}(M_b,\delta_c)\le 1$, where $\mu=1$ is normalized to the maximal line mass for which initial hydrostatic equilibrium is possible. For $\mu<\mu_{cr}$, KHI dominates the stream evolution. A turbulent shear layer expands into the background and leads to stream deceleration at a similar rate to the non-gravitating case. However, with gravity, penetration of the shear layer into the stream is halted at roughly half the initial stream radius by stabilizing buoyancy forces, significantly delaying total stream disruption. Streams with $\mu_{cr}<\mu\le 1$ fragment and form round, long-lived clumps by gravitational instability (GI), with typical separations roughly 8 times the stream radius, similar to the case without KHI. When KHI is still somewhat effective, these clumps are below the spherical Jeans mass and are partially confined by external pressure, but they approach the Jeans mass as $\mu\rightarrow 1$ and GI dominates. We discuss potential applications of our results to filaments in the ISM and dense streams feeding galaxies at high redshift.

As part of the effort to meet the needs of the Large Synoptic Survey Telescope Dark Energy Science Collaboration (LSST DESC) for accurate, realistically complex mock galaxy catalogs, we have developed GalSampler, an open-source python package that assists in generating large volumes of synthetic cosmological data. The key idea behind GalSampler is to recast hydrodynamical simulations and semi-analytic models as physically-motivated galaxy libraries. GalSampler populates a new, larger-volume halo catalog with galaxies drawn from the baseline library; by using weighted sampling guided by empirical modeling techniques, GalSampler inherits statistical accuracy from the empirical model and physically-motivated complexity from the baseline library. We have recently used GalSampler to produce the cosmoDC2 extragalactic catalog made for the LSST DESC Data Challenge 2. Using cosmoDC2 as a guiding example, we outline how GalSampler can continue to support ongoing and near-future galaxy surveys such as the Dark Energy Survey (DES), the Dark Energy Spectroscopic Instrument (DESI), WFIRST, and Euclid.

The standard self-similar model of galaxy cluster formation predicts that the X-ray luminosity-temperature ($L_X$-$T_X$) relation of galaxy clusters should have been $L_X\propto T_X^2$ in absence of the baryonic physics, such as radiative cooling and feedback from stars and black holes. However, this baseline relation is predicted without considering the fact that the halo concentration and the characteristic density of clusters increases as their mass decreases, which is a consequence of hierarchical structure formation of the universe. Here, we show that the actual baseline relation should be $L_X\propto T_X^\alpha$, where $\alpha\sim 1.7$, instead of $\alpha=2$, given the mass dependence of the concentration and the fundamental plane relation of galaxy clusters. Numerical simulations show that $\alpha\sim 1.6$, which is consistent with the prediction. We also show that the baseline luminosity-mass ($L_X$-$M_\Delta$) relation should have been $L_X\propto M_\Delta^\beta$, where $\beta\sim 1.1$-1.2, in contrast with the conventional prediction ($\beta=4/3$). In addition, some of the scatter in the $L_X$-$M_\Delta$ relation can be attributed to the scatter in the concentration-mass ($c$-$M$) relation. The confirmation of the shallow slope could be a proof of hierarchical clustering. As an example, we show that the new baseline relations could be checked by studying the temperature or mass dependence of gas mass fraction of clusters. Moreover, the highest-temperature clusters would follow the shallow baseline relations if the influences of cool cores and cluster mergers are properly removed.

We study the effects of Kelvin-Helmholtz instability (KHI) on the cold streams that feed high-redshift galaxies through their hot haloes, generalizing our earlier analyses of a 2D slab to a 3D cylinder, but still limiting our analysis to the adiabatic case with no gravity. We combine analytic modeling and numerical simulations in the linear and non-linear regimes. For subsonic or transonic streams with respect to the halo sound speed, the instability in 3D is qualitatively similar to 2D, but progresses at a faster pace. For supersonic streams, the instability grows much faster in 3D and can be qualitatively different due to azimuthal modes, which introduce a strong dependence on the initial width of the stream-background interface. Using analytic toy models and approximations supported by high-resolution simulations, we apply our idealized hydrodynamical analysis to the astrophysical scenario. The upper limit for the radius of a stream that disintegrates prior to reaching the central galaxy is ~70% larger than the 2D estimate; it is in the range (0.5-5)% of the halo virial radius, decreasing with increasing stream density and velocity. Stream disruption generates a turbulent mixing zone around the stream with velocities at the level of ~20% of the initial stream velocity. KHI can cause significant stream deceleration and energy dissipation in 3d, contrary to 2D estimates. For typical streams, up to (10-50)% of the gravitational energy gained by inflow down the dark-matter halo potential can be dissipated, capable of powering Lyman-alpha blobs if most of it is dissipated into radiation.

Halos, filaments, sheets and voids in the cosmic web can be defined in terms of the eigenvalues of the smoothed shear tensor and a threshold $\lambda_{\rm th}$. Using analytic methods, we construct mean maps centered on these types of structures for Gaussian random fields corresponding to cosmological initial conditions. Each map also requires a choice of shear at the origin; we consider three possibilities. We find characteristic sizes, shapes and other properties of the central objects in these mean maps and explore how these properties change with varying the threshold and smoothing scale, i.e. varying the separation of the cosmic web into different kinds of components. The mean maps become increasingly complex as the threshold $\lambda_{\rm th}$ decreases to zero. We also describe scatter around these mean maps, subtleties which can arise in their construction, and some comparisons between halos in the maps and collapsed halos at final times.