We jointly analyze Bolocam Sunyaev-Zeldovich (SZ) effect and Chandra X-ray data for a set of 45 clusters to derive gas density and temperature profiles without using spectroscopic information. The sample spans the mass and redshift range $3 \times 10^{14} M_{\odot} \le M_{500} \le 25 \times 10^{14} M_{\odot}$ and $0.15\le z \le 0.89$. We define cool-core (CC) and non-cool core (NCC) subsamples based on the central X-ray luminosity, and 17/45 clusters are classified as CC. In general, the profiles derived from our analysis are found to be in good agreement with previous analyses, and profile constraints beyond $r_{500}$ are obtained for 34/45 clusters. In approximately 30% of the CC clusters our analysis shows a central temperature drop with a statistical significance of $>3\sigma$; this modest detection fraction is due mainly to a combination of coarse angular resolution and modest S/N in the SZ data. Most clusters are consistent with an isothermal profile at the largest radii near $r_{500}$, although 9/45 show a significant temperature decrease with increasing radius. The sample mean density profile is in good agreement with previous studies, and shows a minimum intrinsic scatter of approximately 10% near $0.5 \times r_{500}$. The sample mean temperature profile is consistent with isothermal, and has an intrinsic scatter of approximately 50% independent of radius. This scatter is significantly higher compared to earlier X-ray-only studies, which find intrinsic scatters near 10%, likely due to a combination of unaccounted for non-idealities in the SZ noise, projection effects, and sample selection.

We describe Sunyaev-Zel'dovich (SZ) effect measurements and analysis of the intracluster medium (ICM) pressure profiles of a set of 45 massive galaxy clusters imaged using Bolocam at the Caltech Submillimeter Observatory. We have used masses determined from Chandra X-ray observations to scale each cluster's profile by the overdensity radius R500 and the mass-and-redshift-dependent normalization factor P500. We deproject the average pressure profile of our sample into 13 logarithmically spaced radial bins between 0.07R500 and 3.5R500. We find that a generalized Navarro, Frenk, and White (gNFW) profile describes our data with sufficient goodness-of-fit and best-fit parameters (C500, alpha, beta, gamma, P0 = 1.18, 0.86, 3.67, 0.67, 4.29). We also use the X-ray data to define cool-core and disturbed subsamples of clusters, and we constrain the average pressure profiles of each of these subsamples. We find that given the precision of our data the average pressure profiles of disturbed and cool-core clusters are consistent with one another at R>~0.15R500, with cool-core systems showing indications of higher pressure at R<~0.15R500. In addition, for the first time, we place simultaneous constraints on the mass scaling of cluster pressure profiles, their ensemble mean profile, and their radius-dependent intrinsic scatter between 0.1R500 and 2.0R500. The scatter among profiles is minimized at radii between ~0.2R500 and ~0.5R500, with a value of ~20%. The best-fit mass scaling has a power-law slope of 0.49, which is shallower than the nominal prediction of 2/3 from self-similar hydrostatic equilibrium models. These results for the intrinsic scatter and mass scaling are largely consistent with previous analyses, most of which have relied heavily on X-ray derived pressures of clusters at significantly lower masses and redshifts compared to our sample.

We describe in detail our characterization of the compact radio source population in 140 GHz Bolocam observations of a set of 45 massive galaxy clusters. We use a combination of 1.4 and 30 GHz data to select a total of 28 probable cluster-member radio galaxies and also to predict their 140 GHz flux densities. All of these galaxies are steep-spectrum radio sources and they are found preferentially in the cool-core clusters within our sample. In particular, 11 of the 12 brightest cluster member radio sources are associated with cool-core systems. Although none of the individual galaxies are robustly detected in the Bolocam data, the ensemble-average flux density at 140 GHz is consistent with, but slightly lower than, the extrapolation from lower frequencies assuming a constant spectral index. In addition, our data indicate an intrinsic scatter of 30 percent around the power-law extrapolated flux densities at 140 GHz, although our data do not tightly constrain this scatter. For our cluster sample, which is composed of high-mass and moderate-redshift systems, we find that the maximum fractional change in the Sunyaev-Zel'dovich signal integrated over any single cluster due to the presence of these radio sources is 20 percent, and only 1/4 of the clusters show a fractional change of more than 1 percent. The amount of contamination is strongly dependent on cluster morphology, and nearly all of the clusters with more than 1 percent contamination are cool-core systems. This result indicates that radio contamination is not significant compared to current noise levels in 140 GHz images of massive clusters and is in good agreement with the level of radio contamination found in previous results based on lower frequency data or simulations.

X-ray observations of galaxy clusters reveal a large range of morphologies with various degrees of disturbance, showing that the assumptions of hydrostatic equilibrium and spherical shape which are used to determine the cluster mass from X-ray data are not always satisfied. It is therefore important for the understanding of cluster properties as well as for cosmological applications to detect and quantify substructure in X-ray images of galaxy clusters. Two promising methods to do so are power ratios and center shifts. Since these estimators can be heavily affected by Poisson noise and X-ray background, we performed an extensive analysis of their statistical properties using a large sample of simulated X-ray observations of clusters from hydrodynamical simulations. We quantify the measurement bias and error in detail and give ranges where morphological analysis is feasible. A new, computationally fast method to correct for the Poisson bias and the X-ray background contribution in power ratio and center shift measurements is presented and tested for typical XMM-Newton observational data sets. We studied the morphology of 121 simulated cluster images and establish structure boundaries to divide samples into relaxed, mildly disturbed and disturbed clusters. In addition, we present a new morphology estimator - the peak of the 0.3-1 r500 P3/P0 profile to better identify merging clusters. The analysis methods were applied to a sample of 80 galaxy clusters observed with XMM-Newton. We give structure parameters (P3/P0 in r500, w and P3/P0_max) for all 80 observed clusters. Using our definition of the P3/P0 (w) substructure boundary, we find 41% (47%) of our observed clusters to be disturbed.

Measuring the intrinsic shape and orientation of dark matter (DM) and intracluster (IC) gas in galaxy clusters is crucial to constraining their formation and evolution, and for enhancing the use of clusters as more precise cosmological probes. Extending our previous works, we present for the first time results from a triaxial joint analysis of the galaxy cluster Abell 1835, by means of X-ray, strong lensing (SL) and Sunyaev Zel'dovich (SZ) data. We parametrically reconstruct the full three-dimensional structure (triaxial shape and principal axis orientation) of both the DM and the IC gas, and the level of non-thermal pressure of the IC gas. We find that the intermediate-major and minor-major axis ratios of the DM are 0.71+/-0.08 and 0.59+/-0.05, respectively, and the major axis of the DM halo is inclined with respect to the line of sight at 18.3+/-5.2 deg. We present the first observational measurement of the non-thermal pressure out to R_{200}, which has been evaluated to be a few percent of the total energy budget in the internal regions, while reaching approximately 20% in the outer volumes. We discuss the implications of our method for the viability of the CDM scenario, focusing on the concentration parameter C and the inner slope of the DM gamma in order to test the cold dark matter (CDM) paradigm for structure formation: we measure gamma=1.01+/-0.06 and C=4.32+/-0.44, values which are close to the predictions of the CDM model. The combination of X-ray/SL data at high spatial resolution, capable of resolving the cluster core, with the SZ data, which are more sensitive to the cluster outer volume, allows us to characterize the level and the gradient of the gas entropy distribution and non-thermal pressure out to R_{200}, breaking the degeneracy among the physical models describing the thermal history of the ICM.

We present initial results from our ongoing program to image the Sunyaev-Zel'dovich (SZ) effect in galaxy clusters at 143 GHz using Bolocam; five clusters and one blank field are described in this manuscript. The images have a resolution of 58 arcsec and a radius of 6-7 arcmin, which is approximately r500 - 2r500 for these clusters. The beam-smoothed RMS is ~10 uK_CMB in these images; with this sensitivity we are able to detect SZ signal to beyond r500 in binned radial profiles. We have fit our images to beta and Nagai models, fixing spherical symmetry or allowing for ellipticity in the plane of the sky, and we find that the best-fit parameter values are in general consistent with those obtained from other X-ray and SZ data. Our data show no clear preference for the Nagai model or the beta model due to the limited spatial dynamic range of our images. However, our data show a definitive preference for elliptical models over spherical models. The weighted mean ellipticity of the five clusters is 0.27 +- 0.03, consistent with results from X-ray data. Additionally, we obtain model-independent estimates of Y500, the integrated SZ y-parameter over the cluster face to a radius of r500, with systematics-dominated uncertainties of ~10%. Our Y500 values, which are free from the biases associated with model-derived Y500 values, scale with cluster mass in a way that is consistent with both self-similar predictions and expectations of a 10% intrinsic scatter.

We study the substructure statistics of a representative sample of galaxy clusters by means of two currently popular substructure characterisation methods, power ratios and centroid shifts. We use the 31 clusters from the REXCESS sample, compiled from the southern ROSAT All-Sky cluster survey REFLEX with a morphologically unbiased selection in X-ray luminosity and redshift, all of which have been reobserved with XMM-Newton. We investigate the uncertainties of the substructure parameters and examine the dependence of the results on projection effects, finding that the uncertainties of the parameters can be quite substantial. Thus while the quantification of the dynamical state of individual clusters with these parameters should be treated with extreme caution, these substructure measures provide powerful statistical tools to characterise trends of properties in large cluster samples. The centre shift parameter, w, is found to be more sensitive in general. For the REXCESS sample neither the occurence of substructure nor the presence of cool cores depends on cluster mass. There is a significant anti-correlation between the existence of substantial substructure and cool cores. The simulated clusters show on average larger substructure parameters than the observed clusters, a trend that is traced to the fact that cool regions are more pronounced in the simulated clusters, leading to stronger substructure measures in merging clusters and clusters with offset cores. Moreover, the frequency of cool regions is higher in the simulations than in the observations, implying that the description of the physical processes shaping cluster formation in the simulations requires further improvement.

We present a method to recover mass profiles of galaxy clusters by combining data on thermal Sunyaev-Zeldovich (tSZ) and X-ray imaging, thereby avoiding to use any information on X-ray spectroscopy. This method, which represents a development of the geometrical deprojection technique presented in Ameglio et al. (2007), implements the solution of the hydrostatic equilibrium equation. In order to quantify the efficiency of our mass reconstructions, we apply our technique to a set of hydrodynamical simulations of galaxy clusters. We propose two versions of our method of mass reconstruction. Method 1 is completely model-independent, while Method 2 assumes instead the analytic mass profile proposed by Navarro et al. (1997) (NFW). We find that the main source of bias in recovering the mass profiles is due to deviations from hydrostatic equilibrium, which cause an underestimate of the mass of about 10 per cent at r_500 and up to 20 per cent at the virial radius. Method 1 provides a reconstructed mass which is biased low by about 10 per cent, with a 20 per cent scatter, with respect to the true mass profiles. Method 2 proves to be more stable, reducing the scatter to 10 per cent, but with a larger bias of 20 per cent, mainly induced by the deviations from equilibrium in the outskirts. To better understand the results of Method 2, we check how well it allows to recover the relation between mass and concentration parameter. When analyzing the 3D mass profiles we find that including in the fit the inner 5 per cent of the virial radius biases high the halo concentration. Also, at a fixed mass, hotter clusters tend to have larger concentration. Our procedure recovers the concentration parameter essentially unbiased but with a scatter of about 50 per cent.

Galaxy clusters exhibit regular scaling relations among their bulk properties. These relations establish vital links between halo mass and cluster observables. Precision cosmology studies that depend on these links benefit from a better understanding of scatter in the mass-observable scaling relations. Here we study the role of merger processes in introducing scatter into the $M$-$T_{\rm X}$ relation, using a sample of 121 galaxy clusters simulated with radiative cooling and supernova feedback, along with three statistics previously proposed to measure X-ray surface brightness substructure. These are the centroid variation ($w$), the axial ratio ($\eta$), and the power ratios ($P_{20}$ and $P_{30}$). We find that in this set of simulated clusters, each substructure measure is correlated with a cluster's departures $\delta \ln T_{\rm X}$ and $\delta \ln M$ from the mean $M$-$T_{\rm X}$ relation, both for emission-weighted temperatures $T_{\rm EW}$ and for spectroscopic-like temperatures $T_{\rm SL}$, in the sense that clusters with more substructure tend to be cooler at a given halo mass. In all cases, a three-parameter fit to the $M$-$T_{\rm X}$ relation that includes substructure information has less scatter than a two-parameter fit to the basic $M$-$T_{\rm X}$ relation.

We present two non-parametric deprojection methods aimed at recovering the three-dimensional density and temperature profiles of galaxy clusters from spatially resolved thermal Sunyaev-Zeldovich (tSZ) and X-ray surface brightness maps, thus avoiding the use of X-ray spectroscopic data. In both methods, clusters are assumed to be spherically symmetric and modeled with an onion-skin structure. The first method follows a direct geometrical approach. The second method is based on the maximization of a single joint (tSZ and X-ray) likelihood function, which allows one to fit simultaneously the two signals by following a Monte Carlo Markov Chain approach. These techniques are tested against a set of cosmological simulations of clusters, with and without instrumental noise. We project each cluster along the three orthogonal directions defined by the principal axes of the momentum of inertia tensor. This enables us to check any bias in the deprojection associated to the cluster elongation along the line of sight. After averaging over all the three projection directions, we find an overall good reconstruction, with a small (<~10 per cent) overestimate of the gas density profile. This turns into a comparable overestimate of the gas mass within the virial radius, which we ascribe to the presence of residual gas clumping. Apart from this small bias the reconstruction has an intrinsic scatter of about 5 per cent, which is dominated by gas clumpiness. Cluster elongation along the line of sight biases the deprojected temperature profile upwards at r<~0.2r_vir and downwards at larger radii. A comparable bias is also found in the deprojected temperature profile. Overall, this turns into a systematic underestimate of the gas mass, up to 10 percent. (Abridged)

The angular-diameter distance D_A of a galaxy cluster can be measuread by combining its X-ray emission with the cosmic microwave background fluctution due to the Sunyaev-Zeldovich effect. The application of this distance indicator usually assumes that the cluster is spherically symmetric, the gas is distributed according to the isothermal beta-model, and the X-ray temperature is an unbiased measure of the electron temperature. We test these assumptions with galaxy clusters extracted from an extended set of cosmological N-body/hydrodynamical simulations of a LCDM concordance cosmology, which include the effect of radiative cooling, star formation and energy feedback from supernovae. We find that, due to the steep temperature gradients which are present in the central regions of simulated clusters, the assumption of isothermal gas leads to a significant underestimate of D_A. This bias is efficiently corrected by using the polytropic version of the beta-model to account for the presence of temperature gradients. In this case, once irregular clusters are removed, the correct value of D_A is recovered with a ~ 5 per cent accuracy on average, with a ~ 20 per cent intrinsic scatter due to cluster asphericity. This result is valid when using either the electron temperature or a spectroscopic-like temperature. When using instead the emission-weighted definition for the temperature of the simulated clusters, D_A is biased low by \~ 20 per cent. We discuss the implications of our results for an accurate determination of the Hubble constant H_0 and of the density parameter Omega_m. We find that H_0 can be potentially recovered with exquisite precision, while the resulting estimate of Omega_m, which is unbiased, has typical errors Delta(Omega_m) ~ 0.05.