We study the origins of 122 ultra-diffuse galaxies (UDGs) in the {\sc RomulusC} zoom-in cosmological simulation of a galaxy cluster (M$_{200} = 1.15\times10^{14}$ M$_{\odot}$), one of the only such simulations capable of resolving the evolution and structure of dwarf galaxies (M$_{\star} < 10^9$ M$_{\odot}$). We find broad agreement with observed cluster UDGs and predict that they are not separate from the overall cluster dwarf population. UDGs in cluster environments form primarily from dwarf galaxies that experienced early cluster in-fall and subsequent quenching due to ram pressure. The ensuing dimming of these dwarf galaxies due to passive stellar evolution results in a population of very low surface brightness galaxies that are otherwise typical dwarfs. UDGs and non-UDGs alike are affected by tidal interactions with the cluster potential. Tidal stripping of dark matter, as well as mass loss from stellar evolution, results in the adiabatic expansion of stars, particularly in the lowest mass dwarfs. High mass dwarf galaxies show signatures of tidal heating while low mass dwarfs that survive until $z=0$ typically have not experienced such impulsive interactions. There is little difference between UDGs and non-UDGs in terms of their dark matter halos, stellar morphology, colors, and location within the cluster. In most respects cluster UDG and non-UDGs alike are similar to isolated dwarf galaxies, except for the fact that they are typically quenched.

We present a novel population-based Bayesian inference approach to model the average and population variance of spatial distribution of a set of observables from ensemble analysis of low signal-to-noise ratio measurements. The method consists of (1) inferring the average profile using Gaussian Processes and (2) computing the covariance of the profile observables given a set of independent variables. Our model is computationally efficient and capable of inferring average profiles of a large population size from noisy measurements, without stacking and binning data nor parameterizing the shape of the mean profile. We demonstrate the performance of our method using dark matter, gas and stellar profiles extracted from hydrodynamical cosmological simulations of galaxy formation. Population Profile Estimator (PoPE) is publicly available in a GitHub repository. Our new method should be useful for measuring the spatial distribution and internal structure of a variety of astrophysical systems using large astronomical surveys.

The shape of dark matter haloes plays a critical role in constraining cosmology with upcoming large-scale structure surveys. In this paper, we study the correlations between the triaxial shapes and formation histories in dark matter haloes in the MultiDark Planck 2 N-body cosmological simulation. We find that halo ellipticity is strongly correlated with halo properties that serve as proxies of halo formation history, such as halo concentration and the peak-centroid offset. In particular, the correlation between halo ellipticity and halo concentration is nearly independent of the halo density peak height. We present a simple model for the correlation between halo ellipticity and concentration using conditional abundance matching, and provide fitting formulae for the multi-dimensional distributions of triaxial halo shape as a function of halo peak height. We apply our halo shape model to gauge the effects of halo ellipticity and orientation bias on the excess surface mass density profiles in cluster-size halos. Our model should be useful for exploring the impact of triaxial halo shape on cosmological constraints in upcoming weak lensing surveys of galaxy clusters.

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.

We use the Romulus25 cosmological simulation volume to identify the largest-ever simulated sample of field ultra-diffuse galaxies (UDGs). At z=0, we find that isolated UDGs have average star formation rates, colors, and virial masses for their stellar masses and environment. Although they tend to have slightly elevated HI masses, UDGs are consistent with the general isolated dwarf galaxy population and make up ~20% of all field galaxies with 10$^7$<M$_\star$/M$_\odot$<10$^{9}$. The properties of our UDGs are consistent with existing observations of field UDGs, but we predict that many isolated UDGs have been missed by current surveys. Despite their isolation at z=0, the UDGs in our sample are the products of major mergers at z>1. Mergers are no more common in UDG than non-UDG progenitors, but mergers that create UDGs tend to happen earlier, produce a temporary boost in spin, and cause star formation to migrate to the outskirts of galaxies, resulting in lower central star formation rates. The centers of the galaxies fade as their central stellar populations age, but their global star formation rates are maintained through bursts of star formation at larger radii. This formation channel is unique relative to other proposals for UDG formation in isolated galaxies, or for UDGs in groups and clusters, demonstrating that UDGs can potentially be formed through multiple mechanisms.

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.

CMB-HD is a proposed ultra-deep (0.5 uk-arcmin), high-resolution (15 arcseconds) millimeter-wave survey over half the sky that would answer many outstanding questions in both fundamental physics of the Universe and astrophysics. This survey would be delivered in 7.5 years of observing 20,000 square degrees, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. Each telescope would field 800,000 detectors (200,000 pixels), for a total of 1.6 million detectors.

We present the analysis of deep optical imaging of the galaxy cluster Abell 133 with the IMACS instrument on Magellan. Our multi-band photometry enables stellar mass measurements in the cluster member galaxies down to a mass limit of $M_\star=3\times10^8\,M_\odot$ ($\approx 0.1$ of the Large Magellanic Cloud stellar mass). We observe a clear difference in the spatial distribution of large and dwarf galaxies within the cluster. Modeling these galaxies populations separately, we can confidently track the distribution of stellar mass locked in the galaxies to the cluster's virial radius. The extended envelope of the cluster's brightest galaxy can be tracked to $\sim 200$ kpc. The central galaxy contributes $\sim 1/3$ of the the total cluster stellar mass within the radius $r_{500c}$.

In recent years, observations of the Sunyaev-Zeldovich (SZ) effect have had significant cosmological implications and have begun to serve as a powerful and independent probe of the warm and hot gas that pervades the Universe. As a few pioneering studies have already shown, SZ observations both complement X-ray observations -- the traditional tool for studying the intra-cluster medium -- and bring unique capabilities for probing astrophysical processes at high redshifts and out to the low-density regions in the outskirts of galaxy clusters. Advances in SZ observations have largely been driven by developments in centimetre-, millimetre-, and submillimetre-wave instrumentation on ground-based facilities, with notable exceptions including results from the Planck satellite. Here we review the utility of the thermal, kinematic, relativistic, non-thermal, and polarised SZ effects for studies of galaxy clusters and other large scale structures, incorporating the many advances over the past two decades that have impacted SZ theory, simulations, and observations. We also discuss observational results, techniques, and challenges, and aim to give an overview and perspective on emerging opportunities, with the goal of highlighting some of the exciting new directions in this field.

The intracluster medium (ICM) is a multi-phase environment, dynamically regulated by Active Galactic Nuclei (AGN), the motions of galaxies through it, and mergers with other clusters. AGN as a central heating source are key to preventing runaway cooling flows, but their role in heating cores in a cosmological context is still poorly understood. The activity of the AGN is strongly linked to star formation, especially in the Brightest Cluster Galaxy (BCG), likely because both rely on cold phase gas. A self-consistent model for AGN and star formation in galaxy clusters thus requires cosmological context, higher resolution, and a careful modeling of cooling and heating balance. In this paper, we use the high-resolution hydrodynamical cosmological simulation of the RomulusC galaxy cluster to study in detail the role of AGN and a major, head-on merger in shaping the cluster core. The unprecedented resolution of the RomulusC simulation captures the multiphase structure of the ICM. The realistic large-scale outflows launched by very small-scale thermal injections, the improved modeling of turbulent diffusion and mixing, and the particle nature of the simulation allow us to carefully separate different heating channels. We show that AGN activity, while efficient at regulating star formation, is incapable of destroying a CC. Instead, that process is facilitated by a head-on, 1:8 mass ratio merger. The merger generates bulk and turbulent motions, which in turn mix high entropy gas generated by AGN and merger driven shocks, turbulent dissipation and sloshing of the ICM by infalling substructures. While central cooling times remain shorter than the Hubble time, restoring a CC is made more difficult by the reduced precipitation rates at larger radii, emphasizing that the AGN-ICM connection is truly a multi-scale problem.

Cross-correlation analysis of the thermal Sunyaev-Zel'dovich (tSZ) effect and weak gravitational lensing (WL) provides a powerful probe of cosmology and astrophysics of the intra-cluster medium. We present the measurement of the cross-correlation of tSZ and WL from Planck and Subaru Hyper-Suprime Cam. The combination enables us to study cluster astrophysics at high redshift. We use the tSZ-WL cross-correlation and the tSZ auto-power spectrum measurements to place a tight constraint on the hydrostatic mass bias, which is a measure of the degree of non-thermal pressure support in galaxy clusters. With the prior on cosmological parameters derived from the analysis of the cosmic microwave background anisotropies by Planck and taking into account foreground contributions both in the tSZ auto-power spectrum and the tSZ-WL cross-correlation, the hydrostatic mass bias is estimated to be $26.9^{+8.9}_{-4.4} \%$ ($68\%$ C.L.), which is consistent with recent measurements by mass calibration techniques.

Upcoming multi-wavelength astronomical surveys will soon discover all massive galaxy clusters and provide unprecedented constraints on cosmology and cluster astrophysics. In this paper, we investigate the constraining power of the multi-band cluster surveys, through a joint analysis of three observables associated with clusters of galaxies, including thermal Sunyaev-Zel'dovich (tSZ) effect in cosmic microwave background (CMB), X-ray emission of ionized gas, and gravitational weak lensing effect of background galaxies by the cluster's gravitational potential. We develop a theoretical framework to predict and interpret two-point correlation statistics among the three observables using a semi-analytic model of intracluster medium (ICM) and halo-based approach. In this work, we show that the auto- and cross-angular power spectra in tSZ, X-ray and lensing statistics from upcoming missions (eROSITA, CMB-S4, and LSST) can help break the degeneracy between cosmology and ICM physics. These correlation statistics are less sensitive to selection biases, and are able to probe ICM physics in distant, faint and small clusters that are otherwise difficult to be detected individually. We show that the correlation statistics are able to provide cosmological constraints comparable to the conventional cluster abundance measurements, while constraining cluster astrophysics at the same time. Our results indicate that the correlation statistics can significantly enhance the scientific returns of upcoming multi-wavelength cluster surveys.

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.

We present a machine learning approach for estimating galaxy cluster masses, trained using both Chandra and eROSITA mock X-ray observations of 2,041 clusters from the Magneticum simulations. We train a random forest regressor, an ensemble learning method based on decision tree regression, to predict cluster masses using an input feature set. The feature set uses core-excised X-ray luminosity and a variety of morphological parameters, including surface brightness concentration, smoothness, asymmetry, power ratios, and ellipticity. The regressor is cross-validated and calibrated on a training sample of 1,615 clusters (80% of sample), and then results are reported as applied to a test sample of 426 clusters (20% of sample). This procedure is performed for two different mock observation series in an effort to bracket the potential enhancement in mass predictions that can be made possible by including dynamical state information. The first series is computed from idealized Chandra-like mock cluster observations, with high spatial resolution, long exposure time (1 Ms), and the absence of background. The second series is computed from realistic-condition eROSITA mocks with lower spatial resolution, short exposures (2 ks), instrument effects, and background photons modeled. We report a 20% reduction in the mass estimation scatter when either series is used in our random forest model compared to a standard regression model that only employs core-excised luminosity. The morphological parameters that hold the highest feature importance are smoothness, asymmetry, and surface brightness concentration. Hence, these parameters, which encode the dynamical state of the cluster, can be used to make more accurate predictions of cluster masses in upcoming surveys, offering a crucial step forward for cosmological analyses.

This Science White Paper, prepared in response to the ESA Voyage 2050 call for long-term mission planning, aims to describe the various science possibilities that can be realized with an L-class space observatory that is dedicated to the study of the interactions of cosmic microwave background (CMB) photons with the cosmic web. Our aim is specifically to use the CMB as a backlight -- and survey the gas, total mass, and stellar content of the entire observable Universe by means of analyzing the spatial and spectral distortions imprinted on it. These distortions result from two major processes that impact on CMB photons: scattering by electrons (Sunyaev-Zeldovich effect in diverse forms, Rayleigh scattering, resonant scattering) and deflection by gravitational potential (lensing effect). Even though the list of topics collected in this White Paper is not exhaustive, it helps to illustrate the exceptional diversity of major scientific questions that can be addressed by a space mission that will reach an angular resolution of 1.5 arcmin (goal 1 arcmin), have an average sensitivity better than 1 uK-arcmin, and span the microwave frequency range from roughly 50 GHz to 1 THz. The current paper also highlights the synergy of our BACKLIGHT mission concept with several upcoming and proposed ground-based CMB experiments.

This paper discusses the science case for a sensitive spectro-polarimetric survey of the microwave sky. Such a survey would provide a tomographic and dynamic census of the three-dimensional distribution of hot gas, velocity flows, early metals, dust, and mass distribution in the entire Hubble volume, exploit CMB temperature and polarisation anisotropies down to fundamental limits, and track energy injection and absorption into the radiation background across cosmic times by measuring spectral distortions of the CMB blackbody emission. In addition to its exceptional capability for cosmology and fundamental physics, such a survey would provide an unprecedented view of microwave emissions at sub-arcminute to few-arcminute angular resolution in hundreds of frequency channels, a data set that would be of immense legacy value for many branches of astrophysics. We propose that this survey be carried-out with a large space mission featuring a broad-band polarised imager and a moderate resolution spectro-imager at the focus of a 3.5m aperture telescope actively cooled to about 8K, complemented with absolutely-calibrated Fourier Transform Spectrometer modules observing at degree-scale angular resolution in the 10-2000 GHz frequency range. We propose two observing modes: a survey mode to map the entire sky as well as a few selected wide fields, and an observatory mode for deeper observations of regions of specific interest.

Recent statistical X-ray measurements of the intracluster medium (ICM) indicate that gas temperature profiles in the outskirts of galaxy clusters deviate from self-similar evolution. Using a mass-limited sample of galaxy clusters from cosmological hydrodynamical simulations, we show that the departure from self-similarity can be explained by non-thermal gas motions driven by mergers and accretion. Contrary to previous claims, gaseous substructures only play a minor role in the temperature evolution in cluster outskirts. A careful choice of halo overdensity definition in self-similar scaling mitigates these departures. Our work highlights the importance of non-thermal gas motions in ICM evolution and the use of galaxy clusters as cosmological probes.

The majority of the ordinary matter in the local Universe has been heated by strong structure formation shocks and resides in a largely unexplored hot, diffuse, X-ray emitting plasma that permeates the halos of galaxies, galaxy groups and clusters, and the cosmic web. We propose a next-generation "Cosmic Web Explorer" that will permit a complete and exhaustive understanding of these unseen baryons. This will be the first mission capable to reach the accretion shocks located several times farther than the virial radii of galaxy clusters, and reveal the out-of-equilibrium parts of the intra-cluster medium which are live witnesses to the physics of cosmic accretion. It will also enable a view of the thermal, kinematic, and chemical composition of the circum-galactic medium in galaxies with masses similar to the Milky Way, at the same level of detail that $Athena$ will unravel for the virialized regions of massive galaxy clusters, delivering a transformative understanding of the evolution of those galaxies in which most of the stars and metals in the Universe were formed. Finally, the proposed X-ray satellite will connect the dots of the large-scale structure by mapping, at high spectral resolution, as much as 100% of the diffuse gas hotter than $10^6$ K that fills the filaments of the cosmic web at low redshifts, down to an over-density of 1, both in emission and in absorption against the ubiquitous cosmic X-ray background, surveying up to half of the sky over 5 years in orbit. This requires a large effective area (~10 m$^2$ at 1 keV) over a large field of view (~1 deg$^2$), a megapixel cryogenic microcalorimeter array providing integral field spectroscopy with a resolving power $E/\Delta E$ = 2000 at 0.6 keV and a spatial resolution of 5" in the soft X-ray band, and a Low-Earth orbit that ensures a stable instrumental background and high sensitivity to faint, extended emission.

We provide an overview of the science case, instrument configuration and project plan for the next-generation ground-based cosmic microwave background experiment CMB-S4, for consideration by the 2020 Decadal Survey.

We present an ultraviolet quasar absorption line analysis of metal lines associated with three strong intervening H I absorbers (with $N$(H I) $>$ 10$^{16.5}$ cm$^{-2}$) detected in the outskirts of Sunyaev-Zel'dovich (SZ) effect-selected galaxy clusters ($z_{\rm cl} \sim 0.4 - 0.5$), within clustocentric impact parameters of $\rho_{cl}$ $\sim$ $(1.6 - 4.7)r_{500}$. Discovered in a recent set of targeted far-UV $HST$/COS spectroscopic observations, these absorbers have the highest H I column densities ever observed in the outskirts of galaxy clusters, and are also rich in metal absorption lines. Photoionization models yield single phase solutions for the three absorbers with gas densities of $n_{H} \sim 10^{-3} - 10^{-4}$ cm$^{-3}$ and metallicities of [X/H] $>$ -1.0 (from one-tenth solar to near-solar). The widths of detected absorption lines suggest gas temperatures of $T \sim 10^4$ K. The inferred densities (temperatures) are significantly higher (lower) compared to the X-ray emitting intracluster medium in cluster cores. The absorbers are tracing a cool phase of the intracluster gas in the cluster outskirts, either associated with gas stripped from cluster galaxies via outflows, tidal streams or ram-pressure forces, or denser regions within the intracluster medium that were uniformly chemically enriched from an earlier epoch of enhanced supernova and AGN feedback.

We present the science case, reference design, and project plan for the Stage-4 ground-based cosmic microwave background experiment CMB-S4.

A millimeter-wave survey over half the sky, that spans frequencies in the range of 30 to 350 GHz, and that is both an order of magnitude deeper and of higher-resolution than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. By providing such a deep, high-resolution millimeter-wave survey (about 0.5 uK-arcmin noise and 15 arcsecond resolution at 150 GHz), CMB-HD will enable major advances. It will allow 1) the use of gravitational lensing of the primordial microwave background to map the distribution of matter on small scales (k~10/hMpc), which probes dark matter particle properties. It will also allow 2) measurements of the thermal and kinetic Sunyaev-Zel'dovich effects on small scales to map the gas density and gas pressure profiles of halos over a wide field, which probes galaxy evolution and cluster astrophysics. In addition, CMB-HD would allow us to cross critical thresholds in fundamental physics: 3) ruling out or detecting any new, light (< 0.1eV), thermal particles, which could potentially be the dark matter, and 4) testing a wide class of multi-field models that could explain an epoch of inflation in the early Universe. Such a survey would also 5) monitor the transient sky by mapping the full observing region every few days, which opens a new window on gamma-ray bursts, novae, fast radio bursts, and variable active galactic nuclei. Moreover, CMB-HD would 6) provide a census of planets, dwarf planets, and asteroids in the outer Solar System, and 7) enable the detection of exo-Oort clouds around other solar systems, shedding light on planet formation. CMB-HD will deliver this survey in 5 years of observing half the sky, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. The telescopes will field about 2.4 million detectors (600,000 pixels) in total.

We have imaged the Sunyaev-Zel'dovich (SZ) effect signals at 140 and 270 GHz towards ten galaxy clusters with Bolocam and AzTEC/ASTE. We also used Planck data to constrain the signal at large angular scales, Herschel-SPIRE images to subtract the brightest galaxies that comprise the cosmic infrared background (CIB), Chandra imaging to map the electron temperature $T_e$ of the intra-cluster medium (ICM), and HST imaging to derive models of each galaxy cluster's mass density. The galaxy clusters gravitationally lens the background CIB, which produced an on-average reduction in brightness towards the galaxy clusters' centers after the brightest galaxies were subtracted. We corrected for this deficit, which was between 5-25% of the 270 GHz SZ effect signal within $R_{2500}$. Using the SZ effect measurements, along with the X-ray constraint on $T_e$, we measured each galaxy cluster's average line of sight (LOS) velocity $v_z$ within $R_{2500}$, with a median per-cluster uncertainty of +-700 km/s. We found an ensemble-mean <$v_z$> of 430+-210 km/s, and an intrinsic cluster-to-cluster scatter $\sigma_{int}$ of 470+-340 km/s. We also obtained maps of $v_z$ over each galaxy cluster's face with an angular resolution of 70". All four galaxy clusters previously identified as having a merger oriented along the LOS showed an excess variance in these maps at a significance of 2-4$\sigma$, indicating an internal $v_z$ rms of $\gtrsim$1000 km/s. None of the six galaxy clusters previously identified as relaxed or plane of sky mergers showed any such excess variance.

We present a machine-learning approach for estimating galaxy cluster masses from Chandra mock images. We utilize a Convolutional Neural Network (CNN), a deep machine learning tool commonly used in image recognition tasks. The CNN is trained and tested on our sample of 7,896 Chandra X-ray mock observations, which are based on 329 massive clusters from the IllustrisTNG simulation. Our CNN learns from a low resolution spatial distribution of photon counts and does not use spectral information. Despite our simplifying assumption to neglect spectral information, the resulting mass values estimated by the CNN exhibit small bias in comparison to the true masses of the simulated clusters (-0.02 dex) and reproduce the cluster masses with low intrinsic scatter, 8% in our best fold and 12% averaging over all. In contrast, a more standard core-excised luminosity method achieves 15-18% scatter. We interpret the results with an approach inspired by Google DeepDream and find that the CNN ignores the central regions of clusters, which are known to have high scatter with mass.

Quasar absorption-line studies in the ultraviolet (UV) can uniquely probe the nature of the multiphase cool-warm (10^4 < T < 10^6 K) gas in and around galaxy clusters, promising to provide unprecedented insights into 1) interactions between the circumgalactic medium (CGM) associated with infalling galaxies and the hot (T > 10^6 K) X-ray emitting intracluster medium (ICM), 2) the stripping of metal-rich gas from the CGM, and 3) a multiphase structure of the ICM with a wide range of temperatures and metallicities. In this work, we present results from a high-resolution simulation of a ~10^14 solar mass galaxy cluster to study the physical properties and observable signatures of this cool-warm gas in galaxy clusters. We show that the ICM becomes increasingly multiphased at large radii, with the cool-warm gas becoming dominant in cluster outskirts. The diffuse cool-warm gas also exhibits a wider range of metallicity than the hot X-ray emitting gas. We make predictions for the covering fractions of key absorption-line tracers, both in the ICM and in the CGM of cluster galaxies, typically observed with the Cosmic Origins Spectrograph aboard the Hubble Space Telescope (HST). We further extract synthetic spectra to demonstrate the feasibility of detecting and characterizing the thermal, kinematic, and chemical composition of the cool-warm gas using H I, O VI, and C IV lines, and we predict an enhanced population of broad Ly-alpha absorbers tracing the warm gas. Lastly, we discuss future prospects of probing the multiphase structure of the ICM beyond HST.

Understanding the processes that drive galaxy formation and shape the observed properties of galaxies is one of the most interesting and challenging frontier problems of modern astrophysics. We now know that the evolution of galaxies is critically shaped by the energy injection from accreting supermassive black holes (SMBHs). However, it is unclear how exactly the physics of this feedback process affects galaxy formation and evolution. In particular, a major challenge is unraveling how the energy released near the SMBHs is distributed over nine orders of magnitude in distance throughout galaxies and their immediate environments. The best place to study the impact of SMBH feedback is in the hot atmospheres of massive galaxies, groups, and galaxy clusters, which host the most massive black holes in the Universe, and where we can directly image the impact of black holes on their surroundings. We identify critical questions and potential measurements that will likely transform our understanding of the physics of SMBH feedback and how it shapes galaxies, through detailed measurements of (i) the thermodynamic and velocity fluctuations in the intracluster medium (ICM) as well as (ii) the composition of the bubbles inflated by SMBHs in the centers of galaxy clusters, and their influence on the cluster gas and galaxy growth, using the next generation of high spectral and spatial resolution X-ray and microwave telescopes.

We use a statistical sample of galaxy clusters from a large cosmological $N$-body$+$hydrodynamics simulation to examine the relation between morphology, or shape, of the X-ray emitting intracluster medium (ICM) and the mass accretion history of the galaxy clusters. We find that the mass accretion rate (MAR) of a cluster is correlated with the ellipticity of the ICM. The correlation is largely driven by material accreted in the last $\sim 4.5$~Gyr, indicating a characteristic time-scale for relaxation of cluster gas. Furthermore, we find that the ellipticity of the outer regions ($R\sim R_{\rm 500c}$) of the ICM is correlated with the overall MAR of clusters, while ellipticity of the inner regions ($\lesssim 0.5 R_{\rm 500c}$) is sensitive to recent major mergers with mass ratios of $\geq 1:3$. Finally, we examine the impact of variations in cluster mass accretion history on the X-ray observable-mass scaling relations. We show that there is a {\it continuous\/} anti-correlation between the residuals in the $T_x-M$ relation and cluster MARs, within which merging and relaxed clusters occupy extremes of the distribution rather than form two peaks in a bi-modal distribution, as was often assumed previously. Our results indicate the systematic uncertainties in the X-ray observable-mass relations can be mitigated by using the information encoded in the apparent ICM ellipticity.

(abridged) We present the Significant Cosmic Holes in Universe (SCHU) method for identifying cosmic voids and loops of filaments in cosmological datasets and assigning their statistical significance using techniques from topological data analysis. Persistent homology is used to find different dimensional holes. For dark matter halo catalogs and galaxy surveys, the 0-, 1-, and 2-dimensional holes can be identified with clusters, loops of filaments, and voids. The procedure overlays halos/galaxies on a 3D grid, and a distance-to-measure (DTM) function is calculated at each point of the grid. A filtration is generated over the lower-level sets of the DTM across increasing threshold values. The filtered simplicial complex can be used to summarize the birth/death times of the different dimension homology group generators (i.e., the holes). Persistence diagrams are produced from the dimension and birth/death times of each homology group generator. Using the persistence diagrams and bootstrap sampling, we explain how $p$-values can be assigned to each homology group generator. The homology group generators on a persistence diagram are not, in general, uniquely located back in the original dataset volume so we propose a method for finding a representation of the homology group generators. This method provides a novel, statistically rigorous approach for locating informative generators in cosmological datasets, which may be useful for providing complementary cosmological constraints on the effects of, for example, the sum of the neutrino masses. The method is tested on a Voronoi foam simulation, and then applied to a subset of the SDSS galaxy survey and a cosmological simulation. Lastly, we calculate Betti functions for two of the MassiveNuS simulations and discuss implications for using the persistent homology of the density field to help break degeneracy in the cosmological parameters.

Galaxy clusters are massive dark matter-dominated systems filled with X-ray emitting, optically thin plasma. Their large size and relative simplicity (at least as astrophysical objects go) make them a unique laboratory to measure some of the interesting plasma properties that are inaccessible by other means but fundamentally important for understanding and modeling many astrophysical phenomena -- from solar flares to black hole accretion to galaxy formation and the emergence of the cosmological Large Scale Structure. While every cluster astrophysicist is eagerly anticipating the direct gas velocity measurements from the forthcoming microcalorimeters onboard XRISM, Athena and future missions such as Lynx, a number of those plasma properties can best be probed by high-resolution X-ray imaging of galaxy clusters. Chandra has obtained some trailblazing results, but only grazed the surface of such studies. In this white paper, we discuss why we need arcsecond-resolution, high collecting area, low relative background X-ray imagers (with modest spectral resolution), such as the proposed AXIS and the imaging detector of Lynx.

The past decade has seen an explosion of discoveries and new insights into the diffuse gas within galaxies, galaxy clusters, and the filaments composing the Cosmic Web. A new decade will bring fresh opportunities to further this progress towards developing a comprehensive view of the composition, thermal state, and physical processes of diffuse gas in the Universe. Ultraviolet (UV) spectroscopy, probing diffuse 10^4-10^6 K gas at high spectral resolution, is uniquely poised to (1) witness environmental galaxy quenching processes in action, such as strangulation and tidal- and ram-pressure stripping, (2) directly account for the baryon content of galaxy clusters in the cold-warm (T<10^6 K) gas, (3) determine the phase structure and kinematics of gas participating in the equilibrium-regulating exchange of energy at the cores of galaxy clusters, and (4) map cold streams and filaments of the Cosmic Web that feed galaxies and clusters. With a substantial UV undertaking beyond the Hubble Space Telescope, all of the above would be achievable over the entire epoch of galaxy cluster formation. Such capabilities, coupled with already-planned advancements at other wavelengths, will transform extragalactic astronomy by revealing the dominant formation and growth mechanisms of gaseous halos over the mass spectrum, settling the debate between early- and late-time metal enrichment scenarios, and revealing how the ecosystems in which galaxies reside ultimately facilitate their demise.

Clusters of galaxies, the largest collapsed structures in the Universe, are located at the intersection of extended filaments of baryons and dark matter. Cosmological accretion onto clusters through large scale filaments adds material at cluster outskirts. Kinetic energy in the form of bulk motions and turbulence due to this accretion provides a form of pressure support against gravity, supplemental to thermal pressure. Significant amount of non-thermal pressure support could bias cluster masses derived assuming hydrostatic equilibrium, the primary proxy for cluster cosmology studies. Sensitive measurements of Doppler broadening and shift of astrophysical lines, and the relative fluctuations in thermodynamical quantities (e.g., density, pressure, and entropy) are primary diagnostic tools. Forthcoming planned and proposed X-ray (with large etendue, throughput, and high spectral resolution) and SZ observatories will provide crucial information on the assembly and virialisation processes of clusters, involving turbulent eddies cascading at various spatial scales and larger gas bulk motions in their external regions to the depth or their potential wells.

Sunyaev-Zeldovich (SZ) effects were first proposed in the 1970s as tools to identify the X-ray emitting hot gas inside massive clusters of galaxies and obtain their velocities relative to the cosmic microwave background (CMB). Yet it is only within the last decade that they have begun to significantly impact astronomical research. Thanks to the rapid developments in CMB instrumentation, measurement of the dominant thermal signature of the SZ effects has become a routine tool to find and characterize large samples of galaxy clusters and to seek deeper understanding of several important astrophysical processes via high-resolution imaging studies of many targets. With the notable exception of the Planck satellite and a few combinations of ground-based observatories, much of this "SZ revolution" has happened in the photometric mode, where observations are made at one or two frequencies in the millimeter regime to maximize the cluster detection significance and minimize the foregrounds. Still, there is much more to learn from detailed and systematic analyses of the SZ spectra across multiple wavelengths, specifically in the submillimeter (>300 GHz) domain. The goal of this Science White Paper is to highlight this particular aspect of SZ research, point out what new and potentially groundbreaking insights can be obtained from these studies, and emphasize why the coming decade can be a golden era for SZ spectral measurements.

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.

Achieving a precise understanding of galaxy formation in a cosmological context is one of the great challenges in theoretical astrophysics, due to the vast range of spatial scales involved in the relevant physical processes. Observations in the millimeter bands, particularly those using the cosmic microwave background (CMB) radiation as a "backlight", provide a unique probe of the thermodynamics of these processes, with the capability to directly measure the density, pressure, and temperature of ionized gas. Moreover, these observations have uniquely high sensitivity into the outskirts of the halos of galaxies and clusters, including systems at high redshift. In the next decade, the combination of large spectroscopic and photometric optical galaxy surveys and wide-field, low-noise CMB surveys will transform our understanding of galaxy formation via these probes.

Opening up a new window of millimeter-wave observations that span frequency bands in the range of 30 to 500 GHz, survey half the sky, and are both an order of magnitude deeper (about 0.5 uK-arcmin) and of higher-resolution (about 10 arcseconds) than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. In particular, such a survey would allow for major advances in measuring the distribution of dark matter and gas on small-scales, and yield needed insight on 1.) dark matter particle properties, 2.) the evolution of gas and galaxies, 3.) new light particle species, 4.) the epoch of inflation, and 5.) the census of bodies orbiting in the outer Solar System.

The Sunyaev-Zeldovich (SZ) effect was first predicted nearly five decades ago, but has only recently become a mature tool for performing high resolution studies of the warm and hot ionized gas in and between galaxies, groups, and clusters. Galaxy groups and clusters are powerful probes of cosmology, and they also serve as hosts for roughly half of the galaxies in the Universe. In this white paper, we outline the advances in our understanding of thermodynamic and kinematic properties of the warm-hot universe that can come in the next decade through spatially and spectrally resolved measurements of the SZ effects. Many of these advances will be enabled through new/upcoming millimeter/submillimeter (mm/submm) instrumentation on existing facilities, but truly transformative advances will require construction of new facilities with larger fields of view and broad spectral coverage of the mm/submm bands.

The total mass of a galaxy cluster is one of its most fundamental properties. Together with the redshift, the mass links observation and theory, allowing us to use the cluster population to test models of structure formation and to constrain cosmological parameters. Building on the rich heritage from X-ray surveys, new results from Sunyaev-Zeldovich and optical surveys have stimulated a resurgence of interest in cluster cosmology. These studies have generally found fewer clusters than predicted by the baseline Planck LCDM model, prompting a renewed effort on the part of the community to obtain a definitive measure of the true cluster mass scale. Here we review recent progress on this front. Our theoretical understanding continues to advance, with numerical simulations being the cornerstone of this effort. On the observational side, new, sophisticated techniques are being deployed in individual mass measurements and to account for selection biases in cluster surveys. We summarise the state of the art in cluster mass estimation methods and the systematic uncertainties and biases inherent in each approach, which are now well identified and understood, and explore how current uncertainties propagate into the cosmological parameter analysis. We discuss the prospects for improvements to the measurement of the mass scale using upcoming multi-wavelength data, and the future use of the cluster population as a cosmological probe.

In recent years, machine learning (ML) methods have remarkably improved how cosmologists can interpret data. The next decade will bring new opportunities for data-driven cosmological discovery, but will also present new challenges for adopting ML methodologies and understanding the results. ML could transform our field, but this transformation will require the astronomy community to both foster and promote interdisciplinary research endeavors.

The detailed velocity structure of the diffuse X-ray emitting intra-cluster medium (ICM) remains one of the last missing key ingredients in understanding the microphysical properties of these hot baryons and constraining our models of the growth and evolution of structure on the largest scales in the Universe. Direct measurements of the gas velocities from the widths and shifts of X-ray emission lines were recently provided for the central region of the Perseus Cluster of galaxies by $Hitomi$, and upcoming high-resolution X-ray microcalorimeters onboard $XRISM$ and $Athena$ are expected to extend these studies to many more systems. In the mean time, several other direct and indirect methods have been proposed for estimating the velocity structure in the ICM, ranging from resonant scattering to X-ray surface brightness fluctuation analysis, the kinematic Sunyaev-Zeldovich effect, or using optical line emitting nebulae in the brightest cluster galaxies as tracers of the motions of the ambient plasma. Here, we review and compare the existing estimates of the velocities of the hot baryons, as well as the various overlapping physical processes that drive motions in the ICM, and discuss the implications of these measurements for constraining the viscosity and identifying the source of turbulence in clusters of galaxies.

We use a complete sample of 38 richest maxBCG clusters to study the ICM-galaxy scaling relations and the halo mass selection properties of the maxBCG algorithm, based on X-ray and optical observations. The clusters are selected from the two largest bins of optical richness in the Planck stacking work with the maxBCG richness $N_{200} \geq 78$. We analyze their Chandra and XMM-Newton data to derive the X-ray properties of the ICM. We then use the distribution of $P(X|N)$, $X=T_X,\ L_X,\ Y_X$, to study the mass selection $P(M|N)$ of maxBCG. Compared with previous works based on the whole richness sample, a significant fraction of blended systems with boosted richness is skewed into this richest sample. Parts of the blended halos are picked apart by the redMaPPer, an updated red-sequence cluster finding algorithm with lower mass scatter. Moreover, all the optical blended halos are resolved as individual X-ray halos, following the established $L_X-T_X$ and $L_X-Y_X$ relations. We further discuss that the discrepancy between ICM-galaxy scaling relations, especially for future blind stacking, can come from several factors, including miscentering, projection, contamination of low mass systems, mass bias and covariance bias. We also evaluate the fractions of relaxed and cool core clusters in our sample. Both are smaller than those from SZ or X-ray selected samples. Moreover, disturbed clusters show a higher level of mass bias than relaxed clusters.

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.

As the largest virialized structures in the universe, galaxy clusters continue to grow and accrete matter from the cosmic web. Due to the low gas density in the outskirts of clusters, measurements are very challenging, requiring extremely sensitive telescopes across the entire electromagnetic spectrum. Observations using X-rays, the Sunyaev-Zeldovich effect, and weak lensing and galaxy distributions from the optical band, have over the last decade helped to unravel this exciting new frontier of cluster astrophysics, where the infall and virialization of matter takes place. Here, we review the current state of the art in our observational and theoretical understanding of cluster outskirts, and discuss future prospects for exploration using newly planned and proposed observatories.

We present the first results from RomulusC, the highest resolution cosmological hydrodynamic simulation of a galaxy cluster run to date. RomulusC, a zoom-in simulation of a halo with $z=0$ mass $10^{14}$ M$_{\odot}$, is run with the same sub-grid physics and resolution as Romulus25 (Tremmel et al. 2017). With unprecedented mass and spatial resolution, RomulusC represents a unique opportunity to study the evolution of galaxies in dense environments down to dwarf masses. We demonstrate that RomulusC results in an intracluster medium (ICM) consistent with observations. The star formation history and stellar mass of the brightest cluster galaxy (BCG) is consistent with observations and abundance matching results, indicating that our sub-grid models, optimized only to reproduce observations of field dwarf and Milky Way mass galaxies, are able to produce reasonable galaxy masses and star formation histories in much higher mass systems. Feedback from supermassive black holes (SMBHs) regulates star formation by driving large-scale, collimated outflows that coexist with a low entropy core. We find that non-BCG cluster member galaxies are substantially quenched compared to the field down to dwarf galaxy masses and, at low masses, quenching is seen to have no dependence on mass or distance from the cluster center. This enhanced quenched population extends beyond $R_{200}$ and is in place at high redshift. Similarly, we predict that SMBH activity is significantly suppressed within clusters outside of the BCG, but show how the effect could be lost when only focusing on the brightest AGN in the most massive galaxies.

Filaments of the cosmic web have long been associated with the threadlike structures seen in galaxy redshift surveys. However, despite their baryon content being dominated by hot gas, these filaments have been an elusive target for X-ray observations. Recently, detections of filaments in very deep (2.4 Msec) observations with Chandra were reported around Abell 133 (z=0.0559). To verify these claims, we conducted a multi-object spectrographic campaign on the Baade 6.5m telescope around Abell 133; this resulted in a catalog of ${\sim}3000$ new redshift measurements, of which 254 are of galaxies near the cluster. We investigate the kinematic state of Abell 133 and identify the physical locations of filamentary structure in the galaxy distribution. Contrary to previous studies, we see no evidence that Abell 133 is dynamically disturbed; we reject the hypothesis that there is a kinematically distinct subgroup (p=0.28) and find no velocity offset between the central galaxy and the cluster ($\textrm{Z}_\textrm{score}=0.041^{+0.111}_{-0.106}$). The spatial distribution of galaxies traces the X-ray filaments, as confirmed by angular cross correlation with a significance of ${\sim}5\sigma$. A similar agreement is found in the angular density distribution, where two X-ray structures have corresponding galaxy enhancements. We also identify filaments in the large-scale structure of galaxies; these filaments approach the cluster from the direction the X-ray structures are seen. While more members between $\textrm{R}_{200}$ and $2\times\textrm{R}_{200}$ are required to clarify which large scale filaments connect to the X-ray gas, we argue that this is compelling evidence that the X-ray emission is indeed associated with cosmic filaments.

The diffuse hot medium inside clusters of galaxies typically exhibits turbulent motions whose amplitude increases with radius, as revealed by cosmological hydrodynamical simulations. However, its physical origin remains unclear. It could either be due to an excess injection of turbulence at large radii, or faster turbulence dissipation at small radii. We investigate this by studying the time evolution of turbulence in the intracluster medium (ICM) after major mergers, using the Omega500 non-radiative hydrodynamical cosmological simulations. By applying a novel wavelet analysis to study the radial dependence of the ICM turbulence spectrum, we discover that faster turbulence dissipation in the inner high density regions leads to the increasing turbulence amplitude with radius. We also find that the ICM turbulence at all radii decays in two phases after a major merger: an early fast decay phase followed by a slow secular decay phase. The buoyancy effects resulting from the ICM density stratification becomes increasingly important during turbulence decay, as revealed by a decreasing turbulence Froude number $Fr \sim \mathcal{O}(1)$. Our results indicate that the stronger density stratification and smaller eddy turn-over time are the likely causes of the faster turbulence dissipation rate in the inner regions of the cluster.

Gas motions in galaxy clusters play important roles in determining the properties of the intracluster medium (ICM) and in the constraint of cosmological parameters via X-ray and Sunyaev-Zel'dovich effect observations of galaxy clusters. The Hitomi measurements of gas motions in the core of the Perseus Cluster have provided insights into the physics in galaxy clusters. The XARM mission, equipped with the Resolve X-ray micro-calorimeter, will continue Hitomi's legacy by measuring ICM motions through Doppler shifting and broadening of emission lines in a larger number of galaxy clusters, and at larger radii. In this work, we investigate how well we can measure bulk and turbulent gas motions in the ICM with XARM, by analyzing mock XARM simulations of galaxy clusters extracted from cosmological hydrodynamic simulations. We assess how photon counts, spectral fitting methods, multiphase ICM structure, deprojections, and region selection affect the measurements of gas motions. We first show that XARM is capable of recovering the underlying spherically averaged turbulent and bulk velocity profiles for dynamically relaxed clusters to within $\sim 50\%$ with a reasonable amount of photon counts in the X-ray emission lines. We also find that there are considerable azimuthal variations in the ICM velocities, where the velocities measured in a single azimuthal direction can significantly deviate from the true value even in dynamically relaxed systems. Such variation must be taken into account when interpreting data and developing observing strategies. We will discuss the prospect of using the upcoming XARM mission to measure non-thermal pressure and to correct for the hydrostatic mass bias of galaxy clusters. Our results are broadly applicable for future X-ray missions, such as Athena and Lynx.

Gravitational lensing is a powerful probe of the mass distribution of galaxy clusters and cosmology. However, accurate measurements of the cluster mass profiles are limited by uncertainties in cluster astrophysics. In this work, we present a physically motivated model of baryonic effects on the cluster mass profiles, which self-consistently takes into account the impact of baryons on the concentration as well as mass accretion histories of galaxy clusters. We calibrate this model using the Omega500 hydrodynamical cosmological simulations of galaxy clusters with varying baryonic physics. Our model will enable us to simultaneously constrain cluster mass, concentration, and cosmological parameters using stacked weak lensing measurements from upcoming optical cluster surveys.

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.

Recent detections of the cross-correlation of the thermal Sunyaev-Zel'dovich (tSZ) effect and weak gravitational lensing (WL) enable unique studies of cluster astrophysics and cosmology. In this work, we present constraints on the amplitude of the non-thermal pressure fraction in galaxy clusters, $\alpha_0$, and the amplitude of the matter power spectrum, $\sigma_8$, using measurements of the tSZ power spectrum from Planck, and the tSZ-WL cross-correlation from Planck and the Red Cluster Sequence Lensing Survey. We fit the data to a semi-analytic model with the covariance matrix using $N$-body simulations. We find that the tSZ power spectrum alone prefers $\sigma_8 \sim 0.85$ and a large fraction of non-thermal pressure ($\alpha_0 \sim 0.2$-$0.3$). The tSZ-WL cross-correlation on the other hand prefers a significantly lower $\sigma_8 \sim 0.6$, and low $\alpha_0 \sim 0.05$. We show that this tension can be mitigated by allowing for a steep slope in the stellar-mass-halo-mass relation, which would cause a reduction of the gas in low-mass halos. In such a model, the combined data prefer $\sigma_8 \sim 0.7$ and $\alpha_0 \sim 0.2$, consistent with predictions from hydrodynamical simulations.

We present the "Galaxy Cluster Merger Catalog." This catalog provides an extensive suite of mock observations and related data for N-body and hydrodynamical simulations of galaxy cluster mergers and clusters from cosmological simulations. These mock observations consist of projections of a number of important observable quantities in several different wavebands as well as along different lines of sight through each simulation domain. The web interface to the catalog consists of easily browseable images over epoch and projection direction, as well as download links for the raw data and a JS9 interface for interactive data exploration. The data is presented within a consistent format so that comparison between simulations is straightforward. All of the data products are provided in the standard FITS file format. Data is being stored on the yt Hub (http://hub.yt), which allows for remote access and analysis using a Jupyter notebook server. Future versions of the catalog will include simulations from a number of research groups and a variety of research topics related to the study of interactions of galaxy clusters with each other and with their member galaxies. The catalog is located at http://gcmc.hub.yt.

The Hitomi X-ray satellite has provided the first direct measurements of the plasma velocity dispersion in a galaxy cluster. It finds a relatively "quiescent" gas with a line-of-sight velocity dispersion ~ 160 km/s, at 30 kpc to 60 kpc from the cluster center. This is surprising given the presence of jets and X-ray cavities that indicates on-going activity and feedback from the active galactic nucleus (AGN) at the cluster center. Using a set of mock Hitomi observations generated from a suite of state-of-the-art cosmological cluster simulations, and an isolated but higher resolution simulation of gas physics in the cluster core, including the effects of cooling and AGN feedback, we examine the likelihood of Hitomi detecting a cluster with the observed velocities. As long as the Perseus has not experienced a major merger in the last few gigayears, and AGN feedback is operating in a "gentle" mode, we reproduce the level of gas motions observed by Hitomi. The frequent mechanical AGN feedback generates net line-of-sight velocity dispersions ~100-200 km/s, bracketing the values measured in the Perseus core. The large-scale velocity shear observed across the core, on the other hand, is generated mainly by cosmic accretion such as mergers. We discuss the implications of these results for AGN feedback physics and cluster cosmology and progress that needs to be made in both simulations and observations, including a Hitomi re-flight and calorimeter-based instruments with higher spatial resolution.

We present integral field spectroscopy of galaxy cluster Abell 3827, using ALMA and VLT/MUSE. It reveals an unusual configuration of strong gravitational lensing in the cluster core, with at least seven lensed images of a single background spiral galaxy. Lens modelling based on HST imaging had suggested that the dark matter associated with one of the cluster's central galaxies may be offset. The new spectroscopic data enable better subtraction of foreground light, and better identification of multiple background images. The inferred distribution of dark matter is consistent with being centered on the galaxies, as expected by LCDM. Each galaxy's dark matter also appears to be symmetric. Whilst we do not find an offset between mass and light (suggestive of self-interacting dark matter) as previously reported, the numerical simulations that have been performed to calibrate Abell 3827 indicate that offsets and asymmetry are still worth looking for in collisions with particular geometries. Meanwhile, ALMA proves exceptionally useful for strong lens image identifications.

We report on the detection of three strong HI absorbers originating in the outskirts (i.e., impact parameter, $\rho_{\rm cl} \approx (1.6-4.7) r_{500}$) of three massive ($M_{500}\sim3\times10^{14} M_{\odot}$) clusters of galaxies at redshift $z_{\rm cl} \approx 0.46$, in the $Hubble Space Telescope$ Cosmic Origins Spectrograph ($HST$/COS) spectra of 3 background UV-bright quasars. These clusters were discovered by the 2500 deg$^2$ South Pole Telescope Sunyaev$-$Zel'dovich (SZ) effect survey. All three COS spectra show partial Lyman limit absorber with $N(HI) > 10^{16.5} \ \rm cm^{-2}$ near the photometric redshifts ($|\Delta z/(1+z)| \approx 0.03$) of the clusters. The compound probability of random occurrence of all three absorbers is $<0.02$%, indicating that the absorbers are most likely related to the targeted clusters. We find that the outskirts of these SZ-selected clusters are remarkably rich in cool gas compared to existing observations of other clusters in the literature. The effective Doppler parameters of the Lyman series lines, obtained using single cloud curve-of-growth (COG) analysis, suggest a non-thermal/turbulent velocity of a few $\times10 \ \rm km s^{-1}$ in the absorbing gas. We emphasize the need for uniform galaxy surveys around these fields and for more UV observations of QSO-cluster pairs in general in order to improve the statistics and gain further insights into the unexplored territory of the largest collapsed cosmic structures.

We present a quantitative study of the X-ray morphology of galaxy clusters, as a function of their detection method and redshift. We analyze two separate samples of galaxy clusters: a sample of 36 clusters at 0.35 < z < 0.9 selected in the X-ray with the ROSAT PSPC 400 deg2 survey, and a sample of 90 clusters at 0.25 < z < 1.2 selected via the Sunyaev-Zel'dovich (SZ) effect with the South Pole Telescope. Clusters from both samples have similar-quality Chandra observations, which allow us to quantify their X-ray morphologies via two distinct methods: centroid shifts and photon asymmetry. The latter technique provides nearly unbiased morphology estimates for clusters spanning a broad range of redshift and data quality. We further compare the X-ray morphologies of X-ray- and SZ-selected clusters with those of simulated clusters. We do not find a statistically significant difference in the measured X-ray morphology of X-ray and SZ-selected clusters over the redshift range probed by these samples, suggesting that the two are probing similar populations of clusters. We find that the X-ray morphologies of simulated clusters are statistically indistinguishable from those of X-ray- or SZ-selected clusters, implying that the most important physics for dictating the large-scale gas morphology (outside of the core) is well-approximated in these simulations. Finally, we find no statistically significant redshift evolution in the X-ray morphology (both for observed and simulated clusters), over the range z ~ 0.3 to z ~ 1, seemingly in contradiction with the redshift-dependent halo merger rate predicted by simulations.

We present a $Chandra$ study of the hot intragroup medium (hIGM) of the galaxy group NCG2563. The $Chandra$ mosaic observations, with a total exposure time of ~430 ks, allow the gas density to be detected beyond $R_{200}$ and the gas temperature out to 0.75 $R_{200}$. This represents the first observational measurement of the physical properties of a poor groups beyond $R_{500}$. By capitalizing on the exquisite spatial resolution of $Chandra$ that is capable to remove unrelated emission from point sources and substructures, we are able to radially constrain the inhomogeneities of gas ("clumpiness"), gas fraction, temperature and entropy distribution. Although there is some uncertainty in the measurements, we find evidences of gas clumping in the virialization region, with clumping factor of about 2 - 3 at $R_{200}$. The gas clumping-corrected gas fraction is significantly lower than the cosmological baryon budget. These results may indicate a larger impact of the gas inhomogeneities with respect to the prediction from hydrodynamic numerical simulations, and we discuss possible explanations for our findings.

Future data from galaxy redshift surveys, combined with high-resolutions maps of the cosmic microwave background, will enable measurements of the pairwise kinematic Sunyaev-Zel'dovich (kSZ) signal with unprecedented statistical significance. This signal probes the matter-velocity correlation function, scaled by the average optical depth ($\tau$) of the galaxy groups and clusters in the sample, and is thus of fundamental importance for cosmology. However, in order to translate pairwise kSZ measurements into cosmological constraints, external constraints on $\tau$ are necessary. In this work, we present a new model for the intra-cluster medium, which takes into account star-formation, feedback, non-thermal pressure, and gas cooling. Our semi-analytic model is computationally efficient and can reproduce results of recent hydrodynamical simulations of galaxy cluster formation. By calibrating the model using recent X-ray measurements of gas density profiles of clusters and $M_{\mathrm{gas}}-M$ relations of groups and clusters, we show that the integrated $\tau$ is constrained to better than 6% (with a confidence level of 95%), which sets the tightest constraints on the external prior on the optical depth of groups and clusters reported in the literature to date. Our strong prior on $\tau$ should significantly improve cosmological constraints derived from future pairwise kSZ measurements and shed new insights into the nature of dark energy, modified gravity, and neutrino mass.

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.

Primordial dark matter (DM) haloes are the smallest gravitationally bound DM structures from which the first stars, black holes, and galaxies form and grow in the early universe. However, their structures are sensitive to the free streaming scale of DM, which in turn depends on the nature of DM particles. In this work, we test the hypothesis that the slope of the central cusps in primordial DM haloes near the free streaming scale depends on the nature of merging process. By combining and analysing data from a cosmological simulation with the cutoff in the small-scale matter power spectrum as well as a suite of controlled, high-resolution simulations of binary mergers, we find that (1) the primordial DM haloes form preferentially through major mergers in radial orbits; (2) their central DM density profile is more susceptible to a merging process compared to that of galaxy and cluster-size DM haloes; (3) consecutive major mergers drive the central density slope to approach the {\it universal} form characterized by the NFW profile, which is shown to be robust to the impacts of mergers and serves an attractor solution for the density structure of DM haloes. Our work highlights the importance of dynamical processes on the structure formation during the Dark Ages.

Context. Galaxy clusters are continuously growing through the accretion of matter in their outskirts. This process induces inhomogeneities in the gas density distribution (clumping) which need to be taken into account to recover the physical properties of the intracluster medium (ICM) at large radii. Aims. We studied the thermodynamic properties in the outskirts (R > R500) of the massive galaxy cluster Abell 2142 by combining the Sunyaev Zel'dovich (SZ) effect with the X-ray signal. Methods. We combined the SZ pressure profile measured by Planck with the XMM-Newton gas density profile to recover radial profiles of temperature, entropy and hydrostatic mass out to 2R500. We used a method that is insensitive to clumping to recover the gas density, and we compared the results with traditional X-ray measurement techniques. Results. When taking clumping into account, our joint SZ/X-ray entropy profile is consistent with the predictions from pure gravitational collapse, whereas a significant entropy flattening is found when the effect of clumping is neglected. The hydrostatic mass profile recovered using joint X-ray/SZ data agrees with that obtained from spectroscopic X-ray measurements and with mass reconstructions obtained through weak lensing and galaxy kinematics. Conclusions. We found that clumping can explain the entropy flattening observed by Suzaku in the outskirts of several clusters. When using a method insensitive to clumping for the reconstruction of the gas density, the thermodynamic properties of Abell 2142 are compatible with the assumption that the thermal gas pressure sustains gravity and that the entropy is injected at accretion shocks, with no need to evoke more exotic physics. Our results highlight the need for X-ray observations with sufficient spatial resolution, and large collecting area, to understand the processes at work in cluster outer regions.

The thermal Sunyaev-Zel'dovich (tSZ) effect signal is widely recognized as a robust mass proxy of galaxy clusters with small intrinsic scatter. However, recent observational calibration of the tSZ scaling relation using weak lensing (WL) mass exhibits considerably larger scatter than the intrinsic scatter predicted from numerical simulations. This raises a question as to whether we can realize the full statistical power of ongoing and upcoming tSZ-WL observations of galaxy clusters. In this work, we investigate the origin of observed scatter in the tSZ-WL scaling relation, using mock maps of galaxy clusters extracted from cosmological hydrodynamic simulations. We show that the inferred intrinsic scatter from mock tSZ-WL analyses is considerably larger than the intrinsic scatter measured in simulations, and comparable to the scatter in the observed tSZ-WL relation. We show that this enhanced scatter originates from the combination of the projection of correlated structures along the line of sight and the uncertainty in the cluster radius associated with WL mass estimates, causing the amplitude of the scatter to depend on the covariance between tSZ and WL signals. We present a statistical model to recover the unbiased cluster scaling relation and cosmological parameter by taking into account the covariance in the tSZ-WL mass relation from multi-wavelength cluster surveys.

Using deep Hubble and Spitzer observations Oesch et al. (2016) have identified a bright ($M_{\rm UV}\approx -22$) star forming galaxy candidate at $z \approx 11$. The presence of GN-$z11$ implies a number density $\sim 10^{-6}\,{\rm Mpc^{-3}}$, roughly an order of magnitude higher than the expected value based on extrapolations from lower redshift. Using the unprecedented volume and high resolution of the BlueTides cosmological hydrodynamical simulation, we study the population of luminous rare objects at $z > 10$. The luminosity function in BlueTides implies an enhanced number of massive galaxies, consistent with the observation of GN-$z11$. We find about 30 galaxies at $M_{\rm UV}\approx -22$ at $z = 11$ in the BlueTides volume, including a few objects about 1.5 magnitudes brighter. The probability of observing GN-$z11$ in the volume probed by Oesch et al. (2016) is $\sim 13$ per cent. The predicted properties of the rare bright galaxies at $z = 11$ in BlueTides closely match those inferred from the observations of GN-$z11$. BlueTides predicts a negligible contribution from faint AGN in the observed SED. The enormous increase in volume surveyed by WFIRST will provide observations of $\sim1000$ galaxies with $M_{\rm UV} < -22$ beyond $z = 11$ out to $z = 13.5$.

Non-thermal pressure in galaxy clusters leads to underestimation of the mass of galaxy clusters based on hydrostatic equilibrium with thermal gas pressure. This occurs even for dynamically relaxed clusters that are used for calibrating the mass-observable scaling relations. We show that the analytical model for non-thermal pressure developed in Shi & Komatsu 2014 can correct for this so-called 'hydrostatic mass bias', if most of the non-thermal pressure comes from bulk and turbulent motions of gas in the intracluster medium. Our correction works for the sample average irrespective of the mass estimation method, or the dynamical state of the clusters. This makes it possible to correct for the bias in the hydrostatic mass estimates from X-ray surface brightness and the Sunyaev-Zel'dovich observations that will be available for clusters in a wide range of redshifts and dynamical states.

We present the analysis of a local (z = 0.04 - 0.2) sample of 31 galaxy clusters with the aim of measuring the density of the X-ray emitting gas in cluster outskirts. We compare our results with numerical simulations to set constraints on the azimuthal symmetry and gas clumping in the outer regions of galaxy clusters. We exploit the large field-of-view and low instrumental background of ROSAT/PSPC to trace the density of the intracluster gas out to the virial radius. We perform a stacking of the density profiles to detect a signal beyond r200 and measure the typical density and scatter in cluster outskirts. We also compute the azimuthal scatter of the profiles with respect to the mean value to look for deviations from spherical symmetry. Finally, we compare our average density and scatter profiles with the results of numerical simulations. As opposed to some recent Suzaku results, and confirming previous evidence from ROSAT and Chandra, we observe a steepening of the density profiles beyond \sim r500. Comparing our density profiles with simulations, we find that non-radiative runs predict too steep density profiles, whereas runs including additional physics and/or treating gas clumping are in better agreement with the observed gas distribution. We report for the first time the high-confidence detection of a systematic difference between cool-core and non-cool core clusters beyond \sim 0.3r200, which we explain by a different distribution of the gas in the two classes. Beyond \sim r500, galaxy clusters deviate significantly from spherical symmetry, with only little differences between relaxed and disturbed systems. We find good agreement between the observed and predicted scatter profiles, but only when the 1% densest clumps are filtered out in the simulations. [Abridged]

The analysis of X-ray and Sunyaev-Zel\'dovich measurements of the intracluster medium (ICM) assumes that electrons are in thermal equilibrium with ions in the plasma. However, electron-ion equilibration timescales can be comparable to the Hubble time in the low density galaxy cluster outskirts, leading to differences between the electron and ion temperatures. This temperature difference can lead to systematic biases in cluster mass estimates and mass-observable scaling relations. To quantify the impact of non-equilibrium electrons on the ICM profiles in cluster outskirts, we use a high resolution cosmological simulation with a two-temperature model assuming the Spitzer equilibration timescale for the electrons. First, we show how the radial profile of this temperature bias depends on both the mass and mass accretion rate of the cluster; the bias is most pronounced in the most massive and most rapidly accreting clusters. For the most extreme case in our sample, we find that the bias is of order 10% at half of the cluster virial radius and increases to 40% at the edge of the cluster. We also find that gas in filaments is less susceptible to the non-equilibrium effect, leading to azimuthal variations at large cluster-centric radii. By analyzing mock Chandra observations of simulated clusters, we show that such azimuthal variations can be probed with deep X-ray observations. Finally, the mass-dependent temperature bias introduces biases in hydrostatic mass and cluster temperature, which has implications for cluster-based cosmological inferences. We provide a mass-dependent model for the temperature bias profile which can be useful for modeling the effect of electron-ion equilibration in galaxy clusters.

The cosmic accretion of both dark matter and baryons into halos typically is measured using some evolving virial relation, but recent work suggests that most halo growth at late cosmic time (z < 2) is not physical but is rather the by-product of an evolving virial radius ("pseudo-evolution"). Using Omega25, a suite of cosmological simulations that incorporate both dark matter and gas dynamics with differing treatments of gas cooling, star formation, and thermal feedback, we systematically explore the physics that governs cosmic accretion into halos and their galaxies. Physically meaningful cosmic accretion of both dark matter and baryons occurs at z > 1 across our halo mass range: M_200m = 10^{11-14} M_sun. However, dark matter, because it is dissipationless, is deposited (in a time-average sense) at > R_200m(z) in a shell-like manner, such that dark-matter mass and density experience little-to-no physical growth at any radius within a halo at z < 1. In contrast, gas, because it is able to cool radiatively, experiences significant accretion at all radii, at a rate that roughly tracks the accretion rate at R_200m, at all redshifts. Infalling gas starts to decouple from dark matter at ~2 R_200m and continues to accrete to smaller radii until the onset of strong angular-momentum support at ~ 0.1 R_200m. Thus, while the growth of dark matter is subject to pseudo-evolution, the growth of baryons is not. The fact that the accretion rate of gas on galactic scales tracks the accretion rate near R_200m provides insight into the tight relations between the masses/sizes of galaxies and those of their host halos across cosmic time.

The Sunyaev-Zel'dovich effect (SZE) observable-mass (Y-M) scaling relation is a promising technique for obtaining mass estimates for large samples of galaxy clusters and holds a key to studying the nature of dark matter and dark energy. However, cosmological inference based on SZE cluster surveys is limited by our incomplete knowledge of bias, scatter, and evolution in the Y-M relation. In this work, we investigate the effects of galaxy cluster mergers on the scaling relation using the Omega500 high-resolution cosmological hydrodynamic simulation. We show that the non-thermal pressure associated with merger-induced gas motions contributes significantly to the bias, scatter, and evolution of the scaling relation. After the merger, the kinetic energy of merging systems is slowly converted into thermal energy through dissipation of turbulent gas motions, which causes the thermal SZE signal to increase steadily with time. This post-merger evolution is one of the primary source of bias and scatter in the Y-M relation. However, we show that when the missing non-thermal energy is accounted for, the resulting relation exhibits little to no redshift evolution and the scatter around the scaling relation is ~20-30 % smaller than that of the thermal SZE signal alone. Our work opens up a possibility to further improve the current robust mass proxy, Y, by accounting for the missing non-thermal pressure component. We discuss future prospect of measuring internal gas motions in galaxy clusters and its implication for cluster-based cosmological tests.

Galaxy clusters exhibit remarkable self-similar behavior which allows us to establish simple scaling relationships between observable quantities and cluster masses, making galaxy clusters useful cosmological probes. Recent X-ray observations suggest that self-similarity may be broken in the outskirts of galaxy clusters. In this work, we analyze a mass-limited sample of massive galaxy clusters from the Omega500 cosmological hydrodynamic simulation to investigate the self-similarity of the diffuse X-ray emitting intracluster medium (ICM) in the outskirts of galaxy clusters. We find that the self-similarity of the outer ICM profiles is better preserved if they are normalized with respect to the mean density of the universe, while the inner profiles are more self-similar when normalized using the critical density. However, the outer ICM profiles as well as the location of accretion shock around clusters are sensitive to their mass accretion rate, which causes the apparent breaking of self-similarity in cluster outskirts. We also find that the collisional gas does not follow the distribution of collisionless dark matter perfectly in the infall regions of galaxy clusters, leading to 10% departures in the gas-to-dark matter density ratio from the cosmic mean value. Our results have a number implications for interpreting observations of galaxy clusters in X-ray and through the Sunyaev-Zel'dovich effect and their application to cluster cosmology.

Galaxy cluster Abell 3827 hosts the stellar remnants of four almost equally bright elliptical galaxies within a core of radius 10kpc. Such corrugation of the stellar distribution is very rare, and suggests recent formation by several simultaneous mergers. We map the distribution of associated dark matter, using new Hubble Space Telescope imaging and VLT/MUSE integral field spectroscopy of a gravitationally lensed system threaded through the cluster core. We find that each of the central galaxies retains a dark matter halo, but that (at least) one of these is spatially offset from its stars. The best-constrained offset is 1.62+/-0.48kpc, where the 68% confidence limit includes both statistical error and systematic biases in mass modelling. Such offsets are not seen in field galaxies, but are predicted during the long infall to a cluster, if dark matter self-interactions generate an extra drag force. With such a small physical separation, it is difficult to definitively rule out astrophysical effects operating exclusively in dense cluster core environments - but if interpreted solely as evidence for self-interacting dark matter, this offset implies a cross-section sigma/m=(1.7+/-0.7)x10^{-4}cm^2/g x (t/10^9yrs)^{-2}, where t is the infall duration.

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.

Turbulent gas motion inside galaxy clusters provides a non-negligible non-thermal pressure support to the intracluster gas. If not corrected, it leads to a systematic bias in the estimation of cluster masses from X-ray and Sunyaev-Zel'dovich (SZ) observations assuming hydrostatic equilibrium, and affects interpretation of measurements of the SZ power spectrum and observations of cluster outskirts from ongoing and upcoming large cluster surveys. Recently, Shi & Komatsu developed an analytical model for predicting the radius, mass, and redshift dependence of the non-thermal pressure contributed by the kinetic random motions of intracluster gas sourced by the cluster mass growth. In this paper, we compare the predictions of this analytical model to a state-of-the-art cosmological hydrodynamics simulation. As different mass growth histories result in different non-thermal pressure, we perform the comparison on 65 simulated galaxy clusters on a cluster-by-cluster basis. We find an excellent agreement between the modelled and simulated non-thermal pressure profiles. Our results open up the possibility of using the analytical model to correct the systematic bias in the mass estimation of galaxy clusters. We also discuss tests of the physical picture underlying the evolution of intracluster non-thermal gas motions, as well as a way to further improve the analytical modeling, which may help achieve a unified understanding of non-thermal phenomena in galaxy clusters.

The study of galaxy cluster outskirts has emerged as one of the new frontiers in extragalactic astrophysics and cosmology with the advent of new observations in X-ray and microwave. However, the thermodynamic properties and chemical enrichment of this diffuse and azimuthally asymmetric component of the intracluster medium (ICM) are still not well understood. This work, for the first time, systematically explores potential observational biases in these regions. To assess X-ray measurements of galaxy cluster properties at large radii ($>{R}_{500c}$), we use mock Chandra analyses of cosmological galaxy cluster simulations. The pipeline is identical to that used for Chandra observations, but the biases discussed in this paper are relevant for all X-ray observations outside of ${R}_{500c}$. We find the following from our analysis: (1) filament regions can contribute as much as $50\%$ at $R_{200c}$ to the emission measure; (2) X-ray temperatures and metal abundances from model fitted mock X-ray spectra in a multi-temperature ICM respectively vary to the level of $10\%$ and $50\%$; (3) resulting density profiles vary to within $10\%$ out to $R_{200c}$, and gas mass, total mass, and baryon fractions all vary to within a few percent; (4) the bias from a metal abundance extrapolated a factor of five higher than the true metal abundance results in total mass measurements biased high by $20\%$ and total gas measurements biased low by $10\%$; and (5) differences in projection and dynamical state of a cluster can lead to gas density slope measurements that differ by a factor of $15\%$ and $30\%$, respectively. The presented results can partially account for some of the recent gas profile measurements in cluster outskirts by, e.g., Suzaku. Our findings are pertinent to future X-ray cosmological constraints from cluster outskirts.

Cosmological constraints from X-ray and microwave observations of galaxy clusters are subjected to systematic uncertainties. Non-thermal pressure support due to internal gas motions in galaxy clusters is one of the major sources of astrophysical uncertainties. Using a mass-limited sample of galaxy clusters from a high-resolution hydrodynamical cosmological simulation, we characterize the non-thermal pressure fraction profile and study its dependence on redshift, mass, and mass accretion rate. We find that the non-thermal pressure fraction profile is universal across redshift when galaxy cluster radii are defined with respect to the mean matter density of the universe instead of the commonly used critical density. We also find that the non-thermal pressure is predominantly radial, and the gas velocity anisotropy profile exhibits strong universality when galaxy cluster radii are defined with respect to the mean matter density of the universe. However, we find that the non-thermal pressure fraction is strongly dependent on the mass accretion rate of the galaxy cluster. We provide fitting formulae for the universal non-thermal pressure fraction and velocity anisotropy profiles of gas in galaxy clusters, which should be useful in modeling astrophysical uncertainties pertinent to using galaxy clusters as cosmological probes.

Analyses of cosmological hydrodynamic simulations of galaxy clusters suggest that X-ray masses can be underestimated by 10% to 30%. The largest bias originates by both violation of hydrostatic equilibrium and an additional temperature bias caused by inhomogeneities in the X-ray emitting intra-cluster medium (ICM). To elucidate on this large dispersion among theoretical predictions, we evaluate the degree of temperature structures in cluster sets simulated either with smoothed-particle-hydrodynamics (SPH) and adaptive-mesh-refinement (AMR) codes. We find that the SPH simulations produce larger temperature variations connected to the persistence of both substructures and their stripped cold gas. This difference is more evident in no-radiative simulations, while it is reduced in the presence of radiative cooling. We also find that the temperature variation in radiative cluster simulations is generally in agreement with the observed one in the central regions of clusters. Around R_500 the temperature inhomogeneities of the SPH simulations can generate twice the typical hydrostatic-equilibrium mass bias of the AMR sample. We emphasize that a detailed understanding of the physical processes responsible for the complex thermal structure in ICM requires improved resolution and high sensitivity observations in order to extend the analysis to higher temperature systems and larger cluster-centric radii.

Exploring the ICM power spectrum can help us to probe the physics of galaxy clusters. Using high-resolution 3D plasma simulations, we study the statistics of the velocity field and its relation with the thermodynamic perturbations. The normalization of the ICM spectrum (density, entropy, or pressure) is linearly tied to the level of large-scale motions, which excite both gravity and sound waves due to stratification. For low 3D Mach number M~0.25, gravity waves mainly drive entropy perturbations, traced by preferentially tangential turbulence. For M>0.5, sound waves start to significantly contribute, passing the leading role to compressive pressure fluctuations, associated with isotropic (or slightly radial) turbulence. Density and temperature fluctuations are then characterized by the dominant process: isobaric (low M), adiabatic (high M), or isothermal (strong conduction). Most clusters reside in the intermediate regime, showing a mixture of gravity and sound waves, hence drifting towards isotropic velocities. Remarkably, regardless of the regime, the variance of density perturbations is comparable to the 1D Mach number. This linear relation allows to easily convert between gas motions and ICM perturbations, which can be exploited by Chandra, XMM data and by the forthcoming Astro-H. At intermediate and small scales (10-100 kpc), the turbulent velocities develop a Kolmogorov cascade. The thermodynamic perturbations act as effective tracers of the velocity field, broadly consistent with the Kolmogorov-Obukhov-Corrsin advection theory. Thermal conduction acts to damp the gas fluctuations, washing out the filamentary structures and steepening the spectrum, while leaving unaltered the velocity cascade. The ratio of the velocity and density spectrum thus inverts the downtrend shown by the non-diffusive models, allowing to probe the presence of significant conductivity in the ICM.

We address the problem of evaluating the power spectrum of the velocity field of the ICM using only information on the plasma density fluctuations, which can be measured today by Chandra and XMM-Newton observatories. We argue that for relaxed clusters there is a linear relation between the rms density and velocity fluctuations across a range of scales, from the largest ones, where motions are dominated by buoyancy, down to small, turbulent scales: $(\delta\rho_k/\rho)^2 = \eta_1^2 (V_{1,k}/c_s)^2$, where $\delta\rho_k/\rho$ is the spectral amplitude of the density perturbations at wave number $k$, $V_{1,k}^2=V_k^2/3$ is the mean square component of the velocity field, $c_s$ is the sound speed, and $\eta_1$ is a dimensionless constant of order unity. Using cosmological simulations of relaxed galaxy clusters, we calibrate this relation and find $\eta_1\approx 1 \pm 0.3$. We argue that this value is set at large scales by buoyancy physics, while at small scales the density and velocity power spectra are proportional because the former are a passive scalar advected by the latter. This opens an interesting possibility to use gas density power spectra as a proxy for the velocity power spectra in relaxed clusters, across a wide range of scales.

The use of galaxy clusters as cosmological probes hinges on our ability to measure their masses accurately and with high precision. Hydrostatic mass is one of the most common methods for estimating the masses of individual galaxy clusters, which suffer from biases due to departures from hydrostatic equilibrium. Using a large, mass-limited sample of massive galaxy clusters from a high-resolution hydrodynamical cosmological simulation, in this work we show that in addition to turbulent and bulk gas velocities, acceleration of gas introduces biases in the hydrostatic mass estimate of galaxy clusters. In unrelaxed clusters, the acceleration bias is comparable to the bias due to non-thermal pressure associated with merger-induced turbulent and bulk gas motions. In relaxed clusters, the mean mass bias due to acceleration is small (<3%), but the scatter in the mass bias can be reduced by accounting for gas acceleration. Additionally, this acceleration bias is greater in the outskirts of higher redshift clusters where mergers are more frequent and clusters are accreting more rapidly. Since gas acceleration cannot be observed directly, it introduces an irreducible bias for hydrostatic mass estimates. This acceleration bias places limits on how well we can recover cluster masses from future X-ray and microwave observations. We discuss implications for cluster mass estimates based on X-ray, Sunyaev-Zeldovich effect, and gravitational lensing observations and their impact on cluster cosmology.

We develop a statistical method to measure the interaction cross-section of Dark Matter, exploiting the continuous minor merger events in which small substructures fall into galaxy clusters. We find that by taking the ratio of the distances between the galaxies and Dark Matter, and galaxies and gas in accreting sub-halos, we form a quantity that can be statistically averaged over a large sample of systems whilst removing any inherent line-of-sight projections. In order to interpret this ratio as a cross-section of Dark Matter we derive an analytical description of sub-halo infall which encompasses; the force of the main cluster potential, the drag on a gas sub-halo, a model for Dark Matter self-interactions and the resulting sub-halo drag, the force on the gas and galaxies due to the Dark Matter sub-halo potential, and finally the buoyancy on the gas and Dark Matter. We create mock observations from cosmological simulations of structure formation and find that collisionless Dark Matter becomes physically separated from X-ray gas by up to 20h^-1 kpc. Adding realistic levels of noise, we are able to predict achievable constraints from observational data. Current archival data should be able to detect a difference in the dynamical behaviour of Dark Matter and standard model particles at 6 sigma, and measure the total interaction cross-section sigma/m with 68% confidence limits of +/- 1cm2g^-1. We note that this method is not restricted by the limited number of major merging events and is easily extended to large samples of clusters from future surveys which could potentially push statistical errors to 0.1cm^2g^-1.

The lensing power spectrum from cosmic microwave background (CMB) temperature maps will be measured with unprecedented precision with upcoming experiments, including upgrades to ACT and SPT. Achieving significant improvements in cosmological parameter constraints, such as percent level errors on sigma_8 and an uncertainty on the total neutrino mass of approximately 50 meV, requires percent level measurements of the CMB lensing power. This necessitates tight control of systematic biases. We study several types of biases to the temperature-based lensing reconstruction signal from foreground sources such as radio and infrared galaxies and the thermal Sunyaev-Zel'dovich effect from galaxy clusters. These foregrounds bias the CMB lensing signal due to their non-Gaussian nature. Using simulations as well as some analytical models we find that these sources can substantially impact the measured signal if left untreated. However, these biases can be brought to the percent level if one masks galaxies with fluxes at 150 GHz above 1 mJy and galaxy clusters with masses above M_vir = 10^14 M_sun. To achieve such percent level bias, we find that only modes up to a maximum multipole of l_max ~ 2500 should be included in the lensing reconstruction. We also discuss ways to minimize additional bias induced by such aggressive foreground masking by, for example, exploring a two-step masking and in-painting algorithm.

In the hierarchical structure formation model, clusters of galaxies form through a sequence of mergers and continuous mass accretion, which generate significant random gas motions especially in their outskirts where material is actively accreting. Non-thermal pressure provided by the internal gas motions affects the thermodynamic structure of the X-ray emitting intracluster plasma and introduces biases in the physical interpretation of X-ray and Sunyaev-Zeldovich effect observations. However, we know very little about the nature of gas motions in galaxy clusters. The ASTRO-H X-ray mission, scheduled to launch in 2015, will have a calorimeter capable of measuring gas motions in galaxy clusters at the level of <100 km/s. In this work, we predict the level of merger-induced gas motions expected in the Lambda-CDM model using hydrodynamical simulations of galaxy cluster formation. We show that the gas velocity dispersion is larger in more massive clusters, but exhibits a large scatter. We show that systems with large gas motions are morphologically disturbed, while early forming, relaxed groups show a smaller level of gas motions. By analyzing mock ASTRO-H observations of simulated clusters, we show that such observations can accurately measure the gas velocity dispersion out to the outskirts of nearby relaxed galaxy clusters. ASTRO-H analysis of merging clusters, on the other hand, requires multi-component spectral fitting and enables unique studies of substructures in galaxy clusters by measuring both the peculiar velocities and the velocity dispersion of gas within individual sub-clusters.

Mass estimates of galaxy clusters from X-ray and Sunyeav-Zel'dovich observations assume the intracluster gas is in hydrostatic equilibrium with their gravitational potential. However, since galaxy clusters are dynamically active objects whose dynamical states can deviate significantly from the equilibrium configuration, the departure from the hydrostatic equilibrium assumption is one of the largest sources of systematic uncertainties in cluster cosmology. In the literature there has been two methods for computing the hydrostatic mass bias based on the Euler and the modified Jeans equations, respectively, and there has been some confusion about the validity of these two methods. The word "Jeans" was a misnomer, which incorrectly implies that the gas is collisionless. To avoid further confusion, we instead refer these methods as "summation" and "averaging" methods respectively. In this work, we show that these two methods for computing the hydrostatic mass bias are equivalent by demonstrating that the equation used in the second method can be derived from taking spatial averages of the Euler equation. Specifically, we identify the correspondences of individual terms in these two methods mathematically and show that these correspondences are valid to within a few percent level using hydrodynamical simulations of galaxy cluster formation. In addition, we compute the mass bias associated with the acceleration of gas and show that its contribution is small in the virialized regions in the interior of galaxy cluster, but becomes non-negligible in the cluster outskirts of massive clusters where materials are being actively accreted. We discuss implications for understanding the origin of hydrostatic mass biases and for characterizing them using hydrodynamical simulations as well as observations.

We present a non-parametric method to measure inhomogeneities in the intracluster medium (ICM) from X-ray observations of galaxy clusters. Analysing mock Chandra X-ray observations of simulated clusters, we show that our new method enables the accurate recovery of the 3D gas density and gas clumping factor profiles out to large radii of galaxy clusters. We then apply this method to Chandra X-ray observations of Abell 1835 and present the first determination of the gas clumping factor from the X-ray cluster data. We find that the gas clumping factor in Abell 1835 increases with radius and reaches ~2-3 at r=R_{200}. This is in good agreement with the predictions of hydrodynamical simulations, but it is significantly below the values inferred from recent Suzaku observations. We further show that the radially increasing gas clumping factor causes flattening of the derived entropy profile of the ICM and affects physical interpretation of the cluster gas structure, especially at the large cluster-centric radii. Our new technique should be useful for improving our understanding of the cluster structure and to advance the use of galaxy clusters as cosmological probes, by helping to exploit rich data sets provided by Chandra and XMM-Newton X-ray space telescopes.

The Sunyaev-Zeldovich (SZ) effect is a promising tool to study physical properties of the hot X-ray emitting intracluster medium (ICM) in galaxy clusters. To date, most SZ observations have been interpreted in combination with X-ray follow-up measurements in order to determine the ICM temperature and estimate the cluster mass. Future high-resolution, multifrequency SZ observations promise to enable detailed studies of the ICM structures, by measuring the ICM temperature through the temperature-dependent relativistic corrections. In this work we develop a non-parametric method to derive three-dimensional physical quantities, including temperature, pressure, total mass, and peculiar velocities, of galaxy clusters from SZ observations alone. We test the performance of this method using hydrodynamical simulations of galaxy clusters, in order to assess systematic uncertainties in the reconstructed physical parameters. In particular, we analyze mock Cerro Chajnantor Atacama Telescope (CCAT) SZ observations, taking into account various sources of systematic uncertainties associated with instrumental effects and astrophysical foregrounds. We show that our method enables accurate reconstruction of the three-dimensional ICM profiles, while retaining full information about the gas distribution. We discuss the application of this technique for ongoing and future multifrequency SZ observations.

X-ray observations of galaxy clusters provide emission measure weighted spectra, arising from a range of density and temperature fluctuations in the intra-cluster medium (ICM). This is fitted to a single temperature plasma emission model to provide an estimate of the gas density and temperature, which are sensitive to the gas inhomogeneities. Therefore, X-ray observations yield a potentially biased estimate of the thermal gas pressure, P_X. At the same time Sunyaev-Zeldovich (SZ) observations directly measure the integrated gas pressure, P_SZ. If the X-ray pressure profiles are strongly biased with respect to the SZ, then one has the possibility to probe the gas inhomogeneities, even at scales unresolved by the current generation of telescopes. At the same time, a weak bias has implications for the interchangeable use of mass proxies like Y_SZ and Y_X as cosmological probes. In this paper we investigate the dependence of the bias, defined as b_P(r)=P_X(r)/P_SZ(r)-1, on the characteristics of fluctuations in the ICM taking into account the correlation between temperature and density fluctuations. We made a simple prediction of the irreducible bias in idealised X-ray vs SZ observations using multi-temperature plasma emission model. We also provide a simple fitting form to estimate the bias given the distribution of fluctuations. Analysing a sample of 16 clusters extracted from hydrodynamical simulations, we find that the median value of bias is within +/-3% within R_500, it decreases to -5% at R_500<r<1.5R_500 and then rises back to ~0% at r >~ 2R_500. The scatter of b_P(r) between individual relaxed clusters is small -- at the level of <0.03 within R_500, but turns significantly larger (0.25) and highly skewed at r>~ 1.5R_500. For any relaxed cluster we find |b_P(r)| < 15% within R_500, across different implementations of input physics in the simulations. [abridged]

Future high resolution, high sensitivity Sunyaev-Zeldovich (SZ) observations of individual clusters will provide an exciting opportunity to answer specific questions about the dynamical state of the intra-cluster medium (ICM). In this paper we develop a new method that clearly shows the connection of the SZ signal with the underlying cluster model. We include relativistic temperature and kinematic corrections in the single-scattering approximation, allowing studies of hot clusters. In our approach, particular moments of the temperature and velocity field along the line-of-sight determine the precise spectral shape and morphology of the SZ signal. We illustrate how to apply our method to different cluster models, highlighting parameter degeneracies and instrumental effects that are important for interpreting future high-resolution SZ data. Our analysis shows that line-of-sight temperature variations can introduce significant biases in the derived SZ temperature and peculiar velocity. We furthermore discuss how the position of the SZ null is affected by the cluster's temperature and velocity structure. Our computations indicate that the SZ signal around the null alone is rather insensitive to different cluster models and that high frequency channels add a large leverage in this respect. We also apply our method to recent high sensitivity SZ data of the Bullet cluster, showing how the results can be linked to line-of-sight variations in the electron temperature. The tools developed here as part of SZpack should be useful for analyzing high-resolution SZ data and computing SZ maps from simulated clusters.

We perform a detailed study of the bispectrum of the Sunyaev-Zel'dovich effect. Using an analytical model for the pressure profiles of the intracluster medium, we demonstrate the SZ bispectrum to be a sensitive probe of the amplitude of the matter power spectrum parameter sigma_8. We find that the bispectrum amplitude scales as B_SZ ~ sigma_8^{11-12}, compared to that of the power spectrum, which scales as A_tSZ ~ sigma_8^{7-9}. We show that the SZ bispectrum is principally sourced by massive clusters at redshifts around z~0.4, which have been well-studied observationally. This is in contrast to the SZ power spectrum, which receives a significant contribution from less-well understood low-mass and high-redshift groups and clusters. Therefore, the amplitude of the bispectrum at l~3000 is less sensitive to astrophysical uncertainties than the SZ power spectrum. We show that current high resolution CMB experiments should be able to detect the SZ bispectrum amplitude with high significance, in part due to the low contamination from extra-galactic foregrounds. A combination of the SZ bispectrum and the power spectrum can sharpen the measurements of thermal and kinetic SZ components and help distinguish cosmological and astrophysical information from high-resolution CMB maps.

We present a new method to identify and characterize the structure of the intracluster medium (ICM) in simulated galaxy clusters. The method uses the median of gas properties, such as density and pressure, which we show to be very robust to the presence of gas inhomogeneities. In particular, we show that the radial profiles of median gas properties are smooth and do not exhibit fluctuations at locations of massive clumps in contrast to mean and mode properties. It is shown that distribution of gas properties in a given radial shell can be well described by a log-normal PDF and a tail. The former corresponds to a nearly hydrostatic bulk component, accounting for ~99% of the volume, while the tail corresponds to high density inhomogeneities. We show that this results in a simple and robust separation of the diffuse and clumpy components of the ICM. The FWHM of the density distribution grows with radius and varies from ~0.15 dex in cluster centre to ~0.5 dex at 2r_500 in relaxed clusters. The small scatter in the width between relaxed clusters suggests that the degree of inhomogeneity is a robust characteristic of the ICM. It broadly agrees with the amplitude of density perturbations in the Coma cluster. We discuss the origin of ICM density variations in spherical shells and show that less than 20% of the width can be attributed to the triaxiality of the cluster gravitational potential. As a link to X-ray observations of real clusters we evaluated the ICM clumping factor with and without high density inhomogeneities. We argue that these two cases represent upper and lower limits on the departure of the observed X-ray emissivity from the median value. We find that the typical value of the clumping factor in the bulk component of relaxed clusters varies from ~1.1-1.2 at r_500 up to ~1.3-1.4 at r_200, in broad agreement with recent observations.

Simulations of cluster formation have demonstrated that condensation of baryons into central galaxies during cluster formation can drive the shape of the gas distribution in galaxy clusters significantly rounder, even at radii as large as half of the virial radius. However, such simulations generally predict stellar fractions within cluster virial radii that are ~2 to 3 times larger than the stellar masses deduced from observations. In this work we compare ellipticity profiles of clusters simulated with and without baryonic cooling to the cluster ellipticity profiles derived from Chandra and ROSAT observations in an effort to constrain the fraction of gas that cools and condenses into the central galaxies within clusters. We find that the observed ellipticity profiles are fairly constant with radius, with an average ellipticity of 0.18 +/- 0.05. The observed ellipticity profiles are in good agreement with the predictions of non-radiative simulations. On the other hand, the ellipticity profiles of the clusters in simulations that include radiative cooling, star formation, and supernova feedback (but no AGN feedback) deviate significantly from the observed ellipticity profiles at all radii. The simulations with cooling overpredict (underpredict) ellipticity in the inner (outer) regions of galaxy clusters. By comparing the simulations with and without cooling, we show that the cooling of gas via cooling flows in the central regions of simulated clusters causes the gas distribution to be more oblate in the central regions, but makes the outer gas distribution more spherical. We find that late-time gas cooling and star formation are responsible for the significantly oblate gas distributions in cluster cores, but the gas shapes outside of cluster cores are set primarily by baryon dissipation at high redshift z > 2.

We present a detailed investigation of the impact of astrophysical processes on the shape and amplitude of the kinetic Sunyaev-Zel'dovich (kSZ) power spectrum from the post-reionization epoch. This is achieved by constructing a new model of the kSZ power spectrum which we calibrate to the results of hydrodynamic simulations. By construction, our method accounts for all relevant density and velocity modes and so is unaffected by the limited box size of our simulations. We find that radiative cooling and star-formation can reduce the amplitude of the kSZ power spectrum by up to 33%, or 1 uK^2 at ell = 3000. This is driven by a decrease in the mean gas density in groups and clusters due to the conversion of gas into stars. Variations in the redshifts at which helium reionization occurs can effect the amplitude by a similar fraction, while current constraints on cosmological parameters (namely sigma_8) translate to a further +-15% uncertainty on the kSZ power spectrum. We demonstrate how the models presented in this work can be constrained -- reducing the astrophysical uncertainty on the kSZ signal -- by measuring the redshift dependence of the signal via kSZ tomography. Finally, we discuss how the results of this work can help constrain the duration of reionization via measurements of the kinetic SZ signal sourced by inhomogeneous (or patchy) reionization.

We present an estimate of the projected two-point correlation function (2PCF) of quasars in the Sloan Digital Sky Survey (SDSS) over the full range of one- and two-halo scales, 0.02-120 Mpc/h. This was achieved by combining data from SDSS DR7 on large scales and Hennawi et al. (2006; with appropriate statistical corrections) on small scales. Our combined clustering sample is the largest spectroscopic quasar clustering sample to date, containing ~48,000 quasars in the redshift range 0.4<z<2.5 with median redshift 1.4. We interpret these precise 2PCF measurements within the halo occupation distribution (HOD) framework and constrain the occupation functions of central and satellite quasars in dark matter halos. In order to explain the small-scale clustering, the HOD modeling requires that a small fraction of z~1.4 quasars, fsat=(7.4+/-1.4) 10^(-4), be satellites in dark matter halos. At z~1.4, the median masses of the host halos of central and satellite quasars are constrained to be Mcen=(4.1+0.3/-0.4) 10^12 Msun/h and Msat=(3.6+0.8/-1.0) 10^14 Msun/h, respectively. To investigate the redshift evolution of the quasar-halo relationship, we also perform HOD modeling of the projected 2PCF measured by Shen et al. (2007) for SDSS quasars with median redshift 3.2. We find tentative evidence for an increase in the mass scale of quasar host halos---the inferred median mass of halos hosting central quasars at z~3.2 is Mcen=(14.1+5.8/-6.9) 10^12 Msun/h. The cutoff profiles of the mean occupation functions of central quasars reveal that quasar luminosity is more tightly correlated with halo mass at higher redshifts. The average quasar duty cycle around the median host halo mass is inferred to be fq=(7.3+0.6/-1.5) 10^(-4) at z~1.4 and fq=(8.6+20.4/-7.2) 10^(-2) at z~3.2. We discuss the implications of our results for quasar evolution and quasar-galaxy co-evolution.

New generation ground and space-based CMB experiments have ushered in discoveries of massive galaxy clusters via the Sunyaev-Zeldovich (SZ) effect, providing a new window for studying cluster astrophysics and cosmology. Many of the newly discovered, SZ-selected clusters contain hot intracluster plasma (kTe > 10 keV) and exhibit disturbed morphology, indicative of frequent mergers with large peculiar velocity (v > 1000 km s^{-1}). It is well-known that for the interpretation of the SZ signal from hot, moving galaxy clusters, relativistic corrections must be taken into account, and in this work, we present a fast and accurate method for computing these effects. Our approach is based on an alternative derivation of the Boltzmann collision term which provides new physical insight into the sources of different kinematic corrections in the scattering problem. This allows us to obtain a clean separation of kinematic and scattering terms which differs from previous works. We briefly mention additional complications connected with kinematic effects that should be considered when interpreting future SZ data for individual clusters. One of the main outcomes of this work is SZpack, a numerical library which allows very fast and precise (<~0.001% at frequencies h nu <~ 20kT_g) computation of the SZ signals up to high electron temperature (kT_e ~ 25 keV) and large peculiar velocity (v/c ~ 0.01). The accuracy is well beyond the current and future precision of SZ observations and practically eliminates uncertainties related to more expensive numerical evaluation of the Boltzmann collision term. Our new approach should therefore be useful for analyzing future high-resolution, multi-frequency SZ observations as well as computing the predicted SZ effect signals from numerical simulations.

We analyze 2.8-yr data of 1-100 GeV photons for clusters of galaxies, collected with the Large Area Telescope onboard the Fermi satellite. By analyzing 49 nearby massive clusters located at high Galactic latitudes, we find no excess gamma-ray emission towards directions of the galaxy clusters. Using flux upper limits, we show that the Fornax cluster provides the most stringent constraints on the dark matter annihilation cross section. Stacking a large sample of nearby clusters does not help improve the limit for most dark matter models. This suggests that a detailed modeling of the Fornax cluster is important for setting robust limits on the dark matter annihilation cross section based on clusters. We therefore perform the detailed mass modeling and predict the expected dark matter annihilation signals from the Fornax cluster, by taking into account effects of dark matter contraction and substructures. By modeling the mass distribution of baryons (stars and gas) around a central bright elliptical galaxy, NGC 1399, and using a modified contraction model motivated by numerical simulations, we show that the dark matter contraction boosts the annihilation signatures by a factor of 4. For dark matter masses around 10 GeV, the upper limit obtained on the annihilation cross section times relative velocity is <\sigma v> <~ (2-3)x10^{-25} cm^3 s^{-1}, which is within a factor of 10 from the value required to explain the dark matter relic density. This effect is more robust than the annihilation boost due to substructure, and it is more important unless the mass of the smallest subhalos is much smaller than that of the Sun.

In this work, we examine the effects of mergers on the hydrostatic mass estimate of galaxy clusters using high-resolution Eulerian cosmological simulations. We utilize merger trees to isolate the last merger for each cluster in our sample and follow the time evolution of the hydrostatic mass bias as the systems relax. We find that during a merger, a shock propagates outward from the parent cluster, resulting in an overestimate in the hydrostatic mass bias. After the merger, as a cluster relaxes, the bias in hydrostatic mass estimate decreases but remains at a level of -5-10% with 15-20% scatter within r500. We also investigate the post-merger evolution of the pressure support from bulk motions, a dominant cause of this residual mass bias. At r500, the contribution from random motions peaks at 30% of the total pressure during the merger and quickly decays to \sim 10-15% as a cluster relaxes. Additionally, we use a measure of the random motion pressure to correct the hydrostatic mass estimate. We discover that 4 Gyr after mergers, the direct effects of the merger event on the hydrostatic mass bias have become negligible. Thereafter, the mass bias is primarily due to residual bulk motions in the gas which are not accounted for in the hydrostatic equilibrium equation. We present a hydrostatic mass bias correction method that can recover the unbiased cluster mass for relaxed clusters with 9% scatter at r500 and 11% scatter in the outskirts, within r200.

Using a fully cosmological hydrodynamic simulation that self-consistently incorporates the growth and feedback of supermassive black holes and the physics of galaxy formation, we examine the effects of environmental factors (e.g., local gas density, black hole feedback) on the halo occupation distribution of low luminosity active galactic nuclei (AGN). We decompose the mean occupation function into central and satellite contribution and compute the conditional luminosity functions (CLF). The CLF of the central AGN follows a log-normal distribution with the mean increasing and scatter decreasing with increasing redshifts. We analyze the light curves of individual AGN and show that the peak luminosity of the AGN has a tighter correlation with halo mass compared to instantaneous luminosity. We also compute the CLF of satellite AGN at a given central AGN luminosity. We do not see any significant correlation between the number of satellites with the luminosity of the central AGN at a fixed halo mass. We also show that for a sample of AGN with luminosity above 10^42 ergs/s the mean occupation function can be modeled as a softened step function for central AGN and a power law for the satellite population. The radial distribution of AGN inside halos follows a power law at all redshifts with a mean index of -2.33 +/- 0.08. Incorporating the environmental dependence of supermassive black hole accretion and feedback, our formalism provides a theoretical tool for interpreting current and future measurements of AGN clustering.

The condensation of gas and stars in the inner regions of dark matter halos leads to a more concentrated dark matter distribution. While this effect is based on simple gravitational physics, the question of its validity in hierarchical galaxy formation has led to an active debate in the literature. We use a collection of several state-of-the-art cosmological hydrodynamic simulations to study the halo contraction effect in systems ranging from dwarf galaxies to clusters of galaxies, at high and low redshift. The simulations are run by different groups with different codes and include hierarchical merging, gas cooling, star formation, and stellar feedback. We show that in all our cases the inner dark matter density increases relative to the matching simulation without baryon dissipation, at least by a factor of several. The strength of the contraction effect varies from system to system and cannot be reduced to a simple prescription. We present a revised analytical model that describes the contracted mass profile to an rms accuracy of about 10%. The model can be used to effectively bracket the response of the dark matter halo to baryon dissipation. The halo contraction effect is real and must be included in modeling of the mass distribution of galaxies and galaxy clusters.

We present analysis of the three-dimensional shape of intracluster gas in clusters formed in cosmological simulations of the Lambda-CDM cosmology and compare it to the shape of dark matter distribution and the shape of the overall isopotential surfaces. We find that in simulations with radiative cooling, star formation and stellar feedback (CSF), intracluster gas outside the cluster core is more spherical compared to non-radiative (NR) simulations, while in the core the gas in the CSF runs is more triaxial and has a distinctly oblate shape. The latter reflects the ongoing cooling of gas, which settles into a thick oblate ellipsoid as it loses thermal energy. The shape of the gas in the inner regions of clusters can therefore be a useful diagnostic of gas cooling. We find that gas traces the shape of the underlying potential rather well outside the core, as expected in hydrostatic equilibrium. At smaller radii, however, the gas and potential shapes differ significantly. In the CSF runs, the difference reflects the fact that gas is partly rotationally supported. Interestingly, we find that in NR simulations the difference between gas and potential shape at small radii is due to random gas motions, which make the gas distribution more spherical than the equipotential surfaces. Finally, we use mock Chandra X-ray maps to show that the differences in shapes observed in three-dimensional distribution of gas are discernible in the ellipticity of X-ray isophotes. Contrasting the ellipticities measured in simulated clusters against observations can therefore constrain the amount of cooling of the intracluster medium and the presence of random gas motions in cluster cores.

Recent Suzaku X-ray observations revealed that the observed entropy profile of the intracluster medium (ICM) deviates significantly from the prediction of hydrodynamical simulations of galaxy clusters. In this work, we show that gas clumping introduces significant biases in X-ray measurements of the ICM profiles in the outskirts of galaxy clusters. Using hydrodynamical simulations of galaxy cluster formation in a concordance Lambda-CDM model, we demonstrate that gas clumping leads to an overestimate of the observed gas density and causes flattening of the entropy profile. Our results suggest that gas clumping must be taken into account when interpreting X-ray measurements of cluster outskirts.

We investigate the halo occupation distribution (HOD) of black holes within a hydrodynamic cosmological simulation that directly follows black hole growth. Similar to the HOD of galaxies/subhalos, we find that the black hole occupation number can be described by the form N_BH proportional to 1+ (M_Host)^alpha where alpha evolves mildly with redshift indicating that a given mass halo (M_Host) at low redshift tends to host fewer BHs than at high redshift (as expected as a result of galaxy and BH mergers). We further divide the occupation number into contributions from black holes residing in central and satellite galaxies within a halo. The distribution of M_BH within halos tends to consist of a single massive BH (distributed about a peak mass strongly correlated with M_Host), and a collection of relatively low-mass secondary BHs, with weaker correlation with M_Host. We also examine the spatial distribution of BHs within their host halos, and find they typically follow a power-law radial distribution (i.e. much more centrally concentrated than the subhalo distribution). Finally, we characterize the host mass for which BH growth is feedback dominated (e.g. star formation quenched). We show that halos with M_Host > 3 * 10^12 M_sun have primary BHs that are feedback dominated by z~3 with lower mass halos becoming increasingly more affected at lower redshift.

The unique properties of dark matter are revealed during collisions between clusters of galaxies, like the bullet cluster (1E 0657-56) and baby bullet (MACSJ0025-12). These systems provide evidence for an additional, invisible mass in the separation between the distribution of their total mass, measured via gravitational lensing, and their ordinary 'baryonic' matter, measured via its X-ray emission. Unfortunately, the information available from these systems is limited by their rarity. Constraints on the properties of dark matter, such as its interaction cross-section, are therefore restricted by uncertainties in the individual systems' impact velocity, impact parameter and orientation with respect to the line of sight. Here we develop a complementary, statistical measurement in which every piece of substructure falling into every massive cluster is treated as a bullet. We define 'bulleticity' as the mean separation between dark matter and ordinary matter, and we measure the signal in hydrodynamical simulations. The phase space of substructure orbits also exhibits symmetries that provide an equivalent control test. Any detection of bulleticity in real data would indicate a difference in the interaction cross-sections of baryonic and dark matter that may rule out hypotheses of non-particulate dark matter that are otherwise able to model individual systems. A subsequent measurement of bulleticity could constrain the dark matter cross-section. Even with conservative estimates, the existing HST archive should yield an independent constraint tighter than that from the bullet cluster. This technique is then trivially extendable to, and benefits enormously from larger, future surveys.

We present cosmological simulations of galaxy clusters, with focus on the cluster outskirts. We show that large-scale cosmic accretion and mergers produce significant internal gas motions and inhomogeneous gas distribution ("clumpiness") in the intracluster medium (ICM) and introduce biases in measurements of the ICM profiles and the cluster mass. We also show that non-thermal pressure provided by the gas motions is one of the dominant sources of theoretical uncertainties in cosmic microwave background secondary anisotropies. We briefly discuss implications for cluster cosmology and future prospects for understanding the physics of cluster outskirts using computer simulations and multi-wavelength cluster surveys.

We use an analytic model to investigate the theoretical uncertainty on the thermal Sunyaev-Zel'dovich (SZ) power spectrum due to astrophysical uncertainties in the thermal structure of the intracluster medium. Our model accounts for star formation and energy feedback (from supernovae and active galactic nuclei) as well as radially dependent non-thermal pressure support due to random gas motions, the latter calibrated by recent hydrodynamical simulations. We compare the model against X-ray observations of low redshift clusters, finding excellent agreement with observed pressure profiles. Varying the levels of feedback and non-thermal pressure support can significantly change both the amplitude and shape of the thermal SZ power spectrum. Increasing the feedback suppresses power at small angular scales, shifting the peak of the power spectrum to lower ell. On the other hand, increasing the non-thermal pressure support has the opposite effect, significantly reducing power at large angular scales. In general, including non-thermal pressure at the level measured in simulations has a large effect on the power spectrum, reducing the amplitude by 50% at angular scales of a few arcminutes compared to a model without a non-thermal component. Our results demonstrate that measurements of the shape of the power spectrum can reveal useful information on important physical processes in groups and clusters, especially at high-redshift where there exists little observational data. Comparing with the recent South Pole Telescope measurements of the small-scale cosmic microwave background power spectrum, we find our model reduces the tension between the values of sigma_8 measured from the SZ power spectrum and from cluster abundances.

We investigate effects of baryon dissipation on the dark matter virial scaling relation between total mass and velocity dispersion and the velocity bias of galaxies in groups and clusters using self-consistent cosmological simulations. We show that the baryon dissipation increases the velocity dispersion of dark matter within the virial radius by 5% - 10%. The effect is mainly driven by the change in density and gravitational potential in inner regions of cluster, and it is larger in lower mass systems where gas cooling and star formation are more efficient. We also show that the galaxy velocity bias depends on how galaxies are selected. Galaxies selected based on their stellar mass exhibit no velocity bias, while galaxies selected based on their total mass show positive bias of ~10%, consistent with previous results based on collisionless dark matter- only simulations. We further find that observational estimates of galaxy velocity dispersion are unbiased with respect to the velocity dispersion of dark matter, provided galaxies are selected using their stellar masses and and their velocity dispersions are computed with more than twenty most massive galaxies. Velocity dispersions estimated with fewer galaxies, on the other hand, can lead to significant underestimate of dynamical masses. Results presented in this paper should be useful in interpretating high-redshift groups and clusters as well as cosmological constraints derived from upcoming optical cluster surveys.

Recent observations with Galaxy Evolution Explorer (GALEX) show strong unexpected UV excess in the spectrum of brightest cluster galaxies (BCGs). It is believed that the excess UV signal is produced by old and evolved core-He burning stars, and the UV flux strength could be greatly enhanced if the progenitor stars have high value of He abundance. In this work, we propose that sedimentation process can greatly enhance the He abundance in BCGs. Our model predicts that the UV flux strength is stronger in more massive, low-redshift, and dynamically relaxed BCGs. These predictions are testable with the current generation of GALEX+SDSS observations.

We present analysis of bulk and random gas motions in the intracluster medium using high-resolution Eulerian cosmological simulations of sixteen simulated clusters, including both very relaxed and unrelaxed systems and spanning a virial mass range of 5*10^13 - 2*10^15 Msun/h. We investigate effects of the residual subsonic gas motions on the hydrostatic estimates of mass profiles and concentrations of galaxy clusters. In agreement with previous studies we find that the gas motions contribute up to ~ 5%-15% of the total pressure support in relaxed clusters with contribution increasing with cluster-centric radius. The fractional pressure support is higher in unrelaxed systems. This contribution would not be accounted for in hydrostatic estimates of the total mass profile and would lead to systematic underestimate of mass. We demonstrate that total mass can be recovered accurately if pressure due to gas motions measured in simulations is explicitly taken into account in the equation of hydrostatic equilibrium. Given that the underestimate of mass is increasing at larger radii, where gas is less relaxed and contribution of gas motions to pressure is larger, the total density profile derived from hydrostatic analysis is more concentrated than the true profile. This may at least partially explain some high values of concentrations of clusters estimated from hydrostatic analysis of X-ray data.

We present high-resolution cosmological hydrodynamic simulations of three galaxy clusters employing a two-temperature model for the intracluster medium. We show that electron temperatures in cluster outskirts are significantly lower than the mean gas temperature, because Coulomb collisions are insufficient to keep electrons and ions in thermal equilibrium. This deviation is larger in more massive and less relaxed systems, ranging from 5% in relaxed clusters to 30% for clusters undergoing major mergers. The presence of non-equilibrium electrons leads to significant suppression of the SZE signal at large cluster-centric radius. The suppression of the electron pressure also leads to an underestimate of the hydrostatic mass. Merger-driven, internal shocks may also generate significant populations of non-equilibrium electrons in the cluster core, leading to a 5% bias on the integrated SZ mass proxy during cluster mergers.

A moderate investment of observing time with the International X-ray Observatory to study high-redshift galaxy clusters detected in future large-scale surveys, will provide cosmological measurements of fundamental importance. IXO observations, combined with lensing follow-up, will measure the perturbation growth factor from z=0-2 with an accuracy comparable to, or possibly better than, that expected from observations of cosmic shear with JDEM, and redshift-space distortions with EUCLID. The growth of structure data derived from clusters will significantly improve our knowledge of the dark energy equation of state and will aid in constraining non-GR models for cosmic acceleration. IXO observations of the largest, dynamically relaxed clusters will provide a powerful, independent measurement of the cosmological expansion history using the apparent f_gas(z) trend. Systematic and statistical errors from this technique are competitive with SNIa and BAO studies, making the test extremely useful for improving the accuracy and reliability of the geometric cosmological measurements planned for LSST and JDEM. Only by employing a range of powerful, independent approaches, including those discussed here, can robust answers to puzzles as profound as the origin of cosmic acceleration be expected.

We are learning much about how structure forms, in particular how clusters as nodes in the cosmic web evolve and accrete matter, and about the physical processes within these objects. In the next decade, the study of clusters will enable us to tackle important questions regarding the nature of Dark Matter and Dark Energy, how clusters co-evolve with super-massive black holes at their centers, and to advance our knowledge about fundamental plasma astrophysics. This science white paper outlines the key questions and research opportunities in cluster astrophysics that are emerging in the coming decade and beyond, and serves as an overview to other cluster related white papers.

Groups and clusters of galaxies occupy a special position in the hierarchy of large-scale cosmic structures because they are the largest and the most massive (from ~10^13 Msun to over 10^15 Msun) objects in the universe that have had time to undergo gravitational collapse. The large masses of clusters imply that their contents have been accreted from regions of ~8-40 comoving Mpc in size and should thus be representative of the mean matter content of the universe. During the next decade sensitive multi-wavelength observations should be able to map the radial distributions of all main mass components (stars, cold, warm, and hot gas and total mass) at z<~ 1 out to the virial radius. At the same time, comparative studies of real and simulated cluster samples sould allow us to use clusters as veritable astrophysical laboratories for studying galaxy formation, as well as testing our theoretical models of structure formation and underlying assumptions about fundamental physics governing the universe.

We investigate the utility of a new, self-similar pressure profile for fitting Sunyaev-Zel'dovich (SZ) effect observations of galaxy clusters. Current SZ imaging instruments - such as the Sunyaev-Zel'dovich Array (SZA) - are capable of probing clusters over a large range in physical scale. A model is therefore required that can accurately describe a cluster's pressure profile over a broad range of radii, from the core of the cluster out to a significant fraction of the virial radius. In the analysis presented here, we fit a radial pressure profile derived from simulations and detailed X-ray analysis of relaxed clusters to SZA observations of three clusters with exceptionally high quality X-ray data: A1835, A1914, and CL J1226.9+3332. From the joint analysis of the SZ and X-ray data, we derive physical properties such as gas mass, total mass, gas fraction and the intrinsic, integrated Compton y-parameter. We find that parameters derived from the joint fit to the SZ and X-ray data agree well with a detailed, independent X-ray-only analysis of the same clusters. In particular, we find that, when combined with X-ray imaging data, this new pressure profile yields an independent electron radial temperature profile that is in good agreement with spectroscopic X-ray measurements.

Chandra observations of large samples of galaxy clusters detected in X-rays by ROSAT provide a new, robust determination of the cluster mass functions at low and high redshifts. Statistical and systematic errors are now sufficiently small, and the redshift leverage sufficiently large for the mass function evolution to be used as a useful growth of structure based dark energy probe. In this paper, we present cosmological parameter constraints obtained from Chandra observations of 36 clusters with <z>=0.55 derived from 400deg^2 ROSAT serendipitous survey and 49 brightest z=~0.05 clusters detected in the All-Sky Survey. Evolution of the mass function between these redshifts requires Omega_Lambda>0 with a ~5sigma significance, and constrains the dark energy equation of state parameter to w0=-1.14+-0.21, assuming constant w and flat universe. Cluster information also significantly improves constraints when combined with other methods. Fitting our cluster data jointly with the latest supernovae, WMAP, and baryonic acoustic oscillations measurements, we obtain w0=-0.991+-0.045 (stat) +-0.039 (sys), a factor of 1.5 reduction in statistical uncertainties, and nearly a factor of 2 improvement in systematics compared to constraints that can be obtained without clusters. The joint analysis of these four datasets puts a conservative upper limit on the masses of light neutrinos, Sum m_nu<0.33 eV at 95% CL. We also present updated measurements of Omega_M*h and sigma_8 from the low-redshift cluster mass function.

We discuss the measurements of the galaxy cluster mass functions at z=~0.05 and z=~0.5 using high-quality Chandra observations of samples derived from the ROSAT PSPC All-Sky and 400deg^2 surveys. We provide a full reference for the data analysis procedures, present updated calibration of relations between the total cluster mass and its X-ray indicators (T_X, Mgas, and Y_X) based on a subsample of low-z relaxed clusters, and present a first measurement of the evolving L_X-Mtot relation (with Mtot estimated from Y_X) obtained from a well-defined statistically complete cluster sample and with appropriate corrections for the Malmquist bias applied. Finally, we present the derived cluster mass functions, estimate the systematic uncertainties in this measurement, and discuss the calculation of the likelihood function. We confidently measure the evolution in the cluster comoving number density at a fixed mass threshold, e.g., by a factor of 5.0 +- 1.2 at M_500=2.5e14 h^-1 Msun between z=0 and 0.5. This evolution reflects the growth of density perturbations and can be used for the cosmological constraints complementing those from the distance-redshift relation.

The uniformity of the helium-to-hydrogen abundance ratio in X-ray emitting intracluster medium (ICM) is one of the commonly adopted assumptions in X-ray analyses of galaxy clusters and cosmological constraints derived from these measurements. In this work, we investigate the effect of He sedimentation on X-ray measurements of galaxy clusters in order to assess this assumption and associated systematic uncertainties. By solving a set of flow equations for a H-He plasma, we show that the helium-to-hydrogen mass ratio is significantly enhanced in the inner regions of clusters. The effect of He sedimentation, if not accounted for, introduces systematic biases in observable properties of clusters derived using X-ray observations. We show that these biases also introduce an apparent evolution in the observed gas mass fractions of X-ray luminous, dynamically relaxed clusters and hence biases in observational constraints on the dark energy equation of state parameter, w, derived from the cluster distance-redshift relation. The Hubble parameter derived from the combination of X-ray and Sunyaev-Zel'dovich effect (SZE) measurements is affected by the He sedimentation process as well. Future measurements aiming to constrain w or H_0 to better than 10% may need to take into account the effect of He sedimentation. We propose that the evolution of gas mass fraction in the inner regions of clusters should provide unique observational diagnostics of the He sedimentation process.

The merging galaxy cluster 1E 0657-56, known as the "bullet cluster," is one of the hottest clusters known. The X-ray emitting plasma exhibits bow-shock-like temperature and density jumps. The segregation of this plasma from the peaks of the mass distribution determined by gravitational lensing has been interpreted as a direct proof of collisionless dark matter. If the high shock speed inferred from the shock jump conditions equals the relative speed of the merging CDM halos, however, this merger is predicted to be such a rare event in a LCDM universe that observing it presents a possible conflict with the LCDM model. We examined this question using high resolution, 2D simulations of gas dynamics in cluster collisions to analyze the relative motion of the clusters, the bow shock, and the contact discontinuity, and relate these to the X-ray data for the bullet cluster. We find that the velocity of the fluid shock need not equal the relative velocity of the CDM components. An illustrative simulation finds that the present relative velocity of the CDM halos is 16% lower than that of the shock. While this conclusion is sensitive to the detailed initial mass and gas density profiles of the colliding clusters, such a decrease of the inferred halo relative velocity would significantly increase the likelihood of finding 1E 0657-56 in a LCDM universe. (Conference proceedings based on a poster at Bash Symposium 2007)

Cosmic rays produced in cluster accretion and merger shocks provide pressure to the intracluster medium (ICM) and affect the mass estimates of galaxy clusters. Although direct evidence for cosmic-ray ions in the ICM is still lacking, they produce gamma-ray emission through the decay of neutral pions produced in their collisions with ICM nucleons. We investigate the capability of the Gamma-ray Large Area Space Telescope (GLAST) and imaging atmospheric Cerenkov telescopes (IACTs) for constraining the cosmic-ray pressure contribution to the ICM. We show that GLAST can be used to place stringent upper limits, a few per cent for individual nearby rich clusters, on the ratio of pressures of the cosmic rays and thermal gas. We further show that it is possible to place tight (<~10%) constraints for distant (z <~ 0.25) clusters in the case of hard spectrum, by stacking signals from samples of known clusters. The GLAST limits could be made more precise with the constraint on the cosmic-ray spectrum potentially provided by IACTs. Future gamma-ray observations of clusters can constrain the evolution of cosmic-ray energy density, which would have important implications for cosmological tests with upcoming X-ray and Sunyaev-Zel'dovich effect cluster surveys.

We present Sunyaev-Zel'dovich Effect (SZE) scaling relations for 38 massive galaxy clusters at redshifts 0.14<z0.89, observed with both the Cchandra X-ray Observatory and the centimeter-wave SZE imaging system at the BIMA and OVRO interferometric arrays. An isothermal beta-model with central 100 kpc excluded from the X-ray data is used to model the intracluster medium and to measure global cluster properties. For each cluster, we measure the X-ray spectroscopic temperature, SZE gas mass, total mass and integrated Compton-y parameters within r_2500. Our measurements are in agreement with the expectations based on a simple self-similar model of cluster formation and evolution. We compare the cluster properties derived from our SZE observations with and without Chandra spatial and spectral information and find them to be in good agreement. We compare our results with cosmological numerical simulations, and find that simulations that include radiative cooling, star formation and feedback match well both the slope and normalization of our SZE scaling relations.

We present detailed comparisons of the intracluster medium (ICM) in cosmological Eulerian cluster simulations with deep Chandra observations of nearby relaxed clusters. To assess the impact of galaxy formation, we compare two sets of simulations, one performed in the non-radiative regime and another with radiative cooling and several physical processes critical to various aspects of galaxy formation: star formation, metal enrichment and stellar feedback. We show that the observed ICM properties outside cluster cores are well-reproduced in the simulations that include cooling and star formation, while the non-radiative simulations predict an overall shape of the ICM profiles inconsistent with observations. In particular, we find that the ICM entropy in our runs with cooling is enhanced to the observed levels at radii as large as half of the virial radius. We also find that outside cluster cores entropy scaling with the mean ICM temperature in both simulations and Chandra observations is consistent with being self-similar within current error bars. We find that the pressure profiles of simulated clusters are also close to self-similar and exhibit little cluster-to-cluster scatter. The X-ray observable-total mass relations for our simulated sample agree with the Chandra measurements to \~10%-20% in normalization. We show that this systematic difference could be caused by the subsonic gas motions, unaccounted for in X-ray hydrostatic mass estimates. The much improved agreement of simulations and observations in the ICM profiles and scaling relations is encouraging and the existence of tight relations of X-ray observables, such as Yx, and total cluster mass and the simple redshift evolution of these relations hold promise for the use of clusters as cosmological probes.

X-ray observations of galaxy cluster merger shocks can be used to constrain nonthermal processes in the intracluster medium (ICM). The presence of nonthermal pressure components in the ICM, as well as the shock acceleration of particles and their escape, all affect shock jump conditions in distinct ways. Therefore, these processes can be constrained using X-ray surface brightness and temperature maps of merger shock fronts. Here we use these observations to place constraints on particle acceleration efficiency in intermediate Mach number (M ~ 2-3) shocks and explore the potential to constrain the contribution of nonthermal components (e.g., cosmic rays, magnetic field, and turbulence) to ICM pressure in cluster outskirts. We model the hydrodynamic jump conditions in merger shocks discovered in the galaxy clusters A520 (M ~ 2) and 1E 0657-56 (M ~ 3) using a multifluid model comprised of a thermal plasma, a nonthermal plasma, and a magnetic field. Based on the published X-ray spectroscopic data alone, we find that the fractional contribution of cosmic rays accelerated in these shocks is lower than about 10% of the shock downstream pressure. Current observations do not constrain the fractional contribution of nonthermal components to the pressure of the undisturbed shock upstream. Future X-ray observations, however, have the potential to either detect particle acceleration in these shocks through its effect on the shock dynamics, or to place a lower limit on the nonthermal pressure contributions in the undisturbed ICM. We briefly discuss implications for models of particle acceleration in collisionless shocks and the estimates of galaxy cluster masses derived from X-ray and Sunyaev-Zel'dovich effect observations.

The Sunyaev-Zel'dovich (SZ) effect is a direct probe of thermal energy content of the Universe, induced in the cosmic microwave background (CMB) sky through scattering of CMB photons off hot electrons in the intracluster medium (ICM). We report a 9-sigma detection of the SZ signal in the CMB maps of Wilkinson Microwave Anisotropy Probe (WMAP) 3yr data, through study of a sample of 193 massive galaxy clusters with observed X-ray temperatures greater than 3 keV. For the first time, we make a model-independent measurement of the pressure profile in the outskirts of the ICM, and show that it closely follows the profiles obtained by X-ray observations and numerical simulations. We find that our measurements of the SZ effect would account for only half of the thermal energy of the cluster, if all the cluster baryons were in the hot ICM phase. Our measurements indicate that a significant fraction (35 +/- 8 %) of baryonic mass is missing from the hot ICM, and thus must have cooled to form galaxies, intracluster stars, or an unknown cold phase of the ICM. There does not seem to be enough mass in the form of stars or cold gas in the cluster galaxies or intracluster space, signaling the need for a yet-unknown baryonic component (at 3-sigma level), or otherwise new astrophysical processes in the ICM.

Shock waves driven in the intergalactic medium during the merging of galaxy clusters have been observed in X-ray imaging and spectroscopy. Fluid motions inferred from the shock strength and morphology can be compared to the cold dark matter (CDM) distribution inferred from gravitational lensing. A detailed reconstruction of the CDM kinematics, however, must take into account the nontrivial response of the fluid intracluster medium to the collisionless CDM motions. We have carried out two-dimensional simulations of gas dynamics in cluster collisions. We analyze the relative motion of the clusters, the bow shock wave, and the contact discontinuity and relate these to X-ray data. We focus on the "bullet cluster," 1E 0657-56, a near head-on collision of unequal-mass clusters, for which the gas density and temperature jumps across the prominent bow shock imply a high shock velocity 4,700 km/s. The velocity of the fluid shock has been widely interpreted as the relative velocity of the CDM components. This need not be the case, however. An illustrative simulation finds that the present relative velocity of the CDM halos is 16% lower than that of the shock. While this conclusion is sensitive to the detailed initial mass and gas density profile of the colliding clusters, such a decrease of the inferred halo relative velocity would increase the likelihood of finding 1E 0657-56 in a LambdaCDM universe.

X-ray observations of galaxy clusters potentially provide powerful cosmological probes if systematics due to our incomplete knowledge of the intracluster medium (ICM) physics are understood and controlled. In this paper, we study the effects of galaxy formation on the properties of the ICM and X-ray observable-mass relations using high-resolution self-consistent cosmological simulations of galaxy clusters and comparing their results with recent Chandra X-ray observations. We show that despite complexities of their formation and uncertainties in their modeling, clusters of galaxies both in observations and numerical simulations are remarkably regular outside of their cores, which holds great promise for their use as cosmological probes.

Using ray tracing for simple analytic profiles, we demonstrate that the lensing cross section for producing giant arcs has distinct contributions due to arcs formed through image distortion only, and arcs form from the merging of two or three images. We investigate the dependence of each of these contributions on halo ellipticity and on the slope of the density profile, and demonstrate that at fixed Einstein radius, the lensing cross section increases as the halo profile becomes steeper. We then compare simulations with and without baryonic cooling of the same cluster for a sample of six clusters, and demonstrate that cooling can increase the overall abundance of giant arcs by factors of a few. The net boost to the lensing probability for individual clusters is mass dependent, and can lower the effective low mass limit of lensing clusters. This last effect can potentially increase the number of lensing clusters by an extra 50%. While the magnitude of these effects may be overestimated due to the well known overcooling problem in simulations, it is evident that baryonic cooling has a non-negligible impact on the expected abundance of giant arcs, and hence cosmological constraints from giant arc abundances may be subject to large systematic errors.

X-ray observations of galaxy clusters potentially provide powerful cosmological probes if systematics due to our incomplete knowledge of the intracluster medium (ICM) physics are understood and controlled. In this paper, we present mock Chandra analyses of cosmological cluster simulations and assess X-ray measurements of galaxy cluster properties using a model and procedure essentially identical to that used in real data analysis. We show that reconstruction of three-dimensional ICM density and temperature profiles is excellent for relaxed clusters, but still reasonably accurate for unrelaxed systems. The total ICM mass is measured quite accurately (<6%) in all clusters, while the hydrostatic estimate of the gravitationally bound mass is biased low by about 5%-20% through the virial region, primarily due to additional pressure support provided by subsonic bulk motions in the ICM, ubiquitous in our simulations even in relaxed systems. Gas fraction determinations are therefore biased high; the bias increases toward cluster outskirts and depends sensitively on its dynamical state, but we do not observe significant trends of the bias with cluster mass or redshift. We also find that different average ICM temperatures, such as the X-ray spectroscopic Tspec and gas-mass-weighted Tmg, are related to each other by a constant factor with a relatively small object-to-object scatter and no systematic trend with mass, redshift or the dynamical state of clusters. We briefly discuss direct applications of our results for different cluster-based cosmological tests.

This presentation is a Moriond version of our recent paper (Kravtsov, Vikhlinin & Nagai astro-ph/0603205) where we discussed X-ray proxies for the total cluster mass, including the spectral temperature (Tx), gas mass measured within r500 (Mg), and the new proxy, Yx, which is a simple product of Tx and Mg. We use mock Chandra images constructed for a sample of clusters simulated with high resolution in the concordance LambdaCDM cosmology. The simulated clusters exhibit tight correlations between the considered observables and total mass. The normalizations of the M500-Tx, Mg-Tx, and M500-Yx relations agree to better than =~ 10-15% with the current observational measurements of these relations. Our results show that Yx is the best mass proxy with a remarkably low scatter of only =~ 5-7% in M500 for a fixed Yx, at both low and high redshifts and regardless of whether clusters are relaxed or not. In addition, we show that redshift evolution of the Yx-M500 relation is close to the self-similar prediction, which makes Yx a very attractive mass indicator for measurements of the cluster mass function from X-ray selected samples.

We study the effects of galaxy formation on the Sunyaev-Zel'dovich effect (SZE) observable-mass relations using high-resolution cosmological simulations. The simulations of eleven individual clusters spanning a decade in mass are performed with the shock-capturing Eulerian adaptive mesh refinement N-body+gasdynamics ART code. To assess the impact of galaxy formation, we compare two sets of simulations performed in an adiabatic regime and those with several physical processes critical to various aspects of galaxy formation: radiative cooling, star formation, stellar feedback and metal enrichment. We show that a SZE signal integrated to a sufficiently large fraction of cluster volume correlates strongly with its enclosed mass, independent of details of gas physics and dynamical state of the cluster. The slope and redshift evolution of the SZE flux-mass relation are also insensitive to processes of galaxy formation and are well characterized by a simple self-similar cluster model. Its normalization, on the other hand, is significantly affected by gas cooling and associated star formation. Our simulations show that inclusion of these processes suppresses the normalization by ~30-40%. The effect is due to a decrease in gas mass fraction, which is offset slightly by an increase in gas mass-weighted temperature. Gas cooling and star formation also cause an increase in total mass and modify the normalization by a few percent. Finally, the comparison with recent observations of the SZE scaling relations highlights the importance of galaxy formation in theoretical modelling of clusters and shows that current generation of simulations produce clusters with gross properties quite similar to their observed counterparts.

In this paper we examine how well galaxies and intra-cluster gas trace the gravitational potential of clusters. Utilizing mass profiles derived from gravitational lensing and X-ray observations, coupled with measured galaxy velocities, we solve for the velocity anisotropy parameter using the anisotropic Jeans equation. This is done for five clusters, three at low redshift: A2199, A496 and A576 and two at high redshifts: A2390 and MS1358. With X-ray temperature profiles obtained from Chandra and ASCA/ROSAT data, we estimate the ratio of energy in the galaxies compared to the X-ray gas. We find that none of these clusters is strictly in hydro-static equilibrium. We compare the properties of our sample with clusters that form in high-resolution cosmological N-body simulations that include baryonic physics. Simulations and data show considerable scatter both these profiles. We demonstrate the future feasibility and potential for directly comparing the orbital structure of clusters inferred from multi-wavelength observations with high resolution simulated clusters.

We present gas mass fractions of 38 massive galaxy clusters spanning redshifts from 0.14 to 0.89, derived from Chandra X-ray data and OVRO/BIMA interferometric Sunyaev-Zel'dovich Effect measurements. We use three models for the gas distribution: (1) an isothermal beta-model fit jointly to the X-ray data at radii beyond 100 kpc and to all of the SZE data,(2) a non-isothermal double beta-model fit jointly to all of the X-ray and SZE data, and (3) an isothermal beta-model fit only to the SZE spatial data. We show that the simple isothermal model well characterizes the intracluster medium (ICM) outside of the cluster core in clusters with a wide range of morphological properties. The X-ray and SZE determinations of mean gas mass fractions for the 100 kpc-cut isothermal beta-model are fgas(X-ray)=0.110 +0.003-0.003 +0.006-0.018 and fgas(SZE)=0.116 +0.005-0.005 +0.009-0.026, where uncertainties are statistical followed by systematic at 68% confidence. For the non-isothermal double beta-model, fgas(X-ray)=0.119 +0.003-0.003 +0.007-0.014 and fgas(SZE)=0.121 +0.005-0.005 +0.009-0.016. For the SZE-only model, fgas(SZE)=0.120 +0.009-0.009 +0.009-0.027. Our results indicate that the ratio of the gas mass fraction within r2500 to the cosmic baryon fraction is 0.68 +0.10-0.16 where the range includes statistical and systematic uncertainties. By assuming that cluster gas mass fractions are independent of redshift, we find that the results are in agreement with standard LambdaCDM cosmology and are inconsistent with a flat matter dominated universe.

We present comparison of X-ray proxies for the total cluster mass, including the spectral temperature (Tx), gas mass measured within r500 (Mg), and the new proxy, Yx, which is a simple product of Tx and Mg and is related to the total thermal energy of the ICM. We use mock Chandra images constructed for a sample of clusters simulated with the eulerian N-body+gasdynamics adaptive mesh refinement ART code in the concordance LCDM cosmology. The simulations achieve high spatial and mass resolution and include radiative cooling, star formation, and other processes accompanying galaxy formation. Our analysis shows that simulated clusters exhibit a high degree of regularity and tight correlations between the considered observables and total mass. The normalizations of the M-Tx, Mg-Tx, and M-Yx relations agree to better than 10-15% with the current observational measurements of these relations. Our results show that Yx is the best mass proxy with a remarkably low scatter of only ~5-7% in M500 for a fixed Yx, at both low and high redshifts and regardless of whether clusters are relaxed or not. In addition, we show that redshift evolution of the Yx-M500 relation is close to the self-similar prediction, which makes Yx a very attractive mass indicator for measurements of the cluster mass function from X-ray selected samples.

Observational evidence indicates a mismatch between the shapes of collisionless dark matter (DM) halos and those of observed systems. Using hydrodynamical cosmological simulations we investigate the effect of baryonic dissipation on halo shapes. We show that dissipational simulations produce significantly rounder halos than those formed in equivalent dissipationless simulations. Gas cooling causes an average increase in halo principal axis ratios of ~ 0.2-0.4 in the inner regions and a systematic shift that persists out to the virial radius, alleviating any tension between theory and observations. Although the magnitude of the effect may be overestimated due to overcooling, cluster formation simulations designed to reproduce the observed fraction of cold baryons still produce substantially rounder halos. Subhalos also exhibit a trend of increased axis ratios in dissipational simulations. Moreover, we demonstrate that subhalos are generally rounder than corresponding field halos even in dissipationless simulations. Lastly, we analyze a series of binary, equal-mass merger simulations of disk galaxies. Collisionless mergers reveal a strong correlation between DM halo shape and stellar remnant morphology. In dissipational mergers, the combination of strong gas inflows and star formation leads to an increase of the DM axis ratios in the remnant. All of these results highlight the vital role of baryonic processes in comparing theory with observations and warn against over-interpreting discrepancies with collisionless simulations on small scales.

We study the effects of dissipation on the baryon fractions in clusters using high-resolution cosmological simulations of nine clusters that resolve formation of cluster galaxies. The simulations of each cluster are performed with the shock-capturing eulerian adaptive mesh refinement N-body+gasdynamics ART code with and without radiative cooling. We show that dissipation and associated galaxy formation increase the total baryon fractions within radii as large as the virial radius. The effect is the strongest within cluster cores, where the simulations with cooling have baryon fractions larger than the universal value, while the fraction of baryons in adiabatic simulations are smaller than universal. At larger radii (r >~ r_500) the cumulative baryon fractions in simulations with cooling are close to, while those in the adiabatic runs remain below than, the universal value. The gas fractions in simulations with dissipation are reduced by ~20-40% at r<0.3r_vir and ~10% at larger radii compared to the adiabatic runs, because a fraction of gas is converted into stars. There is an indication that gas fractions within different radii increase with increasing cluster mass as f_gas ~ M_vir^0.2. We find that the total baryon fraction within the virial radius does not evolve with time in both adiabatic simulations and in simulations with cooling. Finally, to evaluate systematic uncertainties in the baryon fraction in cosmological simulations we present a comparison of gas fractions in our adiabatic simulations to re-simulations of the same objects with the entropy-conserving SPH code Gadget. The cumulative gas fraction profiles in the two sets of simulations on average agree to <~3% at r/r_vir>0.2, but differ systematically by up to 10% at small radii.

We study motions of galaxies in galaxy clusters formed in the concordance LCDM cosmology. We use high-resolution cosmological simulations that follow dynamics of dark matter and gas and include various physical processes critical for galaxy formation: gas cooling, heating and star formation. Analysing motions of galaxies and the properties of intracluster gas in the sample of eight simulated clusters at z=0, we study velocity dispersion profiles of the dark matter, gas, and galaxies. We measure the mean velocity of galaxy motions and gas sound speed as a function of radius and calculate the average Mach number of galaxy motions. The simulations show that galaxies, on average, move supersonically with the average Mach number of ~1.4, approximately independent of the cluster-centric radius. The supersonic motions of galaxies may potentially provide an important source of heating for the intracluster gas by driving weak shocks and via dynamical friction, although these heating processes appear to be inefficient in our simulations. We also find that galaxies move faster than the dark matter particles in clusters. The magnitude of the velocity bias, b_v~1.1, is, however, smaller than the bias estimated for subhalos in dissipationless simulations. Interestingly, we find velocity bias in the tangential component of the velocity dispersion, but not in the radial component. Finally, we find significant random bulk motions of gas. The typical gas velocities are of order ~20-30% of the gas sound speed. These random motions provide about 10% of the total pressure support in our simulated clusters. The non-thermal pressure support, if neglected, will bias measurements of the total mass in the hydrostatic analyses of the X-ray cluster observations.

The cooling of gas in the centers of dark matter halos is expected to lead to a more concentrated dark matter distribution. The response of dark matter to the condensation of baryons is usually calculated using the model of adiabatic contraction, which assumes spherical symmetry and circular orbits. In contrast, halos in the hierarchical structure formation scenarios grow via multiple violent mergers and accretion along filaments, and particle orbits in the halos are highly eccentric. We study the effects of the cooling of gas in the inner regions of halos using high-resolution cosmological simulations which include gas dynamics, radiative cooling, and star formation. We find that the dissipation of gas indeed increases the density of dark matter and steepens its radial profile in the inner regions of halos compared to the case without cooling. For the first time, we test the adiabatic contraction model in cosmological simulations and find that the standard model systematically overpredicts the increase of dark matter density in the inner 5% of the virial radius. We show that the model can be improved by a simple modification of the assumed invariant from M(r)r to M(r_av)r, where r and r_av are the current and orbit-averaged particle positions. This modification approximately accounts for orbital eccentricities of particles and reproduces simulation profiles to within 10-20%. We present analytical fitting functions that accurately describe the transformation of the dark matter profile in the modified model and can be used for interpretation of observations.

We study the radial distribution of subhalos and galaxies using high-resolution cosmological simulations of galaxy clusters formed in the concordance LCDM cosmology. In agreement with previous studies, we find that the radial distribution of subhalos is significantly less concentrated than that of the dark matter, when subhalos are selected using their present-day gravitationally bound mass. We show that the difference in the radial distribution is not a numerical artifact and is due to tidal stripping. The subhalos in the cluster core lose more than 70% of their initial mass since accretion, while the average tidal mass loss for halos near the virial radius is ~30%. This introduces a radial bias in the spatial distribution of subhalos when they are selected using their tidally truncated mass. We demonstrate that the radial bias disappears almost entirely if subhalos are selected using their mass or circular velocity at the accretion epoch. The comparisons of the results of dissipationless simulations to the observed distribution of galaxies in clusters are therefore sensitive to the selection criteria used to select subhalo samples. Using the simulations that include cooling and starformation, we show that the radial distribution of subhalos is in reasonable agreement with the observed radial distribution of galaxies in clusters for 0.1<R/R200<2.0, if subhalos are selected using the stellar mass of galaxies. The radial bias is minimized in this case because the stars are located in the centers of dark matter subhalos and are tightly bound. The stellar mass of an object is therefore approximately conserved as the dark matter is stripped from the outer regions. Nevertheless, the concentration of the radial distribution of galaxies is systematically lower than that of the dark matter.

We analyze the effect of dissipation on the shapes of dark matter (DM) halos using high-resolution cosmological gasdynamics simulations of clusters and galaxies in the LCDM cosmology. We find that halos formed in simulations with gas cooling are significantly more spherical than corresponding halos formed in adiabatic simulations. Gas cooling results in an average increase of the principle axis ratios of halos by ~ 0.2-0.4 in the inner regions. The systematic difference decreases slowly with radius but persists almost to the virial radius. We argue that the differences in simulations with and without cooling arise both during periods of quiescent evolution, when gas cools and condenses toward the center, and during major mergers. We perform a series of high-resolution N-body simulations to study the shapes of remnants in major mergers of DM halos and halos with embedded stellar disks. In the DM halo-only mergers, the shape of the remnants depends only on the orbital angular momentum of the encounter and not on the internal structure of the halos. However, significant shape changes in the DM distribution may result if stellar disks are included. In this case the shape of the DM halos is correlated with the morphology of the stellar remnants.

We study the effects of radiative cooling, star formation and stellar feedback on the properties and evolution of galaxy clusters using high-resolution Adaptive Mesh Refinement N-body+gasdynamics simulations of clusters forming in the LCDM universe. Cooling leads to the condensation of gas in the inner regions of clusters, which in turn leads to steepening of the dark matter profile. The cooling gas is replaced by the higher-entropy gas from the outer regions, which raises the entropy and temperature of gas in the cluster core. The magnitude of these effects is likely overestimated in the current simulations because they suffer from the overcooling problem: a much larger fraction of baryons is in the form of cold gas and stars than is observed. We find that the thermal stellar feedback alone does not remedy this problem. Additional ad-hoc preheating can lower the amount of cold gas but a simple uniform preheating results in incorrect star formation history, as it delays the bulk of star formation until z<1. Our analysis shows that the overcooling in a cluster as a whole is really the overcooling in the central galaxy and its progenitors at high redshifts. This indicates that an additional heating mechanism that can continuously heat the gas in the cluster core is required to reproduce the observed cluster properties. Energy injection by the Active Galactic Nuclei, which may provide such heating, may thus be an important missing ingredient in the current theoretical models of cluster formation.

Recently, high-resolution Chandra observations revealed the existence of very sharp features in the X-ray surface brightness and temperature maps of several clusters (Vikhlinin et. al., 2001). These features, called ``cold fronts'', are characterized by an increase in surface brightness by a factor >2 over 10-50 kpc, accompanied by a drop in temperature of a similar magnitude. The existence of such sharp gradients can be used to put interesting constraints on the physics of the intracluster medium (ICM), if their mechanism and longevity are well understood. Here, we present results of a search for cold fronts in high-resolution simulations of galaxy clusters in cold dark matter (CDM) models. We show that sharp gradients with properties similar to those of observed cold fronts naturally arise in cluster mergers when the shocks heat gas surrounding the merging sub-cluster, while its dense core remains relatively cold. The compression induced by supersonic motions and shock heating during the merger enhance the amplitude of gas density and temperature gradients across the front. Our results indicate that cold fronts are non-equilibrium transient phenomena and can be observed for a period of less than a billion years. We show that the velocity and density fields of gas surrounding the cold front can be very irregular which would complicate analyses aiming to put constraints on the physical conditions of the intracluster medium in the vicinity of the front.

Galaxy clusters are potentially powerful probes of the large-scale velocity field in the Universe because their peculiar velocity can be estimated directly via the kinematic Sunyaev-Zeldovich effect (kSZ). Using high-resolution cosmological simulations of an evolving cluster of galaxies, we evaluate how well the average velocity obtained via a kSZ measurement reflects the actual cluster peculiar velocity. We find that the internal velocities in the intracluster gas are comparable to the overall cluster peculiar velocity, 20 to 30% of the sound speed even when a cluster is relatively relaxed. Nevertheless, the velocity averaged over the kSZ map inside a circular aperture matched to the cluster virial region provides an unbiased estimate of a cluster's radial peculiar velocity with a dispersion of 50 to 100 km/s, depending on the line of sight and dynamical state of the cluster. This dispersion puts a lower limit on the accuracy with which cluster peculiar velocity can be measured. Although the dispersion of the average is relatively small, the velocity distribution is broad; regions of low signal must be treated with care to avoid bias. We discuss the extent to which systematic errors might be modelled, and the resulting limitations on using galaxy clusters as cosmological velocity tracers.

We report the results of our continued study of arcminute scale anisotropy in the Cosmic Microwave Background (CMB) with the Berkeley-Illinois-Maryland Array. The survey consists of ten independent fields selected for low infrared dust emission and lack of bright radio point sources. With observations from the Very Large Array (VLA) at 4.8 GHz, we have identified point sources which could act as contaminants in estimates of the CMB power spectrum and removed them in the analysis. Modeling the observed power spectrum with a single flat band power with average multipole of \ell_{eff}=6864, we find dT=14.2^{+4.8}_{-6.0} \mu K at 68% confidence. The signal in the visibility data exceeds the expected contribution from instrumental noise with 96.5% confidence. We have also divided the data into two bins corresponding to different spatial resolutions in the power spectrum. We find dT_1=16.6^{+5.3}_{-5.9} \mu K at 68% confidence for CMB flat band power described by an average multipole of \ell_{eff}=5237 and dT_2<26.5 \mu K at 95% confidence for \ell_{eff}=8748.

We present Sunyaev-Zel'dovich Effect (SZE) imaging observations of three distant (z > 0.8) and highly X-ray luminous clusters of galaxies, Cl1226+33, Cl0152-13 and MS1054-03. Two of the clusters, Cl1226+33 and Cl0152-13, were recently discovered in deep ROSAT x-ray images. Their high X-ray luminosity suggests that they are massive systems which, if confirmed, would provide strong constraints on the cosmological parameters of structure formation models. Our Sunyaev-Zel'dovich Effect data provide confirmation that they are massive clusters similar to the well studied cluster MS1054-03. Assuming the clusters have the same gas mass fraction as that derived from SZE measurements of eighteen known massive clusters, we are able to infer their mass and electron temperature from the SZE data. The derived electron temperatures are 9.8, 8.7, and 10.4 keV, respectively, and we infer total masses of ~2 x 10^14 h^-1 Msun within a radius of 65 arcsec (340 h^{-1} kpc) for all three clusters. For Cl0152-13 and MS1054-03 we find good agreement between our SZE derived temperatures and those inferred from X-ray spectroscopy. No X-ray derived temperatures are available for Cl1226+33, and thus the SZE data provide the first confirmation that it is indeed a massive system. The demonstrated ability to determine cluster temperatures and masses from SZE observations without access to X-ray data illustrates the power of using deep SZE surveys to probe the distant universe.

Constraints on the clustered mass density \Omega_m of the universe derived from the mean intracluster gas fraction of X-ray clusters may be biased by a single-phase assumption for the thermodynamic structure of the intracluster medium (ICM). We propose a descriptive model for multiphase structure in which a spherically symmetric ICM contains isobaric density perturbations with radially dependent variance \sigma^2(r)=\sigma_c^2 (1+r^2/r_c^2)^{-\eps}. The model extends the work of Gunn & Thomas (1996) which assumed radially independent density fluctuations thoughout the ICM. Fixing the X-ray emission profile and emission weighted temperature, we explore two independently observable signatures of the model in the {\sigma_c,\eps} space. For bremsstrahlung dominated emission, the central Sunyaev--Zeldovich (SZ) decrement in the multiphase case is increased over the single-phase case and multiphase X-ray spectra in the range 0.1-20 keV are flatter in the continuum and exhibit stronger low energy emission lines than their single-phase counterpart. We quantify these effects for a fiducial 10^8 K cluster and demonstrate how the combination of SZ and X-ray spectroscopy can be used to identify a preferred location in the {\sigma_c,\eps} plane. From these parameters, the correct value of mean ICM gas fraction in the multiphase model results, allowing an unbiased estimate of \Omega_m to be recovered. The consistency of recent determinations of the Hubble constant from SZ and X-ray observations with values determined by other methods suggests that biases in ICM gas fractions are small, \ltsim 20%.