The non-detection of dark matter (DM) particles in increasingly stringent laboratory searches has encouraged alternative gravity theories where gravity is sourced only from visible matter. Here, we consider whether such theories can pass a two-dimensional test posed by gravitational lensing -- to reproduce a particularly detailed Einstein ring in the core of the galaxy cluster Abell 3827. We find that when we require the lensing mass distribution to strictly follow the shape (ellipticity and position angle) of the light distribution of cluster member galaxies, intracluster stars, and the X-ray emitting intracluster medium, we cannot reproduce the Einstein ring, despite allowing the mass-to-light ratios of these visible components to freely vary with radius to mimic alternative gravity theories. Alternatively, we show that the detailed features of the Einstein ring are accurately reproduced by allowing a smooth, freely oriented DM halo in the lens model, with relatively small contributions from the visible components at a level consistent with their observed brightnesses. This dominant DM component is constrained to have the same orientation as the light from the intracluster stars, indicating that the intracluster stars trace the gravitational potential of this component. The Einstein ring of Abell 3827 therefore presents a new challenge for alternative gravity theories: not only must such theories find agreement between the total lensing mass and visible mass, but they must also find agreement between the projected sky distribution of the lensing mass and that of the visible matter, a more stringent test than has hitherto been posed by lensing data.

We examine the massive colliding cluster El Gordo, one of the most massive clusters at high redshift. We use a free-form lensing reconstruction method that avoids making assumptions about the mass distribution. We use data from the RELICS program and identify new multiply lensed system candidates. The new set of constraints and free-form method provides a new independent mass estimate of this intriguing colliding cluster. Our results are found to be consistent with earlier parametric models, indirectly confirming the assumptions made in earlier work. By fitting a double gNFW profile to the lens model, and extrapolating to the virial radius, we infer a total mass for the cluster of $M_{200c}=(1.08^{+0.65}_{-0.12})\times10^{15}$M$_{\odot}$. We estimate the uncertainty in the mass due to errors in the photometric redshifts, and discuss the uncertainty in the inferred virial mass due to the extrapolation from the lens model. We also find in our lens map a mass overdensity corresponding to the large cometary tail of hot gas, reinforcing its interpretation as a large tidal feature predicted by hydrodynamical simulations that mimic El Gordo. Finally, we discuss the observed relation between the plasma and the mass map, finding that the peak in the projected mass map may be associated with a large concentration of colder gas, exhibiting possible star formation. El Gordo is one of the first clusters that will be observed with JWST, which is expected to unveil new high redshift lensed galaxies around this interesting cluster, and provide a more accurate estimation of its mass.

We study the dynamics of Abell 370 (A370), a highly massive Hubble Frontier Fields galaxy cluster, using self-consistent three-dimensional N-body/hydrodynamical simulations. Our simulations are constrained by X-ray, optical spectroscopic and gravitational lensing, and Sunyaev-Zel'dovich (SZ) effect observations. Analyzing archival Chandra observations of A370 and comparing the X-ray morphology to the latest gravitational lensing mass reconstruction, we find offsets of ~30 kpc and ~100 kpc between the two X-ray surface brightness peaks and their nearest mass surface density peaks, suggesting that it is a merging system, in agreement with previous studies. Based on our dedicated binary cluster merger simulations, we find that initial conditions of the two progenitors with virial masses of 1.7 x 10^(15) M_sun and 1.6 x 10^(15) M_sun, an infall velocity of 3500 km/s, and an impact parameter of 100 kpc can explain the positions and the offsets between the peaks of the X-ray emission and mass surface density, the amplitude of the integrated SZ signal, and the observed relative line-of-sight velocity. Our simulations suggest that A370 is a major merger after the second core passage in the infalling phase, just before the third core passage. In this phase, the gas has not settled down in the gravitational potential well of the cluster, which explains why A370 does not follow closely the galaxy cluster scaling relations.

We present a systematic study of gas density perturbations in cool cores of high-mass galaxy clusters. We select 12 relaxed clusters from the Cluster Lensing And Supernova survey with Hubble (CLASH) sample and analyze their cool core features observed with the Chandra X-ray Observatory. We focus on the X-ray residual image characteristics after subtracting their global profile of the X-ray surface brightness distribution. We find that all the galaxy clusters in our sample have, at least, both one positive and one negative excess regions in the X-ray residual image, indicating the presence of gas density perturbations. We identify and characterize the locally perturbed regions using our detection algorithm, and extract X-ray spectra of the intracluster medium (ICM). The ICM temperature in the positive excess region is lower than that in the negative excess region, whereas the ICM in both regions is in pressure equilibrium in a systematic manner. These results indicate that gas sloshing in cool cores takes place in more than 80% of relaxed clusters (95% CL). We confirm this physical picture by analyzing synthetic X-ray observations of a cool-core cluster from a hydrodynamic simulation, finding that our detection algorithm can accurately extract both the positive and negative excess regions and can reproduce the temperature difference between them. Our findings support the picture that the gas density perturbations are induced by gas sloshing, and a large fraction of cool-core clusters have undergone gas sloshing, indicating that gas sloshing may be capable of suppressing runaway cooling of the ICM.

We present 90, 140, and 268GHz sub-arcminute resolution imaging of the Sunyaev-Zel'dovich effect (SZE) in MACSJ0717.5+3745. Our 90GHz SZE data result in a sensitive, 34uJy/bm map at 13" resolution using MUSTANG. Our 140 and 268GHz SZE imaging, with resolutions of 58" and 31" and sensitivities of 1.8 and 3.3mJy/beam respectively, was obtained using Bolocam. We compare these maps to a 2-dimensional pressure map derived from Chandra X-ray observations. Our MUSTANG data confirm previous indications from Chandra of a pressure enhancement due to shock-heated, >20keV gas immediately adjacent to extended radio emission seen in low-frequency radio maps. The MUSTANG data also detect pressure substructure that is not well-constrained by the X-ray data in the remnant core of a merging subcluster. We find that the small-scale pressure enhancements in the MUSTANG data amount to ~2% of the total pressure measured in the 140GHz Bolocam observations. The X-ray template also fails on larger scales to accurately describe the Bolocam data, particularly at the location of a subcluster known to have a high line of sight optical velocity (~3200km/s). Our Bolocam data are adequately described when we add an additional component - not described by a thermal SZE spectrum - coincident with this subcluster. Using flux densities extracted from our model fits, and marginalizing over the temperature constraints for the region, we fit a thermal+kinetic SZE spectrum to our data and find the subcluster has a best-fit line of sight proper velocity of 3600+3440/-2160km/s. This agrees with the optical velocity estimates for the subcluster. The probability of velocity<0 given our measurements is 2.1%. Repeating this analysis using flux densities measured non-parametrically results in a 3.4% probability of a velocity<=0. We note that this tantalizing result for the kinetic SZE is on resolved, subcluster scales.

We use hydrodynamical/N-body simulations to interpret the newly discovered Bullet-cluster-like merging cluster, ZwCl 0008.8+5215 (ZwCl 0008 hereafter), where a dramatic collision is apparent from multi-wavelength observations. We have been able to find a self-consistent solution for the radio, X-ray, and lensing phenomena by projecting an off-axis, binary cluster encounter viewed just after first core passage. A pair radio relics traces well the leading and trailing shock fronts that our simulation predict, providing constraints on the collision parameters. We can also account for the observed distinctive comet-like X-ray morphology and the positions of the X-ray peaks relative to the two lensing mass centroids and the two shock front locations. Relative to the Bullet cluster, the total mass is about 70% lower, ($1.2\pm0.1) \times 10^{15}$ Msun, with a correspondingly lower infall velocity, $1800\pm300$ km/s, and an impact parameter of $400\pm100$ kpc. As a result, the gas component of the infalling cluster is not trailing significantly behind the associated dark matter as in the case of the Bullet cluster. The degree of agreement we find between all the observables provides strong evidence that dark matter is effectively collisionless on large scales calling into question other claims and theories that advocate modified gravity.

The colliding cluster, CIZA J2242.8+5301, displays a spectacular, almost 2 Mpc long shock front with a radio based Mach number M ~ 5, that is puzzlingly large compared with the X-ray estimate of M ~ 2.5. The extent to which the X-ray temperature jump is diluted by cooler unshocked gas projected through the cluster currently lacks quantification. Thus, here we apply our self-consistent N-body/hydro-dynamical code (based on FLASH) to model this binary cluster encounter. We can account for the location of the shock front and also the elongated X-ray emission by tidal stretching of the gas and dark matter between the two cluster centers. The required total mass is $8.9 \times 10^{14}$ Msun with a 1.3:1 mass ratio favoring the southern cluster component. The relative velocity we derive is $\simeq 2500$ km/s initially between the two main cluster components, with an impact parameter of 120 kpc. This solution implies that the shock temperature jump derived from the low angular resolution X-ray satellite SUZAKU is underestimated by a factor of two, due to cool gas in projection, bringing the observed X-ray and radio estimates into agreement. We propose that the complex southern relics in CIZA J2242.8+5301, have been broken up as the southerly moving "back" shocked gas impacts the gas still falling in along the collision axis. Finally, we use our model to generate Compton-y maps to estimate the reduction in radio flux caused by the thermal Sunyaev-Zel'dovich (SZ) effect. At 30 GHz, this amounts to $\Delta S_n = -0.072$ mJy/arcmin$^2$ and $\Delta S_s = -0.075$ mJy/arcmin$^2$ at the locations of the northern and southern shock fronts respectively. Our model estimate agrees with previous empirical estimates that have inferred the measured radio spectra can be significantly affected by the SZ effect, with implications for charged particle acceleration models of the radio relics.

We report Chandra X-ray observations and optical weak-lensing measurements from Subaru/Suprime-Cam images of the double galaxy cluster Abell 2465 (z=0.245). The X-ray brightness data are fit to a beta-model to obtain the radial gas density profiles of the northeast (NE) and southwest (SW) sub-components, which are seen to differ in structure. We determine core radii, central temperatures, the gas masses within $r_{500c}$, and the total masses for the broader NE and sharper SW components assuming hydrostatic equilibrium. The central entropy of the NE clump is about two times higher than the SW. Along with its structural properties, this suggests that it has undergone merging on its own. The weak-lensing analysis gives virial masses for each substructure, which compare well with earlier dynamical results. The derived outer mass contours of the SW sub-component from weak lensing are more irregular and extended than those of the NE. Although there is a weak enhancement and small offsets between X-ray gas and mass centers from weak lensing, the lack of large amounts of gas between the two sub-clusters indicates that Abell 2465 is in a pre-merger state. A dynamical model that is consistent with the observed cluster data, based on the FLASH program and the radial infall model, is constructed, where the subclusters currently separated by ~1.2Mpc are approaching each other at ~2000km/s and will meet in ~0.4Gyr.

We present a joint analysis of Chandra X-ray observations, Bolocam thermal Sunyaev-Zel'dovich (SZ) effect observations, Hubble Space Telescope (HST) strong lensing data, and HST and Subaru Suprime-Cam weak lensing data. The multiwavelength dataset is used to constrain parametric models for the distribution of dark and baryonic matter in a sample of six massive galaxy clusters selected from the Cluster Lensing And Supernova survey with Hubble (CLASH). For five of the six clusters, the multiwavelength dataset is well described by a relatively simple model that assumes spherical symmetry, hydrostatic equilibrium, and entirely thermal pressure support. The joint analysis yields considerably better constraints on the total mass and concentration of the cluster compared to analysis of any one dataset individually. The subsample of five galaxy clusters is used to place an upper limit on the fraction of pressure support in the intracluster medium (ICM) due to nonthermal processes, such as turbulence and bulk flow of the gas. We constrain the nonthermal pressure fraction at r500c to be less than 0.11 at 95 percent confidence. This is in tension with state-of-the-art hydrodynamical simulations, which predict a nonthermal pressure fraction of approximately 0.25 at r500c for clusters of similar mass and redshift. This tension may be explained by the sample selection and/or our assumption of spherical symmetry.

Radio observations suggest that 3C 75, located in the dumbbell shaped galaxy NGC 1128 at the center of Abell 400, hosts two colliding jets. Motivated by this source, we perform three-dimensional hydrodynamical simulations using a modified version of the GPU-accelerated Adaptive-MEsh-Refinement hydrodynamical parallel code ($\mathit{GAMER}$) to study colliding extragalactic jets. We find that colliding jets can be cast into two categories: 1) bouncing jets, in which case the jets bounce off each other keeping their identities, and 2) merging jets, when only one jet emerges from the collision. Under some conditions the interaction causes the jets to break up into oscillating filaments of opposite helicity, with consequences for their downstream stability. When one jet is significantly faster than the other and the impact parameter is small, the jets merge; the faster jet takes over the slower one. In the case of merging jets, the oscillations of the filaments, in projection, may show a feature which resembles a double helix, similar to the radio image of 3C 75. Thus we interpret the morphology of 3C 75 as a consequence of the collision of two jets with distinctly different speeds at a small impact parameter, with the faster jet breaking up into two oscillating filaments.

The Yuan-Tseh Lee Array for Microwave Background Anisotropy (AMiBA) is a co-planar interferometer array operating at a wavelength of 3mm to measure the Sunyaev-Zeldovich effect (SZE) of galaxy clusters. In the first phase of operation -- with a compact 7-element array with 0.6m antennas (AMiBA-7) -- we observed six clusters at angular scales from 5\arcmin to 23\arcmin. Here, we describe the expansion of AMiBA to a 13-element array with 1.2m antennas (AMiBA-13), its subsequent commissioning, and our cluster SZE observing program. The most important changes compared to AMiBA-7 are (1) array re-configuration with baselines ranging from 1.4m to 4.8m covering angular scales from 2\arcmin to 11.5\arcmin, (2) thirteen new lightweight carbon-fiber-reinforced plastic (CFRP) 1.2m reflectors, and (3) additional correlators and six new receivers. From the AMiBA-13 SZE observing program, we present here maps of a subset of twelve clusters. In highlights, we combine AMiBA-7 and AMiBA-13 observations of Abell 1689 and perform a joint fitting assuming a generalized NFW pressure profile. Our cylindrically integrated Compton-y values for this cluster are consistent with the BIMA/OVRA, SZA, and Planck results. We report the first targeted SZE detection towards the optically selected galaxy cluster RCS J1447+0828, and we demonstrate the ability of AMiBA SZE data to serve as a proxy for the total cluster mass. Finally, we show that our AMiBA-SZE derived cluster masses are consistent with recent lensing mass measurements in the literature.

We measure the SZ signal toward a set of 47 clusters with a median mass of $9.5 \times 10^{14}$ M$_{\odot}$ and a median redshift of 0.40 using data from Planck and the ground-based Bolocam receiver. When Planck XMM-like masses are used to set the scale radius $\theta_{\textrm{s}}$, we find consistency between the integrated SZ signal, $Y_{\textrm{5R500}}$, derived from Bolocam and Planck based on gNFW model fits using A10 shape parameters, with an average ratio of $1.069 \pm 0.030$ (allowing for the $\simeq 5$% Bolocam flux calibration uncertainty). We also perform a joint fit to the Bolocam and Planck data using a modified A10 model with the outer logarithmic slope $\beta$ allowed to vary, finding $\beta = 6.13 \pm 0.16 \pm 0.76$ (measurement error followed by intrinsic scatter). In addition, we find that the value of $\beta$ scales with mass and redshift according to $\beta \propto M^{0.077 \pm 0.026} \times (1+z)^{-0.06 \pm 0.09}$. This mass scaling is in good agreement with recent simulations. We do not observe the strong trend of $\beta$ with redshift seen in simulations, though we conclude that this is most likely due to our sample selection. Finally, we use Bolocam measurements of $Y_{500}$ to test the accuracy of the Planck completeness estimate. We find consistency, with the actual number of Planck detections falling approximately $1 \sigma$ below the expectation from Bolocam. We translate this small difference into a constraint on the the effective mass bias for the Planck cluster cosmology results, with $(1-b) = 0.93 \pm 0.06$.

We present a weak-lensing analysis of the merging {\em Frontier Fields} (FF) cluster Abell~2744 using new Subaru/Suprime-Cam imaging. The wide-field lensing mass distribution reveals this cluster is comprised of four distinct substructures. Simultaneously modeling the two-dimensional reduced shear field using a combination of a Navarro--Frenk--White (NFW) model for the main core and truncated NFW models for the subhalos, we determine their masses and locations. The total mass of the system is constrained as $M_\mathrm{200c} = (2.06\pm0.42)\times10^{15}\,M_\odot$. The most massive clump is the southern component with $M_\mathrm{200c} = (7.7\pm3.4)\times10^{14}\,M_\odot$, followed by the western substructure ($M_\mathrm{200c} = (4.5\pm2.0)\times10^{14}\,M_\odot$) and two smaller substructures to the northeast ($M_\mathrm{200c} = (2.8\pm1.6)\times10^{14}\,M_\odot$) and northwest ($M_\mathrm{200c} = (1.9\pm1.2)\times10^{14}\,M_\odot$). The presence of the four substructures supports the picture of multiple mergers. Using a composite of hydrodynamical binary simulations we explain this complicated system without the need for a "slingshot" effect to produce the northwest X-ray interloper, as previously proposed. The locations of the substructures appear to be offset from both the gas ($87^{+34}_{-28}$ arcsec, 90\% CL) and the galaxies ($72^{+34}_{-53}$ arcsec, 90\% CL) in the case of the northwestern and western subhalos. To confirm or refute these findings, high resolution space-based observations extending beyond the current FF limited coverage to the west and northwestern area are essential.

We present scaling relations between the integrated Sunyaev-Zel'dovich Effect (SZE) signal, $Y_{\rm SZ}$, its X-ray analogue, $Y_{\rm X}\equiv M_{\rm gas}T_{\rm X}$, and total mass, $M_{\rm tot}$, for the 45 galaxy clusters in the Bolocam X-ray-SZ (BOXSZ) sample. All parameters are integrated within $r_{2500}$. $Y_{2500}$ values are measured using SZE data collected with Bolocam, operating at 140 GHz at the Caltech Submillimeter Observatory (CSO). The temperature, $T_{\rm X}$, and mass, $M_{\rm gas,2500}$, of the intracluster medium are determined using X-ray data collected with Chandra, and $M_{\rm tot}$ is derived from $M_{\rm gas}$ assuming a constant gas mass fraction. Our analysis accounts for several potential sources of bias, including: selection effects, contamination from radio point sources, and the loss of SZE signal due to noise filtering and beam-smoothing effects. We measure the $Y_{2500}$--$Y_{\rm X}$ scaling to have a power-law index of $0.84\pm0.07$, and a fractional intrinsic scatter in $Y_{2500}$ of $(21\pm7)\%$ at fixed $Y_{\rm X}$, both of which are consistent with previous analyses. We also measure the scaling between $Y_{2500}$ and $M_{2500}$, finding a power-law index of $1.06\pm0.12$ and a fractional intrinsic scatter in $Y_{2500}$ at fixed mass of $(25\pm9)\%$. While recent SZE scaling relations using X-ray mass proxies have found power-law indices consistent with the self-similar prediction of 5/3, our measurement stands apart by differing from the self-similar prediction by approximately 5$\sigma$. Given the good agreement between the measured $Y_{2500}$--$Y_{\rm X}$ scalings, much of this discrepancy appears to be caused by differences in the calibration of the X-ray mass proxies adopted for each particular analysis.

The distinctive cometary X-ray morphology of the recently discovered massive galaxy cluster "El Gordo" (ACT-CT J0102-4915; z=0.87) indicates that an unusually high-speed collision is ongoing between two massive galaxy clusters. A bright X-ray "bullet" leads a "twin-tailed" wake, with the SZ centroid at the end of the Northern tail. We show how the physical properties of this system can be determined using our FLASH-based, N-body/hydrodynamic model, constrained by detailed X-ray, Sunyaev-Zel'dovich (SZ), and Hubble lensing and dynamical data. The X-ray morphology and the location of the two Dark Matter components and the SZ peak are accurately described by a simple binary collision viewed about 480 million years after the first core passage. We derive an impact parameter of ~300 kpc, and a relative initial infall velocity of ~2250 km/sec when separated by the sum of the two virial radii assuming an initial total mass of 2.15x10^(15) Msun and a mass ratio of 1.9. Our model demonstrates that tidally stretched gas accounts for the Northern X-ray tail along the collision axis between the mass peaks, and that the Southern tail lies off axis, comprising compressed and shock heated gas generated as the massive component plunges through the main cluster. The challenge for LCDM will be to find out if this physically extreme event can be plausibly accommodated when combined with the similarly massive, high infall velocity case of the "Bullet cluster" and other such cases being uncovered in the new SZ based surveys.

We present a free form mass reconstruction of the massive lensing cluster MACSJ0416.1-2403 using the latest Hubble Frontier Fields data. Our model independent method finds that the extended lensing pattern is generated by two elongated, closely projected clusters of similar mass. Our lens model identifies new lensed images with which we improve the accuracy of the dark matter distribution. We find that the bimodal mass distribution is nearly coincident with the bimodal X-ray emission, but with the two dark matter peaks lying closer together than the centroids of the X-ray emisison. We show this can be achieved if the collision has occurred close to the plane and such that the cores are deflected around each other. The projected mass profiles of both clusters are well constrained because of the many interior lensed images, leading to surprisingly flat mass profiles of both components in the region 15-100 kpc. We discuss the extent to which this may be generated by tidal forces in our dynamical model which are large during an encounter of this type as the cores "graze" each other. The relative velocity between the two cores is estimated to be about 1200 km/s and mostly along the line of sight so that our model is consistent with the relative redshift difference between the two cD galaxies (dz = 0.04).

We report our analysis of MACS J0717.5+3745 using 140 and 268 GHz Bolocam data collected at the Caltech Submillimeter Observatory. We detect extended Sunyaev-Zel'dovich (SZ) effect signal at high significance in both Bolocam bands, and we employ Herschel-SPIRE observations to subtract the signal from dusty background galaxies in the 268 GHz data. We constrain the two-band SZ surface brightness toward two of the sub-clusters of MACS J0717.5+3745: the main sub-cluster (named C), and a sub-cluster identified in spectroscopic optical data to have a line-of-sight velocity of +3200 km/s (named B). We determine the surface brightness in two separate ways: via fits of parametric models and via direct integration of the images. For both sub-clusters, we find consistent surface brightnesses from both analysis methods. We constrain spectral templates consisting of relativistically corrected thermal and kinetic SZ signals, using a jointly-derived electron temperature from Chandra and XMM-Newton under the assumption that each sub-cluster is isothermal. The data show no evidence for a kinetic SZ signal toward sub-cluster C, but they do indicate a significant kinetic SZ signal toward sub-cluster B. The model-derived surface brightnesses for sub-cluster B yield a best-fit line-of-sight velocity of v_z = +3450 +- 900 km/s, with (1 - Prob[v_z > 0]) = 1.3 x 10^-5 (4.2 sigma away from 0 for a Gaussian distribution). The directly integrated sub-cluster B SZ surface brightnesses provide a best-fit v_z = +2550 +- 1050 km/s, with (1 - Prob[v_z > 0]) = 2.2 x 10^-3 (2.9 sigma).

Since the discovery of the "bullet cluster" several similar cases have been uncovered suggesting relative velocities well beyond the tail of high speed collisions predicted by the concordance LCDM model. However, quantifying such post-merger events with hydrodynamical models requires a wide coverage of possible initial conditions. Here we show that it is simpler to interpret pre-merger cases, such as A1750, where the gas between the colliding clusters is modestly affected, so that the initial conditions are clear. We analyze publicly available Chandra data confirming a significant increase in the projected X-ray temperature between the two cluster centers in A1750 consistent with our expectations for a merging cluster. We model this system with a self-consistent hydrodynamical simulation of dark matter and gas using the FLASH code. Our simulations reproduce well the X-ray data, and the measured redshift difference between the two clusters in the phase before the first core passage viewed at an intermediate projection angle. The deprojected initial relative velocity derived using our model is 1460 km/sec which is considerably higher than the predicted mean impact velocity for simulated massive haloes derived by recent LCDM cosmological simulations, but it is within the allowed range. Our simulations demonstrate that such systems can be identified using a multi-wavelength approach and numerical simulations, for which the statistical distribution of relative impact velocities may provide a definitive examination of a broad range of dark matter scenarios.

We show that the fast moving component of the "bullet cluster" (1E0657-56) can induce potentially resolvable redshift differences between multiply-lensed images of background galaxies. The moving cluster effect can be expressed as the scalar product of the lensing deflection angle with the tangential velocity of the mass components, and it is maximal for clusters colliding in the plane of the sky with velocities boosted by their mutual gravity. The bullet cluster is likely to be the best candidate for the first measurement of this effect due to the large collision velocity and because the lensing deflection and the cluster fields can be calculated in advance. We derive the deflection field using multiply-lensed background galaxies detected with the Hubble Space Telescope. The velocity field is modeled using self-consistent N-body/hydrodynamical simulations constrained by the observed X-ray and gravitational lensing features of this system. We predict that the triply-lensed images of systems "G" and "H" straddling the critical curve of the bullet component will show the largest frequency shifts up to ~0.5 km/sec. This is within the range of the Atacama Large Millimeter/sub-millimeter Array (ALMA) for molecular emission, and is near the resolution limit of the new generation high-throughput optical-IR spectrographs. A detection of this effect measures the tangential motion of the subclusters directly, thereby clarifying the tension with LCDM, which is inferred from gas motion less directly. This method may be extended to smaller redshift differences using the Ly-alpha forest towards QSOs lensed by more typical clusters of galaxies. More generally, the tangential component of the peculiar velocities of clusters derived by our method complements the radial component determined by the kinematic SZ effect, providing a full 3-dimensional description of velocities.

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 present a new way to solve the platform deformation problem of co-planar interferometers. The platform of a co-planar interferometer can be deformed due to driving forces and gravity. A deformed platform will induce extra components into the geometric delay of each baseline, and change the phases of observed visibilities. The reconstructed images will also be diluted due to the errors of the phases. The platform deformations of The Yuan-Tseh Lee Array for Microwave Background Anisotropy (AMiBA) were modelled based on photogrammetry data with about 20 mount pointing positions. We then used the differential optical pointing error between two optical telescopes to fit the model parameters in the entire horizontal coordinate space. With the platform deformation model, we can predict the errors of the geometric phase delays due to platform deformation with given azimuth and elevation of the targets and calibrators. After correcting the phases of the radio point sources in the AMiBA interferometric data, we recover 50% - 70% flux loss due to phase errors. This allows us to restore more than 90% of a source flux. The method outlined in this work is not only applicable to the correction of deformation for other co-planar telescopes but also to single dish telescopes with deformation problems. This work also forms the basis of the upcoming science results of AMiBA-13.

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.

We test the assumption of strict hydrostatic equilibrium in galaxy cluster MS2137.3-2353 (MS 2137) using the latest CHANDRA X-ray observations and results from a combined strong and weak lensing analysis based on optical observations. We deproject the two-dimensional X-ray surface brightness and mass surface density maps assuming spherical and spheroidal dark matter distributions. We find a significant, 40%-50%, contribution from non-thermal pressure in the core assuming a spherical model. This non-thermal pressure support is similar to what was found by Molnar et al. (2010) using a sample of massive relaxed clusters drawn from high resolution cosmological simulations. We have studied hydrostatic equilibrium in MS 2137 under the assumption of elliptical cluster geometry adopting prolate models for the dark matter density distribution with different axis ratios. Our results suggest that the main effect of ellipticity (compared to spherical models) is to decrease the non-thermal pressure support required for equilibrium at all radii without changing the distribution qualitatively. We find that a prolate model with an axis ratio of 1.25 (axis in the line of sight over perpendicular to it) provides a physically acceptable model implying that MS 2137 is close to hydrostatic equilibrium at about 0.04-0.15 Rvir and have an about 25% contribution from non-thermal pressure at the center. Our results provide further evidence that there is a significant contribution from non-thermal pressure in the core region of even relaxed clusters, i.e., the assumption of hydrostatic equilibrium is not valid in this region, independently of the assumed shape of the cluster.

We derive an accurate mass distribution of the galaxy cluster MACS J1206.2-0847 (z=0.439) from a combined weak-lensing distortion, magnification, and strong-lensing analysis of wide-field Subaru BVRIz' imaging and our recent 16-band Hubble Space Telescope observations taken as part of the Cluster Lensing And Supernova survey with Hubble (CLASH) program. We find good agreement in the regions of overlap between several weak and strong lensing mass reconstructions using a wide variety of modeling methods, ensuring consistency. The Subaru data reveal the presence of a surrounding large scale structure with the major axis running approximately north-west south-east (NW-SE), aligned with the cluster and its brightest galaxy shapes, showing elongation with a \sim 2:1 axis ratio in the plane of the sky. Our full-lensing mass profile exhibits a shallow profile slope dln\Sigma/dlnR\sim -1 at cluster outskirts (R>1Mpc/h), whereas the mass distribution excluding the NW-SE excess regions steepens further out, well described by the Navarro-Frenk-White form. Assuming a spherical halo, we obtain a virial mass M_{vir}=(1.1\pm 0.2\pm 0.1)\times 10^{15} M_{sun}/h and a halo concentration c_{vir} = 6.9\pm 1.0\pm 1.2 (\sim 5.7 when the central 50kpc/h is excluded), which falls in the range 4< <c> <7 of average c(M,z) predictions for relaxed clusters from recent Lambda cold dark matter simulations. Our full lensing results are found to be in agreement with X-ray mass measurements where the data overlap, and when combined with Chandra gas mass measurements, yield a cumulative gas mass fraction of 13.7^{+4.5}_{-3.0}% at 0.7Mpc/h (\approx 1.7r_{2500}), a typical value observed for high mass clusters.

We present Bolocam observations of two galaxy cluster candidates reported as unconfirmed in the Planck early Sunyaev-Zel'dovich (eSZ) sample, PLCKESZ G115.71+17.52 and PLCKESZ G189.84-37.24. We observed each of these candidates with Bolocam at 140 GHz from the Caltech Submm Observatory in October 2011. The resulting images have white noise levels of ~30 {\mu}KCMB-arcmin in their central regions. We find a significant SZ decrement towards PLCKESZ G115.71. This decrement has a false detection probability of 5.3\times10-5, and we therefore confirm PLCKESZ G115.71 as a cluster. The maximum SZ decrement towards PLCKESZ G189.84 corresponds to a false detection probability of 0.027, and it therefore remains as an unconfirmed cluster candidate. In order to make our SZ-derived results more robust, we have also analyzed data from the Wide-field Infrared Survey Explorer (WISE) at the location of each cluster candidate. We find an overdensity of WISE sources consistent with other clusters in the eSZ at the location of PLCKESZ G115.71, providing further evidence that it is a cluster. We do not find a significant overdensity of WISE sources at the location of PLCKESZ G189.84.

Galaxy clusters, the most massive collapsed structures, have been routinely used to determine cosmological parameters. When using clusters for cosmology, the crucial assumption is that they are relaxed. However, subarcminute resolution Sunyaev-Zel'dovich (SZ) effect images compared with high resolution X-ray images of some clusters show significant offsets between the two peaks. We have carried out self-consistent N-body/hydrodynamical simulations of merging galaxy clusters using FLASH to study these offsets quantitatively. We have found that significant displacements result between the SZ and X-ray peaks for large relative velocities for all masses used in our simulations as long as the impact parameters were about 100-250 kpc. Our results suggest that the SZ peak coincides with the peak in the pressure times the line-of-sight characteristic length and not the pressure maximum (as it would for clusters in equilibrium). The peak in the X-ray emission, as expected, coincides with the density maximum of the main cluster. As a consequence, the morphology of the SZ signal and therefore the offset between the SZ and X-ray peaks change with viewing angle. As an application, we compare the morphologies of our simulated images to observed SZ and X-ray images and mass surface densities derived from weak lensing observations of the merging galaxy cluster CL0152-1357. We find that a large relative velocity of 4800 km/s is necessary to explain these observations. We conclude that an analysis of the morphologies of multi-frequency observations of merging clusters can be used to put meaningful constraints on the initial parameters of the progenitors.

We investigate the contamination of the Sunyaev--Zel'dovich (SZ) effect for six galaxy clusters, A1689, A1995, A2142, A2163, A2261, and A2390, observed by the Y. T. Lee Array for Microwave Background Anisotropy during 2007. With the range of baselines used, we find that the largest effect (of order 13%-50% of the central SZ flux density) comes from primary anisotropies in the cosmic microwave background and exceeds the thermal noise in all six cases. Contamination from discrete radio sources is estimated to be at a level of (3%-60%) of the central SZ flux density. We use the statistics of these contaminating sources to estimate and correct the errors in the measured SZ effects of these clusters.

Accurate mass determination of clusters of galaxies is crucial if they are to be used as cosmological probes. However, there are some discrepancies between cluster masses determined based on gravitational lensing, and X-ray observations assuming strict hydrostatic equilibirium (i.e., the equilibrium gas pressure is provided entirely by thermal pressure). Cosmological simulations suggest that turbulent gas motions remaining from hierarchical structure formation may provide a significant contribution to the equilibrium pressure in clusters. We analyze a sample of massive clusters of galaxies drawn from high resolution cosmological simulations, and find a significant contribution (20%-45%) from non-thermal pressure near the center of relaxed clusters, and, in accord with previous studies, a minimum contribution at about 0.1 Rvir, growing to about 30%-45% at the virial radius, Rvir. Our results strongly suggest that relaxed clusters should have significant non-thermal support in their core region. As an example, we test the validity of strict hydrostatic equilibirium in the well-studied massive galaxy cluster Abell 1689 using the latest high resolution gravitational lensing and X-ray observations. We find a contribution of about 40% from non-thermal pressure within the core region of A1689, suggesting an alternate explanation for the mass discrepancy: the strict hydrostatic equilibirium is not valid in this region.

Clusters of galaxies have been used extensively to determine cosmological parameters. A major difficulty in making best use of Sunyaev-Zel'dovich (SZ) and X-ray observations of clusters for cosmology is that using X-ray observations it is difficult to measure the temperature distribution and therefore determine the density distribution in individual clusters of galaxies out to the virial radius. Observations with the new generation of SZ instruments are a promising alternative approach. We use clusters of galaxies drawn from high-resolution adaptive mesh refinement (AMR) cosmological simulations to study how well we should be able to constrain the large-scale distribution of the intra-cluster gas (ICG) in individual massive relaxed clusters using AMiBA in its configuration with 13 1.2-m diameter dishes (AMiBA13) along with X-ray observations. We show that non-isothermal beta models provide a good description of the ICG in our simulated relaxed clusters. We use simulated X-ray observations to estimate the quality of constraints on the distribution of gas density, and simulated SZ visibilities (AMiBA13 observations) for constraints on the large-scale temperature distribution of the ICG. We find that AMiBA13 visibilities should constrain the scale radius of the temperature distribution to about 50% accuracy. We conclude that the upgraded AMiBA, AMiBA13, should be a powerful instrument to constrain the large-scale distribution of the ICG.

The most X-ray luminous cluster known, RXJ1347-1145 (z=0.45), has been the object of extensive study across the electromagnetic spectrum. We have imaged the Sunyaev-Zel'dovich Effect (SZE) at 90 GHz (3.3 mm) in RXJ1347-1145 at 10" resolution with the 64-pixel MUSTANG bolometer array on the Green Bank Telescope (GBT), confirming a previously reported strong, localized enhancement of the SZE 20" to the South-East of the center of X-ray emission. This enhancement of the SZE has been interpreted as shock-heated (> 20 keV) gas caused by an ongoing major (low mass-ratio) merger event. Our data support this interpretation. We also detect a pronounced asymmetry in the projected cluster pressure profile, with the pressure just east of the cluster core ~1.6 times higher than just to the west. This is the highest resolution image of the SZE made to date.

We derive an accurate mass distribution of the rich galaxy cluster Cl0024+1654 (z=0.395) based on deep Subaru BR_{c}z' imaging and our recent comprehensive strong lensing analysis of HST/ACS/NIC3 observations. We obtain the weak lensing distortion and magnification of undilted samples of red and blue background galaxies by carefully combining all color and positional information. Unlike previous work, the weak and strong lensing are in excellent agreement where the data overlap. The joint mass profile continuously steepens out to the virial radius with only a minor contribution \sim 10% in the mass from known subcluster at a projected distance of \sim 700kpc/h. The projected mass distribution for the entire cluster is well fitted with a single Navarro-Frenk-White model with a virial mass, M_{vir} = (1.2 \pm 0.2) \times 10^{15} M_{sun}/h, and a concentration, c_{vir} = 9.2^{+1.4}_{-1.2}. This model fit is fully consistent with the depletion of the red background counts, providing independent confirmation. Careful examination and interpretation of X-ray and dynamical data strongly suggest that this cluster system is in a post collision state, which we show is consistent with our well-defined mass profile for a major merger occurring along the line of sight, viewed approximately 2-3Gyr after impact when the gravitational potential has had time to relax in the center, before the gas has recovered and before the outskirts are fully virialized. Finally, our full lensing analysis provides a model-independent constraint of M_{2D}(<r_{vir}) = (1.4 \pm 0.3) \times 10^{15} M_{sun}/h for the projected mass of the whole system, including any currently unbound material beyond the virial radius, which can constrain the sum of the two pre-merger cluster masses when designing simulations to explore this system.

The Sunyaev-Zel'dovich Effect (SZE) has been observed toward six massive galaxy clusters, at redshifts 0.091 \leq z \leq 0.322 in the 86-102 GHz band with the Y. T. Lee Array for Microwave Background Anisotropy (AMiBA). We modify an iterative method, based on the isothermal \beta-models, to derive the electron temperature T_e, total mass M_t, gas mass M_g, and integrated Compton Y within r_2500, from the AMiBA SZE data. Non-isothermal universal temperature profile (UTP) \beta models are also considered in this paper. These results are in good agreement with those deduced from other observations. We also investigate the embedded scaling relations, due to the assumptions that have been made in the method we adopted, between these purely SZE-deduced T_e, M_t, M_g and Y. Our results suggest that cluster properties may be measurable with SZE observations alone. However, the assumptions built into the pure-SZE method bias the results of scaling relation estimations and need further study.

We investigate the scaling relations between the X-ray and the thermal Sunyaev-Zel'dovich Effect (SZE) properties of clusters of galaxies, using data taken during 2007 by the Y.T. Lee Array for Microwave Background Anisotropy (AMiBA) at 94 GHz for the six clusters A1689, A1995, A2142, A2163, A2261, and A2390. The scaling relations relate the integrated Compton-y parameter Y_{2500} to the X-ray derived gas temperature T_{e}, total mass M_{2500}, and bolometric luminosity L_X within r_{2500}. Our results for the power-law index and normalization are both consistent with the self-similar model and other studies in the literature except for the Y_{2500}-L_X relation, for which a physical explanation is given though further investigation may be still needed. Our results not only provide confidence for the AMiBA project but also support our understanding of galaxy clusters.

The Array for Microwave Background Anisotropy (AMiBA) is a radio interferometer for research in cosmology, currently operating 7 0.6m diameter antennas co-mounted on a 6m diameter platform driven by a hexapod mount. AMiBA is currently the largest hexapod telescope. We briefly summarize the hexapod operation with the current pointing error model. We then focus on the upcoming 13-element expansion with its potential difficulties and solutions. Photogrammetry measurements of the platform reveal deformations at a level which can affect the optical pointing and the receiver radio phase. In order to prepare for the 13-element upgrade, two optical telescopes are installed on the platform to correlate optical pointing tests. Being mounted on different locations, the residuals of the two sets of pointing errors show a characteristic phase and amplitude difference as a function of the platform deformation pattern. These results depend on the telescope's azimuth, elevation and polarization position. An analytical model for the deformation is derived in order to separate the local deformation induced error from the real hexapod pointing error. Similarly, we demonstrate that the deformation induced radio phase error can be reliably modeled and calibrated, which allows us to recover the ideal synthesized beam in amplitude and shape of up to 90% or more. The resulting array efficiency and its limits are discussed based on the derived errors.

The Y. T. Lee Array for Microwave Background (AMiBA) has reported the first science results on the detection of galaxy clusters via the Sunyaev Zel'dovich effect. The science objectives required small reflectors in order to sample large scale structures (20') while interferometry provided modest resolutions (2'). With these constraints, we designed for the best sensitivity by utilizing the maximum possible continuum bandwidth matched to the atmospheric window at 86-102GHz, with dual polarizations. A novel wide-band analog correlator was designed that is easily expandable for more interferometer elements. MMIC technology was used throughout as much as possible in order to miniaturize the components and to enhance mass production. These designs will find application in other upcoming astronomy projects. AMiBA is now in operations since 2006, and we are in the process to expand the array from 7 to 13 elements.

Cold dark matter (CDM) hierarchical structure formation models predict the existence of large-scale accretion shocks between the virial and turnaround radii of clusters of galaxies. Kocsis et al. (2005) suggest that the Sunyaev-Zel'dovich (SZ) signal associated with such shocks might be observable with the next generation radio interferometer, ALMA. We study the three--dimensional distribution of accretion shocks around individual clusters of galaxies drawn from adaptive mesh refinement (AMR) and smoothed particle hydrodynamics (SPH) simulations of LCDM (dark energy dominated CDM) models. In relaxed clusters, we find two distinct sets of shocks. One set ("virial shocks"), with Mach numbers of 2.5-4, is located at radii 0.9-1.3 Rvir, where Rvir is the spherical infall estimate of the virial radius, covering about 40-50% of the total surface area around clusters at these radii. Another set of stronger shocks ("external shocks") is located farther out, at about 3 Rvir, with large Mach numbers (~100), covering about 40-60% of the surface area. We simulate SZ surface brightness maps of relaxed massive galaxy clusters drawn from high resolution AMR runs, and conclude that ALMA should be capable of detecting the virial shocks in massive clusters of galaxies. More simulations are needed to improve estimates of astrophysical noise and to determine optimal observational strategies.

The Y.T. Lee Array for Microwave Background Anisotropy (AMiBA) started scientific operation in early 2007. This work describes the optimization of the system performance for the measurements of the Sunyaev-Zel'dovich effect for six massive galaxy clusters at redshifts $0.09 - 0.32$. We achieved a point source sensitivity of $63\pm 7$ mJy with the seven 0.6m dishes in 1 hour of on-source integration in 2-patch differencing observations. We measured and compensated for the delays between the antennas of our platform-mounted interferometer. Beam switching was used to cancel instrumental instabilities and ground pick up. Total power and phase stability were good on time scales of hours, and the system was shown to integrate down on equivalent timescales of 300 hours per baseline/correlation, or about 10 hours for the entire array. While the broadband correlator leads to good sensitivity, the small number of lags in the correlator resulted in poorly measured bandpass response. We corrected for this by using external calibrators (Jupiter and Saturn). Using Jupiter as the flux standard, we measured the disk brightness temperature of Saturn to be $149^{+5}_{-12}$ K.

AMiBA is the largest hexapod astronomical telescope in current operation. We present a description of this novel hexapod mount with its main mechanical components -- the support cone, universal joints, jack screws, and platform -- and outline the control system with the pointing model and the operating modes that are supported. The AMiBA hexapod mount performance is verified based on optical pointing tests and platform photogrammetry measurements. The photogrammetry results show that the deformations in the inner part of the platform are less than 120 micron rms. This is negligible for optical pointing corrections, radio alignment and radio phase errors for the currently operational 7-element compact configuration. The optical pointing error in azimuth and elevation is successively reduced by a series of corrections to about 0.4 arcmin rms which meets our goal for the 7-element target specifications.

We present observations, analysis and results for the first-year operation of AMiBA, an interferometric experiment designed to study cosmology via the measurement of Cosmic Microwave Background (CMB). AMiBA is the first CMB interferometer operating at 3 mm to have reported successful results, currently with seven close-packed antennas of 60-cm diameter giving a synthesized resolution of around 6'. During 2007 AMiBA detected the Sunyaev-Zel'dovich effects (SZE) of six galaxy clusters at redshift 0.091 <= z <= 0.322. An observing strategy with on-off-source switching is used to minimize the effects from electronic offset and ground pickup. Planets were used to test the observational capability of AMiBA and to calibrate the conversion from correlator time-lag data to visibilities. The detailed formalism for data analysis is given. We summarize our early tests including observations of planets and quasars, and present images, visibility profiles, the estimated central coordinates, sizes, and SZE amplitudes of the galaxy clusters. Science results are summarized. We also discuss possible systematic effects in the results.

The Yuan-Tseh Lee Array for Microwave Background Anisotropy (AMiBA) is the first interferometer dedicated to studying the cosmic microwave background (CMB) radiation at 3mm wavelength. The choice of 3mm was made to minimize the contributions from foreground synchrotron radiation and Galactic dust emission. The initial configuration of seven 0.6m telescopes mounted on a 6-m hexapod platform was dedicated in October 2006 on Mauna Loa, Hawaii. Scientific operations began with the detection of a number of clusters of galaxies via the thermal Sunyaev-Zel'dovich effect. We compare our data with Subaru weak lensing data in order to study the structure of dark matter. We also compare our data with X-ray data in order to derive the Hubble constant.

We present a multiwavelength analysis of a sample of four hot (T_X>8keV) X-ray galaxy clusters (A1689, A2261, A2142, and A2390) using joint AMiBA Sunyaev-Zel'dovich effect (SZE) and Subaru weak lensing observations, combined with published X-ray temperatures, to examine the distribution of mass and the intracluster medium (ICM) in massive cluster environments. Our observations show that A2261 is very similar to A1689 in terms of lensing properties. Many tangential arcs are visible around A2261, with an effective Einstein radius \sim 40 arcsec (at z \sim 1.5), which when combined with our weak lensing measurements implies a mass profile well fitted by an NFW model with a high concentration c_{vir} \sim 10, similar to A1689 and to other massive clusters. The cluster A2142 shows complex mass substructure, and displays a shallower profile (c_{vir} \sim 5), consistent with detailed X-ray observations which imply recent interaction. The AMiBA map of A2142 exhibits an SZE feature associated with mass substructure lying ahead of the sharp north-west edge of the X-ray core suggesting a pressure increase in the ICM. For A2390 we obtain highly elliptical mass and ICM distributions at all radii, consistent with other X-ray and strong lensing work. Our cluster gas fraction measurements, free from the hydrostatic equilibrium assumption, are overall in good agreement with published X-ray and SZE observations, with the sample-averaged gas fraction of <f_{gas}> = 0.133 \pm 0.027, for our sample <M_{vir}> = (1.2 \pm 0.1) \times 10^{15} M_{sun} h^{-1}. When compared to the cosmic baryon fraction f_b = \Omega_b/\Omega_m constrained by the WMAP 5-year data, this indicates <f_{gas}>/f_b = 0.78 \pm 0.16, i.e., (22 \pm 16)% of the baryons are missing from the hot phase of clusters.

We describe methods used to validate data from the Y.T. Lee Array for Microwave Background Anisotropy (AMiBA), an interferometric array designed to measure the Sunyaev-Zel'dovich effect and the anisotropy of the Cosmic Microwave Background (CMB). We perform several statistical tests on data from pointed galaxy cluster observations taken in 2007 and noise data from long-term blank sky observations and measurements with the feeds covered by the absorbers. We apply power spectrum analysis, cross power spectrum analysis among different outputs with different time lags in our analog correlator, and sample variance law tests to noise data. We find that (1) there is no time variation of electronic offsets on the time scale of our two-patch observations (~10 minutes); (2) noise is correlated by less than 10% between different lags; and (3) the variance of noise scales with the inverse of time. To test the Gaussianity of the data, we apply Kolmogorov-Smirnov (K-S) tests to cluster data, and find that a 5% significance level efficiently detects data sets with known hardware problems without rejecting an excess of acceptable data. We also calculate third- and fourth-order moments and cumulants for the noise residual visibilities and find that about 95% of our data are within the 99% confidence regions of Gaussianity.

Bright clusters of galaxies can be seen out to cosmological distances, and thus they can be used to derive cosmological parameters. Although the continuum X-ray emission from the intra-cluster gas is optically thin, the optical depth of resonant lines of ions of heavy elements can be larger than unity. In this Letter we study the feasibility of deriving distances to clusters of galaxies by determining the spatial distribution of the intra-cluster gas from X-ray imaging and the optical depth from resonant emission lines (the XREL method). We solve the radiative transfer problem for line scattering in the hot intra-cluster gas using Monte Carlo simulations. We discuss the spatial and spectral resolutions needed to use the XREL method for accurate determination of distances, and hence cosmological parameters, and show that accurate distances will be obtained by applying this technique with the next generation of high resolution X-ray spectrometers.

We report on preliminary results of our X-ray survey of the most massive clusters currently identified from the Deep Lens Survey (DLS). The DLS cluster sample is selected based on weak lensing shear, which makes it possible for the first time to study clusters in a baryon-independent way. In this article we present X-ray properties of a subset of the shear-selected cluster sample.

We demonstrate that constraints on cosmological parameters from the distribution of clusters as a function of redshift (dN/dz) are complementary to accurate angular diameter distance (D_A) measurements to clusters, and their combination significantly tightens constraints on the energy density content of the Universe. The number counts can be obtained from X-ray and/or SZ (Sunyaev-Zel'dovich effect) surveys, and the angular diameter distances can be determined from deep observations of the intra-cluster gas using their thermal bremsstrahlung X-ray emission and the SZ effect. We combine constraints from simulated cluster number counts expected from a 12 deg^2 SZ cluster survey and constraints from simulated angular diameter distance measurements based on the X-ray/SZ method assuming a statistical accuracy of 10% in the angular diameter distance determination of 100 clusters with redshifts less than 1.5. We find that Omega_m can be determined within about 25%, Omega_Lambda within 20%, and w within 16%. We show that combined dN/dz + D_A constraints can be used to constrain the different energy densities in the Universe even in the presence of a few percent redshift dependent systematic error in D_A. We also address the question of how best to select clusters of galaxies for accurate diameter distance determinations. We show that the joint dN/dz + D_A constraints on cosmological parameters for a fixed target accuracy in the energy density parameters are optimized by selecting clusters with redshift upper cut--offs in the range 0.5 < z < 1.

We propose two new methods for measuring tangential peculiar velocities of rich clusters of galaxies. Our first method is based on weak gravitational lensing and takes advantage of the differing images of background galaxies caused by moving and stationary gravitational potentials. Our second method is based on measuring relative frequency shifts between multiple images of a single strongly lensed background galaxy. We illustrate this method using the example of galaxy cluster CL 0024+1654.

We have searched for unresolved X-ray sources in the vicinity of two rich clusters of galaxies: Abell 1995 (A1995) and MS 0451.6-0305 (MS0451), using the Chandra X-ray observatory. We detected significantly more unresolved sources around A1995 than expected based on the number of X-ray sources to the same flux limit detected in deep \CHANDRA observations of blank fields. Previous studies have also found excess X-ray sources in the vicinity of several nearby clusters of galaxies using \ROSAT$ $, and recently in more distant ($z ~ 0.5$) clusters (RXJ0030 and 3C295) using Chandra. In contrast, we detect only 14 unresolved X-ray sources near MS0451, which is consistent with the number expected from a cluster-free background. We determine the luminosity functions of the extra sources under the assumption that they are at the distance of their respective clusters. The characteristic luminosity of the extra sources around A1995 must be an order of magnitude fainter than that of the extra sources around RXJ0030 and 3C295. The apparent lack of extra sources around MS0451 is consistent with its greater distance and the same characteristic luminosity as the A1995 sources. Hardness ratios suggest that, on average, the extra sources in A1995 may have harder spectra than those of RXJ0030 and 3C295. These results indicate that different classes of objects may dominate in different clusters, perhaps depending on the formation history and/or dynamical state of the accompanying cluster.

We discuss how the space of possible cosmological parameters is constrained by the angular diameter distance function, D_A(z), as measured using the SZ/X-ray method which combines Sunyaev-Zel'dovich (SZ) effect and X-ray brightness data for clusters of galaxies. New X-ray satellites, and ground-based interferometers dedicated to SZ observations, should soon lead to D_A(z) measurements limited by systematic rather than random error. We analyze the systematic and random error budgets to make a realistic estimate of the accuracy achievable in the determination of (Omega_m,Lambda,h), the density parameters of matter and cosmological constant, and the dimensionless Hubble constant, using D_A(z) derived from the SZ/X-ray method, and the position of the first ``Doppler'' peak in the cosmic microwave background fluctuations. We briefly study the effect of systematic errors. We find that Omega_m, Lambda, and w are affected, but h is not by systematic errors which grow with redshift. With as few as 70 clusters, each providing a measurement of D_A(z) with a 7% random and 5% systematic error, Omega_m can be constrained to +/-0.2, Lambda to +/-0.2, and h to +/-0.11 (all at 3 sigma). We also estimate constraints for the alternative three-parameter set (Omega_m,w,h), where w is the equation of state parameter. The measurement of D_A(z) provides constraints complementary to those from the number density of clusters in redshift space. A sample of 70 clusters (D_A measured with the same accuracy as before) combined with cluster evolution results (or a known matter density), can constrain w within +/-0.45 (at 3 sigma). Studies of X-ray and SZ properties of clusters of galaxies promise an independent and powerful test for cosmological parameters.

The formation and evolution of clusters of galaxies are sensitive to the underlying cosmological model. Constraints on cosmological parameters of cold dark matter models have been derived from mass, temperature and luminosity functions of clusters. We study the importance of including cluster formation history and a correction for collapsed fraction of objects in determining the cluster X-ray temperature function. We find that both effects are important. We compare temperature functions obtained by using a power law approximation for the mass variance normalized to X-ray clusters to those obtained by using a COBE normalized full CDM treatment. We conclude that the temperature function could be a powerful test on the average density of the Universe if we could find the correct way of interpreting the data.

We estimate the contributions to the cosmic microwave background radiation (CMBR) power spectrum from the static and kinematic Sunyaev-Zel'dovich (SZ) effects, and from the moving cluster of galaxies (MCG) effect. We conclude, in agreement with other studies, that at sufficiently small scales secondary fluctuations caused by clusters provide important contributions to the CMBR. At $\ell \gtrsim 3000$, these secondary fluctuations become important relative to lensed primordial fluctuations. Gravitational lensing at small angular scales has been proposed as a way to break the ``geometric degeneracy'' in determining fundamental cosmological parameters. We show that this method requires the separation of the static SZ effect, but the kinematic SZ effect and the MCG effect are less important. The power spectrum of secondary fluctuations caused by clusters of galaxies, if separated from the spectrum of lensed primordial fluctuations, might provide an independent constraint on several important cosmological parameters.

We investigated the effect of inverse Compton scattering in mildly relativistic static and moving plasmas with low optical depth using Monte Carlo simulations, and calculated the Sunyaev-Zel'dovich effect in the cosmic background radiation. Our semi-analytic method is based on a separation of photon diffusion in frequency and real space. We use Monte Carlo simulation to derive the intensity and frequency of the scattered photons for a monochromatic incoming radiation. The outgoing spectrum is determined by integrating over the spectrum of the incoming radiation using the intensity to determine the correct weight. This method makes it possible to study the emerging radiation as a function of frequency and direction. As a first application we have studied the effects of finite optical depth and gas infall on the Sunyaev-Zel'dovich effect (not possible with the extended Kompaneets equation) and discuss the parameter range in which the Boltzmann equation and its expansions can be used. For high temperature clusters ($k_B T_e \gtrsim 15$ keV) relativistic corrections based on a fifth order expansion of the extended Kompaneets equation seriously underestimate the Sunyaev-Zel'dovich effect at high frequencies. The contribution from plasma infall is less important for reasonable velocities. We give a convenient analytical expression for the dependence of the cross-over frequency on temperature, optical depth, and gas infall speed. Optical depth effects are often more important than relativistic corrections, and should be taken into account for high-precision work, but are smaller than the typical kinematic effect from cluster radial velocities.