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.

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.

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

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.

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

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.

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.

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.

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.

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 use ROSAT All Sky Survey (RASS) broadband X-ray images and the optical clusters identified from SDSS DR7 to estimate the X-ray luminosities around $\sim 65,000$ candidate clusters with masses $\ga 10^{13}\msunh$ based on an Optical to X-ray (OTX) code we develop. We obtain a catalogue with X-ray luminosity for each cluster. This catalog contains 817 clusters (473 at redshift $z\le 0.12$) with $S/N> 3$ in X-ray detection. We find about $65\%$ of these X-ray clusters have their most massive member located near the X-ray flux peak; for the rest $35\%$, the most massive galaxy is separated from the X-ray peak, with the separation following a distribution expected from a NFW profile. We investigate a number of correlations between the optical and X-ray properties of these X-ray clusters, and find that: the cluster X-ray luminosity is correlated with the stellar mass (luminosity) of the clusters, as well as with the stellar mass (luminosity) of the central galaxy and the mass of the halo, but the scatter in these correlations is large. Comparing the properties of X-ray clusters of similar halo masses but having different X-ray luminosities, we find that massive halos with masses $\ga 10^{14}\msunh$ contain a larger fraction of red satellite galaxies when they are brighter in X-ray. ... A cluster catalog containing the optical properties of member galaxies and the X-ray luminosity is available at {\it http://gax.shao.ac.cn/data/Group.html}.

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.

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.

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.

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]

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 use the ROSAT all sky survey X-ray cluster catalogs and the optical SDSS DR7 galaxy and group catalogs to cross-identify X-ray clusters with their optical counterparts, resulting in a sample of 201 X-ray clusters in the sky coverage of SDSS DR7. We investigate various correlations between the optical and X-ray properties of these X-ray clusters, and find that the following optical properties are correlated with the X-ray luminosity: the central galaxy luminosity, the central galaxy mass, the characteristic group luminosity ($\propto \Lx^{0.43}$), the group stellar mass ($\propto \Lx^{0.46}$), with typical 1-$\sigma$ scatter of $\sim 0.67$ in $\log \Lx$. Using the observed number distribution of X-ray clusters, we obtain an unbiased scaling relation between the X-ray luminosity, the central galaxy stellar mass and the characteristic satellite stellar mass as ${\log L_X} = -0.26 + 2.90 [\log (M_{\ast, c} + 0.26 M_{\rm sat}) -12.0]$ (and in terms of luminosities, as ${\log L_X} = -0.15 + 2.38 [\log (L_{c} + 0.72 L_{\rm sat}) -12.0]$). We find that the systematic difference between different halo mass estimations, e.g., using the ranking of characteristic group stellar mass or using the X-ray luminosity scaling relation can be used to constrain cosmology. Comparing the properties of groups of similar stellar mass (or optical luminosities) and redshift that are X-ray luminous or under-luminous, we find that X-ray luminous groups have more faint satellite galaxies and higher red fraction in their satellites. The cross-identified X-ray clusters together with their optical properties are provided in Appendix B.

We propose a simple estimator for the gravitational potential of cluster-size halos using the temperature and density profiles of the intracluster gas based on the assumptions of hydrostatic equilibrium and spherical symmetry. Using high resolution cosmological simulations of galaxy clusters, we show that the scaling relation between this estimator and the gravitational potential has a small intrinsic scatter of ~8%-15%, and it is insensitive to baryon physics outside the cluster core. The slope and the normalization of the scaling relation vary weakly with redshift, and they are relatively independent of the choice of radial range used and the dynamical states of the clusters. The results presented here provide a possible way for using the cluster potential function as an alternative to the cluster mass function in constraining cosmology using 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 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.

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.

The distribution of galaxies in position and velocity around the centers of galaxy clusters encodes important information about cluster mass and structure. Using the maxBCG galaxy cluster catalog identified from imaging data obtained in the Sloan Digital Sky Survey, we study the BCG-galaxy velocity correlation function. By modeling its non-Gaussianity, we measure the mean and scatter in velocity dispersion at fixed richness. The mean velocity dispersion increases from 202+/-10 km/s for small groups to more than 854+/-102 km/s for large clusters. We show the scatter to be at most 40.5+/-3.5%, declining to 14.9+/-9.4% in the richest bins. We test our methods in the C4 cluster catalog, a spectroscopic cluster catalog produced from the Sloan Digital Sky Survey DR2 spectroscopic sample, and in mock galaxy catalogs constructed from N-body simulations. Our methods are robust, measuring the scatter to well within one-sigma of the true value, and the mean to within 10%, in the mock catalogs. By convolving the scatter in velocity dispersion at fixed richness with the observed richness space density function, we measure the velocity dispersion function of the maxBCG galaxy clusters. Although velocity dispersion and richness do not form a true mass-observable relation, the relationship between velocity dispersion and mass is theoretically well characterized and has low scatter. Thus our results provide a key link between theory and observations up to the velocity bias between dark matter and galaxies.

Our previous analyses of radio Doppler and ranging data from distant spacecraft in the solar system indicated that an apparent anomalous acceleration is acting on Pioneer 10 and 11, with a magnitude a_P ~ 8 x 10^{-8} cm/s^2, directed towards the Sun (anderson,moriond). Much effort has been expended looking for possible systematic origins of the residuals, but none has been found. A detailed investigation of effects both external to and internal to the spacecraft, as well as those due to modeling and computational techniques, is provided. We also discuss the methods, theoretical models, and experimental techniques used to detect and study small forces acting on interplanetary spacecraft. These include the methods of radio Doppler data collection, data editing, and data reduction. There is now further data for the Pioneer 10 orbit determination. The extended Pioneer 10 data set spans 3 January 1987 to 22 July 1998. [For Pioneer 11 the shorter span goes from 5 January 1987 to the time of loss of coherent data on 1 October 1990.] With these data sets and more detailed studies of all the systematics, we now give a result, of a_P = (8.74 +/- 1.33) x 10^{-8} cm/s^2. (Annual/diurnal variations on top of a_P, that leave a_P unchanged, are also reported and discussed.)

We conclude that Katz's proposal (anisotropic heat reflection off of the back of the spacecraft high-gain antennae, the heat coming from the RTGs) does not provide enough power and so can not explain the Pioneer anomaly.

We conclude that Murphy's proposal (radiation of the power of the main-bus electrical systems from the rear of the craft) can not explain the anomalous Pioneer acceleration.

Recently we reported that radio Doppler data generated by NASA's Deep Space Network (DSN) from the Pioneer 10 and 11 spacecraft indicate an apparent anomalous, constant, spacecraft acceleration with a magnitude $\sim 8.5\times 10^{-8}$ cm s$^{-2}$, directed towards the Sun (gr-qc/9808081). Analysis of similar Doppler and ranging data from the Galileo and Ulysses spacecraft yielded ambiguous results for the anomalous acceleration, but it was useful in that it ruled out the possibility of a systematic error in the DSN Doppler system that could easily have been mistaken as a spacecraft acceleration. Here we present some new results, including a critique suggestions that the anomalous acceleration could be caused by collimated thermal emission. Based partially on a further data for the Pioneer 10 orbit determination, the data now spans January 1987 to July 1998, our best estimate of the average Pioneer 10 acceleration directed towards the Sun is $\sim 7.5 \times 10^{-8}$ cm s$^{-2}$.

Radio metric data from the Pioneer 10/11, Galileo, and Ulysses spacecraft indicate an apparent anomalous, constant, acceleration acting on the spacecraft with a magnitude $\sim 8.5\times 10^{-8}$ cm/s$^2$, directed towards the Sun. Two independent codes and physical strategies have been used to analyze the data. A number of potential causes have been ruled out. We discuss future kinematic tests and possible origins of the signal.

The year 1996 will mark the initiation of a number of new space missions to the planet Mars from which we expect to obtain a rich set of data, including spacecraft radio tracking data. Anticipating these events, we have analyzed the feasibility of testing a violation of the strong equivalence principle (SEP) with Earth-Mars ranging. Using analytic and numerical methods, we have demonstrated that ranging data can provide a useful estimate of the SEP parameter $\eta$. Two estimates of the predicted accuracy are quoted, one based on conventional covariance analysis, and the other based on ``modified worst case'' analysis, which assumes that systematic errors dominate the experiment. If future Mars missions provide ranging measurements with an accuracy of $\sigma$ meters, after ten years of ranging the expected accuracy for the parameter $\eta$ will be of order $\sigma_\eta\approx (1-12)\times 10^{-4}\sigma $. In addition, these ranging measurements will provide a significantly improved determination of the mass of the Jupiter system, independent of the test of the {\small SEP} polarization effect.

We review the results from {\it Mariner 10} regarding Mercury's gravity field and the results from radar ranging regarding topography. We discuss the implications of improving these results, including a determination of the polar component, as well as the opportunity to perform relativistic gravity tests with a future {\it Mercury Orbiter}. With a spacecraft placed in orbit with periherm at 400 km altitude, apherm at 16,800 km, period 13.45 hr and latitude of periherm at +30 deg, one can expect a significant improvement in our knowledge of Mercury's gravity field and geophysical properties. The 2000 Plus mission that evolved during the European Space Agency (ESA) {\it Mercury Orbiter} assessment study can provide a global gravity field complete through the 25th degree and order in spherical harmonics. If after completion of the main mission, the periherm could be lowered to 200 km altitude, the gravity field could be extended to 50th degree and order. We discuss the possibility that a search for a Hermean ionosphere could be performed during the mission phases featuring Earth occultations. Because of its relatively large eccentricity and close proximity to the Sun, Mercury's orbital motion provides one of the best solar-system tests of general relativity. Consequently, we emphasize the number of feasible relativistic gravity tests that can be performed within the context of the parameterized post-Newtonian formalism - a useful framework for testing modern gravitational theories. We pointed out that current results on relativistic precession of Mercury's perihelion are uncertain by 0.5 %, and we discuss the expected improvement using {\it Mercury Orbiter}. We discuss the importance of {\it Mercury Orbiter} for setting limits on a possible time variation in the

We know that the generally accepted theories of gravity and quantum mechanics are fundamentally incompatible. Thus, when we try to combine these theories, we must beware of physical pitfalls. Modern theories of quantum gravity are trying to overcome these problems. Any ideas must confront the present agreement with general relativity, but yet be free to wonder about not understood phenomena, such as the dark matter problem and the anomalous spacecraft data which we announce here. This all has led some ``intrepid" theorists to consider a new gravitational regime, that of antimatter. Even more ``daring" experimentalists are attempting, or considering attempting, the measurement of the gravitational force on antimatter, including low-energy antiprotons and, perhaps most enticing, antihydrogen.