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

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 present the design, simulation, and planned fabrication process of a flat high resistivity silicon gradient index (GRIN) lens for millimeter and submillimeter wavelengths with very low absorption losses. The gradient index is created by subwavelength holes whose size increases with the radius of the lens. The effective refractive index created by the subwavelength holes is constant over a very wide bandwidth, allowing the fabrication of achromatic lenses up to submillimeter wavelengths. The designed GRIN lens was successfully simulated and shows an expected efficiency better than that of a classic silicon plano-concave spherical lens with approximately the same thickness and focal length. Deep reactive ion etching (DRIE) and wafer-bonding of several patterned wafers will be used to realize our first GRIN lens prototype.

We present the on-going activity to characterize the geometrical properties of the gas and dark matter haloes using multi-wavelength observations of galaxy clusters. The role of the SZ signal in describing the gas distribution is discussed for the pilot case of the CLASH object MACS J1206.2-0847. Preliminary images of the NIKA2 and ALMA exposures are presented.

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

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

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

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

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

We reconstruct the two-dimensional (2D) matter distributions in 20 high-mass galaxy clusters selected from the CLASH survey by using the new approach of performing a joint weak lensing analysis of 2D shear and azimuthally averaged magnification measurements. This combination allows for a complete analysis of the field, effectively breaking the mass-sheet degeneracy. In a Bayesian framework, we simultaneously constrain the mass profile and morphology of each individual cluster assuming an elliptical Navarro-Frenk-White halo characterized by the mass, concentration, projected axis ratio, and position angle of the projected major axis.. We find that spherical mass estimates of the clusters from azimuthally averaged weak-lensing measurements in previous work are in excellent agreement with our results from a full 2D analysis. Combining all 20 clusters in our sample, we detect the elliptical shape of weak-lensing halos at the $5\sigma$ significance level within a scale of 2Mpc$/h$. The median projected axis ratio is $0.67\pm 0.07$ at a virial mass of $M_\mathrm{vir}=(15.2\pm 2.8) \times 10^{14} M_\odot$, which is in agreement with theoretical predictions of the standard collisionless cold dark matter model. We also study misalignment statistics of the brightest cluster galaxy, X-ray, thermal Sunyaev-Zel'dovich effect, and strong-lensing morphologies with respect to the weak-lensing signal. Among the three baryonic tracers studied here, we find that the X-ray morphology is best aligned with the weak-lensing mass distribution, with a median misalignment angle of $21\pm 7$ degrees. We also conduct a stacked quadrupole shear analysis assuming that the X-ray major axis is aligned with that of the projected mass distribution. This yields a consistent axis ratio of $0.67\pm 0.10$, suggesting again a tight alignment between the intracluster gas and dark matter.

Silicon optics with wide bandwidth anti-reflection (AR) coatings, made of multi-layer textured silicon surfaces, are developed for millimeter and submillimeter wavelengths. Single and double layer AR coatings were designed for an optimal transmission centered on 250 GHz, and fabricated using the DRIE (Deep Reaction Ion Etching) technique. Tests of high resistivity silicon wafers with single-layer coatings between 75 GHz and 330 GHz are presented and compared with the simulations.

Although high-resistivity, low-loss silicon is an excellent material for THz transmission optics, its high refractive index necessitates antireflection treatment. We fabricated a wide-bandwidth, two-layer antireflection treatment by cutting subwavelength structures into the silicon surface using multi-depth deep reactive ion etching (DRIE). A wafer with this treatment on both sides has <-20 dB (<1%) reflectance over 190-310 GHz. We also demonstrated that bonding wafers introduces no reflection features above the -20 dB level, reproducing previous work. Together these developments immediately enable construction of wide-bandwidth silicon vacuum windows and represent two important steps toward gradient-index silicon optics with integral broadband antireflection treatment.

The $\Lambda$CDM model of structure formation makes strong predictions on concentration and shape of DM (dark matter) halos, which are determined by mass accretion processes. Comparison between predicted shapes and observations provides a geometric test of the $\Lambda$CDM model. Accurate and precise measurements needs a full three-dimensional analysis of the cluster mass distribution. We accomplish this with a multi-probe 3D analysis of the X-ray regular CLASH (Cluster Lensing And Supernova survey with Hubble) clusters combining strong and weak lensing, X-ray photometry and spectroscopy, and the Sunyaev-Zel'dovich effect. The cluster shapes and concentrations are consistent with $\Lambda$CDM predictions. The CLASH clusters are randomly oriented, as expected given the sample selection criteria. Shapes agree with numerical results for DM-only halos, which hints at baryonic physics being not so effective in making halos rounder.

We perform a three-dimensional triaxial analysis of 16 X-ray regular and 4 high-magnification galaxy clusters selected from the CLASH survey by combining two-dimensional weak-lensing and central strong-lensing constraints. In a Bayesian framework, we constrain the intrinsic structure and geometry of each individual cluster assuming a triaxial Navarro-Frenk-White halo with arbitrary orientations, characterized by the mass $M_{200\mathrm{c}}$, halo concentration $C_{200\mathrm{c}}$, and triaxial axis ratios ($q_{\mathrm{a}} \le q_{\mathrm{b}}$), and investigate scaling relations between these halo structural parameters. From triaxial modeling of the X-ray-selected subsample, we find that the halo concentration decreases with increasing cluster mass, with a mean concentration of $C_{200\mathrm{c}} = 4.82\pm0.30$ at the pivot mass $M_{200\mathrm{c}}=10^{15}M_{\odot}h^{-1}$. This is consistent with the result from spherical modeling, $C_{200\mathrm{c}}=4.51\pm 0.14$. Independently of the priors, the minor-to-major axis ratio $q_{\mathrm{a}}$ of our full sample exhibits a clear deviation from the spherical configuration ($q_{\mathrm{a}}=0.52 \pm 0.04$ at $10^{15}M_{\odot}h^{-1}$ with uniform priors), with a weak dependence on the cluster mass. Combining all 20 clusters, we obtain a joint ensemble constraint on the minor-to-major axis ratio of $q_{\mathrm{a}}=0.652^{+0.162}_{-0.078}$ and a lower bound on the intermediate-to-major axis ratio of $q_{\mathrm{b}}>0.63$ at the $2\sigma$ level from an analysis with uniform priors. Assuming priors on the axis ratios derived from numerical simulations, we constrain the degree of triaxiality for the full sample to be $\mathcal{T}=0.79 \pm 0.03$ at $10^{15}M_{\odot}h^{-1}$, indicating a preference for a prolate geometry of cluster halos. We find no statistical evidence for an orientation bias ($f_{\mathrm{geo}}=0.93 \pm 0.07$) (abridged)

We have selected a sample of eleven massive clusters of galaxies observed by the Hubble Space Telescope in order to study the impact of the dynamical state on the IntraCluster Light (ICL) fraction, the ratio of total integrated ICL to the total galaxy member light. With the exception of the Bullet cluster, the sample is drawn from the Cluster Lensing and Supernova Survey and the Frontier Fields program, containing five relaxed and six merging clusters. The ICL fraction is calculated in three optical filters using the CHEFs IntraCluster Light Estimator, a robust and accurate algorithm free of a priori assumptions. We find that the ICL fraction in the three bands is, on average, higher for the merging clusters, ranging between $\sim7-23\%$, compared with the $\sim 2-11\%$ found for the relaxed systems. We observe a nearly constant value (within the error bars) in the ICL fraction of the regular clusters at the three wavelengths considered, which would indicate that the colors of the ICL and the cluster galaxies are, on average, coincident and, thus, their stellar populations. However, we find a higher ICL fraction in the F606W filter for the merging clusters, consistent with an excess of lower-metallicity/younger stars in the ICL, which could have migrated violently from the outskirts of the infalling galaxies during the merger event.

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.

The V380 Ori NE bipolar outflow was imaged in the SiO and CO J = 1 - 0 lines, and dense cores in L1641 were observed in the 2.0-0.89 mm continuum. The highly collimated SiO jet shows point-symmetric oscillation patterns in both position and velocity, which suggests that the jet axis is precessing and the driving source may belong to a non-coplanar binary system. By considering the position and velocity variabilities together, accurate jet parameters were derived. The protostellar system is viewed nearly edge-on, and the jet has a flow speed of 35 km/s and a precession period of 1600 years. The CO outflow length gives a dynamical timescale of 6300 years, and the protostar must be extremely young. The inferred binary separation of 6-70 au implies that this protobinary system may have been formed through the disk instability process. The continuum spectra of L1641 dense cores indicate that the emission comes from dust, and the fits with modified blackbody functions give emissivity power indices of beta = 0.3-2.2. The emissivity index shows a positive correlation with the molecular line width, but no strong correlation with bolometric luminosity or temperature. V380 Ori NE has a particularly low value of beta = 0.3, which tentatively suggests the presence of millimeter-sized dust grains. Because the dust growth takes millions of years, much longer than the protostellar age, this core may have produced large grains in the starless core stage. HH 34 MMS and HH 147 MMS also have low emissivity indices.

We present accurate photometric redshifts for galaxies observed by the Cluster Lensing and Supernova survey with Hubble (CLASH). CLASH observed 25 massive galaxy cluster cores with the Hubble Space Telescope in 16 filters spanning 0.2 - 1.7 $\mu$m. Photometry in such crowded fields is challenging. Compared to our previously released catalogs, we make several improvements to the photometry, including smaller apertures, ICL subtraction, PSF matching, and empirically measured uncertainties. We further improve the Bayesian Photometric Redshift (BPZ) estimates by adding a redder elliptical template and by inflating the photometric uncertainties of the brightest galaxies. The resulting photometric redshift accuracies are dz/(1+z) $\sim$ 0.8\%, 1.0\%, and 2.0\% for galaxies with I-band F814W AB magnitudes $<$ 18, 20, and 23, respectively. These results are consistent with our expectations. They improve on our previously reported accuracies by a factor of 4 at the bright end and a factor of 2 at the faint end. Our new catalog includes 1257 spectroscopic redshifts, including 382 confirmed cluster members. We also provide stellar mass estimates. Finally, we include lensing magnification estimates of background galaxies based on our public lens models. Our new catalog of all 25 CLASH clusters is available via MAST. The analysis techniques developed here will be useful in other surveys of crowded fields, including the Frontier Fields and surveys carried out with J-PAS and JWST.

We present pressure profiles of galaxy clusters determined from high resolution Sunyaev-Zel'dovich (SZ) effect observations of fourteen clusters, which span the redshift range $ 0.25 < z < 0.89$. The procedure simultaneously fits spherical cluster models to MUSTANG and Bolocam data. In this analysis, we adopt the generalized NFW parameterization of pressure profiles to produce our models. Our constraints on ensemble-average pressure profile parameters, in this study $\gamma$, $C_{500}$, and $P_0$, are consistent with those in previous studies, but for individual clusters we find discrepancies with the X-ray derived pressure profiles from the ACCEPT2 database. We investigate potential sources of these discrepancies, especially cluster geometry, electron temperature of the intracluster medium, and substructure. We find that the ensemble mean profile for all clusters in our sample is described by the parameters: $[\gamma,C_{500},P_0] = [0.3_{-0.1}^{+0.1}, 1.3_{-0.1}^{+0.1}, 8.6_{-2.4}^{+2.4}]$, for cool core clusters: $[\gamma,C_{500},P_0] = [0.6_{-0.1}^{+0.1}, 0.9_{-0.1}^{+0.1}, 3.6_{-1.5}^{+1.5}]$, and for disturbed clusters: $[\gamma,C_{500},P_0] = [0.0_{-0.0}^{+0.1}, 1.5_{-0.2}^{+0.1},13.8_{-1.6}^{+1.6}]$. Four of the fourteen clusters have clear substructure in our SZ observations, while an additional two clusters exhibit potential substructure.

Multi-wavelength techniques can probe the distribution and the physical properties of baryons and dark matter in galaxy clusters from the inner regions out to the peripheries. We present a full three-dimensional analysis combining strong and weak lensing, X-ray surface brightness and temperature, and the Sunyaev-Zel'dovich effect. The method is applied to MACS J1206.2-0847, a remarkably regular, face-on, massive, M_{200}=(1.1+-0.2)*10^{15}M_Sun/h, cluster at z=0.44. The measured concentration, c_{200}=6.3+-1.2, and the triaxial shape are common to halos formed in a LCDM scenario. The gas has settled in and follows the shape of the gravitational potential, which is evidence of pressure equilibrium via the shape theorem. There is no evidence for significant non-thermal pressure and the equilibrium is hydrostatic.

To investigate the relationship between thermal and non-thermal components in merger galaxy clusters, we present deep JVLA and Chandra observations of the HST Frontier Fields cluster MACS J0717.5+3745. The Chandra image shows a complex merger event, with at least four components belonging to different merging subclusters. NW of the cluster, $\sim 0.7$ Mpc from the center, there is a ram-pressure-stripped core that appears to have traversed the densest parts of the cluster after entering the ICM from the direction of a galaxy filament to the SE. We detect a density discontinuity NNE of this core which we speculate is associated with a cold front. Our radio images reveal new details for the complex radio relic and radio halo in this cluster. In addition, we discover several new filamentary radio sources with sizes of 100-300 kpc. A few of these seem to be connected to the main radio relic, while others are either embedded within the radio halo or projected onto it. A narrow-angled-tailed (NAT) radio galaxy, a cluster member, is located at the center of the radio relic. The steep spectrum tails of this AGN leads into the large radio relic where the radio spectrum flattens again. This morphological connection between the NAT radio galaxy and relic provides evidence for re-acceleration (revival) of fossil electrons. The presence of hot $\gtrsim 20$ keV ICM gas detected by Chandra near the relic location provides additional support for this re-acceleration scenario.

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.

Selecting sources with rising flux densities towards longer wavelengths from Herschel/SPIRE maps is an efficient way to produce a catalogue rich in high-redshift (z > 4) dusty star-forming galaxies. The effectiveness of this approach has already been confirmed by spectroscopic follow-up observations, but the previously available catalogues made this way are limited by small survey areas. Here we apply a map-based search method to 274 deg$^2$ of the HerMES Large Mode Survey (HeLMS) and create a catalogue of 477 objects with SPIRE flux densities $S_{500} > S_{350} >S_{250}$ and a $5 \sigma$ cut-off $S_{500} > $ 52 mJy. From this catalogue we determine that the total number of these "red" sources is at least an order of magnitude higher than predicted by galaxy evolution models. These results are in agreement with previous findings in smaller HerMES fields; however, due to our significantly larger sample size we are also able to investigate the shape of the red source counts for the first time. We have obtained spectroscopic redshift measurements for two of our sources using the Atacama Large Millimeter/submillimeter Array (ALMA). The redshifts z = 5.1 and z = 3.8 confirm that with our selection method we can indeed find high-redshift dusty star-forming galaxies.

We report on high-resolution JVLA and Chandra observations of the HST Frontier Cluster MACS J0717.5+3745. MACS J0717.5+3745 offers the largest contiguous magnified area of any known cluster, making it a promising target to search for lensed radio and X-ray sources. With the high-resolution 1.0-6.5 GHz JVLA imaging in A and B configuration, we detect a total of 51 compact radio sources within the area covered by the HST imaging. Within this sample we find 7 lensed sources with amplification factors larger than $2$. None of these sources are identified as multiply-lensed. Based on the radio luminosities, the majority of these sources are likely star forming galaxies with star formation rates of 10-50 M$_\odot$ yr$^{-1}$ located at $1 \lesssim z \lesssim 2$. Two of the lensed radio sources are also detected in the Chandra image of the cluster. These two sources are likely AGN, given their $2-10$ keV X-ray luminosities of $\sim 10^{43-44}$ erg s$^{-1}$. From the derived radio luminosity function, we find evidence for an increase in the number density of radio sources at $0.6<z<2.0$, compared to a $z < 0.3$ sample. Our observations indicate that deep radio imaging of lensing clusters can be used to study star forming galaxies, with star formation rates as low as $\sim10$ M$_{\odot}$ yr$^{-1}$, at the peak of cosmic star formation history.

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

The HST Frontier Fields cluster MACS J1149.6+2223 is one of the most complex merging clusters, believed to consist of four dark matter halos. We present results from deep (365 ks) Chandra observations of the cluster, which reveal the most distant cold front (z=0.544) discovered to date. In the cluster outskirts, we also detect hints of a surface brightness edge that could be the bow shock preceding the cold front. The substructure analysis of the cluster identified several components with large relative radial velocities, thus indicating that at least some collisions occur almost along the line of sight. The inclination of the mergers with respect to the plane of the sky poses significant observational challenges at X-ray wavelengths. MACS J1149.6+2223 possibly hosts a steep-spectrum radio halo. If the steepness of the radio halo is confirmed, then the radio spectrum, combined with the relatively regular ICM morphology, could indicate that MACS J1149.6+2223 is an old merging cluster.

We present Sunyaev-Zel'dovich (SZ) effect measurements from wide-field images towards the galaxy cluster RX J1347.5-1145 obtained from the Caltech Submillimeter Observatory with the Multiwavelength Submillimeter Inductance Camera (MUSIC) at 147, 213, 281, and 337 GHz and with Bolocam at 140 GHz. As part of our analysis, we have used higher frequency data from Herschel-SPIRE and previously published lower frequency radio data to subtract the signal from the brightest dusty star-forming galaxies behind RX J1347.5-1145 and from the AGN in RX J1347.5-1145's BCG. Using these five-band SZ effect images, combined with X-ray spectroscopic measurements of the temperature of the intra-cluster medium (ICM) from Chandra, we constrain the ICM optical depth to be $\tau_e = 7.33^{+0.96}_{-0.97} \times 10^{-3}$ and the ICM line of sight peculiar velocity to be $v_{pec} = -1040^{+870}_{-840}$ km s$^{-1}$. The errors for both quantities are limited by measurement noise rather than calibration uncertainties or astrophysical contamination, and significant improvements are possible with deeper observations. Our best-fit velocity is in good agreement with one previously published SZ effect analysis and in mild tension with the other, although some or all of that tension may be because that measurement samples a much smaller cluster volume. Furthermore, our best-fit optical depth implies a gas mass slightly larger than the Chandra-derived value, implying the cluster is elongated along the line of sight.

Morphology is often used to infer the state of relaxation of galaxy clusters. The regularity, symmetry, and degree to which a cluster is centrally concentrated inform quantitative measures of cluster morphology. The Cluster Lensing and Supernova survey with Hubble Space Telescope (CLASH) used weak and strong lensing to measure the distribution of matter within a sample of 25 clusters, 20 of which were deemed to be relaxed based on their X-ray morphology and alignment of the X-ray emission with the BCG. Towards a quantitative characterization of this important sample of clusters, we present uniformly estimated X-ray morphological statistics for all 25 CLASH clusters. We compare X-ray morphologies of CLASH clusters with those identically measured for a large sample of simulated clusters from the MUSIC-2 simulations, selected by mass. We confirm a threshold in X-ray surface brightness concentration of C>0.4 for cool-core clusters, where C is the ratio of X-ray emission inside 100 kpc/h70 compared to inside 500 kpc/h70. We report and compare morphologies of these clusters inferred from Sunyaev-Zeldovich Effect (SZE) maps of the hot gas and in from projected mass maps based on strong and weak lensing. We find a strong agreement in alignments of the orientation of major axes for the lensing, X-ray, and SZE maps of nearly all of the CLASH clusters at radii of 500 kpc (approximately 0.5R500 for these clusters). We also find a striking alignment of clusters shapes at the 500 kpc scale, as measured with X-ray, SZE, and lensing, with that of the near-infrared stellar light at 10 kpc scales for the 20 "relaxed" clusters. This strong alignment indicates a powerful coupling between the cluster- and galaxy-scale galaxy formation processes.

Merging galaxy clusters leave long-lasting signatures on the baryonic and non-baryonic cluster constituents, including shock fronts, cold fronts, X-ray substructure, radio halos, and offsets between the dark matter and the gas components. Using observations from Chandra, the Jansky Very Large Array, the Giant Metrewave Radio Telescope, and the Hubble Space Telescope, we present a multiwavelength analysis of the merging Frontier Fields cluster MACS J0416.1-2403 (z=0.396), which consists of a NE and a SW subclusters whose cores are separated on the sky by ~250 kpc. We find that the NE subcluster has a compact core and hosts an X-ray cavity, yet it is not a cool core. Approximately 450 kpc south-south west of the SW subcluster, we detect a density discontinuity that corresponds to a compression factor of ~1.5. The discontinuity was most likely caused by the interaction of the SW subcluster with a less massive structure detected in the lensing maps SW of the subcluster's center. For both the NE and the SW subclusters, the dark matter and the gas components are well-aligned, suggesting that MACS J0416.1-2403 is a pre-merging system. The cluster also hosts a radio halo, which is unusual for a pre-merging system. The halo has a 1.4 GHz power of (1.06 +/- 0.09) x 10^{24} W Hz^{-1}, which is somewhat lower than expected based on the X-ray luminosity of the cluster. We suggest that we are either witnessing the birth of a radio halo, or have discovered a rare ultra-steep spectrum halo.

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.

We present profiles of temperature (Tx), gas mass, and hydrostatic mass estimated from new and archival X-ray observations of CLASH clusters. We compare measurements derived from XMM and Chandra observations with one another and compare both to gravitational lensing mass profiles derived with CLASH HST and ground-based lensing data. Radial profiles of Chandra and XMM electron density and enclosed gas mass are nearly identical, indicating that differences in hydrostatic masses inferred from X-ray observations arise from differences in Tx measurements. Encouragingly, cluster Txs are consistent with one another at ~100-200 kpc radii but XMM Tx systematically decline relative to Chandra Tx at larger radii. The angular dependence of the discrepancy suggests additional investigation on systematics such as the XMM point spread function correction, vignetting and off-axis responses. We present the CLASH-X mass-profile comparisons in the form of cosmology-independent and redshift-independent circular-velocity profiles. Ratios of Chandra HSE mass profiles to CLASH lensing profiles show no obvious radial dependence in the 0.3-0.8 Mpc range. However, the mean mass biases inferred from the WL and SaWLens data are different. e.g., the weighted-mean value at 0.5 Mpc is <b> = 0.12 for the WL comparison and <b> = -0.11 for the SaWLens comparison. The ratios of XMM HSE mass profiles to CLASH lensing profiles show a pronounced radial dependence in the 0.3-1.0 Mpc range, with a weighted mean mass bias of value rising to <b>~0.3 at ~1 Mpc for the WL comparison and <b> of 0.25 for SaWLens comparison. The enclosed gas mass profiles from both Chandra and XMM rise to a value 1/8 times the total-mass profiles inferred from lensing at 0.5 Mpc and remain constant outside of that radius, suggesting that [8xMgas] profiles may be an excellent proxy for total-mass profiles at >0.5 Mpc in massive galaxy clusters.

We present a technique to constrain galaxy cluster pressure profiles by jointly fitting Sunyaev-Zel'dovich effect (SZE) data obtained with MUSTANG and Bolocam for the clusters Abell 1835 and MACS0647. Bolocam and MUSTANG probe different angular scales and are thus highly complementary. We find that the addition of the high resolution MUSTANG data can improve constraints on pressure profile parameters relative to those derived solely from Bolocam. In Abell 1835 and MACS0647, we find gNFW inner slopes of $\gamma = 0.36_{-0.21}^{+0.33}$ and $\gamma = 0.38_{-0.25}^{+0.20}$, respectively when $\alpha$ and $\beta$ are constrained to 0.86 and 4.67 respectively. The fitted SZE pressure profiles are in good agreement with X-ray derived pressure profiles.

We present 20 WISE-selected galaxies with bolometric luminosities L_bol > 10^14 L_sun, including five with infrared luminosities L_IR = L(rest 8-1000 micron) > 10^14 L_sun. These "extremely luminous infrared galaxies," or ELIRGs, were discovered using the "W1W2-dropout" selection criteria which requires marginal or non-detections at 3.4 and 4.6 micron (W1 and W2, respectively) but strong detections at 12 and 22 micron in the WISE survey. Their spectral energy distributions are dominated by emission at rest-frame 4-10 micron, suggesting that hot dust with T_d ~ 450K is responsible for the high luminosities. These galaxies are likely powered by highly obscured AGNs, and there is no evidence suggesting these systems are beamed or lensed. We compare this WISE-selected sample with 116 optically selected quasars that reach the same L_bol level, corresponding to the most luminous unobscured quasars in the literature. We find that the rest-frame 5.8 and 7.8 micron luminosities of the WISE-selected ELIRGs can be 30-80% higher than that of the unobscured quasars. The existence of AGNs with L_bol > 10^14 L_sun at z > 3 suggests that these supermassive black holes are born with large mass, or have very rapid mass assembly. For black hole seed masses ~ 10^3 M_sun, either sustained super-Eddington accretion is needed, or the radiative efficiency must be <15%, implying a black hole with slow spin, possibly due to chaotic accretion.

We present high resolution (9$^{\prime \prime}$) imaging of the Sunyaev-Zel'dovich Effect (SZE) toward two massive galaxy clusters, MACS J0647.7+7015 ($z=0.591$) and MACS J1206.2-0847 ($z=0.439$). We compare these 90 GHz measurements, taken with the MUSTANG receiver on the Green Bank Telescope, with generalized Navarro-Frenk-White (gNFW) models derived from Bolocam 140 GHz SZE data as well as maps of the thermal gas derived from {\it Chandra} X-ray observations. For MACS J0647.7+7015, we find a gNFW profile with core slope parameter $\gamma= 0.9$ fits the MUSTANG image with $\chi^{2}_{red}=1.005$ and probability to exceed (PTE) = 0.34. For MACS J1206.2-0847, we find $\gamma=0.7$, $\chi^{2}_{red}=0.993$, and PTE = 0.70. In addition, we find a significant ($>$3-$\sigma$) residual SZE feature in MACS J1206.2-0847 coincident with a group of galaxies identified in VLT data and filamentary structure found in a weak-lensing mass reconstruction. We suggest the detected sub-structure may be the SZE decrement from a low mass foreground group or an infalling group. GMRT measurements at 610 MHz reveal diffuse extended radio emission to the west, which we posit is either an AGN-driven radio lobe, a bubble expanding away from disturbed gas associated with the SZE signal, or a bubble detached and perhaps re-accelerated by sloshing within the cluster. Using the spectroscopic redshifts available, we find evidence for a foreground ($z=0.423$) or infalling group, coincident with the residual SZE feature.

CCAT will be a 25m diameter sub-millimeter telescope capable of operating in the 0.2 to 2.1mm wavelength range. It will be located at an altitude of 5600m on Cerro Chajnantor in northern Chile near the ALMA site. The anticipated first generation instruments include large format (60,000 pixel) kinetic inductance detector (KID) cameras, a large format heterodyne array and a direct detection multi-object spectrometer. The paper describes the architecture of the CCAT software and the development strategy.

Current ground-based submillimeter instruments (e.g. SCUBA-2, SHARC-2 and LABOCA) have hundreds to thousands of detectors, sampled at tens to hundreds of hertz, generating up to hundreds of gigabytes per night. Since noise is correlated between detectors and in time, due to atmospheric signals and temperature oscillations, naive map-making is not applicable. In addition, the size of the data sets makes direct likelihood based inversion techniques intractable. As a result, the data reduction approach for most current submm cameras is to adopt iterative methods in order to separate noise from sky signal, and hence effectively produce astronomical images. We investigate how today's map-makers scale to the next generation of instruments, which will have tens of thousands of detectors sampled at thousands of hertz, leading to data sets of challenging size. We propose strategies for reducing such large data sets.

We present a new determination of the concentration-mass relation for galaxy clusters based on our comprehensive lensing analysis of 19 X-ray selected galaxy clusters from the Cluster Lensing and Supernova Survey with Hubble (CLASH). Our sample spans a redshift range between 0.19 and 0.89. We combine weak lensing constraints from the Hubble Space Telescope (HST) and from ground-based wide field data with strong lensing constraints from HST. The result are reconstructions of the surface-mass density for all CLASH clusters on multi-scale grids. Our derivation of NFW parameters yields virial masses between 0.53 x 10^15 and 1.76 x 10^15 M_sol/h and the halo concentrations are distributed around c_200c ~ 3.7 with a 1-sigma significant negative trend with cluster mass. We find an excellent 4% agreement between our measured concentrations and the expectation from numerical simulations after accounting for the CLASH selection function based on X-ray morphology. The simulations are analyzed in 2D to account for possible biases in the lensing reconstructions due to projection effects. The theoretical concentration-mass (c-M) relation from our X-ray selected set of simulated clusters and the c-M relation derived directly from the CLASH data agree at the 90% confidence level.

We present the results of a numerical study based on the analysis of the MUSIC-2 simulations, aimed at estimating the expected concentration-mass relation for the CLASH cluster sample. We study nearly 1400 halos simulated at high spatial and mass resolution, which were projected along many lines-of-sight each. We study the shape of both their density and surface-density profiles and fit them with a variety of radial functions, including the Navarro-Frenk-White, the generalised Navarro-Frenk-White, and the Einasto density profiles. We derive concentrations and masses from these fits and investigate their distributions as a function of redshift and halo relaxation. We use the X-ray image simulator X-MAS to produce simulated Chandra observations of the halos and we use them to identify objects resembling the X-ray morphologies and masses of the clusters in the CLASH X-ray selected sample. We also derive a concentration-mass relation for strong-lensing clusters. We find that the sample of simulated halos which resemble the X-ray morphology of the CLASH clusters is composed mainly by relaxed halos, but it also contains a significant fraction of un-relaxed systems. For such a sample we measure an average 2D concentration which is ~11% higher than found for the full sample of simulated halos. After accounting for projection and selection effects, the average NFW concentrations of CLASH clusters are expected to be intermediate between those predicted in 3D for relaxed and super-relaxed halos. Matching the simulations to the individual CLASH clusters on the basis of the X-ray morphology, we expect that the NFW concentrations recovered from the lensing analysis of the CLASH clusters are in the range [3-6], with an average value of 3.87 and a standard deviation of 0.61. Simulated halos with X-ray morphologies similar to those of the CLASH clusters are affected by a modest orientation bias.

CCAT will be a 25-meter telescope for sub millimeter astronomy located at 5600m altitude on Cerro Chajnantor in northern Chile. CCAT will combine high sensitivity, a wide field of view, and a broad wavelength range (0.35 to 2.1mm) to provide an unprecedented capability for deep, large-area multicolor submillimeter surveys. It is planned to have a suite of instruments including large format KID cameras, a large heterodyne array and a KID-based direct detection multi-object spectrometer. The remote location drives a desire for fully autonomous observing coupled with data reduction pipelines and fast feedback to principal investigators.

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

On behalf of the WISE Science team, we present the discovery of a class of distant dust-enshrouded galaxies with extremely high luminosity. These galaxies are selected to have extreme red colors in the mid-IR using NASA's Wide-field Infrared Survey Explorer (WISE). They are faint in the optical and near-IR, predominantly at z=2-4, and with IR luminosity > $10^{13}\, L_{Sun}$, making them Hyper-Luminous Infrared Galaxies (HyLIRGs). SEDs incorporating the WISE, Spitzer, and Herschel PACS and SPIRE photometry indicate hot dust dominates the bolometric luminosity, presumably powered by AGN. Preliminary multi-wavelength follow-up suggests that they are different from normal populations in the local M-sigma relation. Their low source density implies that these objects are either intrinsically rare, or a short-lived phase in a more numerous population. If the latter is the case, these hot, dust-enshrouded galaxies may be an early stage in the interplay between AGN and galaxies.

We present a method for selecting $z>4$ dusty, star forming galaxies (DSFGs) using Herschel/SPIRE 250/350/500 $\mu m$ flux densities to search for red sources. We apply this method to 21 deg$^2$ of data from the HerMES survey to produce a catalog of 38 high-$z$ candidates. Follow-up of the first 5 of these sources confirms that this method is efficient at selecting high-$z$ DSFGs, with 4/5 at $z=4.3$ to $6.3$ (and the remaining source at $z=3.4$), and that they are some of the most luminous dusty sources known. Comparison with previous DSFG samples, mostly selected at longer wavelengths (e.g., 850 $\mu m$) and in single-band surveys, shows that our method is much more efficient at selecting high-$z$ DSFGs, in the sense that a much larger fraction are at $z>3$. Correcting for the selection completeness and purity, we find that the number of bright ($S_{500\,\mu m} \ge 30$ mJy), red Herschel sources is $3.3 \pm 0.8$ deg$^{-2}$. This is much higher than the number predicted by current models, suggesting that the DSFG population extends to higher redshifts than previously believed. If the shape of the luminosity function for high-$z$ DSFGs is similar to that at $z\sim2$, rest-frame UV based studies may be missing a significant component of the star formation density at $z=4$ to $6$, even after correction for extinction.

The galaxy cluster MACS J0717.5+3745 (z=0.55) is the largest known cosmic lens, with complex internal structures seen in deep X-ray, Sunyaev-Zel'dovich effect and dynamical observations. We perform a combined weak and strong lensing analysis with wide-field BVRi'z' Subaru/Suprime-Cam observations and 16-band Hubble Space Telescope observations taken as part of the Cluster Lensing And Supernova survey with Hubble (CLASH). We find consistent weak distortion and magnification measurements of background galaxies, and combine these signals to construct an optimally estimated radial mass profile of the cluster and its surrounding large-scale structure out to 5 Mpc/h. We find consistency between strong-lensing and weak-lensing in the region where these independent data overlap, <500 kpc/h. The two-dimensional weak-lensing map reveals a clear filamentary structure traced by distinct mass halos. We model the lensing shear field with 9 halos, including the main cluster, corresponding to mass peaks detected above 2.5\sigma_\kappa. The total mass of the cluster as determined by the different methods is M_{vir}=(2.8\pm0.4) \times 10^15 M_sun. Although this is the most massive cluster known at z>0.5, in terms of extreme value statistics we conclude that the mass of MACS J0717.5+3745 by itself is not in serious tension with LambdaCDM, representing only a ~2{\sigma} departure above the maximum simulated halo mass at this redshift.

The Sunyaev-Zel'dovich (SZ) effect spectrum crosses through a null where dT_CMB = 0 near nu_0 = 217 GHz. In a cluster of galaxies, nu_0 can be shifted from the canonical thermal SZ effect value by corrections to the SZ effect scattering due to the properties of the inter-cluster medium. We have measured the SZ effect in the hot galaxy cluster RX J1347.5-1145 with Z-Spec, an R ~ 300 grating spectrometer sensitive between 185 and 305 GHz. These data comprise a high spectral resolution measurement around the null of the SZ effect and clearly exhibit the transition from negative to positive dT_CMB over the Z-Spec band. The SZ null position is measured to be nu_0 = 225.8 \pm 2.5 (stat.) \pm 1.2 (sys.) GHz, which differs from the canonical null frequency by 3.0 sigma and is evidence for modifications to the canonical thermal SZ effect shape. Assuming the measured shift in nu_0 is due only to relativistic corrections to the SZ spectrum, we place the limit T_e = 17.1 \pm 5.3 keV from the zero-point measurement alone. By simulating the response of the instrument to the sky, we are able to generate likelihood functions in {y_0, T_e, v_pec} space. For v_pec = 0 km/s, we measure the best fitting SZ model to be y_0 = 4.6 (+0.6, -0.9) x 10^-4, T_e,0 = 15.2 (+12,-7.4) keV. When v_pec is allowed to vary, a most probable value of v_pec = +450 \pm 810 km/s is found.

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.

We present the status of MUSIC, the MUltiwavelength Sub/millimeter Inductance Camera, a new instrument for the Caltech Submillimeter Observatory. MUSIC is designed to have a 14', diffraction-limited field-of-view instrumented with 2304 detectors in 576 spatial pixels and four spectral bands at 0.87, 1.04, 1.33, and 1.98 mm. MUSIC will be used to study dusty star-forming galaxies, galaxy clusters via the Sunyaev-Zeldovich effect, and star formation in our own and nearby galaxies. MUSIC uses broadband superconducting phased-array slot-dipole antennas to form beams, lumped-element on-chip bandpass filters to define spectral bands, and microwave kinetic inductance detectors to senseincoming light. The focal plane is fabricated in 8 tiles consisting of 72 spatial pixels each. It is coupled to the telescope via an ambient temperature ellipsoidal mirror and a cold reimaging lens. A cold Lyot stop sits at the image of the primary mirror formed by the ellipsoidal mirror. Dielectric and metal mesh filters are used to block thermal infrared and out-of-band radiation. The instrument uses a pulse tube cooler and 3He/3He/4He closed-cycle cooler to cool the focal plane to below 250 mK. A multilayer shield attenuates Earth's magnetic field. Each focal plane tile is read out by a single pair of coaxes and a HEMT amplifier. The readout system consists of 16 copies of custom-designed ADC/DAC and IF boards coupled to the CASPER ROACH platform. We focus on recent updates on the instrument design and results from the commissioning of the full camera in 2012.

We report the discovery by the Wide-field Infrared Survey Explorer of the z = 2.452 source WISE J181417.29+341224.9, the first hyperluminous source found in the WISE survey. WISE 1814+3412 is also the prototype for an all-sky sample of ~1000 extremely luminous "W1W2-dropouts" (sources faint or undetected by WISE at 3.4 and 4.6 microns and well detected at 12 or 22 microns). The WISE data and a 350 micron detection give a minimum bolometric luminosity of 3.7 x 10^13 Lsun, with ~10^14 Lsun plausible. Followup images reveal four nearby sources: a QSO and two Lyman Break Galaxies (LBGs) at z = 2.45, and an M dwarf star. The brighter LBG dominates the bolometric emission. Gravitational lensing is unlikely given the source locations and their different spectra and colors. The dominant LBG spectrum indicates a star formation rate ~300 Msun/yr, accounting for < 10% of the bolometric luminosity. Strong 22 micron emission relative to 350 microns implies that warm dust contributes significantly to the luminosity, while cooler dust normally associated with starbursts is constrained by an upper limit at 1.1 mm. Radio emission is ~10x above the far-infrared/radio correlation, indicating an active galactic nucleus is present. An obscured AGN combined with starburst and evolved stellar components can account for the observations. If the black hole mass follows the local M_BH-bulge mass relation, the implied Eddington ratio is >~4. WISE 1814+3412 may be a heavily obscured object where the peak AGN activity occurred prior to the peak era of star formation.

We have used the Caltech Submillimeter Observatory (CSO) to follow-up a sample of WISE-selected, hyperluminous galaxies, so called W1W2-dropout galaxies. This is a rare (~ 1000 all-sky) population of galaxies at high redshift (peaks at z=2-3), that are faint or undetected by WISE at 3.4 and 4.6 um, yet are clearly detected at 12 and 22 um. The optical spectra of most of these galaxies show significant AGN activity. We observed 14 high-redshift (z > 1.7) W1W2-dropout galaxies with SHARC-II at 350 to 850 um, with 9 detections; and observed 18 with Bolocam at 1.1 mm, with five detections. Warm Spitzer follow-up of 25 targets at 3.6 and 4.5 um, as well as optical spectra of 12 targets are also presented in the paper. Combining WISE data with observations from warm Spitzer and CSO, we constructed their mid-IR to millimeter spectral energy distributions (SEDs). These SEDs have a consistent shape, showing significantly higher mid-IR to submm ratios than other galaxy templates, suggesting a hotter dust temperature. We estimate their dust temperatures to be 60-120 K using a single-temperature model. Their infrared luminosities are well over 10^{13} Lsun. These SEDs are not well fitted with existing galaxy templates, suggesting they are a new population with very high luminosity and hot dust. They are likely among the most luminous galaxies in the Universe. We argue that they are extreme cases of luminous, hot dust-obscured galaxies (DOGs), possibly representing a short evolutionary phase during galaxy merging and evolution. A better understanding of their long-wavelength properties needs ALMA as well as Herschel data.

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.

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.

Bolocam began collecting science data in 2003 as the long-wavelength imaging camera at the Caltech Submillimeter Observatory. The planets, along with a handful of secondary calibrators, have been used to determine the flux calibration for all of the data collected with Bolocam. Uranus and Neptune stand out as the only two planets that are bright enough to be seen with high signal-to-noise in short integrations without saturating the standard Bolocam readout electronics. By analyzing all of the 143 GHz observations made with Bolocam between 2003 and 2010, we find that the brightness ratio of Uranus to Neptune is 1.027 +- 0.006, with no evidence for any variations over that period. Including previously published results at \simeq 150 GHz, we find a brightness ratio of 1.029 +- 0.006 with no evidence for time variability over the period 1983-2010. Additionally, we find no evidence for time-variability in the brightness ratio of either Uranus or Neptune to the ultracompact HII region G34.3 or the protostellar source NGC 2071IR. Using recently published WMAP results we constrain the absolute 143 GHz brightness of both Uranus and Neptune to ~3%. Finally, we present ~3% absolute 143 GHz peak flux density values for the ultracompact HII regions G34.3 and K3-50A and the protostellar source NGC 2071IR.

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.

The Cluster Lensing And Supernova survey with Hubble (CLASH) is a 524-orbit multi-cycle treasury program to use the gravitational lensing properties of 25 galaxy clusters to accurately constrain their mass distributions. The survey, described in detail in this paper, will definitively establish the degree of concentration of dark matter in the cluster cores, a key prediction of CDM. The CLASH cluster sample is larger and less biased than current samples of space-based imaging studies of clusters to similar depth, as we have minimized lensing-based selection that favors systems with overly dense cores. Specifically, twenty CLASH clusters are solely X-ray selected. The X-ray selected clusters are massive (kT > 5 keV; 5 - 30 x 10^14 M_solar) and, in most cases, dynamically relaxed. Five additional clusters are included for their lensing strength (Einstein radii > 35 arcsec at z_source = 2) to further quantify the lensing bias on concentration, to yield high resolution dark matter maps, and to optimize the likelihood of finding highly magnified high-redshift (z > 7) galaxies. The high magnification, in some cases, provides angular resolutions unobtainable with any current UVOIR facility and can yield z > 7 candidates bright enough for spectroscopic follow-up. A total of 16 broadband filters, spanning the near-UV to near-IR, are employed for each 20-orbit campaign on each cluster. These data are used to measure precise (sigma_phz < 0.02(1+z)) photometric redshifts for dozens of newly discovered multiply-lensed images per cluster. Observations of each cluster are spread over 8 epochs to enable a search, primarily in the parallel fields, for Type Ia supernovae at z > 1 to improve constraints on the time dependence of the dark energy equation of state and the evolution of such supernovae in an epoch when the universe is matter dominated.

Extensive surveys of galaxy clusters motivate us to assess the likelihood of cluster-cluster lensing (CCL), namely, gravitational-lensing of a background cluster by a foreground cluster. We briefly describe the characteristics of CCLs in optical, X-ray and SZ measurements, and calculate their predicted numbers for $\Lambda$CDM parameters and a viable range of cluster mass functions and their uncertainties. The predicted number of CCLs in the strong-lensing regime varies from several ($<10$) to as high as a few dozen, depending mainly on whether lensing triaxiality bias is accounted for, through the c-M relation. A much larger number is predicted when taking into account also CCL in the weak-lensing regime. In addition to few previously suggested CCLs, we report a detection of a possible CCL in A383, where background candidate high-$z$ structures are magnified, as seen in deep Subaru observations.

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 present an analysis of an 8 arcminute diameter map of the area around the galaxy cluster Abell 1835 from jiggle map observations at a wavelength of 1.1 mm using the Bolometric Camera (Bolocam) mounted on the Caltech Submillimeter Observatory (CSO). The data is well described by a model including an extended Sunyaev-Zel'dovich (SZ) signal from the cluster gas plus emission from two bright background submm galaxies magnified by the gravitational lensing of the cluster. The best-fit values for the central Compton value for the cluster and the fluxes of the two main point sources in the field: SMM J140104+0252, and SMM J14009+0252 are found to be $y_{0}=(4.34\pm0.52\pm0.69)\times10^{-4}$, 6.5$\pm{2.0}\pm0.7$ mJy and 11.3$\pm{1.9}\pm1.1$ mJy, where the first error represents the statistical measurement error and the second error represents the estimated systematic error in the result. This measurement assumes the presence of dust emission from the cluster's central cD galaxy of $1.8\pm0.5$ mJy, based on higher frequency observations of Abell 1835. The cluster image represents one of the highest-significance SZ detections of a cluster in the positive region of the thermal SZ spectrum to date. The inferred central intensity is compared to other SZ measurements of Abell 1835 and this collection of results is used to obtain values for $y_{0} = (3.60\pm0.24)\times10^{-4}$ and the cluster peculiar velocity $v_{z} = -226\pm275$ km/s.

We report measurements of the fluctuations in atmospheric emission (atmospheric noise) above Mauna Kea recorded with Bolocam at 143 and 268 GHz from the Caltech Submillimeter Observatory (CSO). The 143 GHz data were collected during a 40 night observing run in late 2003, and the 268 GHz observations were made in early 2004 and early 2005 over a total of 60 nights. Below 0.5 Hz, the data time-streams are dominated by atmospheric noise in all observing conditions. The atmospheric noise data are consistent with a Kolmogorov-Taylor (K-T) turbulence model for a thin wind-driven screen, and the median amplitude of the fluctuations is 280 mK^2 rad^(-5/3) at 143 GHz and 4000 mK^2 rad^(-5/3) at 268 GHz. Comparing our results with previous ACBAR data, we find that the normalization of the power spectrum of the atmospheric noise fluctuations is a factor of 80 larger above Mauna Kea than above the South Pole at millimeter wavelengths. Most of this difference is due to the fact that the atmosphere above the South Pole is much drier than the atmosphere above Mauna Kea. However, the atmosphere above the South Pole is slightly more stable as well: the fractional fluctuations in the column depth of precipitable water vapor are a factor of sqrt(2) smaller at the South Pole compared to Mauna Kea. Based on our atmospheric modeling, we developed several algorithms to remove the atmospheric noise, and the best results were achieved when we described the fluctuations using a low-order polynomial in detector position over the 8 arcmin field of view (FOV). However, even with these algorithms, we were not able to reach photon-background-limited instrument photometer (BLIP) performance at frequencies below 0.5 Hz in any observing conditions.

In this paper we present a dust emission map of the starless core TMC-1C taken at 2100 microns. Along with maps at 160, 450, 850 and 1200 microns, we study the dust emissivity spectral index from the (sub)millimeter spectral energy distribution, and find that it is close to the typically assumed value of beta = 2. We also map the dust temperature and column density in TMC-1C, and find that at the position of the dust peak (A_V ~ 50), the line-of-sight-averaged temperature is ~7 K. Employing simple Monte Carlo modeling, we show that the data are consistent with a constant value for the emissivity spectral index over the whole map of TMC-1C.

We have surveyed two science fields totaling one square degree with Bolocam at 2.1 mm to search for secondary CMB anisotropies caused by the Sunyaev- Zel'dovich effect (SZE). The fields are in the Lynx and Subaru/XMM SDS1 fields. Our survey is sensitive to angular scales with an effective angular multipole of l_eff = 5700 with FWHM_l = 2800 and has an angular resolution of 60 arcseconds FWHM. Our data provide no evidence for anisotropy. We are able to constrain the level of total astronomical anisotropy, modeled as a flat bandpower in C_l, with frequentist 68%, 90%, and 95% CL upper limits of 590, 760, and 830 uKCMB^2. We statistically subtract the known contribution from primary CMB anisotropy, including cosmic variance, to obtain constraints on the SZE anisotropy contribution. Now including flux calibration uncertainty, our frequentist 68%, 90% and 95% CL upper limits on a flat bandpower in C_l are 690, 960, and 1000 uKCMB^2. When we instead employ the analytic spectrum suggested by Komatsu and Seljak (2002), and account for the non-Gaussianity of the SZE anisotropy signal, we obtain upper limits on the average amplitude of their spectrum weighted by our transfer function of 790, 1060, and 1080 uKCMB^2. We obtain a 90% CL upper limit on sigma8, which normalizes the power spectrum of density fluctuations, of 1.57. These are the first constraints on anisotropy and sigma8 from survey data at these angular scales at frequencies near 150 GHz.

We present a large-scale millimeter continuum map of the Ophiuchus molecular cloud. Nearly 11 square degrees, including all of the area in the cloud with visual extinction more than 3 magnitudes, was mapped at 1.1 mm with Bolocam on the Caltech Submillimeter Observatory (CSO). By design, the map also covers the region mapped in the infrared with the Spitzer Space Telescope. We detect 44 definite sources, and a few likely sources are also seen along a filament in the eastern streamer. The map indicates that dense cores in Ophiuchus are very clustered and often found in filaments within the cloud. Most sources are round, as measured at the half power point, but elongated when measured at lower contour levels, suggesting spherical sources lying within filaments. The masses, for an assumed dust temperature of 10 K, range from 0.24 to 3.9 solar masses, with a mean value of 0.96 solar masses. The total mass in distinct cores is 42 solar masses, 0.5 to 2% of the total cloud mass, and the total mass above 4 sigma is about 80 solar masses. The mean densities in the cores are quite high, with an average of 1.6 x 10^6 per cc, suggesting short free-fall times. The core mass distribution can be fitted with a power law with slope of 2.1 plus or minus 0.3 for M>0.5 solar masses, similar to that found in other regions, but slightly shallower than that of some determinations of the local IMF. In agreement with previous studies, our survey shows that dense cores account for a very small fraction of the cloud volume and total mass. They are nearly all confined to regions with visual extinction at least 9 mag, a lower threshold than found previously.

We have completed a 1.1 mm continuum survey of 7.5 sq deg of the Perseus Molecular Cloud using Bolocam at the Caltech Submillimeter Observatory. This represents the largest millimeter or submillimeter continuum map of Perseus to date. Our map covers more than 30,000 31" (FWHM) resolution elements to a 1 sigma RMS of 15 mJy/beam. We detect a total of 122 cores above a 5 sigma point source mass detection limit of 0.18 M_sun, assuming a dust temperature of 10 K, 60 of which are new millimeter or submillimeter detections. The 1.1 mm mass function is consistent with a broken power law of slope -1.3 (0.5 M_sun<M<2.5 M_sun) and -2.6 (M>2.5 M_sun), similar to the local initial mass function slope. No more than 5% of the total cloud mass is contained in discrete 1.1 mm cores, which account for a total mass of 285 M_sun. We suggest an extinction threshold for millimeter cores of Av~5 mag, based on our calculation of the probability of finding a 1.1 mm core as a function of Av. Much of the cloud is devoid of compact millimeter emission; despite the significantly greater area covered compared to previous surveys, only 5-10 of the newly identified sources lie outside previously observed areas. The two-point correlation function confirms that dense cores in the cloud are highly structured, with significant clustering on scales as large as 2e5 AU. These 1.1 mm results, especially when combined with recently acquired c2d Spitzer Legacy data, will provide a census of dense cores and protostars in Perseus and improve our understanding of the earliest stages of star formation in molecular clouds.

We perform a fluctuation analysis of the 1.1mm Bolocam Lockman Hole Survey, which covers 324 square arcmin to a very uniform point source-filtered RMS noise level of 1.4 mJy/beam. The fluctuation analysis has the significant advantage of utilizing all of the available data. We constrain the number counts in the 1-10 mJy range, and derive significantly tighter constraints than in previous work: the power-law index is 2.7 (+0.18, -0.15), while the amplitude is equal to 1595 (+85,-238) sources per mJy per square degree, or N(>1 mJy) = 940 (+50,-140) sources/square degree (95% confidence). Our results agree extremely well with those derived from the extracted source number counts by Laurent et al (2005). Our derived normalization is about 2.5 times smaller than determined by MAMBO at 1.2mm by Greve et al (2004). However, the uncertainty in the normalization for both data sets is dominated by the systematic (i.e., absolute flux calibration) rather than statistical errors; within these uncertainties, our results are in agreement. We estimate that about 7% of the 1.1mm background has been resolved at 1 mJy.

We present results of a new deep 1.1 mm survey using Bolocam, a millimeter-wavelength bolometer array camera designed for mapping large fields at fast scan rates, without chopping. A map, galaxy candidate list, and derived number counts are presented. This survey encompasses 324 arcmin^2 to an rms noise level (filtered for point sources) of 1.4 mJy/beam and includes the entire regions surveyed by the published 8 mJy 850 micron JCMT SCUBA and 1.2 mm IRAM MAMBO surveys. We reduced the data using a custom software pipeline to remove correlated sky and instrument noise via a principal component analysis. Extensive simulations and jackknife tests were performed to confirm the robustness of our source candidates and estimate the effects of false detections, bias, and completeness. In total, 17 source candidates were detected at a significance > 3.0 sigma, with six expected false detections. Nine candidates are new detections, while eight candidates have coincident SCUBA 850 micron and/or MAMBO 1.2 mm detections. From our observed number counts, we estimate the underlying differential number count distribution of submillimeter galaxies and find it to be in general agreement with previous surveys. Modeling the spectral energy distributions of these submillimeter galaxies after observations of dusty nearby galaxies suggests extreme luminosities of L = 1.0-1.6 x 10^13 L_solar and, if powered by star formation, star formation rates of 500-800 M_solar/yr.