We present Gemini-S and {\it Spitzer}-IRAC optical-through-near-IR observations in the field of the SPT2349-56 proto-cluster at $z=4.3$. We detect optical/IR counterparts for only nine of the 14 submillimetre galaxies (SMGs) previously identified by ALMA in the core of SPT2349-56. In addition, we detect four $z\sim4$ Lyman-break galaxies (LBGs) in the 30 arcsec diameter region surrounding this proto-cluster core. Three of the four LBGs are new systems, while one appears to be a counterpart of one of the nine observed SMGs. We identify a candidate brightest cluster galaxy (BCG) with a stellar mass of $(3.2^{+2.5}_{-1.4})\times10^{11}\,{\rm M}_{\odot}$. The stellar masses of the eight other SMGs place them on, above, and below the main sequence of star formation at $z\approx4.5$. The cumulative stellar mass for the SPT2349-56 core is at least $(11.5\pm2.9)\times10^{11}\,{\rm M}_{\odot}$, a sizeable fraction of the stellar mass in local BCGs, and close to the universal baryon fraction (0.16) relative to the virial mass of the core ($10^{13}\,{\rm M}_{\odot}$). As all 14 of these SMGs are destined to quickly merge, we conclude that the proto-cluster core has already developed a significant stellar mass at this early stage, comparable to $z=1$ BCGs. Importantly, we also find that the SPT2349-56 core structure would be difficult to uncover in optical surveys, with none of the ALMA sources being easily identifiable or constrained through $g,r,$ and $i$ colour-selection in deep optical surveys and only a modest overdensity of LBGs over the extended core structure. SPT2349-56 therefore represents a truly dust-obscured phase of a massive cluster core under formation.

X-ray luminosity ($L_X$) originating from high-mass X-ray binaries (HMXBs) is tightly correlated with the host galaxy's star-formation rate (SFR). We explore this connection at sub-galactic scales spanning ${\sim}$7 dex in SFR and ${\sim}$8 dex in specific SFR (sSFR). There is good agreement with established relations down to ${\rm SFR {\simeq} 10^{-3}\,M_\odot \, yr^{-1}}$, below which an excess of X-ray luminosity emerges. This excess likely arises from low mass X-ray binaries. The intrinsic scatter of the $L_X$-SFR relation is constant, not correlated with SFR. Different star formation indicators scale with $L_X$ in different ways, and we attribute the differences to the effect of star formation history. The SFR derived from H$\alpha$ shows the tightest correlation with X-ray luminosity because H$\alpha$ emission probes stellar populations with ages similar to HMXB formation timescales, but the H$\alpha$-based SFR is reliable only for $\rm sSFR{>}10^{-12}\,M_\odot \, yr^{-1}/M_\odot$.

We describe the observations and resultant galaxy cluster catalog from the 2770 deg$^2$ SPTpol Extended Cluster Survey (SPT-ECS). Clusters are identified via the Sunyaev-Zel'dovich (SZ) effect, and confirmed with a combination of archival and targeted follow-up data, making particular use of data from the Dark Energy Survey (DES). With incomplete followup we have confirmed as clusters 244 of 266 candidates at a detection significance $\xi \ge 5$ and an additional 204 systems at $4<\xi<5$. The confirmed sample has a median mass of $M_{500c} \sim {4.4 \times 10^{14} M_\odot h_{70}^{-1}}$, a median redshift of $z=0.49$, and we have identified 44 strong gravitational lenses in the sample thus far. Radio data are used to characterize contamination to the SZ signal; the median contamination for confirmed clusters is predicted to be $\sim$1% of the SZ signal at the $\xi>4$ threshold, and $<4\%$ of clusters have a predicted contamination $>10\% $ of their measured SZ flux. We associate SZ-selected clusters, from both SPT-ECS and the SPT-SZ survey, with clusters from the DES redMaPPer sample, and find an offset distribution between the SZ center and central galaxy in general agreement with previous work, though with a larger fraction of clusters with significant offsets. Adopting a fixed Planck-like cosmology, we measure the optical richness-to-SZ-mass ($\lambda-M$) relation and find it to be 28% shallower than that from a weak-lensing analysis of the DES data---a difference significant at the 4 $\sigma$ level---with the relations intersecting at $\lambda=60$ . The SPT-ECS cluster sample will be particularly useful for studying the evolution of massive clusters and, in combination with DES lensing observations and the SPT-SZ cluster sample, will be an important component of future cosmological analyses.

Here we present new red sequence overdensity measurements for 77 fields in the high-$z$ Clusters Occupied by Bent Radio AGN (COBRA) survey, based on $r$- and $i$-band imaging taken with Lowell Observatory's Discovery Channel Telescope. We observe 38 COBRA fields in $r$-band and 90 COBRA fields in $i$-band. By combining the $r$- and $i$-band photometry with our 3.6$\mu$m and 4.5$\mu$m $Spitzer$ IRAC observations, we identify 39 red sequence cluster candidates that host a strong overdensity of galaxies when measuring the excess of red sequence galaxies relative to a background field. We initially treat the radio host as the cluster center and then determine a new cluster center based on the surface density of red sequence sources. Using our color selection, we identify which COBRA cluster candidates have strong red sequence populations. By removing foreground and background contaminants, we more securely determine which fields include cluster candidates with a higher significance than our single-band observations. Additionally, of the 77 fields we analyze with a redshift estimate, 26 include newly estimated photometric redshifts.

We present 16 new ultrabright $H_{AB}\lesssim25$ galaxy candidates at z~8 identified over the COSMOS/UltraVISTA field. The new search takes advantage of the deepest-available ground-based optical and near-infrared observations, including the DR3 release of UltraVISTA and full-depth Spitzer/IRAC observations from the SMUVS and SPLASH programs. Candidates are selected using Lyman-break criteria, combined with strict optical non-detection and SED-fitting criteria, minimizing contamination by low-redshift galaxies and low-mass stars. HST/WFC3 coverage from the DASH program reveals that one source evident in our ground-based near-IR data has significant substructure and may actually correspond to 3 separate z~8 objects, resulting in a sample of 18 galaxies, 10 of which seem to be fairly robust (with a >97% probability of being at z>7). The UV-continuum slope $\beta$ for the bright z~8 sample is $\beta=-2.2\pm0.6$, bluer but still consistent with that of similarly bright galaxies at z~6 ($\beta=-1.55\pm0.17$) and z~7 ($\beta=-1.75\pm0.18$). Their typical stellar masses are 10$^{9.1^{+0.5}_{-0.4}}M_{\odot}$, with the SFRs of $32^{+44}_{-32}M_{\odot}$/year, specific SFR of $4^{+8}_{-4}$ Gyr$^{-1}$, stellar ages of $\sim22^{+69}_{-22}$\,Myr, and low dust content A$_V=0.15^{+0.30}_{-0.15}$ mag. Using this sample we constrain the bright end of the z~8 UV luminosity function (LF). When combined with recent empty field LF estimates at z~8-9, the resulting z~8 LF can be equally well represented by either a Schechter or a double power-law (DPL) form. Assuming a Schechter parameterization, the best-fit characteristic magnitude is $M^*= -20.95^{+0.30}_{-0.35}$ mag with a very steep faint end slope $\alpha=-2.15^{+0.20}_{-0.19}$. These new candidates include amongst the brightest yet found at these redshifts, 0.5-1.0 mag brighter than found over CANDELS, providing excellent targets for follow-up studies.

The largest Herschel extragalactic surveys, H-ATLAS and HerMES, have selected a sample of "ultrared" dusty, star-forming galaxies (DSFGs) with rising SPIRE flux densities ($S_{500} > S_{350} > S_{250}$; so-called "500 $\mu$m-risers") as an efficient way for identifying DSFGs at higher redshift ($z > 4$). In this paper, we present a large Spitzer follow-up program of 300 Herschel ultrared DSFGs. We have obtained high-resolution ALMA, NOEMA, and SMA data for 63 of them, which allow us to securely identify the Spitzer/IRAC counterparts and classify them as gravitationally lensed or unlensed. Within the 63 ultrared sources with high-resolution data, $\sim$65% appear to be unlensed, and $\sim$27% are resolved into multiple components. We focus on analyzing the unlensed sample by directly performing multi-wavelength spectral energy distribution (SED) modeling to derive their physical properties and compare with the more numerous $z \sim 2$ DSFG population. The ultrared sample has a median redshift of 3.3, stellar mass of 3.7 $\times$ 10$^{11}$ $M_{\odot}$, star formation rate (SFR) of 730 $M_{\odot}$yr$^{-1}$, total dust luminosity of 9.0 $\times$ 10$^{12}$ $L_{\odot}$, dust mass of 2.8 $\times$ 10$^9$ $M_{\odot}$, and V-band extinction of 4.0, which are all higher than those of the ALESS DSFGs. Based on the space density, SFR density, and stellar mass density estimates, we conclude that our ultrared sample cannot account for the majority of the star-forming progenitors of the massive, quiescent galaxies found in infrared surveys. Our sample contains the rarer, intrinsically most dusty, luminous and massive galaxies in the early universe that will help us understand the physical drivers of extreme star formation.

We present catalogs of calibrated photometry and spectroscopic redshifts in the Extended Groth Strip, intended for studies of photometric redshifts (photo-z's). The data includes ugriz photometry from CFHTLS and Y-band photometry from the Subaru Suprime camera, as well as spectroscopic redshifts from the DEEP2, DEEP3 and 3D-HST surveys. These catalogs incorporate corrections to produce effectively matched-aperture photometry across all bands, based upon object size information available in the catalog and Moffat profile point spread function fits. We test this catalog with a simple machine learning-based photometric redshift algorithm based upon Random Forest regression, and find that the corrected aperture photometry leads to significant improvement in photo-z accuracy compared to the original SExtractor catalogs from CFHTLS and Subaru. The deep ugrizY photometry and spectroscopic redshifts are well-suited for empirical tests of photometric redshift algorithms for LSST. The resulting catalogs are publicly available. We include a basic summary of the strategy of the DEEP3 Galaxy Redshift Survey to accompany the recent public release of DEEP3 data.

We present a WFC3 F160W ($H$-band) selected catalog in the CANDELS/GOODS-N field containing photometry from the ultraviolet (UV) to the far-infrared (IR), photometric redshifts and stellar parameters derived from the analysis of the multi-wavelength data. The catalog contains 35,445 sources over the 171 arcmin$^{2}$ of the CANDELS F160W mosaic. The 5$\sigma$ detection limits (within an aperture of radius 0\farcs17) of the mosaic range between $H=27.8$, 28.2 and 28.7 in the wide, intermediate and deep regions, that span approximately 50\%, 15\% and 35\% of the total area. The multi-wavelength photometry includes broad-band data from UV (U band from KPNO and LBC), optical (HST/ACS F435W, F606W, F775W, F814W, and F850LP), near-to-mid IR (HST/WFC3 F105W, F125W, F140W and F160W, Subaru/MOIRCS Ks, CFHT/Megacam K, and \spitzer/IRAC 3.6, 4.5, 5.8, 8.0 $\mu$m) and far IR (\spitzer/MIPS 24$\mu$m, HERSCHEL/PACS 100 and 160$\mu$m, SPIRE 250, 350 and 500$\mu$m) observations. In addition, the catalog also includes optical medium-band data (R$\sim50$) in 25 consecutive bands, $\lambda=500$ to 950~nm, from the SHARDS survey and WFC3 IR spectroscopic observations with the G102 and G141 grisms (R$\sim210$ and 130). The use of higher spectral resolution data to estimate photometric redshifts provides very high, and nearly uniform, precision from $z=0-2.5$. The comparison to 1,485 good quality spectroscopic redshifts up to $z\sim3$ yields $\Delta z$/(1+$z_{\rm spec}$)$=$0.0032 and an outlier fraction of $\eta=$4.3\%. In addition to the multi-band photometry, we release added-value catalogs with emission line fluxes, stellar masses, dust attenuations, UV- and IR- based star formation rates and rest-frame colors.

We derive cosmological constraints using a galaxy cluster sample selected from the 2500~deg$^2$ SPT-SZ survey. The sample spans the redshift range $0.25< z<1.75$ and contains 343 clusters with SZ detection significance $\xi>5$. The sample is supplemented with optical weak gravitational lensing measurements of 32 clusters with $0.29<z<1.13$ (from Magellan and HST) and X-ray measurements of 89 clusters with $0.25<z<1.75$ (from Chandra). We rely on minimal modeling assumptions: i) weak lensing provides an accurate means of measuring halo masses, ii) the mean SZ and X-ray observables are related to the true halo mass through power-law relations in mass and dimensionless Hubble parameter $E(z)$ with a-priori unknown parameters, iii) there is (correlated, lognormal) intrinsic scatter and measurement noise relating these observables to their mean relations. We simultaneously fit for these astrophysical modeling parameters and for cosmology. Assuming a flat $\nu\Lambda$CDM model, in which the sum of neutrino masses is a free parameter, we measure $\Omega_\mathrm{m}=0.276\pm0.047$, $\sigma_8=0.781\pm0.037$, and $\sigma_8(\Omega_\mathrm{m}/0.3)^{0.2}=0.766\pm0.025$. The redshift evolution of the X-ray $Y_\mathrm{X}$-mass and $M_\mathrm{gas}$-mass relations are both consistent with self-similar evolution to within $1\sigma$. The mass-slope of the $Y_\mathrm{X}$-mass relation shows a $2.3\sigma$ deviation from self-similarity. Similarly, the mass-slope of the $M_\mathrm{gas}$-mass relation is steeper than self-similarity at the $2.5\sigma$ level. In a $\nu w$CDM cosmology, we measure the dark energy equation of state parameter $w=-1.55\pm0.41$ from the cluster data. We perform a measurement of the growth of structure since redshift $z\sim1.7$ and find no evidence for tension with the prediction from General Relativity. We provide updated redshift and mass estimates for the SPT sample. (abridged)

We present measurements of the stellar mass fractions ($f_\star$) for a sample of high-redshift ($0.93 \le z \le 1.32$) infrared-selected galaxy clusters from the Massive and Distant Clusters of WISE Survey (MaDCoWS) and compare them to the stellar mass fractions of Sunyaev-Zel'dovich (SZ) effect-selected clusters in a similar mass and redshift range from the South Pole Telescope (SPT)-SZ Survey. We do not find a significant difference in mean $f_\star$ between the two selection methods, though we do find an unexpectedly large range in $f_\star$ for the SZ-selected clusters. In addition, we measure the luminosity function of the MaDCoWS clusters and find $m^*= 19.41\pm0.07$, similar to other studies of clusters at or near our redshift range. Finally, we present SZ detections and masses for seven MaDCoWS clusters and new spectroscopic redshifts for five MaDCoWS clusters. One of these new clusters, MOO J1521+0452 at $z=1.31$, is the most distant MaDCoWS cluster confirmed to date.

In this work, we use measurements of galaxy stellar mass and two-point angular correlation functions to constrain the stellar-to-halo mass ratios (SHMRs) of passive and \np\ galaxies at $z\sim2-3$, as identified in the \emph{Spitzer} Matching Survey of the UltraVISTA ultra-deep Stripes (SMUVS). We adopt a sophisticated halo modeling approach to statistically divide our two populations into central and satellite galaxies. For central galaxies, we find that the normalization of the SHMR is greater for our passive population. Through the modeling of $\Lambda$ cold dark matter halo mass accretion histories, we show that this can only arise if the conversion of baryons into stars was more efficient at higher redshifts and additionally that passive galaxies can be plausibly explained as residing in halos with the highest formation redshifts (i.e., those with the lowest accretion rates) at a given halo mass. At a fixed stellar mass, satellite galaxies occupy host halos with a greater mass than central galaxies, and we find further that the fraction of passive galaxies that are satellites is higher than for the combined population. This, and our derived satellite quenching timescales, combined with earlier estimates from the literature, support dynamical/environmental mechanisms as the dominant process for satellite quenching at $z\lesssim3$.

While the upcoming telescopes will reveal correspondingly fainter, more distant galaxies, a question will persist: what more is there that these telescopes cannot see? One answer is the source-subtracted Cosmic Infrared Background (CIB). The CIB is comprised of the collective light from all sources remaining after known, resolved sources are accounted for. Ever-more-sensitive surveys will identify the brightest of these, allowing them to be removed, and - like peeling layers off an onion - reveal deeper layers of the CIB. In this way it is possible to measure the contributions from populations not accessible to direct telescopic observation. Measurement of fluctuations in the source-subtracted CIB, i.e., the spatial power spectrum of the CIB after subtracting resolved sources, provides a robust means of characterizing its faint, and potentially new, populations. Studies over the past 15 years have revealed source-subtracted CIB fluctuations on scales out to ~100' which cannot be explained by extrapolating from known galaxy populations. Moreover, they appear highly coherent with the unresolved Cosmic X-ray Background, hinting at a significant population of accreting black holes among the CIB sources. Characterizing the source-subtracted CIB with high accuracy, and thereby constraining the nature of the new populations, is feasible with upcoming instruments and would produce critically important cosmological information in the next decade. New coextensive deep and wide-area near-infrared, X-ray, and microwave surveys will bring decisive opportunities to examine, with high fidelity, the spatial spectrum and origin of the CIB fluctuations and their cross-correlations with cosmic microwave and X-ray backgrounds, and determine the formation epochs and the nature of the new sources (stellar nucleosynthetic or accreting black holes).

We present first results from a galaxy population study in the highest redshift galaxy clusters identified in the 2500 deg$^2$ South Pole Telescope Sunyaev Zel'dovich effect (SPT-SZ) survey. The cluster selection is to first order independent of galaxy properties, making the SPT-SZ sample particularly well suited for cluster galaxy population studies. We carry out a 4-band imaging campaign with the {\it Hubble} and {\it Spitzer} Space Telescopes of the five $z\gtrsim 1.4$, S/N$_{SZE}>$5 clusters, that are among the rarest most massive clusters known at this redshift. All five show clear overdensities of red galaxies whose colors agree with the initial cluster redshift estimates. The highest redshift cluster in this sample, SPT-CLJ0459-4947 at $z\sim1.72$, is the most distant $M_{500}>10^{14}~M_{\odot}$ ICM-selected cluster discovered thus far, and is one of only three known clusters in this mass range at $z\gtrsim 1.7$, regardless of selection. Based on UVJ-like photometric classification of quiescent and star-forming galaxies, the passive fraction in the cluster central regions ($r/r_{500}<0.7$) is higher than in the field at the same redshift, with corresponding environmental quenching efficiencies typically in the range $\sim0.5-0.8$ for stellar masses $\log(M/M_{\odot})>10.85$. We have explored the impact of emission from star formation on the selection of this sample, concluding that all five clusters studied here would still have been detected with S/N$_{SZE}>$5, even if they had the same passive fraction as measured in the field. Our results thus point towards an efficient suppression of star formation in massive galaxies in the central regions of the most massive clusters, occurring already earlier than $z\sim1.5$. [Abridged]

We present deep $J$ and $H$-band images in the extended Great Observatories Origins Deep Survey-North (GOODS-N) field covering an area of 0.22 $\rm{deg}^{2}$. The observations were taken using WIRCam on the 3.6-m Canada France Hawaii Telescope (CFHT). Together with the reprocessed $K_{\rm s}$-band image, the $5\sigma$ limiting AB magnitudes (in 2" diameter apertures) are 24.7, 24.2, and 24.4 AB mag in the $J$, $H$, and $K_{\rm s}$ bands, respectively. We also release a multi-band photometry and photometric redshift catalog containing 93598 sources. For non-X-ray sources, we obtained a photometric redshift accuracy $\sigma_{\mathrm{NMAD}}=0.036$ with an outlier fraction $\eta = 7.3\%$. For X-ray sources, which are mainly active galactic nuclei (AGNs), we cross-matched our catalog with the updated 2M-CDFN X-ray catalog from Xue et al. (2016) and found that 658 out of 683 X-ray sources have counterparts. $GALEX$ UV data are included in the photometric redshift computation for the X-ray sources to give $\sigma_{\mathrm{NMAD}} = 0.040$ with $\eta=10.5\%$. Our approach yields more accurate photometric redshift estimates compared to previous works in this field. In particular, by adopting AGN-galaxy hybrid templates, our approach delivers photometric redshifts for the X-ray counterparts with fewer outliers compared to the 3D-HST catalog, which fit these sources with galaxy-only templates.

We present a detailed census of galaxies in and around PC217.96+32.3, a spectroscopically confirmed Coma analog at z=3.78. Diverse galaxy types identified in the field include Lya emitters (LAEs), massive star-forming galaxies, and ultra-massive galaxies (log (Mstar/Msun)>= 11) which may have already halted their star formation. The sky distribution of the star-forming galaxies suggests the presence of a significant overdensity (delta_g=8+/-2), which is spatially offset from the previously confirmed members by 3-4 Mpc to the west. Candidate quiescent and post-starburst galaxies are also found in large excess (a factor of ~ 8-15 higher surface density than the field) although their redshifts are less certain. We estimate that the total enclosed mass traced by star-forming galaxy candidates is roughly comparable to that of PC217.96+32.3 traced by the LAEs. We speculate that the true extent of PC217.96+32.3 may be larger than previously known, a half of which is missed by our LAE selection. Alternatively, the newly discovered overdensity may belong to another Coma progenitor not associated with PC217.96+32.3. Expectations from theory suggest that both scenarios are equally unlikely (<1%), particularly in the cosmic volume probed in our survey. If confirmed as a single structure, its total mass will be well in excess of Coma, making this a singularly large cosmic structure rarely seen even in large cosmological simulations. Finally, we find that the protocluster galaxies follow the same SFR-M_star scaling relation as the field galaxies, suggesting that the environmental effect at z~4 is a subtle one at best for normal star-forming galaxies.

We present multi-wavelength global star formation rate (SFR) estimates for 326 galaxies from the Star Formation Reference Survey (SFRS) in order to determine the mutual scatter and range of validity of different indicators. The widely used empirical SFR recipes based on 1.4 GHz continuum, 8.0 $\mu$m polycyclic aromatic hydrocarbons (PAH), and a combination of far-infrared (FIR) plus ultraviolet (UV) emission are mutually consistent with scatter of $\raise{-0.8ex}\stackrel{\textstyle <}{\sim }$0.3 dex. The scatter is even smaller, $\raise{-0.8ex}\stackrel{\textstyle <}{\sim }$0.24 dex, in the intermediate luminosity range 9.3<log(L(60 $\mu$m/L$_\odot$)<10.7. The data prefer a non-linear relation between 1.4 GHz luminosity and other SFR measures. PAH luminosity underestimates SFR for galaxies with strong UV emission. A bolometric extinction correction to far-ultraviolet luminosity yields SFR within 0.2 dex of the total SFR estimate, but extinction corrections based on UV spectral slope or nuclear Balmer decrement give SFRs that may differ from the total SFR by up to 2 dex. However, for the minority of galaxies with UV luminosity ${>}5\times10^9$ L$_{\odot}$ or with implied far-UV extinction <1 mag, the UV spectral slope gives extinction corrections with 0.22~dex uncertainty.

Optical and near-mid-infrared reverberation mapping data obtained at Universit\"{a}tssternwarte Bochum in Chile and with the Spitzer Space Telescope allow us to explore the geometry of both the H$\alpha$ BLR and the dust torus for the nearby Seyfert 1 galaxy WPVS\,48. On average, the H$\alpha$ variations lag the blue AGN continuum by about 18 days, while the dust emission variations lag by 70 days in the J+K and by 90 days in the L+M bands. The IR echoes are sharp, while the H$\alpha$ echo is smeared. This together favours a bowl-shaped toroidal geometry where the dust sublimation radius is defined by a bowl surface, which is virtually aligned with a single iso-delay surface, thus leading to the sharp IR echoes. The BLR clouds, however, are located inside the bowl and spread over a range of iso-delay surfaces, leading to a smeared echo.

Sagittarius A* (Sgr A*) is the variable radio, near-infrared (NIR), and X-ray source associated with accretion onto the Galactic center black hole. We present an analysis of the most comprehensive NIR variability dataset of Sgr A* to date: eight 24-hour epochs of continuous monitoring of Sgr A* at 4.5 $\mu$m with the IRAC instrument on the Spitzer Space Telescope, 93 epochs of 2.18 $\mu$m data from Naos Conica at the Very Large Telescope, and 30 epochs of 2.12 $\mu$m data from the NIRC2 camera at the Keck Observatory, in total 94,929 measurements. A new approximate Bayesian computation method for fitting the first-order structure function extracts information beyond current Fast Fourier Transformation (FFT) methods of power spectral density (PSD) estimation. With a combined fit of the data of all three observatories, the characteristic coherence timescale of Sgr A* is $\tau_{b} = 243^{+82}_{-57}$ minutes ($90\%$ credible interval). The PSD has no detectable features on timescales down to 8.5 minutes ($95\%$ credible level), which is the ISCO orbital frequency for a dimensionless spin parameter $a = 0.92$. One light curve measured simultaneously at 2.12 and 4.5 $\mu$m during a low flux-density phase gave a spectral index $\alpha_s = 1.6 \pm 0.1$ ($F_\nu \propto \nu^{-\alpha_s}$). This value implies that the Sgr A* NIR color becomes bluer during higher flux-density phases. The probability densities of flux densities of the combined datasets are best fit by log-normal distributions. Based on these distributions, the Sgr A* spectral energy distribution is consistent with synchrotron radiation from a non-thermal electron population from below 20 GHz through the NIR.

The Spitzer Matching Survey of the UltraVISTA Ultra-deep Stripes (SMUVS) has obtained the largest ultra-deep Spitzer maps to date in a single field of the sky. We considered the sample of about 66,000 SMUVS sources at $z=2-6$ to investigate the evolution of dusty and non-dusty galaxies with stellar mass through the analysis of the galaxy stellar mass function (GSMF). We further divide our non-dusty galaxy sample with rest-frame optical colours to isolate red quiescent (`passive') galaxies. At each redshift, we identify a characteristic stellar mass in the GSMF above which dusty galaxies dominate, or are at least as important as non-dusty galaxies. Below that stellar mass, non-dusty galaxies comprise about 80% of all sources, at all redshifts except at $z=4-5$. The percentage of dusty galaxies at $z=4-5$ is unusually high: 30-40% for $M_{*}=10^9 - 10^{10.5} \, \rm M_\odot$ and $>80\%$ at $M_*>10^{11} \, \rm M_\odot$, which indicates that dust obscuration is of major importance in this cosmic period. The overall percentage of massive ($\log_{10} (M_*/M_\odot)>10.6$) galaxies that are quiescent increases with decreasing redshift, reaching $>30\%$ at $z\sim2$. Instead, the quiescent percentage among intermediate-mass galaxies (with $\log_{10} (M_*/M_\odot)=9.7-10.6$) stays roughly constant at a $\sim 10\%$ level. Our results indicate that massive and intermediate-mass galaxies clearly have different evolutionary paths in the young Universe, and are consistent with the scenario of galaxy downsizing.

We report the discovery and spectroscopic confirmation of a very large star-forming Lyman Break galaxy, G6025, at z_spec=3.721+/-0.003. In the rest-frame ~2100A, G6025 subtends ~24 kpc in physical extent when measured from the 1.5-sigma isophote, in agreement with the parametric size measurements which yield the half-light radius of 4.9+/-0.5 kpc and the semi-major axis of 12.5+/-0.1 kpc. G6025 is also very UV-luminous (~5L*(z~4}) and young (~140+/-60 Myr). Despite its unusual size and luminosity, the stellar population parameters and dust reddening (M_star~M*(z~4)$, and E(B-V)=0.18+/-0.05) estimated from the integrated light, are similar to those of smaller galaxies at comparable redshifts. The ground-based morphology and spectroscopy show two dominant components, both located off-center, embedded in more diffuse emission. We speculate that G6025 may be a scaled-up version of chain galaxies seen in deep HST imaging, or alternatively, a nearly equal-mass merger involving two super-L* galaxies in its early stage. G6025 lies close to but not within a known massive protocluster at z=3.78. We find four companions within 6 Mpc from G6025, two of which lie within 1.6 Mpc. While the limited sensitivity of the existing spectroscopy does not allow us to robustly characterize the local environment of G6025, it likely resides in a locally overdense environment. The luminosity, size, and youth of G6025 make it uniquely suited to study the early formation of massive galaxies in the universe.

This paper describes new deep 3.6 and 4.5 micron imaging of three UltraVISTA near-infrared survey stripes within the COSMOS field. The observations were carried out with Spitzer's Infrared Array Camera (IRAC) for the Spitzer Matching Survey of the Ultra-VISTA Deep Stripes (SMUVS). In this work we present our data reduction techniques, and document the resulting mosaics, coverage maps, and catalogs in both IRAC passbands for the three easternmost UltraVISTA survey stripes, covering a combined area of about 0.66 square degrees, of which 0.45 square degrees have at least 20 hr integration time. SMUVS reaches point-source sensitivities of about 25.0 AB mag at both 3.6 and 4.5 microns with a significance of 4-sigma accounting for both survey sensitivity and source confusion. To this limit the SMUVS catalogs contain a total of about 350,000 sources, each of which is detected significantly in at least one IRAC band. Because of its uniform and high sensitivity, relatively large area coverage, and the wide array of ancillary data available in COSMOS, the SMUVS survey will be useful for a large number of cosmological investigations. We will make all images and catalogues described herein publicly available via the Spitzer Science Center.

Sgr A*, the supermassive black hole (SMBH) at the center of our Milky Way Galaxy, is known to be a variable source of X-ray, near-infrared (NIR), and submillimeter (submm) radiation and therefore a prime candidate to study the electromagnetic radiation generated by mass accretion flow onto a black hole and/or a related jet. Disentangling the power source and emission mechanisms of this variability is a central challenge to our understanding of accretion flows around SMBHs. Simultaneous multiwavelength observations of the flux variations and their time correlations can play an important role in obtaining a better understanding of possible emission mechanisms and their origin. This paper presents observations of two flares that both apparently violate the previously established patterns in the relative timing of submm/NIR/X-ray flares from Sgr A*. One of these events provides the first evidence of coeval structure between NIR and submm flux increases, while the second event is the first example of the sequence of submm/X-ray/NIR flux increases all occurring within ~1 hr. Each of these two events appears to upend assumptions that have been the basis of some analytic models of flaring in Sgr A*. However, it cannot be ruled out that these events, even though unusual, were just coincidental. These observations demonstrate that we do not fully understand the origin of the multiwavelength variability of Sgr A*, and show that there is a continued and important need for long-term, coordinated, and precise multiwavelength observations of Sgr A* to characterize the full range of variability behavior.

Galaxy mergers are key events in galaxy evolution, often causing massive starbursts and fueling active galactic nuclei (AGN). In these highly dynamic systems, it is not yet precisely known how much starbursts and AGN respectively contribute to the total luminosity, at what interaction stages they occur, and how long they persist. Here we estimate the fraction of the bolometric infrared (IR) luminosity that can be attributed to AGN by measuring and modeling the full ultraviolet to far-infrared spectral energy distributions (SEDs) in up to 33 broad bands for 24 merging galaxies with the Code for Investigating Galaxy Emission. In addition to a sample of 12 confirmed AGN in late-stage mergers, found in the $Infrared$ $Array$ $Satellite$ Revised Bright Galaxy Sample or Faint Source Catalog, our sample includes a comparison sample of 12 galaxy mergers from the $Spitzer$ Interacting Galaxies Survey, mostly early-stage. We perform identical SED modeling of simulated mergers to validate our methods, and we supplement the SED data with mid-IR spectra of diagnostic lines obtained with $Spitzer$ InfraRed Spectrograph. The estimated AGN contributions to the IR luminosities vary from system to system from 0% up to 91% but are significantly greater in the later-stage, more luminous mergers, consistent with what is known about galaxy evolution and AGN triggering.

We present spectroscopic confirmation of five galaxy clusters at $1.25 < \textit{z} < 1.5$, discovered in the $2500$ deg$^{2}$ South Pole Telescope Sunyaev-Zel'dovich (SPT-SZ) survey. These clusters, taken from a mass-limited sample with a nearly redshift independent selection function, have multi-wavelength follow-up imaging data from the X-ray to near-infrared, and currently form the most homogeneous massive high-redshift cluster sample known. We identify $44$ member galaxies, along with $25$ field galaxies, among the five clusters, and describe the full set of observations and data products from Magellan/LDSS3 multi-object spectroscopy of these cluster fields. We briefly describe the analysis pipeline, and present ensemble analyses of cluster member galaxies that demonstrate the reliability of the measured redshifts. We report $\textit{z} = 1.259, 1.288, 1.316, 1.401$ and $1.474$ for the five clusters from a combination of absorption-line (Ca II H$\&$K doublet - $3968,3934$ {\AA}) and emission-line ([OII] $3727,3729$ {\AA}) spectral features. Moreover, the calculated velocity dispersions yield dynamical cluster masses in good agreement with SZ masses for these clusters. We discuss the velocity and spatial distributions of passive and [OII]-emitting galaxies in these clusters, showing that they are consistent with velocity segregation and biases observed in lower redshift SPT clusters. We identify modest [OII] emission and pronounced CN and H$\delta$ absorption in a stacked spectrum of $28$ passive galaxies with Ca II H$\&$K-derived redshifts. This work increases the number of spectroscopically-confirmed SZ-selected galaxy clusters at $\textit{z} > 1.25$ from three to eight, further demonstrating the efficacy of SZ selection for the highest redshift massive clusters, and enabling detailed study of these systems.

We estimate total mass ($M_{500}$), intracluster medium (ICM) mass ($M_{\mathrm{ICM}}$) and stellar mass ($M_{\star}$) in a Sunyaev-Zel'dovich effect (SZE) selected sample of 91 galaxy clusters with masses $M_{500}\gtrsim2.5\times10^{14}M_{\odot}$ and redshift $0.2 < z < 1.25$ from the 2500 deg$^2$ South Pole Telescope SPT-SZ survey. The total masses $M_{500}$ are estimated from the SZE observable, the ICM masses $M_{\mathrm{ICM}}$ are obtained from the analysis of $Chandra$ X-ray observations, and the stellar masses $M_{\star}$ are derived by fitting spectral energy distribution templates to Dark Energy Survey (DES) $griz$ optical photometry and $WISE$ or $Spitzer$ near-infrared photometry. We study trends in the stellar mass, the ICM mass, the total baryonic mass and the cold baryonic fraction with cluster mass and redshift. We find significant departures from self-similarity in the mass scaling for all quantities, while the redshift trends are all statistically consistent with zero, indicating that the baryon content of clusters at fixed mass has changed remarkably little over the past $\approx9$ Gyr. We compare our results to the mean baryon fraction (and the stellar mass fraction) in the field, finding that these values lie above (below) those in cluster virial regions in all but the most massive clusters at low redshift. Using a simple model of the matter assembly of clusters from infalling groups with lower masses and from infalling material from the low density environment or field surrounding the parent halos, we show that the measured mass trends without strong redshift trends in the stellar mass scaling relation could be explained by a mass and redshift dependent fractional contribution from field material. Similar analyses of the ICM and baryon mass scaling relations provide evidence for the so-called "missing baryons" outside cluster virial regions.

SPHEREx is a proposed NASA MIDEX mission selected for Phase A study. SPHEREx would carry out the first all-sky spectral survey in the near infrared. At the end of its two-year mission, SPHEREx would obtain 0.75-to-5$\mu$m spectra of every 6.2 arcsec pixel on the sky, with spectral resolution R>35 and a 5-$\sigma$ sensitivity AB$>$19 per spectral/spatial resolution element. More details concerning SPHEREx are available at http://spherex.caltech.edu. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. Though these three themes are undoubtedly compelling, they are far from exhausting the scientific output of SPHEREx. Indeed, SPHEREx would create a unique all-sky spectral database including spectra of very large numbers of astronomical and solar system targets, including both extended and diffuse sources. These spectra would enable a wide variety of investigations, and the SPHEREx team is dedicated to making the data available to the community to enable these investigations, which we refer to as Legacy Science. To that end, we have sponsored two workshops for the general scientific community to identify the most interesting Legacy Science themes and to ensure that the SPHEREx data products are responsive to their needs. In February of 2016, some 50 scientists from all fields met in Pasadena to develop these themes and to understand their implications for the SPHEREx mission. The 2016 workshop highlighted many synergies between SPHEREx and other contemporaneous astronomical missions, facilities, and databases. Consequently, in January 2018 we convened a second workshop at the Center for Astrophysics in Cambridge to focus specifically on these synergies. This white paper reports on the results of the 2018 SPHEREx workshop.

Massive galaxy clusters are now found as early as 3 billion years after the Big Bang, containing stars that formed at even earlier epochs. The high-redshift progenitors of these galaxy clusters, termed 'protoclusters', are identified in cosmological simulations with the highest dark matter overdensities. While their observational signatures are less well defined compared to virialized clusters with a substantial hot intra-cluster medium (ICM), protoclusters are expected to contain extremely massive galaxies that can be observed as luminous starbursts. Recent claimed detections of protoclusters hosting such starbursts do not support the kind of rapid cluster core formation expected in simulations because these structures contain only a handful of starbursting galaxies spread throughout a broad structure, with poor evidence for eventual collapse into a protocluster. Here we report that the source SPT2349-56 consists of at least 14 gas-rich galaxies all lying at z = 4.31 based on sensitive observations of carbon monoxide and ionized carbon. We demonstrate that each of these galaxies is forming stars between 50 and 1000 times faster than our own Milky Way, and all are located within a projected region only $\sim$ 130 kiloparsecs in diameter. This galaxy surface density is more than 10 times the average blank field value (integrated over all redshifts) and $>$1000 times the average field volume density. The velocity dispersion ($\sim$ 410 km s$^{-1}$) of these galaxies and enormous gas and star formation densities suggest that this system represents a galaxy cluster core at an advanced stage of formation when the Universe was only 1.4 billion years old. A comparison with other known protoclusters at high redshifts shows that SPT2349-56 is a uniquely massive and dense system that could be building one of the most massive structures in the Universe today.

The oxygen-bearing molecular ions OH+, H2O+, and H3O+ are key species that probe the ionization rate of (partially) molecular gas that is ionized by X-rays and cosmic rays permeating the interstellar medium. We report Herschel far-infrared and submillimeter spectroscopic observations of OH+ in Mrk 231, showing both ground-state P-Cygni profiles, and excited line profiles with blueshifted absorption wings extending up to ~1000 km s^{-1}. In addition, OH+ probes an excited component peaking at central velocities, likely arising from the torus probed by the OH centimeter-wave megamaser. Four lines of H2O+ are also detected at systemic velocities, but H3O+ is undetected. Based on our earlier OH studies, we estimate an abundance ratio of OH/OH+~5-10 for the outflowing components and ~20 for the torus, and an OH+ abundance relative to H nuclei of ~>10^{-7}. We also find high OH+/H2O+ and OH+/H3O+ ratios, both are ~>4 in the torus and ~>10-20 in the outflowing gas components. Chemical models indicate that these high OH+ abundances relative to OH, H2O+, and H3O+ are characteristic of gas with a high ionization rate per unit density, \zeta/n_H~(1-5)x10^{-17} cm^3 s^{-1} and ~(1-2)x10^{-16} cm^3 s^{-1} for the above components, respectively, and an ionization rate of \zeta~(0.5-2)x10^{-12} s^{-1}. X-rays appear to be unable to explain the inferred ionization rate, and thus we suggest that low-energy (10-400 MeV) cosmic-rays are primarily responsible for the ionization with \dot{M}_{CR}~0.01 M_{sun} yr^{-1} and \dot{E}_{CR}~10^{44} erg s^{-1}, the latter corresponding to 1% of the AGN luminosity and similar to the energetics of the molecular outflow. We suggest that cosmic-rays accelerated in the forward shock associated with the molecular outflow are responsible for the ionization, as they diffuse through the outflowing molecular phase downstream.

The \emph{Spitzer} Matching Survey of the UltraVISTA ultra-deep Stripes (SMUVS) provides unparalleled depth at $3.6$ and $4.5$~$\mu$m over $\sim0.66$~deg$^2$ of the COSMOS field, allowing precise photometric determinations of redshift and stellar mass. From this unique dataset we can connect galaxy samples, selected by stellar mass, to their host dark matter halos for $1.5<z<5.0$, filling in a large hitherto unexplored region of the parameter space. To interpret the observed galaxy clustering we utilize a phenomenological halo model, combined with a novel method to account for uncertainties arising from the use of photometric redshifts. We find that the satellite fraction decreases with increasing redshift and that the clustering amplitude (e.g., comoving correlation length / large-scale bias) displays monotonic trends with redshift and stellar mass. Applying $\Lambda$CDM halo mass accretion histories and cumulative abundance arguments for the evolution of stellar mass content we propose pathways for the coevolution of dark matter and stellar mass assembly. Additionally, we are able to estimate that the halo mass at which the ratio of stellar to halo mass is maximized is $10^{12.5_{-0.08}^{+0.10}}$~M$_{\odot}$ at $z\sim2.5$. This peak halo mass is here inferred for the first time from stellar mass-selected clustering measurements at $z\gtrsim2$, and implies mild evolution of this quantity for $z\lesssim3$, consistent with constraints from abundance-matching techniques.

We present activity demographics and host-galaxy properties of infrared-selected galaxies in the local Universe, using the representative Star Formation Reference Survey (SFRS). Our classification scheme is based on a combination of optical emission-line diagrams (BPT) and IR-color diagnostics. Using the weights assigned to the SFRS galaxies based on its parent sample, a far-infrared-selected sample comprises 71\% H\,\textsc{ii} galaxies, 13\% Seyferts, 3\% Transition Objects (TOs), and 13\% Low-Ionization Nuclear Emission-Line Regions (LINERs). For the SFRS H\,\textsc{ii} galaxies we derive nuclear star-formation rates and gas-phase metallicities. We measure host-galaxy metallicities for all galaxies with available long-slit spectroscopy and abundance gradients for a subset of 12 face-on galaxies. The majority of H\,\textsc{ii} galaxies show a narrow range of metallicities, close to solar, and flat metallicity profiles. Based on their host-galaxy and nuclear properties, the dominant ionizing source in the far-infrared selected TOs is star-forming activity. LINERs are found mostly in massive hosts (median of $10^{10.5}$ M$_{\odot} $), median $L(60\mu m) = 10^{9}$ L$_{\odot}$, median dust temperatures of $ F60/F100 = 0.36 $, and median $L_{\textrm{H}\alpha}$ surface density of $ 10^{40.2} $ erg s$ ^{-1} $kpc$ ^{-2} $, indicating older stellar populations as their main ionizing source rather than AGN activity.

According to the current understanding of cosmic structure formation, the precursors of the most massive structures in the Universe began to form shortly after the Big Bang, in regions corresponding to the largest fluctuations in the cosmic density field. Observing these structures during their period of active growth and assembly - the first few hundred million years of the Universe - is challenging because it requires surveys that are sensitive enough to detect the distant galaxies that act as signposts for these structures and wide enough to capture the rarest objects. As a result, very few such objects have been detected so far. Here we report observations of a far-infrared-luminous object at redshift 6.900 (less than 800 Myr after the Big Bang) that was discovered in a wide-field survey. High-resolution imaging reveals this source to be a pair of extremely massive star-forming galaxies. The larger of these galaxies is forming stars at a rate of 2900 solar masses per year, contains 270 billion solar masses of gas and 2.5 billion solar masses of dust, and is more massive than any other known object at a redshift of more than 6. Its rapid star formation is probably triggered by its companion galaxy at a projected separation of just 8 kiloparsecs. This merging companion hosts 35 billion solar masses of stars and has a star-formation rate of 540 solar masses per year, but has an order of magnitude less gas and dust than its neighbor and physical conditions akin to those observed in lower-metallicity galaxies in the nearby Universe. These objects suggest the presence of a dark-matter halo with a mass of more than 400 billion solar masses, making it among the rarest dark-matter haloes that should exist in the Universe at this epoch.

Although the final observations of the Spitzer Warm Mission are currently scheduled for March 2019, it can continue operations through the end of the decade with no loss of photometric precision. As we will show, there is a strong science case for extending the current Warm Mission to December 2020. Spitzer has already made major impacts in the fields of exoplanets (including microlensing events), characterizing near Earth objects, enhancing our knowledge of nearby stars and brown dwarfs, understanding the properties and structure of our Milky Way galaxy, and deep wide-field extragalactic surveys to study galaxy birth and evolution. By extending Spitzer through 2020, it can continue to make ground-breaking discoveries in those fields, and provide crucial support to the NASA flagship missions JWST and WFIRST, as well as the upcoming TESS mission, and it will complement ground-based observations by LSST and the new large telescopes of the next decade. This scientific program addresses NASA's Science Mission Directive's objectives in astrophysics, which include discovering how the universe works, exploring how it began and evolved, and searching for life on planets around other stars.

Luminous radio sources are thought to lie in galaxy clusters or proto-clusters. The complete sample of 64 high-redshift 3C sources at $1<z<2.5$ has been mapped with the Spitzer Space Telescope. The IRAC $3.6\mu$m and $4.5\mu$m 5-$\sigma$ detection limit of $4\mu$Jy (22.4 AB mag) allows us to search for the brightest candidate cluster member galaxies associated with the 3C sources. To remove the contamination of foreground stars and galaxies along the lines of sight toward the 3C sources we apply color cuts: removed sources satisfy either the IRAC1/2 cut $[3.6] - [4.5] < -0.3$ or the Pan-STARRS/IRAC cut $i - [4.5] < 0.5$ if detected by Pan-STARRS. For both selection methods, about half of the 3C radio sources show significant overdensities ($>3\sigma$) within 30" ($\sim$250 kpc) projected distance from the radio source compared to the surrounding galaxy densities measured in the $50''-120''$ annulus. The Pan-STARRS/IRAC cut reveals higher average overdensities than the IRAC1/2 cut and also a decline of overdensities at $z>1.4$. To infer the nature of the cluster members, we rerun the analysis using a stronger IRAC1/2 cut $[3.6] - [4.5] < -0.1$ which removes $1<z<1.4$ passive ellipticals but not star-forming galaxies. For the strong cut, the overdensities, on average, completely disappear at $1<z<1.4$. We therefore suggest that the 4.5$\mu$m detected cluster member galaxies are mainly passive ellipticals.

We use observations made with the Giant Metrewave Radio Telescope (GMRT) to probe the neutral hydrogen (HI) gas content of field galaxies in the VIMOS VLT Deep Survey (VVDS) 14h field at $z \approx 0.32$. Because the HI emission from individual galaxies is too faint to detect at this redshift, we use an HI spectral stacking technique using the known optical positions and redshifts of the 165 galaxies in our sample to co-add their HI spectra and thus obtain the average HI mass of the galaxies. Stacked HI measurements of 165 galaxies show that 95 per cent of the neutral gas is found in blue, star-forming galaxies. Among these galaxies, those having lower stellar mass are more gas-rich than more massive ones. We apply a volume correction to our HI measurement to evaluate the HI gas density at $z \approx 0.32$ as $\Omega_{HI}=(0.50\pm0.18) \times 10^{-3}$ in units of the cosmic critical density. This value is in good agreement with previous results at z < 0.4, suggesting no evolution in the neutral hydrogen gas density over the last $\sim 4$ Gyr. However the $z \approx 0.32$ gas density is lower than that at $z \sim 5$ by at least a factor of two.

We study a large galaxy sample from the Spitzer Matching Survey of the UltraVISTA ultra-deep Stripes (SMUVS) to search for sources with enhanced 3.6 micron fluxes indicative of strong Halpha emission at z=3.9-4.9. We find that the percentage of "Halpha excess" sources reaches 37-40% for galaxies with stellar masses log10(M*/Msun) ~ 9-10, and decreases to <20% at log10(M*/Msun) ~ 10.7. At higher stellar masses, however, the trend reverses, although this is likely due to AGN contamination. We derive star formation rates (SFR) and specific SFR (sSFR) from the inferred Halpha equivalent widths (EW) of our "Halpha excess" galaxies. We show, for the first time, that the "Halpha excess" galaxies clearly have a bimodal distribution on the SFR-M* plane: they lie on the main sequence of star formation (with log10(sSFR/yr^{-1})<-8.05) or in a starburst cloud (with log10(sSFR/yr^{-1}) >-7.60). The latter contains ~15% of all the objects in our sample and accounts for >50% of the cosmic SFR density at z=3.9-4.9, for which we derive a robust lower limit of 0.066 Msun yr^{-1} Mpc^{-3}. Finally, we identify an unusual >50sigma overdensity of z=3.9-4.9 galaxies within a 0.20 x 0.20 sq. arcmin region. We conclude that the SMUVS unique combination of area and depth at mid-IR wavelengths provides an unprecedented level of statistics and dynamic range which are fundamental to reveal new aspects of galaxy evolution in the young Universe.

We present 190 galaxy cluster candidates (most at high redshift) based on galaxy overdensity measurements in the \spitzer/IRAC imaging of the fields surrounding 646 bent, double-lobed radio sources drawn from the Clusters Occupied by Bent Radio AGN (COBRA) Survey. The COBRA sources were chosen as objects in the VLA FIRST survey that lack optical counterparts in the Sloan Digital Sky Survey (SDSS) to a limit of $m_r=22$, making them likely to lie at high redshift. This is confirmed by our observations: the redshift distribution of COBRA sources with estimated redshifts peaks near $z=1$, and extends out to $z\approx3$. Cluster candidates were identified by comparing our target fields to a background field and searching for statistically significant ($\ge2\sigma$) excesses in the galaxy number counts surrounding the radio sources; 190 fields satisfy the $\ge2\sigma$ limit. We find that 530 fields (82.0\%) have a net positive excess of galaxies surrounding the radio source. Many of the fields with positive excesses but below the $2\sigma$ cutoff are likely to be galaxy groups. Forty-one COBRA sources are quasars with known spectroscopic redshifts, which may be tracers of some of the most distant clusters known.

We report the discovery and constrain the physical conditions of the interstellar medium of the highest-redshift millimeter-selected dusty star-forming galaxy (DSFG) to date, SPT-S J031132-5823.4 (hereafter SPT0311-58), at $z=6.900 +/- 0.002$. SPT0311-58 was discovered via its 1.4mm thermal dust continuum emission in the South Pole Telescope (SPT)-SZ survey. The spectroscopic redshift was determined through an ALMA 3mm frequency scan that detected CO(6-5), CO(7-6) and [CI](2-1), and subsequently confirmed by detections of CO(3-2) with ATCA and [CII] with APEX. We constrain the properties of the ISM in SPT0311-58 with a radiative transfer analysis of the dust continuum photometry and the CO and [CI] line emission. This allows us to determine the gas content without ad hoc assumptions about gas mass scaling factors. SPT0311-58 is extremely massive, with an intrinsic gas mass of $M_{\rm gas} = 3.3 \pm 1.9 \times10^{11}\,M_{\odot}$. Its large mass and intense star formation is very rare for a source well into the Epoch of Reionization.

This paper uses radial colour profiles to infer the distributions of dust, gas and star formation in z=0.4-1.4 star-forming main sequence galaxies. We start with the standard UVJ-based method to estimate dust extinction and specific star formation rate (sSFR). By replacing J with I band, a new calibration method suitable for use with ACS+WFC3 data is created (i.e. UVI diagram). Using a multi-wavelength multi-aperture photometry catalogue based on CANDELS, UVI colour profiles of 1328 galaxies are stacked in stellar mass and redshift bins. The resulting colour gradients, covering a radial range of 0.2--2.0 effective radii, increase strongly with galaxy mass and with global $A_V$. Colour gradient directions are nearly parallel to the Calzetti extinction vector, indicating that dust plays a more important role than stellar population variations. With our calibration, the resulting $A_V$ profiles fall much more slowly than stellar mass profiles over the measured radial range. sSFR gradients are nearly flat without central quenching signatures, except for $M_*>10^{10.5} M_{\odot}$, where central declines of 20--25 per cent are observed. Both sets of profiles agree well with previous radial sSFR and (continuum) $A_V$ measurements. They are also consistent with the sSFR profiles and, if assuming a radially constant gas-to-dust ratio, gas profiles in recent hydrodynamic models. We finally discuss the striking findings that SFR scales with stellar mass density in the inner parts of galaxies, and that dust content is high in the outer parts despite low stellar-mass surface densities there.

We present a 0.4-8$\mu$m multi-wavelength photometric catalog in the Extended Groth Strip (EGS) field. This catalog is built on the Hubble Space Telescope (HST) WFC3 and ACS data from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS), and it incorporates the existing HST data from the All-wavelength Extended Groth strip International Survey (AEGIS) and the 3D-HST program. The catalog is based on detections in the F160W band reaching a depth of F160W=26.62 AB (90% completeness, point-sources). It includes the photometry for 41457 objects over an area of $\approx 206$ arcmin$^2$ in the following bands: HST ACS F606W and F814W; HST WFC3 F125W, F140W and F160W; CFHT/Megacam $u^*$, $g'$, $r'$, $i'$ and $z'$; CFHT/WIRCAM $J$, $H$ and $K_\mathrm{S}$; Mayall/NEWFIRM $J1$, $J2$, $J3$, $H1$, $H2$, $K$; Spitzer IRAC $3.6\mu$m, $4.5\mu$m, $5.8\mu$m and $8.0\mu$m. We are also releasing value-added catalogs that provide robust photometric redshifts and stellar mass measurements. The catalogs are publicly available through the CANDELS repository.

We report on the energetics of molecular outflows in 14 local Ultraluminous Infrared Galaxies (ULIRGs) that show unambiguous outflow signatures (P-Cygni profiles or high-velocity absorption wings) in the far-infrared lines of OH measured with the Herschel/PACS spectrometer. Detection of both ground-state (at 119 and 79 um) and one or more radiatively-excited (at 65 and 84 um) lines allows us to model the nuclear gas (<~300 pc) as well as the more extended components using spherically symmetric radiative transfer models. The highest molecular outflow velocities are found in buried sources, in which slower but massive expansion of the nuclear gas is also observed. With the exception of a few outliers, the outflows have momentum fluxes of (2-5)xL_IR/c and mechanical luminosities of (0.1-0.3)% of L_IR. The moderate momentum boosts in these sources (<~3) suggest that the outflows are mostly momentum-driven by the combined effects of AGN and nuclear starbursts, as a result of radiation pressure, winds, and supernovae remnants. In some sources (~20%), however, powerful (10^{10.5-11} Lsun) AGN feedback and (partially) energy-conserving phases are required, with momentum boosts in the range 3-20. These outflows appear to be stochastic strong-AGN feedback events that occur throughout the merging process. In a few sources, the outflow activity in the innermost regions has subsided in the last ~1 Myr. While OH traces the molecular outflows at sub-kpc scales, comparison of the masses traced by OH with those previously inferred from tracers of more extended outflowing gas suggests that most mass is loaded (with loading factors of Mdot/SFR=1-10) from the central galactic cores (a few x 100 pc). Outflow depletion timescales are <10^8 yr, shorter than the gas consumption timescales by factors of 1.1-15, and are anti-correlated with the AGN luminosity.

We present the detailed characterization of two extremely red submillimeter galaxies (SMGs), ASXDF1100.053.1 and 231.1, with the Atacama Large Millimeter/submillimeter Array (ALMA) and the Jansky Very Large Array (VLA). These SMGs were selected originally using AzTEC at 1100 micron, and are observed by Herschel to be faint at 100--500 micron. Their (sub)millimeter colors are as red as -- or redder -- than known z>~5 SMGs; indeed, ASXDF1100.053.1 is redder than HFLS 3, which lies at z=6.3. They are also faint and red in the near-/mid-infrared: ~1 microJy at IRAC 4.5 micron and <0.2 microJy in the Ks filter. These SMGs are also faint in the radio waveband, where F_6GHz=4.5 microJy for ASXDF1100.053.1 and F_1.4GHz=28 microJy for ASXDF1100.231.1, suggestive of z=6.5^{+1.4}_{-1.1} and z=4.1^{+0.6}_{-0.7} for ASXDF1100.053.1 and 231.1, respectively. ASXDF1100.231.1 has a flux excess in the 3.6-micron filter, probably due to H$\alpha$ emission at z=4--5. Derived properties of ASXDF1100.053.1 for z=5.5--7.5 and 231.1 for z=3.5--5.5 are as follows: their infrared luminosities are [6.5-7.4]x10^{12} and [4.2-4.5]x10^{12} L_sun; their stellar masses are [0.9-2]x10^{11} and [0.4-3]x10^{10} M_sun; their circularized half-light radii in the ALMA maps are ~1 and <~0.2 kpc (~2--3 kpc for 90% of the total flux). Lastly, their surface infrared luminosity densities, Sigma_IR, are ~1x10^{12} and >~1.5x10^{13} L_sun kpc^{-2}, similar to values seen for local (U)LIRGs. These data suggest that ASXDF1100.053.1 and 231.1 are compact SMGs at z>~4 and can plausibly evolve into z>~3 compact quiescent galaxies.

We present a multi-wavelength photometric catalog in the COSMOS field as part of the observations by the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). The catalog is based on Hubble Space Telescope Wide Field Camera 3 (HST/WFC3) and Advanced Camera for Surveys (ACS) observations of the COSMOS field (centered at RA: $10^h00^m28^s$, Dec:$+02^{\circ}12^{\prime}21^{\prime\prime}$). The final catalog has 38671 sources with photometric data in forty two bands from UV to the infrared ($\rm \sim 0.3-8\,\mu m$). This includes broad-band photometry from the HST, CFHT, Subaru, VISTA and Spitzer Space Telescope in the visible, near infrared and infrared bands along with intermediate and narrow-band photometry from Subaru and medium band data from Mayall NEWFIRM. Source detection was conducted in the WFC3 F160W band (at $\rm 1.6\,\mu m$) and photometry is generated using the Template FITting algorithm. We further present a catalog of the physical properties of sources as identified in the HST F160W band and measured from the multi-band photometry by fitting the observed spectral energy distributions of sources against templates.

We describe a sub-galactic main sequence (SGMS) relating star formation rate surface density ($\Sigma_{\textrm{SFR}}$) and stellar mass density ($\Sigma_{\star}$) for distinct regions within star forming galaxies, including their nuclei. We use a sample of 246 nearby star-forming galaxies from the "Star Formation Reference Survey" and demonstrate that the SGMS holds down to $ \sim $1 kpc scales with a slope of $\alpha=0.91$ and a dispersion of 0.31 dex, similar to the well-known main sequence (MS) measured for globally integrated star formation rates (SFRs) and stellar masses. The SGMS slope depends on galaxy morphology, with late-type galaxies (Sc$-$Irr) having $\alpha = 0.97$ and early-type spirals (Sa$-$Sbc) having $\alpha = 0.81$. The SGMS constructed from sub-regions of individual galaxies has on average the same characteristics as the composite SGMS from all galaxies. The SGMS for galaxy nuclei shows a dispersion similar to that seen for other sub-regions. Sampling a limited range of SFR$-$M$_{\star} $ space may produce either sub-linearity or super-linearity of the SGMS slope. For nearly all galaxies, both SFR and stellar mass peak in the nucleus, indicating that circumnuclear clusters are among the most actively star-forming regions in the galaxy and the most massive. The nuclear SFR also correlates with total galaxy mass, forming a distinct sequence from the standard MS of star-formation. The nuclear main sequence will be useful for studying bulge growth and for characterizing feedback processes connecting AGN and star formation.

We present Chandra ACIS-S and ATCA radio continuum observations of the strongly lensed dusty, star-forming galaxy SPT-S J034640-5204.9 (hereafter SPT0346-52) at $z$ = 5.656. This galaxy has also been observed with ALMA, HST, Spitzer, Herschel, APEX, and the VLT. Previous observations indicate that if the infrared (IR) emission is driven by star formation, then the inferred lensing-corrected star formation rate ($\sim$ 4500 $M_{\sun}$ yr$^{-1}$) and star formation rate surface density $\Sigma_{\rm SFR}$ ($\sim$ 2000 $M_{\sun} {yr^{-1}} {kpc^{-2}}$) are both exceptionally high. It remained unclear from the previous data, however, whether a central active galactic nucleus (AGN) contributes appreciably to the IR luminosity. The {\it Chandra} upper limit shows that SPT0346-52 is consistent with being star-formation dominated in the X-ray, and any AGN contribution to the IR emission is negligible. The ATCA radio continuum upper limits are also consistent with the FIR-to-radio correlation for star-forming galaxies with no indication of an additional AGN contribution. The observed prodigious intrinsic IR luminosity of (3.6 $\pm$ 0.3) $\times$ 10$^{13}$ $L_{\sun}$ originates almost solely from vigorous star formation activity. With an intrinsic source size of 0.61 $\pm$ 0.03 kpc, SPT0346-52 is confirmed to have one of the highest $\Sigma_{SFR}$ of any known galaxy. This high $\Sigma_{SFR}$, which approaches the Eddington limit for a radiation pressure supported starburst, may be explained by a combination of very high star formation efficiency and gas fraction.

We present the results of SPT-GMOS, a spectroscopic survey with the Gemini Multi-Object Spectrograph (GMOS) on Gemini South. The targets of SPT-GMOS are galaxy clusters identified in the SPT-SZ survey, a millimeter-wave survey of 2500 squ. deg. of the southern sky using the South Pole Telescope (SPT). Multi-object spectroscopic observations of 62 SPT-selected galaxy clusters were performed between January 2011 and December 2015, yielding spectra with radial velocity measurements for 2595 sources. We identify 2243 of these sources as galaxies, and 352 as stars. Of the galaxies, we identify 1579 as members of SPT-SZ galaxy clusters. The primary goal of these observations was to obtain spectra of cluster member galaxies to estimate cluster redshifts and velocity dispersions. We describe the full spectroscopic dataset and resulting data products, including galaxy redshifts, cluster redshifts and velocity dispersions, and measurements of several well-known spectral indices for each galaxy: the equivalent width, W, of [O II] 3727,3729 and H-delta, and the 4000A break strength, D4000. We use the spectral indices to classify galaxies by spectral type (i.e., passive, post-starburst, star-forming), and we match the spectra against photometric catalogs to characterize spectroscopically-observed cluster members as a function of brightness (relative to m*). Finally, we report several new measurements of redshifts for ten bright, strongly-lensed background galaxies in the cores of eight galaxy clusters. Combining the SPT-GMOS dataset with previous spectroscopic follow-up of SPT-SZ galaxy clusters results in spectroscopic measurements for >100 clusters, or ~20% of the full SPT-SZ sample.

The deep, wide-area (~800-900 arcmin**2) near-infrared/WFC3/IR + Spitzer/IRAC observations over the CANDELS fields have been a remarkable resource for constraining the bright end of high redshift UV luminosity functions (LFs). However, the lack of HST 1.05-micron observations over the CANDELS fields has made it difficult to identify z~9-10 sources robustly, since such data are needed to confirm the presence of an abrupt Lyman break at 1.2 microns. We report here on the successful identification of many such z~9-10 sources from a new HST program (z9-CANDELS) that targets the highest-probability z~9-10 galaxy candidates with observations at 1.05 microns, to search for a robust Lyman-break at 1.2 microns. The potential z~9-10 candidates are preselected from the full HST, Spitzer/IRAC S-CANDELS observations, and the deepest-available ground-based optical+near-infrared observations. We identified 15 credible z~9-10 galaxies over the CANDELS fields. Nine of these galaxies lie at z~9 and 5 are new identifications. Our targeted follow-up strategy has proven to be very efficient in making use of scarce HST time to secure a reliable sample of z~9-10 galaxies. Through extensive simulations, we replicate the selection process for our sample (both the preselection and follow-up) and use it to improve current estimates for the volume density of bright z~9 and z~10 galaxies. The volume densities we find are 5(-2)(+3)x and 8(-3)(+9)x lower, respectively, than found at z~8. When compared with the best-fit evolution (i.e., dlog_{10} rho(UV)/dz=-0.29+/-0.02) in the UV luminosities densities from z~8 to z~4 integrated to 0.3L*(z=3) (-20 mag), these luminosity densities are 2.6(-0.9)(+1.5)x and 2.2(-1.1)(+2.0)x lower, respectively, than the extrapolated trends. Our new results are broadly consistent with the "accelerated evolution" scenario at z>8, as seen in many theoretical models.

We report the orbital distribution of the Trans-Neptunian objects (TNOs) discovered during the High Ecliptic Latitude (HiLat) extension of the Canada-France Ecliptic Plane Survey (CFEPS), conducted from June 2006 to July 2009. The HiLat component was designed to address one of the shortcomings of ecliptic surveys (like CFEPS), their lack of sensitivity to high-inclination objects. We searched 701~deg$^2$ of sky ranging from 12$^\circ$ to 85$^\circ$ ecliptic latitude and discovered \lKBO TNOs, with inclinations between 15$^\circ$ to 104$^\circ$. This survey places a very strong constraint on the inclination distribution of the hot component of the classical Kuiper Belt, ruling out any possibility of a large intrinsic fraction of highly inclined orbits. Using the parameterization of \citet{2001AJ....121.2804B}, the HiLat sample combined with CFEPS imposes a width $14^\circ \le \sigma \le 15.5^\circ$, with a best match for $\sigma = 14.5^\circ$. HiLat discovered the first retrograde TNO, 2008~KV$_{42}$, with an almost polar orbit with inclination 104$^\circ$, and (418993), a scattering object with perihelion in the region of Saturn's influence, with $a \sim 400$~AU and $i = 68^\circ$.

Stellar streams have become central to studies of the interaction histories of nearby galaxies. To characterize the most prominent parts of the stellar stream around the well-known nearby (d = 17 Mpc) edge-on disk galaxy NGC 5907, we have obtained and analyzed new, deep gri Subaru/Suprime-Cam and 3.6 micron Spitzer/Infrared Array Camera (IRAC) observations. Combining the near-infrared 3.6 micron data with visible-light images allows us to use a long wavelength baseline to estimate the metallicity and age of the stellar population along a ~60 kpc long segment of the stream. We have fitted the stellar spectral energy distribution (SED) with a single-burst stellar population synthesis model and we use it to distinguish between the proposed satellite accretion and minor/major merger formation models of the stellar stream around this galaxy. We conclude that a massive minor merger (stellar mass ratio of at least 1:8) can best account for the metallicity of -0.3 inferred along the brightest parts of the stream.

We present the first results from an extended survey of the Andromeda galaxy (M31) using 41.1 h of observations by Spitzer-IRAC at 3.6 and 4.5$\,\mu$m. This survey extends previous observations to the outer disc and halo, covering total lengths of 4$.\!\!^\circ$4 and 6$.\!\!^\circ$6 along the minor and major axes, respectively. We have produced surface brightness profiles by combining the integrated light from background-corrected maps with stellar counts from a new catalogue of point sources. Using auxiliary catalogues we have carried out a statistical analysis in colour-magnitude space to discriminate M31 objects from foreground Milky Way stars and background galaxies. The catalogue includes 426,529 sources, of which 66 per cent have been assigned probability values to identify M31 objects with magnitude depths of [3.6]$\,=\,$19.0$\,\pm\,$0.2, [4.5]$\,=\,$18.7$\,\pm\,$0.2. We discuss applications of our data for constraining the stellar mass and characterising point sources in the outer radii.

SPHEREx is a proposed SMEX mission selected for Phase A. SPHEREx will carry out the first all-sky spectral survey and provide for every 6.2" pixel a spectra between 0.75 and 4.18 $\mu$m [with R$\sim$41.4] and 4.18 and 5.00 $\mu$m [with R$\sim$135]. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. It is readily apparent, however, that many other questions in astrophysics and planetary sciences could be addressed with the SPHEREx data. The SPHEREx team convened a community workshop in February 2016, with the intent of enlisting the aid of a larger group of scientists in defining these questions. This paper summarizes the rich and varied menu of investigations that was laid out. It includes studies of the composition of main belt and Trojan/Greek asteroids; mapping the zodiacal light with unprecedented spatial and spectral resolution; identifying and studying very low-metallicity stars; improving stellar parameters in order to better characterize transiting exoplanets; studying aliphatic and aromatic carbon-bearing molecules in the interstellar medium; mapping star formation rates in nearby galaxies; determining the redshift of clusters of galaxies; identifying high redshift quasars over the full sky; and providing a NIR spectrum for most eROSITA X-ray sources. All of these investigations, and others not listed here, can be carried out with the nominal all-sky spectra to be produced by SPHEREx. In addition, the workshop defined enhanced data products and user tools which would facilitate some of these scientific studies. Finally, the workshop noted the high degrees of synergy between SPHEREx and a number of other current or forthcoming programs, including JWST, WFIRST, Euclid, GAIA, K2/Kepler, TESS, eROSITA and LSST.

We present galaxy stellar mass functions (GSMFs) at $z=$ 4-8 from a rest-frame ultraviolet (UV) selected sample of $\sim$4500 galaxies, found via photometric redshifts over an area of $\sim$280 arcmin$^2$ in the CANDELS/GOODS fields and the Hubble Ultra Deep Field. The deepest Spitzer/IRAC data yet-to-date and the relatively large volume allow us to place a better constraint at both the low- and high-mass ends of the GSMFs compared to previous space-based studies from pre-CANDELS observations. Supplemented by a stacking analysis, we find a linear correlation between the rest-frame UV absolute magnitude at 1500 \AA\ ($M_{\rm UV}$) and logarithmic stellar mass ($\log M_*$) that holds for galaxies with $\log(M_*/M_{\odot}) \lesssim 10$. We use simulations to validate our method of measuring the slope of the $\log M_*$-$M_{\rm UV}$ relation, finding that the bias is minimized with a hybrid technique combining photometry of individual bright galaxies with stacked photometry for faint galaxies. The resultant measured slopes do not significantly evolve over $z=$ 4-8, while the normalization of the trend exhibits a weak evolution toward lower masses at higher redshift. We combine the $\log M_*$-$M_{\rm UV}$ distribution with observed rest-frame UV luminosity functions at each redshift to derive the GSMFs, finding that the low-mass-end slope becomes steeper with increasing redshift from $\alpha=-1.55^{+0.08}_{-0.07}$ at $z=4$ to $\alpha=-2.25^{+0.72}_{-0.35}$ at $z=8$. The inferred stellar mass density, when integrated over $M_*=10^8$-$10^{13} M_{\odot}$, increases by a factor of $10^{+30}_{-2}$ between $z=7$ and $z=4$ and is in good agreement with the time integral of the cosmic star formation rate density.

[abridged] We quantify the evolution of the stellar mass functions of star-forming and quiescent galaxies as a function of morphology from $z\sim 3$ to the present. Our sample consists of ~50,000 galaxies in the CANDELS fields ($\sim880$ $arcmin^2$), which we divide into four main morphological types, i.e. pure bulge dominated systems, pure spiral disk dominated, intermediate 2-component bulge+disk systems and irregular disturbed galaxies. Our main results are: Star-formation: At $z\sim 2$, 80\% of the stellar mass density of star-forming galaxies is in irregular systems. However, by $z\sim 0.5$, irregular objects only dominate at stellar masses below $10^9M\odot$. A majority of the star-forming irregulars present at $z\sim 2$ undergo a gradual transformation from disturbed to normal spiral disk morphologies by $z\sim 1$ without significant interruption to their star-formation. Rejuvenation after a quenching event does not seem to be common except perhaps for the most massive objects. Quenching: We confirm that galaxies reaching a stellar mass of $M_*\sim10^{10.8}M_\odot$ ($M^*$) tend to quench. Also, quenching implies the presence of a bulge: the abundance of massive red disks is negligible at all redshifts over 2~dex in stellar mass. However the dominant quenching mechanism evolves. At $z>2$, the SMF of quiescent galaxies above $M^*$ is dominated by compact spheroids. Quenching at this early epoch destroys the disk and produces a compact remnant unless the star-forming progenitors at even higher redshifts are significantly more dense. At $1<z<2$, the majority of newly quenched galaxies are disks with a significant central bulge. This suggests that mass-quenching at this epoch starts from the inner parts and preserves the disk. At $z<1$, the high mass end of the passive SMF is globally in place and the evolution mostly happens at stellar masses below $10^{10}M_\odot$.

We report on the identification of the old stellar population galaxy candidates at z > 5. We developed a new infrared color selection scheme to isolate galaxies with the strong Balmer breaks at z > 5, and applied it to the ultra-deep and wide infrared survey data from the Spitzer Extended Deep Survey (SEDS) and the UKIRT Infrared Deep Sky Survey. The eight objects satisfying K - [3.6] > 1.3 and K - [3.6] > 2.4 ([3.6] - [4.5]) + 0.6 are selected in the 0.34 deg^2 SEDS Ultra Deep Survey field. Rich multi-wavelength imaging data from optical to far-infrared are also used to reject blending sources and strong nebular line emitters, and we finally obtained the three most likely evolved galaxies at z > 5. Their stacked spectral energy distribution is fitted well with the old stellar population template with M_{*} = (7.5+-1.5) x 10^{10} Msun, star formation rate = 0.9 +- 0.2 Msun yr^{-1}, dust A_V < 1, and age = 0.7+-0.4 Gyr at z = 5.7+-0.6, where the dusty star-forming galaxies at z ~ 2.8 are disfavored because of the faintness in the 24um. The stellar mass density of these evolved galaxy candidates, (6+-4) x 10^4 Msun Mpc^{-3}, is much lower than that of star-forming galaxies, but the non-zero fraction suggests that initial star-formation and quenching have been completed by z ~ 6.

(abridged) We present cosmological constraints obtained from galaxy clusters identified by their Sunyaev-Zel'dovich effect signature in the 2500 square degree South Pole Telescope Sunyaev Zel'dovich survey. We consider the 377 cluster candidates identified at z>0.25 with a detection significance greater than five, corresponding to the 95% purity threshold for the survey. We compute constraints on cosmological models using the measured cluster abundance as a function of mass and redshift. We include additional constraints from multi-wavelength observations, including Chandra X-ray data for 82 clusters and a weak lensing-based prior on the normalization of the mass-observable scaling relations. Assuming a LCDM cosmology, where the species-summed neutrino mass has the minimum allowed value (mnu = 0.06 eV) from neutrino oscillation experiments, we combine the cluster data with a prior on H0 and find sigma_8 = 0.797+-0.031 and Omega_m = 0.289+-0.042, with the parameter combination sigma_8(Omega_m/0.27)^0.3 = 0.784+-0.039. These results are in good agreement with constraints from the CMB from SPT, WMAP, and Planck, as well as with constraints from other cluster datasets. Adding mnu as a free parameter, we find mnu = 0.14+-0.08 eV when combining the SPT cluster data with Planck CMB data and BAO data, consistent with the minimum allowed value. Finally, we consider a cosmology where mnu and N_eff are fixed to the LCDM values, but the dark energy equation of state parameter w is free. Using the SPT cluster data in combination with an H0 prior, we measure w = -1.28+-0.31, a constraint consistent with the LCDM cosmological model and derived from the combination of growth of structure and geometry. When combined with primarily geometrical constraints from Planck CMB, H0, BAO and SNe, adding the SPT cluster data improves the w constraint from the geometrical data alone by 14%, to w = -1.023+-0.042.

We use the Atacama Large Millimeter/submillimeter Array (ALMA) in Cycle 1 to determine spectroscopic redshifts of high-redshift dusty star-forming galaxies (DSFGs) selected by their 1.4mm continuum emission in the South Pole Telescope (SPT) survey. We present ALMA 3mm spectral scans between 84-114GHz for 15 galaxies and targeted ALMA 1mm observations for an additional eight sources. Our observations yield 30 new line detections from CO, [CI] , [NII] , H_2O and NH_3. We further present APEX [CII] and CO mid-J observations for seven sources for which only a single line was detected in spectral-scan data from ALMA Cycle 0 or Cycle 1. We combine the new observations with previously published and new mm/submm line and photometric data of the SPT-selected DSFGs to study their redshift distribution. The combined data yield 39 spectroscopic redshifts from molecular lines, a success rate of >85%. Our sample represents the largest data set of its kind today and has the highest spectroscopic completeness among all redshift surveys of high-z DSFGs. The median of the redshift distribution is z=3.9+/-0.4, and the highest-redshift source in our sample is at z=5.8. We discuss how the selection of our sources affects the redshift distribution, focusing on source brightness, selection wavelength, and strong gravitational lensing. We correct for the effect of gravitational lensing and find the redshift distribution for 1.4mm-selected sources with a median redshift of z=3.1+/-0.3. Comparing to redshift distributions selected at shorter wavelengths from the literature, we show that selection wavelength affects the shape of the redshift distribution.

We present Hubble WFC3/IR slitless grism spectra of a remarkably bright $z\gtrsim10$ galaxy candidate, GN-z11, identified initially from CANDELS/GOODS-N imaging data. A significant spectroscopic continuum break is detected at $\lambda=1.47\pm0.01~\mu$m. The new grism data, combined with the photometric data, rule out all plausible lower redshift solutions for this source. The only viable solution is that this continuum break is the Ly$\alpha$ break redshifted to ${z_\mathrm{grism}=11.09^{+0.08}_{-0.12}}$, just $\sim$400 Myr after the Big Bang. This observation extends the current spectroscopic frontier by 150 Myr to well before the Planck (instantaneous) cosmic reionization peak at z~8.8, demonstrating that galaxy build-up was well underway early in the reionization epoch at z>10. GN-z11 is remarkably and unexpectedly luminous for a galaxy at such an early time: its UV luminosity is 3x larger than L* measured at z~6-8. The Spitzer IRAC detections up to 4.5 $\mu$m of this galaxy are consistent with a stellar mass of ${\sim10^{9}~M_\odot}$. This spectroscopic redshift measurement suggests that the James Webb Space Telescope (JWST) will be able to similarly and easily confirm such sources at z>10 and characterize their physical properties through detailed spectroscopy. Furthermore, WFIRST, with its wide-field near-IR imaging, would find large numbers of similar galaxies and contribute greatly to JWST's spectroscopy, if it is launched early enough to overlap with JWST.

The Decadal IRAC Bootes Survey is a mid-IR variability survey of the ~9 sq. deg. of the NDWFS Bootes Field and extends the time baseline of its predecessor, the Spitzer Deep, Wide-Field Survey (SDWFS), from 4 to 10 years. The Spitzer Space Telescope visited the field five times between 2004 and 2014 at 3.6 and 4.5 microns. We provide the difference image analysis photometry for a half a million mostly extragalactic sources. In mid-IR color-color plane, sources with quasar colors constitute the largest variability class (75%), 16% of the variable objects have stellar colors and the remaining 9% have the colors of galaxies. Adding the fifth epoch doubles the number of variable active galactic nuclei (AGNs) for the same false positive rates as in SDWFS, or increases the number of sources by 20% while decreasing the false positive rates by factors of 2-3 for the same variability amplitude. We quantify the ensemble mid-IR variability of ~1500 spectroscopically confirmed AGNs using single power-law structure functions (SFs), which we find to be steeper (index $\gamma=0.45$) than in the optical ($\gamma=0.3$), leading to much lower amplitudes at short time-lags. This provides evidence for large emission regions, smoothing out any fast UV/optical variations, as the origin of infrared quasar variability. The mid-IR AGN SF slope $\gamma$ seems to be uncorrelated with both the luminosity and rest-frame wavelength, while the amplitude shows an anti-correlation with the luminosity and a correlation with the rest-frame wavelength.

(abridged) The majority of star-forming galaxies follow a simple empirical correlation in the star formation rate (SFR) versus stellar mass ($M_*$) plane, usually referred to as the star formation Main Sequence (MS). Here we combine a set of hydro-dynamical simulations of interacting galactic disks with state-of-the-art radiative transfer codes to analyze how the evolution of mergers is reflected upon the properties of the MS. We present \textsc{Chiburst}, a Markov Chain Monte Carlo (MCMC) Spectral Energy Distribution (SED) code that fits the multi-wavelength, broad-band photometry of galaxies and derives stellar masses, star formation rates, and geometrical properties of the dust distribution. We apply this tool to the SEDs of simulated mergers and compare the derived results with the reference output from the simulations. Our results indicate that changes in the SEDs of mergers as they approach coalescence and depart from the MS are related to an evolution of dust geometry in scales larger than a few hundred parsecs. This is reflected in a correlation between the specific star formation rate (sSFR), and the compactness parameter $\mathcal{C}$, that parametrizes this geometry and hence the evolution of dust temperature ($T_{\rm{dust}}$) with time. As mergers approach coalescence, they depart from the MS and increase their compactness, which implies that moderate outliers of the MS are consistent with late-type mergers. By further applying our method to real observations of Luminous Infrared Galaxies (LIRGs), we show that the merger scenario is unable to explain these extreme outliers of the MS. Only by significantly increasing the gas fraction in the simulations are we able to reproduce the SEDs of LIRGs.

We introduce a new color-selection technique to identify high-redshift, massive galaxies that are systematically missed by Lyman-break selection. The new selection is based on the H_{160} and IRAC 4.5um bands, specifically H - [4.5] > 2.25 mag. These galaxies, dubbed "HIEROs", include two major populations that can be separated with an additional J - H color. The populations are massive and dusty star-forming galaxies at z > 3 (JH-blue) and extremely dusty galaxies at z < 3 (JH-red). The 350 arcmin^2 of the GOODS-N and GOODS-S fields with the deepest HST/WFC3 and IRAC data contain 285 HIEROs down to [4.5] < 24 mag. We focus here primarily on JH-blue (z > 3) HIEROs, which have a median photometric redshift z ~4.4 and stellar massM_{*}~10^{10.6} Msun, and are much fainter in the rest-frame UV than similarly massive Lyman-break galaxies (LBGs). Their star formation rates (SFRs) reaches ~240 Msun yr^{-1} leading to a specific SFR, sSFR ~4.2 Gyr^{-1}, suggesting that the sSFRs for massive galaxies continue to grow at z > 2 but at a lower growth rate than from z=0 to z=2. With a median half-light radius of 2 kpc, including ~20% as compact as quiescent galaxies at similar redshifts, JH-blue HIEROs represent perfect star-forming progenitors of the most massive (M_{*} > 10^{11.2} Msun) compact quiescent galaxies at z ~ 3 and have the right number density. HIEROs make up ~60% of all galaxies with M_{*} > 10^{10.5} Msun identified at z > 3 from their photometric redshifts. This is five times more than LBGs with nearly no overlap between the two populations. While HIEROs make up 15-25% of the total SFR density at z ~ 4-5, they completely dominate the SFR density taking place in M_{*} >10^{10.5} Msun galaxies, and are therefore crucial to understanding the very early phase of massive galaxy formation.

Spatially resolved kinematics have been used to determine the dynamical status of star-forming galaxies with ambiguous morphologies, and constrain the importance of galaxy interactions during the assembly of galaxies. However, measuring the importance of interactions or galaxy merger rates requires knowledge of the systematics in kinematic diagnostics and the visible time with merger indicators. We analyze the dynamics of star-forming gas in a set of binary merger hydrodynamic simulations with stellar mass ratios of 1:1 and 1:4. We find that the evolution of kinematic asymmetries traced by star-forming gas mirrors morphological asymmetries derived from mock optical images, in which both merger indicators show the largest deviation from isolated disks during strong interaction phases. Based on a series of simulations with various initial disk orientations, orbital parameters, gas fractions, and mass ratios, we find that the merger signatures are visible for ~0.2-0.4 Gyr with kinematic merger indicators but can be approximately twice as long for equal-mass mergers of massive gas-rich disk galaxies designed to be analogs of z~2-3 submillimeter galaxies. Merger signatures are most apparent after the second passage and before the black holes coalescence, but in some cases they persist up to several hundred Myr after coalescence. About 20-60% of the simulated galaxies are not identified as mergers during the strong interaction phase, implying that galaxies undergoing violent merging process do not necessarily exhibit highly asymmetric kinematics in their star-forming gas. The lack of identifiable merger signatures in this population can lead to an underestimation of merger abundances in star-forming galaxies, and including them in samples of star-forming disks may bias the measurements of disk properties such as intrinsic velocity dispersion.

Recent observations have shown that the characteristic luminosity of the rest-frame ultraviolet (UV) luminosity function does not significantly evolve at 4 < z < 7 and is approximately M*_UV ~ -21. We investigate this apparent non-evolution by examining a sample of 178 bright, M_UV < -21 galaxies at z=4 to 7, analyzing their stellar populations and host halo masses. Including deep Spitzer/IRAC imaging to constrain the rest-frame optical light, we find that M*_UV galaxies at z=4-7 have similar stellar masses of log(M/Msol)=9.6-9.9 and are thus relatively massive for these high redshifts. However, bright galaxies at z=4-7 are less massive and have younger inferred ages than similarly bright galaxies at z=2-3, even though the two populations have similar star formation rates and levels of dust attenuation. We match the abundances of these bright z=4-7 galaxies to halo mass functions from the Bolshoi Lambda-CDM simulation to estimate the halo masses. We find that the typical halo masses in ~M*_UV galaxies decrease from log(M_h/Msol)=11.9 at z=4 to log(M_h/Msol)=11.4 at z=7. Thus, although we are studying galaxies at a similar mass across multiple redshifts, these galaxies live in lower mass halos at higher redshift. The stellar baryon fraction in units of the cosmic mean Omega_b/Omega_m rises from 5.1% at z=4 to 11.7% at z=7; this evolution is significant at the ~3-sigma level. This rise does not agree with simple expectations of how galaxies grow, and implies that some effect, perhaps a diminishing efficiency of feedback, is allowing a higher fraction of available baryons to be converted into stars at high redshifts.

We probe star formation in the environments of massive $\sim10^{13}\,M_{\odot}$ dark matter halos at redshifts of $z$$\sim$$1$. This star formation is linked to a sub-millimetre clustering signal which we detect in maps of the Planck High Frequency Instrument that are stacked at the positions of a sample of high-redshift ($z$$>$$2$) strongly-lensed dusty star-forming galaxies (DSFGs) selected from the South Pole Telescope (SPT) 2500 deg$^2$ survey. The clustering signal has sub-millimetre colours which are consistent with the mean redshift of the foreground lensing halos ($z$$\sim$$1$). We report a mean excess of star formation rate (SFR) compared to the field, of $(2700\pm700)\,M_{\odot}\,{yr}^{-1}$ from all galaxies contributing to this clustering signal within a radius of 3.5' from the SPT DSFGs. The magnitude of the Planck excess is in broad agreement with predictions of a current model of the cosmic infrared background. The model predicts that 80$\%$ of the excess emission measured by Planck originates from galaxies lying in the neighbouring halos of the lensing halo. Using Herschel maps of the same fields, we find a clear excess, relative to the field, of individual sources which contribute to the Planck excess. The mean excess SFR compared to the field is measured to be ($370\pm40)$$\,M_{\odot}\,{yr}^{-1}$ per resolved, clustered source. Our findings suggest that the environments around these massive $z$$\sim$$1$ lensing halos host intense star formation out to about $2\,$Mpc. The flux enhancement due to clustering should also be considered when measuring flux densities of galaxies in Planck data.

We study the stellar, Brightest Cluster Galaxy (BCG) and intracluster medium (ICM) masses of 14 South Pole Telescope (SPT) selected galaxy clusters with median redshift $z=0.9$ and median mass $M_{500}=6\times10^{14}M_{\odot}$. We estimate stellar masses for each cluster and BCG using six photometric bands spanning the range from the ultraviolet to the near-infrared observed with the VLT, HST and Spitzer. The ICM masses are derived from Chandra and XMM-Newton X-ray observations, and the virial masses are derived from the SPT Sunyaev-Zel'dovich Effect signature. At $z=0.9$ the BCG mass $M_{\star}^{\textrm{BCG}}$ constitutes $0.12\pm0.01$% of the halo mass for a $6\times10^{14}M_{\odot}$ cluster, and this fraction falls as $M_{500}^{-0.58\pm0.07}$. The cluster stellar mass function has a characteristic mass $M_{0}=10^{11.0\pm0.1}M_{\odot}$, and the number of galaxies per unit mass in clusters is larger than in the field by a factor $1.65\pm0.2$. Both results are consistent with measurements on group scales and at lower redshift. We combine our SPT sample with previously published samples at low redshift that we correct to a common initial mass function and for systematic differences in virial masses. We then explore mass and redshift trends in the stellar fraction (fstar), the ICM fraction (fICM), the cold baryon fraction (fc) and the baryon fraction (fb). At a pivot mass of $6\times10^{14}M_{\odot}$ and redshift $z=0.9$, the characteristic values are fstar=$1.1\pm0.1$%, fICM=$9.6\pm0.5$%, fc=$10.4\pm1.2$% and fb=$10.7\pm0.6$%. These fractions all vary with cluster mass at high significance, indicating that higher mass clusters have lower fstar and fc and higher fICM and fb. When accounting for a 15% systematic virial mass uncertainty, there is no statistically significant redshift trend at fixed mass in these baryon fractions. (abridged)

We present Spitzer/Infrared Spectrograph (IRS) 5-21 micron spectroscopic maps towards 12 regions in the Andromeda galaxy (M31). These regions include the nucleus, bulge, an active region in the star-forming ring, and 9 other regions chosen to cover a range of mid-to-far-infrared colours. In line with previous results, PAH feature ratios (6.2 micron and 7.7 micron features compared to the 11.2 micron feature) measured from our extracted M31 spectra, except the nucleus, strongly correlate. The equivalent widths of the main PAH features, as a function of metallicity and radiation hardness, are consistent with those observed for other nearby spiral and starburst galaxies. Reprocessed data from the ISOCAM instrument on the Infrared Space Observatory agree with the IRS data; early reports of suppressed 6-8 micron features and enhanced 11.3 micron feature intensity and FWHM apparently resulted from background-subtraction problems. The nucleus does not show any PAH emission but does show strong silicate emission at 9.7 micron. Furthermore, different spectral features (11.3 micron PAH emission, silicate emission and [NeIII] 15.5 micron line emission) have distinct spatial distributions in the nuclear region: the silicate emission is strongest towards the stellar nucleus, while the PAH emission peaks 15 arcsec north of the nucleus. The PAH feature ratios at this position are atypical with strong emission at 11.2 microns and 15-20 microns but weak emission at 6--8 microns. The nucleus itself is dominated by stellar light giving rise to a strong blue continuum and silicate emission.

The spatial fluctuations of the extragalactic background light trace the total emission from all stars and galaxies in the Universe. A multi-wavelength study can be used to measure the integrated emission from first galaxies during reionization when the Universe was about 500 million years old. Here we report arcminute-scale spatial fluctuations in one of the deepest sky surveys with the Hubble Space Telescope in five wavebands between 0.6 and 1.6 $\mu$m. We model-fit the angular power spectra of intensity fluctuation measurements to find the ultraviolet luminosity density of galaxies at $z$ > 8 to be $\log \rho_{\rm UV} = 27.4^{+0.2}_{-1.2}$ erg s$^{-1}$ Hz$^{-1}$ Mpc$^{-3}$ $(1\sigma)$. This level of integrated light emission allows for a significant surface density of fainter primeval galaxies that are below the point source detection level in current surveys.

To understand cosmic mass assembly in the Universe at early epochs, we primarily rely on measurements of stellar mass and star formation rate of distant galaxies. In this paper, we present stellar masses and star formation rates of six high-redshift ($2.8\leq z \leq 5.7$) dusty, star-forming galaxies (DSFGs) that are strongly gravitationally lensed by foreground galaxies. These sources were first discovered by the South Pole Telescope (SPT) at millimeter wavelengths and all have spectroscopic redshifts and robust lens models derived from ALMA observations. We have conducted follow-up observations, obtaining multi-wavelength imaging data, using {\it HST}, {\it Spitzer}, {\it Herschel} and the Atacama Pathfinder EXperiment (APEX). We use the high-resolution {\it HST}/WFC3 images to disentangle the background source from the foreground lens in {\it Spitzer}/IRAC data. The detections and upper limits provide important constraints on the spectral energy distributions (SEDs) for these DSFGs, yielding stellar masses, IR luminosities, and star formation rates (SFRs). The SED fits of six SPT sources show that the intrinsic stellar masses span a range more than one order of magnitude with a median value $\sim$ 5 $\times 10^{10}M_{\Sun}$. The intrinsic IR luminosities range from 4$\times 10^{12}L_{\Sun}$ to 4$\times 10^{13}L_{\Sun}$. They all have prodigious intrinsic star formation rates of 510 to 4800 $M_{\Sun} {\rm yr}^{-1}$. Compared to the star-forming main sequence (MS), these six DSFGs have specific SFRs that all lie above the MS, including two galaxies that are a factor of 10 higher than the MS. Our results suggest that we are witnessing the ongoing strong starburst events which may be driven by major mergers.

We present a multi-wavelength study of 90 brightest cluster galaxies (BCGs) in a sample of galaxy clusters selected via the Sunyaev Zel'dovich effect by the South Pole Telescope, utilizing data from various ground- and space-based facilities. We infer the star formation rate (SFR) for the BCG in each cluster, based on the UV and IR continuum luminosity, as well as the [O II] emission line luminosity in cases where spectroscopy is available, finding 7 systems with SFR > 100 Msun/yr. We find that the BCG SFR exceeds 10 Msun/yr in 31 of 90 (34%) cases at 0.25 < z < 1.25, compared to ~1-5% at z ~ 0 from the literature. At z > 1, this fraction increases to 92(+6)(-31)%, implying a steady decrease in the BCG SFR over the past ~9 Gyr. At low-z, we find that the specific star formation rate in BCGs is declining more slowly with time than for field or cluster galaxies, most likely due to the replenishing fuel from the cooling ICM in relaxed, cool core clusters. At z > 0.6, the correlation between cluster central entropy and BCG star formation - which is well established at z ~ 0 - is not present. Instead, we find that the most star-forming BCGs at high-z are found in the cores of dynamically unrelaxed clusters. We investigate the rest-frame near-UV morphology of a subsample of the most star-forming BCGs using data from the Hubble Space Telescope, finding complex, highly asymmetric UV morphologies on scales as large as ~50-60 kpc. The high fraction of star-forming BCGs hosted in unrelaxed, non-cool core clusters at early times suggests that the dominant mode of fueling star formation in BCGs may have recently transitioned from galaxy-galaxy interactions to ICM cooling.

Clusters of galaxies are expected to gravitationally lens the cosmic microwave background (CMB) and thereby generate a distinct signal in the CMB on arcminute scales. Measurements of this effect can be used to constrain the masses of galaxy clusters with CMB data alone. Here we present a measurement of lensing of the CMB by galaxy clusters using data from the South Pole Telescope (SPT). We develop a maximum likelihood approach to extract the CMB cluster lensing signal and validate the method on mock data. We quantify the effects on our analysis of several potential sources of systematic error and find that they generally act to reduce the best-fit cluster mass. It is estimated that this bias to lower cluster mass is roughly $0.85\sigma$ in units of the statistical error bar, although this estimate should be viewed as an upper limit. We apply our maximum likelihood technique to 513 clusters selected via their SZ signatures in SPT data, and rule out the null hypothesis of no lensing at $3.1\sigma$. The lensing-derived mass estimate for the full cluster sample is consistent with that inferred from the SZ flux: $M_{200,\mathrm{lens}} = 0.83_{-0.37}^{+0.38}\, M_{200,\mathrm{SZ}}$ (68% C.L., statistical error only).

We present a robust measurement and analysis of the rest-frame ultraviolet (UV) luminosity function at z=4-8. We use deep Hubble Space Telescope imaging over the CANDELS/GOODS fields, the Hubble Ultra Deep Field and the Year 1 Hubble Frontier Field deep parallel observations. These surveys provides an effective volume of 0.6-1.2 x 10^6 Mpc^3 over this epoch, allowing us to perform a robust search for faint (M_UV=-18) and bright (M_UV < -21) galaxies. We select candidate galaxies using a well-tested photometric redshift technique with careful screening of contaminants, finding a sample of 7446 galaxies at 3.5<z<8.5, with >1000 galaxies at z~6-8. We measure the luminosity function using a Markov Chain Monte Carlo analysis to measure robust uncertainties. At the faint end our results agree with previous studies, yet we find a higher abundance of UV-bright galaxies at z>6, with M* ~ -21 at z>5, different than that inferred based on previous trends at lower redshift. At z=8, a single power-law provides an equally good fit to the UV luminosity function, while at z=6 and 7, an exponential cutoff at the bright-end is moderately preferred. We compare to semi-analytical models, and find that the lack of evolution in M* is consistent with models where the impact of dust attenuation on the bright-end of the luminosity function decreases at higher redshift. We measure the evolution of the cosmic star-formation rate density, correcting for dust attenuation, and find that it declines as (1+z)^(-4.3 +/- 0.5) at z>4, consistent with observations at z>9. Our observations are consistent with a reionization history that starts at z>10, completes at z>6, and reaches a midpoint (x_HII = 0.5) at 6.7<z<9.4. Finally, our observations predict that the abundance of bright z=9 galaxies is likely higher than previous constraints, though consistent with recent estimates of bright z~10 galaxies. [abridged]

The Spitzer-Cosmic Assembly Deep Near-Infrared Extragalactic Legacy Survey (S-CANDELS; PI G. Fazio) is a Cycle 8 Exploration Program designed to detect galaxies at very high redshifts (z > 5). To mitigate the effects of cosmic variance and also to take advantage of deep coextensive coverage in multiple bands by the Hubble Space Telescope Multi-Cycle Treasury Program CANDELS, S-CANDELS was carried out within five widely separated extragalactic fields: the UKIDSS Ultra-Deep Survey, the Extended Chandra Deep Field South, COSMOS, the HST Deep Field North, and the Extended Groth Strip. S-CANDELS builds upon the existing coverage of these fields from the Spitzer Extended Deep Survey (SEDS) by increasing the integration time from 12 hours to a total of 50 hours but within a smaller area, 0.16 square degrees. The additional depth significantly increases the survey completeness at faint magnitudes. This paper describes the S-CANDELS survey design, processing, and publicly-available data products. We present IRAC dual-band 3.6+4.5 micron catalogs reaching to a depth of 26.5 AB mag. Deep IRAC counts for the roughly 135,000 galaxies detected by S-CANDELS are consistent with models based on known galaxy populations. The increase in depth beyond earlier Spitzer/IRAC surveys does not reveal a significant additional contribution from discrete sources to the diffuse Cosmic Infrared Background (CIB). Thus it remains true that only roughly half of the estimated CIB flux from COBE/DIRBE is resolved.

(Abridged) We use 95, 150, and 220GHz observations from the SPT to examine the SZE signatures of a sample of 46 X-ray selected groups and clusters drawn from ~6 deg^2 of the XMM-BCS. These systems extend to redshift z=1.02, have characteristic masses ~3x lower than clusters detected directly in the SPT data and probe the SZE signal to the lowest X-ray luminosities (>10^42 erg s^-1) yet. We develop an analysis tool that combines the SZE information for the full ensemble of X-ray-selected clusters. Using X-ray luminosity as a mass proxy, we extract selection-bias corrected constraints on the SZE significance- and Y_500-mass relations. The SZE significance- mass relation is in good agreement with an extrapolation of the relation obtained from high mass clusters. However, the fit to the Y_500-mass relation at low masses, while in good agreement with the extrapolation from high mass SPT clusters, is in tension at 2.8 sigma with the constraints from the Planck sample. We examine the tension with the Planck relation, discussing sample differences and biases that could contribute. We also present an analysis of the radio galaxy point source population in this ensemble of X-ray selected systems. We find 18 of our systems have 843 MHz SUMSS sources within 2 arcmin of the X-ray centre, and three of these are also detected at significance >4 by SPT. Of these three, two are associated with the group brightest cluster galaxies, and the third is likely an unassociated quasar candidate. We examine the impact of these point sources on our SZE scaling relation analyses and find no evidence of biases. We also examine the impact of dusty galaxies using constraints from the 220 GHz data. The stacked sample provides 2.8$\sigma$ significant evidence of dusty galaxy flux, which would correspond to an average underestimate of the SPT Y_500 signal that is (17+-9) per cent in this sample of low mass systems.

In this paper we perform a comprehensive study of the main sources of random and systematic errors in stellar mass measurement for galaxies using their Spectral Energy Distributions (SEDs). We use mock galaxy catalogs with simulated multi-waveband photometry (from U-band to mid-infrared) and known redshift, stellar mass, age and extinction for individual galaxies. Given different parameters affecting stellar mass measurement (photometric S/N ratios, SED fitting errors, systematic effects, the inherent degeneracies and correlated errors), we formulated different simulated galaxy catalogs to quantify these effects individually. We studied the sensitivity of stellar mass estimates to the codes/methods used, population synthesis models, star formation histories, nebular emission line contributions, photometric uncertainties, extinction and age. For each simulated galaxy, the difference between the input stellar masses and those estimated using different simulation catalogs, $\Delta\log(M)$, was calculated and used to identify the most fundamental parameters affecting stellar masses. We measured different components of the error budget, with the results listed as follows: (1). no significant bias was found among different codes/methods, with all having comparable scatter; (2). A source of error is found to be due to photometric uncertainties and low resolution in age and extinction grids; (3). The median of stellar masses among different methods provides a stable measure of the mass associated with any given galaxy; (4). The deviations in stellar mass strongly correlate with those in age, with a weaker correlation with extinction; (5). the scatter in the stellar masses due to free parameters are quantified, with the sensitivity of the stellar mass to both the population synthesis codes and inclusion of nebular emission lines studied.

We present the discovery of two galaxy overdensities in the Hubble Space Telescope UDF: a proto-cluster, HUDFJ0332.4-2746.6 at $z = 1.84 \pm 0.01$, and a group, HUDFJ0332.5-2747.3 at $z =1.90 \pm 0.01$. Assuming viralization, the velocity dispersion of HUDFJ0332.4-2746.6 implies a mass of $M_{200}= (2.2 \pm 1.8) \times 10^{14} M_{\odot}$, consistent with the lack of extended X-ray emission. Neither overdensity shows evidence of a red sequence. About $50\%$ of their members show interactions and/or disturbed morphologies, which are signatures of merger remnants or disk instability. Most of their ETGs have blue colors and show recent star-formation. These observations reveal for the first time large fractions of spectroscopically confirmed star-forming blue ETGs in proto-clusters at $z\approx 2$. These star-forming ETGs are most likely among the progenitors of the quiescent population in clusters at more recent epochs. Their mass-size relation is consistent with that of passive ETGs in clusters at $z\sim0.7-1.5$. If these galaxies are the progenitors of cluster ETGs at these lower redshifts, their size would evolve according to a similar mass-size relation. It is noteworthy that quiescent ETGs in clusters at $z=1.8-2$ also do not show any significant size evolution over this redshift range, contrary to field ETGs. The ETG fraction is $\lesssim 50\%$, compared to the typical quiescent ETG fraction of $\approx 80\%$ in cluster cores at $z< 1$. The fraction, masses, and colors of the newly discovered ETGs imply that other cluster ETGs will be formed/accreted at a later time.

SPHEREx (Spectro-Photometer for the History of the Universe, Epoch of Reionization, and Ices Explorer) ( http://spherex.caltech.edu ) is a proposed all-sky spectroscopic survey satellite designed to address all three science goals in NASA's Astrophysics Division: probe the origin and destiny of our Universe; explore whether planets around other stars could harbor life; and explore the origin and evolution of galaxies. SPHEREx will scan a series of Linear Variable Filters systematically across the entire sky. The SPHEREx data set will contain R=40 spectra fir 0.75$<\lambda<$4.1$\mu$m and R=150 spectra for 4.1$<\lambda<$4.8$\mu$m for every 6.2 arc second pixel over the entire-sky. In this paper, we detail the extra-galactic and cosmological studies SPHEREx will enable and present detailed systematic effect evaluations. We also outline the Ice and Galaxy Evolution Investigations.

We present a catalog of galaxy clusters selected via their Sunyaev-Zel'dovich (SZ) effect signature from 2500 deg$^2$ of South Pole Telescope (SPT) data. This work represents the complete sample of clusters detected at high significance in the 2500-square-degree SPT-SZ survey, which was completed in 2011. A total of 677 (409) cluster candidates are identified above a signal-to-noise threshold of $\xi$ =4.5 (5.0). Ground- and space-based optical and near-infrared (NIR) imaging confirms overdensities of similarly colored galaxies in the direction of 516 (or 76%) of the $\xi$>4.5 candidates and 387 (or 95%) of the $\xi$>5 candidates; the measured purity is consistent with expectations from simulations. Of these confirmed clusters, 415 were first identified in SPT data, including 251 new discoveries reported in this work. We estimate photometric redshifts for all candidates with identified optical and/or NIR counterparts; we additionally report redshifts derived from spectroscopic observations for 141 of these systems. The mass threshold of the catalog is roughly independent of redshift above $z$~0.25 leading to a sample of massive clusters that extends to high redshift. The median mass of the sample is $M_{\scriptsize 500c}(\rho_\mathrm{crit})$ ~ 3.5 x 10$^{14} M_\odot h^{-1}$, the median redshift is $z_{med}$ =0.55, and the highest-redshift systems are at $z$>1.4. The combination of large redshift extent, clean selection, and high typical mass makes this cluster sample of particular interest for cosmological analyses and studies of cluster formation and evolution.

We present Herschel (PACS and SPIRE) far-infrared (FIR) photometry of a complete sample of z>1 3CR sources, from the Herschel GT project The Herschel Legacy of distant radio-loud AGN (PI: Barthel). Combining these with existing Spitzer photometric data, we perform an infrared (IR) spectral energy distribution (SED) analysis of these landmark objects in extragalactic research to study the star formation in the hosts of some of the brightest active galactic nuclei (AGN) known at any epoch. Accounting for the contribution from an AGN-powered warm dust component to the IR SED, about 40% of our objects undergo episodes of prodigious, ULIRG-strength star formation, with rates of hundreds of solar masses per year, coeval with the growth of the central supermassive black hole. Median SEDs imply that the quasar and radio galaxy hosts have similar FIR properties, in agreement with the orientation-based unification for radio-loud AGN. The star-forming properties of the AGN hosts are similar to those of the general population of equally massive non-AGN galaxies at comparable redshifts, thus there is no strong evidence of universal quenching of star formation (negative feedback) within this sample. Massive galaxies at high redshift may be forming stars prodigiously, regardless of whether their supermassive black holes are accreting or not.

Herschel/PACS observations of 29 local (Ultra-)Luminous Infrared Galaxies, including both starburst and AGN-dominated sources as diagnosed in the mid-infrared/optical, show that the equivalent width of the absorbing OH 65 um Pi_{3/2} J=9/2-7/2 line (W_{eq}(OH65)) with lower level energy E_{low}~300 K, is anticorrelated with the [C ii]158 um line to far-infrared luminosity ratio, and correlated with the far-infrared luminosity per unit gas mass and with the 60-to-100 um far-infrared color. While all sources are in the active L_{IR}/M_{H2}>50 Lsun/Msun mode as derived from previous CO line studies, the OH65 absorption shows a bimodal distribution with a discontinuity at L_{FIR}/M_{H2}~100 Lsun/Msun. In the most buried sources, OH65 probes material partially responsible for the silicate 9.7 um absorption. Combined with observations of the OH 71 um Pi_{1/2} J=7/2-5/2 doublet (E_{low}~415 K), radiative transfer models characterized by the equivalent dust temperature, Tdust, and the continuum optical depth at 100 um, tau_{100}, indicate that strong [C ii]158 um deficits are associated with far-IR thick (tau_{100}>~0.7, N_{H}>~10^{24} cm^{-2}), warm (Tdust>~60 K) structures where the OH 65 um absorption is produced, most likely in circumnuclear disks/tori/cocoons. With their high L_{FIR}/M_{H2} ratios and columns, the presence of these structures is expected to give rise to strong [C ii] deficits. W_{eq}(OH65) probes the fraction of infrared luminosity arising from these compact/warm environments, which is >~30-50% in sources with high W_{eq}({OH65}). Sources with high W_{eq}({OH65}) have surface densities of both L_{IR} and M_{H2} higher than inferred from the half-light (CO or UV/optical) radius, tracing coherent structures that represent the most buried/active stage of (circum)nuclear starburst-AGN co-evolution.

We present a velocity dispersion-based mass calibration of the South Pole Telescope Sunyaev-Zel'dovich effect survey (SPT-SZ) galaxy cluster sample. Using a homogeneously selected sample of 100 cluster candidates from 720 deg2 of the survey along with 63 velocity dispersion ($\sigma_v$) and 16 X-ray Yx measurements of sample clusters, we simultaneously calibrate the mass-observable relation and constrain cosmological parameters. The calibrations using $\sigma_v$ and Yx are consistent at the $0.6\sigma$ level, with the $\sigma_v$ calibration preferring ~16% higher masses. We use the full cluster dataset to measure $\sigma_8(\Omega_ m/0.27)^{0.3}=0.809\pm0.036$. The SPT cluster abundance is lower than preferred by either the WMAP9 or Planck+WMAP9 polarization (WP) data, but assuming the sum of the neutrino masses is $\sum m_\nu=0.06$ eV, we find the datasets to be consistent at the 1.0$\sigma$ level for WMAP9 and 1.5$\sigma$ for Planck+WP. Allowing for larger $\sum m_\nu$ further reconciles the results. When we combine the cluster and Planck+WP datasets with BAO and SNIa, the preferred cluster masses are $1.9\sigma$ higher than the Yx calibration and $0.8\sigma$ higher than the $\sigma_v$ calibration. Given the scale of these shifts (~44% and ~23% in mass, respectively), we execute a goodness of fit test; it reveals no tension, indicating that the best-fit model provides an adequate description of the data. Using the multi-probe dataset, we measure $\Omega_ m=0.299\pm0.009$ and $\sigma_8=0.829\pm0.011$. Within a $\nu$CDM model we find $\sum m_\nu = 0.148\pm0.081$ eV. We present a consistency test of the cosmic growth rate. Allowing both the growth index $\gamma$ and the dark energy equation of state parameter $w$ to vary, we find $\gamma=0.73\pm0.28$ and $w=-1.007\pm0.065$, demonstrating that the expansion and the growth histories are consistent with a LCDM model ($\gamma=0.55; \,w=-1$).

One of the most challenging aspects of studying galaxies in the z>~7 universe is the infrequent confirmation of their redshifts through spectroscopy, a phenomenon thought to occur from the increasing opacity of the intergalactic medium to Lya photons at z>6.5. The resulting redshift uncertainties inhibit the efficient search for [C II] in z~7 galaxies with sub-mm instruments such as ALMA, given their limited scan speed for faint lines. One means by which to improve the precision of the inferred redshifts is to exploit the potential impact of strong nebular emission lines on the colors of z~4-8 galaxies as observed by Spitzer/IRAC. At z~6.8, galaxies exhibit IRAC colors as blue as [3.6]-[4.5] ~-1, likely due to the contribution of [O III]+Hb to the 3.6 mum flux combined with the absence of line contamination in the 4.5 mum band. In this paper we explore the use of extremely blue [3.6]-[4.5] colors to identify galaxies in the narrow redshift window z~6.6-6.9. When combined with an I-dropout criterion, we demonstrate that we can plausibly select a relatively clean sample of z~6.8 galaxies. Through a systematic application of this selection technique to our catalogs from all five CANDELS fields, we identify 20 probable z~6.6-6.9 galaxies. We estimate that our criteria select the ~50% strongest line emitters at z~6.8 and from the IRAC colors we estimate a typical [O III]+Hb rest-frame equivalent width of 1085A for this sample. The small redshift uncertainties on our sample make it particularly well suited for follow-up studies with facilities such as ALMA.

The galaxy stellar mass function (GSMF) at high-z provides key information on star-formation history and mass assembly in the young Universe. We aimed to use the unique combination of deep optical/NIR/MIR imaging provided by HST, Spitzer and the VLT in the CANDELS-UDS, GOODS-South, and HUDF fields to determine the GSMF over the redshift range 3.5<z<7.5. We utilised the HST WFC3/IR NIR imaging from CANDELS and HUDF09, reaching H~27-28.5 over a total area of 369 arcmin2, in combination with associated deep HST ACS optical data, deep Spitzer IRAC imaging from the SEDS programme, and deep Y and K-band VLT Hawk-I images from the HUGS programme, to select a galaxy sample with high-quality photometric redshifts. These have been calibrated with more than 150 spectroscopic redshifts in the range 3.5<z<7.5, resulting in an overall precision of sigma_z/(1+z)~0.037. We have determined the low-mass end of the high-z GSMF with unprecedented precision, reaching down to masses as low as M*~10^9 Msun at z=4 and ~6x10^9 Msun at z=7. We find that the GSMF at 3.5<z<7.5 depends only slightly on the recipes adopted to measure the stellar masses, namely the photo-z, the SFHs, the nebular contribution or the presence of AGN on the parent sample. The low-mass end of the GSMF is steeper than has been found at lower redshifts, but appears to be unchanged over the redshift range probed here. Our results are very different from previous GSMF estimates based on converting UV galaxy luminosity functions into mass functions via tight M/L relations. Integrating our evolving GSMF over mass, we find that the growth of stellar mass density is barely consistent with the time-integral of the SFR density over cosmic time at z>4. These results confirm the unique synergy of the CANDELS+HUDF, HUGS, and SEDS surveys for the discovery and study of moderate/low-mass galaxies at high redshifts.

We present an analysis of the clustering of high-redshift galaxies in the recently completed 94 deg$^2$ Spitzer-SPT Deep Field survey. Applying flux and color cuts to the mid-infrared photometry efficiently selects galaxies at $z\sim1.5$ in the stellar mass range $10^{10}-10^{11}M_\odot$, making this sample the largest used so far to study such a distant population. We measure the angular correlation function in different flux-limited samples at scales $>6^{\prime \prime}$ (corresponding to physical distances $>0.05$ Mpc) and thereby map the one- and two-halo contributions to the clustering. We fit halo occupation distributions and determine how the central galaxy's stellar mass and satellite occupation depend on the halo mass. We measure a prominent peak in the stellar-to-halo mass ratio at a halo mass of $\log(M_{\rm halo} / M_\odot) = 12.44\pm0.08$, 4.5 times higher than the $z=0$ value. This supports the idea of an evolving mass threshold above which star formation is quenched. We estimate the large-scale bias in the range $b_g=2-4$ and the satellite fraction to be $f_\mathrm{sat}\sim0.2$, showing a clear evolution compared to $z=0$. We also find that, above a given stellar mass limit, the fraction of galaxies that are in similar mass pairs is higher at $z=1.5$ than at $z=0$. In addition, we measure that this fraction mildly increases with the stellar mass limit at $z=1.5$, which is the opposite of the behavior seen at low-redshift.

We are conducting a large survey of distant clusters of galaxies using radio sources with bent jets and lobes as tracers. These radio sources are driven by AGN and achieve their bent morphologies through interaction with the surrounding gas found in clusters of galaxies. Based on low-redshift studies, these types of sources can be used to identify clusters very efficiently. We present initial results from our survey of 653 bent-double radio sources with optical hosts too faint to appear in the SDSS. The sample was observed in the infrared with Spitzer, and it has revealed $\sim$200 distant clusters or proto-clusters in the redshift range $z\sim0.7 - 3.0$. The sample of bent-doubles contains both quasars and radio galaxies enabling us to study both radiative and kinetic mode feedback in cluster and group environments at a wide range of redshifts.

We combine photometry from the UDS, and CANDELS UDS and CANDELS GOODS-S surveys to construct the galaxy stellar mass function probing both the low and high mass end accurately in the redshift range 0.3<z<3. The advantages of using a homogeneous concatenation of these datasets include meaningful measures of environment in the UDS, due to its large area (0.88 deg^2), and the high resolution deep imaging in CANDELS (H_160 > 26.0), affording us robust measures of structural parameters. We construct stellar mass functions for the entire sample as parameterised by the Schechter function, and find that there is a decline in the values of phi and of alpha with higher redshifts, and a nearly constant M* up to z~3. We divide the galaxy stellar mass function by colour, structure, and environment and explore the links between environmental over-density, morphology, and the quenching of star formation. We find that a double Schechter function describes galaxies with high Sersic index (n>2.5), similar to galaxies which are red or passive. The low-mass end of the n>2.5 stellar mass function is dominated by blue galaxies, whereas the high-mass end is dominated by red galaxies. This hints that possible links between morphological evolution and star formation quenching are only present in high-mass galaxies. This is turn suggests that there are strong mass dependent quenching mechanisms. In addition, we find that the number density of high mass systems is elevated in dense environments, suggesting that an environmental process is building up massive galaxies quicker in over densities than in lower densities.

X-ray cavities are key tracers of mechanical (or radio mode) heating arising from the active galactic nuclei (AGN) in brightest cluster galaxies. We report on a survey for X-ray cavities in 83 massive, high-redshift (0.4<z<1.2) clusters of galaxies selected by their Sunyaev-Zel'dovich signature in the South Pole Telescope data. Based on Chandra X-ray images, we find a total of 6 clusters having symmetric pairs of surface brightness depressions consistent with the picture of radio jets inflating X-ray cavities in the intracluster medium. The majority of these detections are of relatively low significance and require deeper follow-up data in order to be confirmed. Further, due to the limitations of Chandra at high redshift, this search misses small (<10 kpc), unresolved X-ray cavities at high (z>0.5) redshift. Despite these limitations, our results suggest that the power generated by AGN feedback in brightest cluster galaxies has remained unchanged for over half of the age of the Universe (>7 Gyrs at z=0.8). On average, the detected X-ray cavities have powers of 0.8-5*10^45 erg/s, enthalpies of 3-6*10^59 erg, and radii of 17 kpc. Integrating over 7 Gyrs, we find that the supermassive black holes in the brightest cluster galaxies may have accreted 10^8 to several 10^9M_sun of material to power these outflows. This level of accretion indicates that significant supermassive black hole growth may occur not only at early times, in the quasar era, but at late times as well. We also find that X-ray cavities at high-redshift may inject an excess heat of 0.1-1.0 keV per particle into the hot intracluster medium above and beyond the energy needed to offset cooling. This value is similar to the energy needed to preheat clusters, break self-similarity, and explain the excess entropy in hot atmospheres.

We present photometric redshifts and associated probability distributions for all detected sources in the Extended Chandra Deep Field South (ECDFS). The work makes use of the most up-to-date data from the Cosmic Assembly Near-IR Deep Legacy Survey (CANDELS) and the Taiwan ECDFS Near-Infrared Survey (TENIS) in addition to other data. We also revisit multi-wavelength counterparts for published X-ray sources from the 4Ms-CDFS and 250ks-ECDFS surveys, finding reliable counterparts for 1207 out of 1259 sources ($\sim 96\%$). Data used for photometric redshifts include intermediate-band photometry deblended using the TFIT method, which is used for the first time in this work. Photometric redshifts for X-ray source counterparts are based on a new library of AGN/galaxy hybrid templates appropriate for the faint X-ray population in the CDFS. Photometric redshift accuracy for normal galaxies is 0.010 and for X-ray sources is 0.014, and outlier fractions are $4\%$ and $5.4\%$ respectively. The results within the CANDELS coverage area are even better as demonstrated both by spectroscopic comparison and by galaxy-pair statistics. Intermediate-band photometry, even if shallow, is valuable when combined with deep broad-band photometry. For best accuracy, templates must include emission lines.

The total infrared (IR) luminosity is very useful for estimating the star formation rate (SFR) of galaxies, but converting the IR luminosity into an SFR relies on assumptions that do not hold for all galaxies. We test the effectiveness of the IR luminosity as an SFR indicator by applying it to synthetic spectral energy distributions generated from three-dimensional hydrodynamical simulations of isolated disc galaxies and galaxy mergers. In general, the SFR inferred from the IR luminosity agrees well with the true instantaneous SFR of the simulated galaxies. However, for the major mergers in which a strong starburst is induced, the SFR inferred from the IR luminosity can overestimate the instantaneous SFR during the post-starburst phase by greater than two orders of magnitude. Even though the instantaneous SFR decreases rapidly after the starburst, the stars that were formed in the starburst remain dust-obscured and thus produce significant IR luminosity. Consequently, use of the IR luminosity as an SFR indicator may cause one to conclude that post-starburst galaxies are still star-forming, whereas in reality, star formation was recently quenched.

In a LCDM universe, disk galaxies' outer regions are the last to form. Characterizing their contents is critical for understanding the ongoing process of disk formation, but observing outer disk stellar populations is challenging due to their low surface brightness. We present extremely deep 3.6 micron observations (Spitzer IRAC) of NGC 4625, a galaxy known for its radially extended ultraviolet emitting stellar population. We combine the new imaging with archival UV imaging from the GALEX mission to derive multiwavelength radial profiles for NGC 4625 and compare them to stellar populations models. The colors can be explained by the young stellar population that is responsible for the UV emission and indicate that the current star formation rates in the outermost disk are recent. Extended star formation in NGC 4625 may have been initiated by an interaction with neighboring galaxies NGC 4618 and NGC 4625a, supporting speculation that minor interactions are a common trigger for outer disk star formation and late stage disk growth.

We use data taken as part of HST/WFC3 observations of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) to identify massive and evolved galaxies at 3<z<4.5. This is performed using the strength of the Balmer break feature at rest-frame 3648A, which is a diagnostic of the age of the stellar population in galaxies. Using WFC3 H-band selected catalog for the CANDELS GOODS-S field and deep multi-waveband photometry from optical (HST) to mid-infrared (Spitzer) wavelengths, we identify a population of old and evolved post-starburst galaxies based on the strength of their Balmer breaks (Balmer Break Galaxies- BBGs). The galaxies are also selected to be bright in rest-frame near-IR wavelengths and hence, massive. We identify a total of 16 BBGs. Fitting the spectral energy distribution (SED) of the BBGs show that the candidate galaxies have average estimated ages of ~800 Myr and average stellar masses of ~5x10^10 M_sun, consistent with being old and massive systems. Two of our BBG candidates are also identified by the criteria that is sensitive to star forming galaxies (LBG selection). We find a number density of ~3.2x10^-5 Mpc^-3 for the BBGs corresponding to a mass density of ~2.0x10^6 M_sun/Mpc^3 in the redshift range covering the survey. Given the old age and the passive evolution, it is argued that some of these objects formed the bulk of their mass only a few hundred million years after the Big Bang.

We measure new estimates for the galaxy stellar mass function and star formation rates for samples of galaxies at $z \sim 4,~5,~6~\&~7$ using data in the CANDELS GOODS South field. The deep near-infrared observations allow us to construct the stellar mass function at $z \geq 6$ directly for the first time. We estimate stellar masses for our sample by fitting the observed spectral energy distributions with synthetic stellar populations, including nebular line and continuum emission. The observed UV luminosity functions for the samples are consistent with previous observations, however we find that the observed $M_{UV}$ - M$_{*}$ relation has a shallow slope more consistent with a constant mass to light ratio and a normalisation which evolves with redshift. Our stellar mass functions have steep low-mass slopes ($\alpha \approx -1.9$), steeper than previously observed at these redshifts and closer to that of the UV luminosity function. Integrating our new mass functions, we find the observed stellar mass density evolves from $\log_{10} \rho_{*} = 6.64^{+0.58}_{-0.89}$ at $z \sim 7$ to $7.36\pm0.06$ $\text{M}_{\odot} \text{Mpc}^{-3}$ at $z \sim 4$. Finally, combining the measured UV continuum slopes ($\beta$) with their rest-frame UV luminosities, we calculate dust corrected star-formation rates (SFR) for our sample. We find the specific star-formation rate for a fixed stellar mass increases with redshift whilst the global SFR density falls rapidly over this period. Our new SFR density estimates are higher than previously observed at this redshift.

We present the first detection from the Spitzer Space Telescope of 4.5 micron variability from Sgr A*, the emitting source associated with the Milky Way's central black hole. The >23 hour continuous light curve was obtained with the IRAC instrument in 2013 December. The result characterizes the variability of Sgr A* prior to the closest approach of the G2 object, a putative infalling gas cloud that orbits close to Sgr A*. The high stellar density at the location of Sgr A* produces a background of ~250 mJy at 4.5 microns in each pixel with a large pixel-to-pixel gradient, but the light curve for the highly variable Sgr A* source was successfully measured by modeling and removing the variations due to pointing wobble. The observed flux densities range from the noise level of ~0.7 mJy rms in a 6.4-s measurement to ~10 mJy. Emission was seen above the noise level ~34% of the time. The light curve characteristics, including the flux density distribution and structure function, are consistent with those previously derived at shorter infrared wavelengths. We see no evidence in the light curve for activity attributable to the G2 interaction at the observing epoch, ~100 days before the expected G2 periapsis passage. The IRAC light curve is more than a factor of two longer than any previous infrared observation, improving constraints on the timescale of the break in the power spectral distribution of Sgr A* flux densities. The data favor the longer of the two previously published values for the timescale.

SPT-CLJ2040-4451 -- spectroscopically confirmed at z = 1.478 -- is the highest redshift galaxy cluster yet discovered via the Sunyaev-Zel'dovich effect. SPT-CLJ2040-4451 was a candidate galaxy cluster identified in the first 720 deg^2 of the South Pole Telescope Sunyaev-Zel'dovich (SPT-SZ) survey, and confirmed in follow-up imaging and spectroscopy. From multi-object spectroscopy with Magellan-I/Baade+IMACS we measure spectroscopic redshifts for 15 cluster member galaxies, all of which have strong [O II] 3727 emission. SPT-CLJ2040-4451 has an SZ-measured mass of M_500,SZ = 3.2 +/- 0.8 X 10^14 M_Sun/h_70, corresponding to M_200,SZ = 5.8 +/- 1.4 X 10^14 M_Sun/h_70. The velocity dispersion measured entirely from blue star forming members is sigma_v = 1500 +/- 520 km/s. The prevalence of star forming cluster members (galaxies with > 1.5 M_Sun/yr) implies that this massive, high-redshift cluster is experiencing a phase of active star formation, and supports recent results showing a marked increase in star formation occurring in galaxy clusters at z >1.4. We also compute the probability of finding a cluster as rare as this in the SPT-SZ survey to be >99%, indicating that its discovery is not in tension with the concordance Lambda-CDM cosmological model.

We present optical spectroscopy of galaxies in clusters detected through the Sunyaev-Zel'dovich (SZ) effect with the South Pole Telescope (SPT). We report our own measurements of $61$ spectroscopic cluster redshifts, and $48$ velocity dispersions each calculated with more than $15$ member galaxies. This catalog also includes $19$ dispersions of SPT-observed clusters previously reported in the literature. The majority of the clusters in this paper are SPT-discovered; of these, most have been previously reported in other SPT cluster catalogs, and five are reported here as SPT discoveries for the first time. By performing a resampling analysis of galaxy velocities, we find that unbiased velocity dispersions can be obtained from a relatively small number of member galaxies ($\lesssim 30$), but with increased systematic scatter. We use this analysis to determine statistical confidence intervals that include the effect of membership selection. We fit scaling relations between the observed cluster velocity dispersions and mass estimates from SZ and X-ray observables. In both cases, the results are consistent with the scaling relation between velocity dispersion and mass expected from dark-matter simulations. We measure a $\sim$30% log-normal scatter in dispersion at fixed mass, and a $\sim$10% offset in the normalization of the dispersion-mass relation when compared to the expectation from simulations, which is within the expected level of systematic uncertainty.

We present the results from a VLT/SINFONI and Keck/NIRSPEC near-infrared spectroscopic survey of 16 Lyman-alpha emitters (LAEs) at $z$ = 2.1 - 2.5 in the COSMOS and GOODS-N fields discovered from the HETDEX Pilot Survey. We detect rest-frame optical nebular lines (H$\alpha$ and/or [OIII]$\lambda$5007) for 10 of the LAEs and measure physical properties, including the star formation rate (SFR), gas-phase metallicity, gas-mass fraction, and Ly$\alpha$ velocity offset. We find that LAEs may lie below the mass-metallicity relation for continuum-selected star-forming galaxies at the same redshift. The LAEs all show velocity shifts of Ly$\alpha$ relative to the systemic redshift ranging between +85 and +296 km s$^{-1}$ with a mean of +180 km s$^{-1}$. This value is smaller than measured for continuum-selected star-forming galaxies at similar redshifts. The Ly$\alpha$ velocity offsets show a moderate correlation with the measured star formation rate (2.5$\sigma$), but no significant correlations are seen with the SFR surface density, specific SFR, stellar mass, or dynamical mass ($\lesssim$ 1.5$\sigma$). Exploring the role of dust, kinematics of the interstellar medium (ISM), and geometry on the escape of Ly$\alpha$ photons, we find no signature of selective quenching of resonantly scattered Ly$\alpha$ photons. However, we also find no evidence that a clumpy ISM is enhancing the Ly$\alpha$ equivalent width. Our results suggest that the low metallicity in LAEs may be responsible for yielding an environment with a low neutral hydrogen column density as well as less dust, easing the escape of Ly$\alpha$ photons over that in continuum-selected star-forming galaxies.

Herschel Space Observatory photometry and extensive multiwavelength followup have revealed that the powerful radio galaxy 3C 220.3 at z=0.685 acts as a gravitational lens for a background submillimeter galaxy (SMG) at z=2.221. At an observed wavelength of 1mm, the SMG is lensed into three distinct images. In the observed near infrared, these images are connected by an arc of 1.8" radius forming an Einstein half-ring centered near the radio galaxy. In visible light, only the arc is apparent. 3C 220.3 is the only known instance of strong galaxy-scale lensing by a powerful radio galaxy not located in a galaxy cluster and therefore it offers the potential to probe the dark matter content of the radio galaxy host. Lens modeling rejects a single lens, but two lenses centered on the radio galaxy host A and a companion B, separated by 1.5", provide a fit consistent with all data and reveal faint candidates for the predicted fourth and fifth images. The model does not require an extended common dark matter halo, consistent with the absence of extended bright X-ray emission on our Chandra image. The projected dark matter fractions within the Einstein radii of A (1.02") and B (0.61") are about 0.4 +/- 0.3 and 0.55 +/- 0.3. The mass to i-band light ratios of A and B, M/L ~ 8 +/- 4 Msun/Lsun, appear comparable to those of radio-quiet lensing galaxies at the same redshift in the CASTLES, LSD, and SL2S samples. The lensed SMG is extremely bright with observed f(250um) = 440mJy owing to a magnification factor mu~10. The SMG spectrum shows luminous, narrow CIV 154.9nm emission, revealing that the SMG houses a hidden quasar in addition to a violent starburst. Multicolor image reconstruction of the SMG indicates a bipolar morphology of the emitted ultraviolet (UV) light suggestive of cones through which UV light escapes a dust-enshrouded nucleus.

We report Plateau de Bure Interferometer (PdBI) 1.1 mm continuum imaging towards two extremely red H-[4.5]>4 (AB) galaxies at z>3, which we have previously discovered making use of Spitzer SEDS and Hubble Space Telescope CANDELS ultra-deep images of the Ultra Deep Survey field. One of our objects is detected on the PdBI map with a 4.3 sigma significance, corresponding to Snu(1.1mm)=(0.78 +/- 0.18) mJy. By combining this detection with the Spitzer 8 and 24 micron photometry for this source, and SCUBA2 flux density upper limits, we infer that this galaxy is a composite active galactic nucleus/star-forming system. The infrared (IR)-derived star formation rate is SFR~(200 +/- 100) Msun/yr, which implies that this galaxy is a higher-redshift analogue of the ordinary ultra-luminous infrared galaxies (ULIRGs) more commonly found at z~2-3. In the field of the other target, we find a tentative 3.1 sigma detection on the PdBI 1.1 mm map, but 3.7 arcsec away of our target position, so it likely corresponds to a different object. In spite of the lower significance, the PdBI detection is supported by a close SCUBA2 3.3 sigma detection. No counterpart is found on either the deep SEDS or CANDELS maps, so, if real, the PdBI source could be similar in nature to the sub-millimetre source GN10. We conclude that the analysis of ultra-deep near- and mid-IR images offers an efficient, alternative route to discover new sites of powerful star formation activity at high redshifts.

The Herschel very wide-field surveys have charted hundreds of square degrees in multiple far-IR (FIR) bands. While the Sloan Digital Sky Survey (SDSS) is currently the best resource for optical counterpart identifications over such wide areas, it does not detect a large number of Herschel FIR sources and leaves their nature undetermined. As a test case, we studied seven "SDSS-invisible", very bright 250um sources (S_{250} > 55 mJy) in the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) fields where we have a rich multi-wavelength data set. We took a new approach to decompose the FIR sources, using the near-IR or the optical images directly for position priors. This is an improvement over the previous decomposition efforts where the priors are from mid-IR data that still suffer from the source blending problem in the first place. We found that in most cases the single Herschel sources are made of multiple components that are not necessarily at the same redshifts. Our decomposition succeeded in identifying and extracting their major contributors. We show that these are all ULIRGs at z ~ 1--2 whose high L_IR is mainly due to dust-obscured star formation. Most of them would not be selected as sub-mm galaxies. They all have complicated morphologies indicative of merger or violent instability, and their stellar populations are heterogeneous in terms of stellar masses, ages and formation histories. Their current ULIRG phases are of various degrees of importance in their stellar mass assembly. Our practice provides a promising starting point to develop an automatic routine to reliably study bright Herschel sources.

Activity classification of galaxies based on long-slit and fiber spectroscopy can be strongly influenced by aperture effects. Here we investigate how activity classification for 14 nearby galaxies depends on the proportion of the host galaxy's light that is included in the aperture. We use both observed long-slit spectra and simulated elliptical-aperture spectra of different sizes. The degree of change varies with galaxy morphology and nuclear activity type. Starlight removal techniques can mitigate but not remove the effect of host galaxy contamination in the nuclear aperture. Galaxies with extra-nuclear star formation can show higher [O III] {\lambda}5007/H{\beta} ratios with increasing aperture, in contrast to the naive expectation that integrated light will only dilute the nuclear emission lines. We calculate the mean dispersion for the diagnostic line ratios used in the standard BPT diagrams with respect to the central aperture of spectral extraction to obtain an estimate of the uncertainties resulting from aperture effects.

We present 279 galaxy cluster candidates at $z > 1.3$ selected from the 94 deg$^{2}$ Spitzer South Pole Telescope Deep Field (SSDF) survey. We use a simple algorithm to select candidate high-redshift clusters of galaxies based on Spitzer/IRAC mid-infrared data combined with shallow all-sky optical data. We identify distant cluster candidates in SSDF adopting an overdensity threshold that results in a high purity (80%) cluster sample based on tests in the Spitzer Deep, Wide-Field Survey of the Bo\"otes field. Our simple algorithm detects all three $1.4 < z \leq 1.75$ X-ray detected clusters in the Bo\"otes field. The uniqueness of the SSDF survey resides not just in its area, one of the largest contiguous extragalactic fields observed with Spitzer, but also in its deep, multi-wavelength coverage by the South Pole Telescope (SPT), Herschel/SPIRE and XMM-Newton. This rich dataset will allow direct or stacked measurements of Sunyaev-Zel'dovich effect decrements or X-ray masses for many of the SSDF clusters presented here, and enable systematic study of the most distant clusters on an unprecedented scale. We measure the angular correlation function of our sample and find that these candidates show strong clustering. Employing the COSMOS/UltraVista photometric catalog in order to infer the redshift distribution of our cluster selection, we find that these clusters have a comoving number density $n_c = (0.7^{+6.3}_{-0.6}) \times 10^{-7} h^{3} \mathrm{Mpc}^{-3}$ and a spatial clustering correlation scale length $r_0 = (32 \pm 7) h^{-1} \rm{Mpc}$. Assuming our sample is comprised of dark matter halos above a characteristic minimum mass, $M_{{\rm min}}$, we derive that at $z=1.5$ these clusters reside in halos larger than $M_{{\rm min}} = 1.5^{+0.9}_{-0.7} \times 10^{14} h^{-1} M_{\odot}$. (abridged)

We present the average rest-frame spectrum of high-redshift dusty, star-forming galaxies from 250-770GHz. This spectrum was constructed by stacking ALMA 3mm spectra of 22 such sources discovered by the South Pole Telescope and spanning z=2.0-5.7. In addition to multiple bright spectral features of 12CO, [CI], and H2O, we also detect several faint transitions of 13CO, HCN, HNC, HCO+, and CN, and use the observed line strengths to characterize the typical properties of the interstellar medium of these high-redshift starburst galaxies. We find that the 13CO brightness in these objects is comparable to that of the only other z>2 star-forming galaxy in which 13CO has been observed. We show that the emission from the high-critical density molecules HCN, HNC, HCO+, and CN is consistent with a warm, dense medium with T_kin ~ 55K and n_H2 >~ 10^5.5 cm^-3. High molecular hydrogen densities are required to reproduce the observed line ratios, and we demonstrate that alternatives to purely collisional excitation are unlikely to be significant for the bulk of these systems. We quantify the average emission from several species with no individually detected transitions, and find emission from the hydride CH and the linear molecule CCH for the first time at high redshift, indicating that these molecules may be powerful probes of interstellar chemistry in high-redshift systems. These observations represent the first constraints on many molecular species with rest-frame transitions from 0.4-1.2mm in star-forming systems at high redshift, and will be invaluable in making effective use of ALMA in full science operations.

We present the first systematic comparison of ultraviolet-millimeter spectral energy distributions (SEDs) of observed and simulated interacting galaxies. Our sample is drawn from the Spitzer Interacting Galaxy Survey, and probes a range of galaxy interaction parameters. We use 31 galaxies in 14 systems which have been observed with Herschel, Spitzer, GALEX, and 2MASS. We create a suite of GADGET-3 hydrodynamic simulations of isolated and interacting galaxies with stellar masses comparable to those in our sample of interacting galaxies. Photometry for the simulated systems is then calculated with the SUNRISE radiative transfer code for comparison with the observed systems. For most of the observed systems, one or more of the simulated SEDs match reasonably well. The best matches recover the infrared luminosity and the star formation rate of the observed systems, and the more massive systems preferentially match SEDs from simulations of more massive galaxies. The most morphologically distorted systems in our sample are best matched to simulated SEDs close to coalescence, while less evolved systems match well with SEDs over a wide range of interaction stages, suggesting that an SED alone is insufficient to identify interaction stage except during the most active phases in strongly interacting systems. This result is supported by our finding that the SEDs calculated for simulated systems vary little over the interaction sequence.

We derive physical properties of 10 submillimeter galaxies located in the CANDELS coverage of the GOODS-S field. The galaxies were first identified as submillimeter sources with the LABOCA bolometer and subsequently targeted for 870um continuum observation with ALMA. The high angular resolution of the ALMA imaging allows secure counterparts to be identified in the CANDELS multiband dataset. The CANDELS data provide deep photometric data from UV through near-infrared wavelengths. Using synthetic spectral energy distributions, we derive photometric redshifts, stellar masses, extinction, ages, and the star formation history. The redshift range is z=1.65-4.76, with two of the galaxies located at z>4. Two SMG counterparts have stellar masses 2-3 orders of magnitude lower than the rest. The remaining SMG counterparts have stellar masses around 1x10^11 Msun. The stellar population in the SMGs is typically older than the expected duration of the submillimeter phase, suggesting that the star formation history of submillimeter galaxies is more complex than a single burst. Non-parametric morphology indices suggest that the SMG counterparts are among the most asymmetric systems compared with galaxies of the same stellar mass and redshift. The HST images shows that 3 of the SMGs are associated with on-going mergers. The remaining counterparts are isolated. Estimating the dust and molecular gas mass from the submm fluxes, and comparing with our stellar masses shows that the molecular gas mass fraction of SMGs is ~28% and that the final stellar mass is likely to be (1-2)x10^11 Msun.

[abridged] We present the discovery of four surprisingly bright (H_160 ~ 26 - 27 mag AB) galaxy candidates at z~9-10 in the complete HST CANDELS WFC3/IR GOODS-N imaging data, doubling the number of z~10 galaxy candidates that are known, just ~500 Myr after the Big Bang. Two similarly bright sources are also detected in a systematic re-analysis of the GOODS-S data set. Three of the four galaxies in GOODS-N are significantly detected at 4.5-6.2sigma in the very deep Spitzer/IRAC 4.5 micron data, as is one of the GOODS-S candidates. Furthermore, the brightest of our candidates (at z=10.2+-0.4) is robustly detected also at 3.6 micron (6.9sigma), revealing a flat UV spectral energy distribution with a slope beta=-2.0+-0.2, consistent with demonstrated trends with luminosity at high redshift. The abundance of such luminous candidates suggests that the luminosity function evolves more significantly in phi_* than in L_* at z>~8 with a higher number density of bright sources than previously expected. Despite the discovery of these luminous candidates, the cosmic star formation rate density for galaxies with SFR >0.7 M_sun/yr shows an order-of-magnitude increase in only 170 Myr from z ~ 10 to z ~ 8, consistent with previous results. Based on the IRAC detections, we derive galaxy stellar masses at z~10, finding that these luminous objects are typically 10^9 M_sun. The cosmic stellar mass density at z~10 is log10 rho_* = 4.7^+0.5_-0.8 M_sun Mpc^-3 for galaxies brighter than M_UV~-18. The remarkable brightness, and hence luminosity, of these z~9-10 candidates highlights the opportunity for deep spectroscopy to determine their redshift and nature, demonstrates the value of additional search fields covering a wider area to understand star-formation in the very early universe, and highlights the opportunities for JWST to map the buildup of galaxies at redshifts much earlier than z~10.

The adiabatic evolution of the temperature of the cosmic microwave background (CMB) is a key prediction of standard cosmology. We study deviations from the expected adiabatic evolution of the CMB temperature of the form $T(z) =T_0(1+z)^{1-\alpha}$ using measurements of the spectrum of the Sunyaev Zel'dovich Effect with the South Pole Telescope (SPT). We present a method for using the ratio of the Sunyaev Zel'dovich signal measured at 95 and 150 GHz in the SPT data to constrain the temperature of the CMB. We demonstrate that this approach provides unbiased results using mock observations of clusters from a new set of hydrodynamical simulations. We apply this method to a sample of 158 SPT-selected clusters, spanning the redshift range $0.05 < z < 1.35$, and measure $\alpha = 0.017^{+0.030}_{-0.028}$, consistent with the standard model prediction of $\alpha=0$. In combination with other published results, we constrain $\alpha = 0.011 \pm 0.016$, an improvement of $\sim 20\%$ over published constraints. This measurement also provides a strong constraint on the effective equation of state in models of decaying dark energy $w_\mathrm{eff} = -0.987^{+0.016}_{-0.017}$.

We analyze the star-forming and structural properties of 45 massive (log(M/Msun)>10) compact star-forming galaxies (SFGs) at 2<z<3 to explore whether they are progenitors of compact quiescent galaxies at z~2. The optical/NIR and far-IR Spitzer/Herschel colors indicate that most compact SFGs are heavily obscured. Nearly half (47%) host an X-ray bright AGN. In contrast, only about 10% of other massive galaxies at that time host AGNs. Compact SFGs have centrally-concentrated light profiles and spheroidal morphologies similar to quiescent galaxies, and are thus strikingly different from other SFGs. Most compact SFGs lie either within the SFR-M main sequence (65%) or below (30%), on the expected evolutionary path towards quiescent galaxies. These results show conclusively that galaxies become more compact before they lose their gas and dust, quenching star formation. Using extensive HST photometry from CANDELS and grism spectroscopy from the 3D-HST survey, we model their stellar populations with either exponentially declining (tau) star formation histories (SFHs) or physically-motivated SFHs drawn from semi-analytic models (SAMs). SAMs predict longer formation timescales and older ages ~2 Gyr, which are nearly twice as old as the estimates of the tau models. While both models yield good SED fits, SAM SFHs better match the observed slope and zero point of the SFR-M main sequence. Some low-mass compact SFGs (log(M/Msun)=10-10.6) have younger ages but lower sSFRs than that of more massive galaxies, suggesting that the low-mass galaxies reach the red sequence faster. If the progenitors of compact SFGs are extended SFGs, state-of-the-art SAMs show that mergers and disk instabilities are both able to shrink galaxies, but disk instabilities are more frequent (60% versus 40%) and form more concentrated galaxies. We confirm this result via high-resolution hydrodynamic simulations.

We present the results of a combined galaxy population analysis for the host galaxies of active galactic nuclei (AGN) identified at 0 < z < 1.4 within the SDSS, Bootes and DEEP2 surveys. We identified AGN in a uniform and unbiased manner at X-ray, infrared and radio wavelengths. Supermassive black holes undergoing radiatively-efficient accretion (detected as X-ray and/or infrared AGN) appear to be hosted in a separate and distinct galaxy population than AGN undergoing powerful mechanically dominated accretion (radio AGN). Consistent with some previous studies, radiatively efficient AGN appear to be preferentially hosted in modest star-forming galaxies, with little dependence on AGN or galaxy luminosity. AGN exhibiting radio-emitting jets due to mechanically-dominated accretion are almost exclusively observed in massive, passive galaxies. Crucially, we now provide strong evidence that the observed host-galaxy trends are independent of redshift. In particular, these different accretion-mode AGN have remained as separate galaxy populations throughout the last 9 Gyr. Furthermore, it appears that galaxies hosting AGN have evolved along the same path as galaxies that are not hosting AGN with little evidence for distinctly separate evolution.

Out of several dozen z > 7 candidate galaxies observed spectroscopically, only five have been confirmed via Lyman-alpha emission, at z=7.008, 7.045, 7.109, 7.213 and 7.215. The small fraction of confirmed galaxies may indicate that the neutral fraction in the intergalactic medium (IGM) rises quickly at z > 6.5, as Lyman-alpha is resonantly scattered by neutral gas. However, the small samples and limited depth of previous observations makes these conclusions tentative. Here we report the results of a deep near-infrared spectroscopic survey of 43 z > 6.5 galaxies. We detect only a single galaxy, confirming that some process is making Lyman-alpha difficult to detect. The detected emission line at 1.0343 um is likely to be Lyman-alpha emission, placing this galaxy at a redshift z = 7.51, an epoch 700 million years after the Big Bang. This galaxy's colors are consistent with significant metal content, implying that galaxies become enriched rapidly. We measure a surprisingly high star formation rate of 330 Msol/yr, more than a factor of 100 greater than seen in the Milky Way. Such a galaxy is unexpected in a survey of our size, suggesting that the early universe may harbor more intense sites of star-formation than expected.

We analyze the star formation properties of 16 infrared-selected, spectroscopically confirmed galaxy clusters at $1 < z < 1.5$ from the Spitzer/IRAC Shallow Cluster Survey (ISCS). We present new spectroscopic confirmation for six of these high-redshift clusters, five of which are at $z>1.35$. Using infrared luminosities measured with deep Spitzer/MIPS observations at 24 $\mu$m, along with robust optical+IRAC photometric redshifts and SED-fitted stellar masses, we present the dust-obscured star-forming fractions, star formation rates and specific star formation rates in these clusters as functions of redshift and projected clustercentric radius. We find that $z\sim 1.4$ represents a transition redshift for the ISCS sample, with clear evidence of an unquenched era of cluster star formation at earlier times. Beyond this redshift the fraction of star-forming cluster members increases monotonically toward the cluster centers. Indeed, the specific star formation rate in the cores of these distant clusters is consistent with field values at similar redshifts, indicating that at $z>1.4$ environment-dependent quenching had not yet been established in ISCS clusters. Combining these observations with complementary studies showing a rapid increase in the AGN fraction, a stochastic star formation history, and a major merging episode at the same epoch in this cluster sample, we suggest that the starburst activity is likely merger-driven and that the subsequent quenching is due to feedback from merger-fueled AGN. The totality of the evidence suggests we are witnessing the final quenching period that brings an end to the era of star formation in galaxy clusters and initiates the era of passive evolution.

The Spitzer-South Pole Telescope Deep Field (SSDF) is a wide-area survey using Spitzer's Infrared Array Camera (IRAC) to cover 94 square degrees of extragalactic sky, making it the largest IRAC survey completed to date outside the Milky Way midplane. The SSDF is centered at 23:30,-55:00, in a region that combines observations spanning a broad wavelength range from numerous facilities. These include millimeter imaging from the South Pole Telescope, far-infrared observations from Herschel/SPIRE, X-ray observations from the XMM XXL survey, near-infrared observations from the VISTA Hemisphere Survey, and radio-wavelength imaging from the Australia Telescope Compact Array, in a panchromatic project designed to address major outstanding questions surrounding galaxy clusters and the baryon budget. Here we describe the Spitzer/IRAC observations of the SSDF, including the survey design, observations, processing, source extraction, and publicly available data products. In particular, we present two band-merged catalogs, one for each of the two warm IRAC selection bands. They contain roughly 5.5 and 3.7 million distinct sources, the vast majority of which are galaxies, down to the SSDF 5-sigma sensitivity limits of 19.0 and 18.2 Vega mag (7.0 and 9.4 microJy) at 3.6 and 4.5 microns, respectively.

We use the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS) GOODS-S multi-wavelength catalog to identify counterparts for 20 Lya Emitting (LAE) galaxies at z=2.1. We build several types of stacked Spectral Energy Distributions (SEDs) of these objects. We combine photometry to form average and median flux-stacked SEDs, and postage stamp images to form average and median image-stacked SEDs. We also introduce scaled flux stacks that eliminate the influence of variation in overall brightness. We use the SED fitting code SpeedyMC to constrain the physical properties of individual objects and stacks. Our LAEs at z = 2.1 have stellar masses ranging from 2x10^7 Msun - 8x10^9 Msun (median = 3x10^8 Msun), ages ranging from 4 Myr to 500 Myr (median =100 Myr), and E(B-V) between 0.02 and 0.24 (median = 0.12). We do not observe strong correlations between Lya equivalent width (EW) and stellar mass, age, or E(B-V). The Lya radiative transfer (q) factors of our sample are predominantly close to one and do not correlate strongly with EW or E(B-V), implying that Lya radiative transfer prevents Lya photons from resonantly scattering in dusty regions. The SED parameters of the flux stacks match the average and median values of the individual objects, with the flux-scaled median SED performing best with reduced uncertainties. Median image-stacked SEDs provide a poor representation of the median individual object, and none of the stacking methods captures the large dispersion of LAE properties.

We present observations of SPT-S J053816-5030.8, a gravitationally-lensed dusty star forming galaxy (DSFG) at z = 2.7817, first discovered at millimeter wavelengths by the South Pole Telescope. SPT 0538-50 is typical of the brightest sources found by wide-field millimeter-wavelength surveys, being lensed by an intervening galaxy at moderate redshift (in this instance, at z = 0.441). We present a wide array of multi-wavelength spectroscopic and photometric data on SPT 0538-50, including data from ALMA, Herschel PACS and SPIRE, Hubble, Spitzer, VLT, ATCA, APEX, and the SMA. We use high resolution imaging from HST to de-blend SPT 0538-50, separating DSFG emission from that of the foreground lens. Combined with a source model derived from ALMA imaging (which suggests a magnification factor of 21 +/- 4), we derive the intrinsic properties of SPT 0538-50, including the stellar mass, far-IR luminosity, star formation rate, molecular gas mass, and - using molecular line fluxes - the excitation conditions within the ISM. The derived physical properties argue that we are witnessing compact, merger-driven star formation in SPT 0538-50, similar to local starburst galaxies, and unlike that seen in some other DSFGs at this epoch.

We present first results on the cooling properties derived from Chandra X-ray observations of 83 high-redshift (0.3 < z < 1.2) massive galaxy clusters selected by their Sunyaev-Zel'dovich signature in the South Pole Telescope data. We measure each cluster's central cooling time, central entropy, and mass deposition rate, and compare to local cluster samples. We find no significant evolution from z~0 to z~1 in the distribution of these properties, suggesting that cooling in cluster cores is stable over long periods of time. We also find that the average cool core entropy profile in the inner ~100 kpc has not changed dramatically since z ~ 1, implying that feedback must be providing nearly constant energy injection to maintain the observed "entropy floor" at ~10 keV cm^2. While the cooling properties appear roughly constant over long periods of time, we observe strong evolution in the gas density profile, with the normalized central density (rho_0/rho_crit) increasing by an order of magnitude from z ~ 1 to z ~ 0. When using metrics defined by the inner surface brightness profile of clusters, we find an apparent lack of classical, cuspy, cool-core clusters at z > 0.75, consistent with earlier reports for clusters at z > 0.5 using similar definitions. Our measurements indicate that cool cores have been steadily growing over the 8 Gyr spanned by our sample, consistent with a constant, ~150 Msun/yr cooling flow that is unable to cool below entropies of 10 keV cm^2 and, instead, accumulates in the cluster center. We estimate that cool cores began to assemble in these massive systems at z ~ 1, which represents the first constraints on the onset of cooling in galaxy cluster cores. We investigate several potential biases which could conspire to mimic this cool core evolution and are unable to find a bias that has a similar redshift dependence and a substantial amplitude.

We report deep ALMA observations complemented with associated HST imaging for a luminous (m_uv=25) galaxy, `Himiko', at a redshift z=6.595. The galaxy is remarkable for its high star formation rate, 100 Mo/yr, securely estimated from our deep HST and Spitzer photometry, and the absence of any evidence for strong AGN activity or gravitational lensing magnification. Our ALMA observations probe an order of magnitude deeper than previous IRAM observations, yet fail to detect a 1.2mm dust continuum, indicating a flux <52uJy comparable with or weaker than that of local dwarf irregulars with much lower star formation rates. We likewise provide a strong upper limit for the flux of [CII] 158um, L([CII]) < 5.4x10^7 Lo, a diagnostic of the hot interstellar gas often described as a valuable probe for early galaxies. In fact, our observations indicate Himiko lies off the local L([CII]) - star formation rate scaling relation by a factor of more than 30. Both aspects of our ALMA observations suggest Himiko is an unique object with a very low dust content and perhaps nearly primordial interstellar gas. Our HST images provide unique insight into the morphology of this remarkable source, highlighting an extremely blue core of activity and two less extreme associated clumps. Himiko is undergoing a triple major merger event whose extensive ionized nebula of Lyman alpha emitting gas, discovered in our earlier work with Subaru, is powered by star formation and the dense circum-galactic gas. We are likely witnessing an early massive galaxy during a key period of its mass assembly close to the end of the reionization era.

We present a UV-to-mid infrared multi-wavelength catalog in the CANDELS/GOODS-S field, combining the newly obtained CANDELS HST/WFC3 F105W, F125W, and F160W data with existing public data. The catalog is based on source detection in the WFC3 F160W band. The F160W mosaic includes the data from CANDELS deep and wide observations as well as previous ERS and HUDF09 programs. The mosaic reaches a 5$\sigma$ limiting depth (within an aperture of radius 0.17 arcsec) of 27.4, 28.2, and 29.7 AB for CANDELS wide, deep, and HUDF regions, respectively. The catalog contains 34930 sources with the representative 50% completeness reaching 25.9, 26.6, and 28.1 AB in the F160W band for the three regions. In addition to WFC3 bands, the catalog also includes data from UV (U-band from both CTIO/MOSAIC and VLT/VIMOS), optical (HST/ACS F435W, F606W, F775W, F814W, and F850LP), and infrared (HST/WFC3 F098M, VLT/ISAAC Ks, VLT/HAWK-I Ks, and Spitzer/IRAC 3.6, 4.5, 5.8, 8.0 $\mu$m) observations. The catalog is validated via stellar colors, comparison with other published catalogs, zeropoint offsets determined from the best-fit templates of the spectral energy distribution of spectroscopically observed objects, and the accuracy of photometric redshifts. The catalog is able to detect unreddened star-forming (passive) galaxies with stellar mass of 10^{10}M_\odot at a 50% completeness level to z$\sim$3.4 (2.8), 4.6 (3.2), and 7.0 (4.2) in the three regions. As an example of application, the catalog is used to select both star-forming and passive galaxies at z$\sim$2--4 via the Balmer break. It is also used to study the color--magnitude diagram of galaxies at 0<z<4.

Chandra observations of a complete, flux-limited sample of 38 high-redshift (1<z<2), low-frequency selected (and so unbiased in orientation) 3CRR radio sources are reported. The sample includes 21 quasars (= broad line radio galaxies) and 17 narrow line radio galaxies (NLRGs) with matched 178 MHz radio luminosity (log L_R ~ 44-45). The quasars have high radio core-fraction, high X-ray luminosities (log L_X ~ 45-46) and soft X-ray hardness ratios (HR ~ -0.5) indicating low obscuration. The NLRGs have lower core-fraction, lower apparent X-ray luminosities (log L_X ~ 43-45) and mostly hard X-ray hardness ratios (HR>0) indicating obscuration (log N_H ~ 22-24 cm^-2). These properties and the correlation between obscuration and radio core-fraction are consistent with orientation-dependent obscuration as in Unification models. About half the NLRGs have soft X-ray hardness ratios and/or high [OIII] emission line to X-ray luminosity ratio suggesting obscuration by Compton thick (CT) material so that scattered nuclear or extended X-ray emission dominates (as in NGC1068). The ratios of unobscured to Compton-thin (10^{22}< N_H(int) < 1.5 x 10^{24} cm^-2) to CT (N_H(int) > 1.5 x 10^{24} cm^-2) is 2.5:1.4:1 in this high luminosity, radio-selected sample. The obscured fraction is 0.5, higher than is typically reported for AGN at comparable luminosities from multi-wavelength surveys (0.1-0.3). Assuming random nuclear orientation, the unobscured half-opening angle of the disk/wind/torus structure is ~ 60deg and the obscuring material covers 30deg of which ~ 12deg is Compton thick. The multi-wavelength properties reveal that many NLRGs have intrinsic absorption 10-1000x higher than indicated by their X-ray hardness ratios, and their true L_X values are ~10--100x larger than the hardness-ratio absorption corrections would indicate.

We present the multiwavelength - ultraviolet to mid-infrared - catalog of the UKIDSS Ultra-Deep Survey (UDS) field observed as part of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). Based on publicly available data, the catalog includes: the CANDELS data from the Hubble Space Telescope (near-infrared WFC3 F125W and F160W data and visible ACS F606W and F814W data), u-band data from CFHT/Megacam, B, V, Rc, i' and z' band data from Subaru/Suprime-Cam, Y and Ks band data from VLT/HAWK-I, J, H and K bands data from UKIDSS (Data Release 8), and Spitzer/IRAC data (3.6, 4.5 from SEDS, 5.8 and 8.0um from SpUDS). The present catalog is F160W-selected and contains 35932 sources over an area of 201.7 square arcmin and includes radio and X-ray detected sources and spectroscopic redshifts available for 210 sources.

We report the detection of CO(1-0) line emission in the bright, lensed star-forming galaxies SPT-S 233227-5358.5 (z=2.73) and SPT-S 053816-5030.8 (z=2.78), using the Australia Telescope Compact Array (ATCA). Both galaxies were discovered in a large-area millimeter survey with the South Pole Telescope (SPT) and found to be gravitationally lensed by intervening structures. The measured CO intensities imply galaxies with molecular gas masses of (3.2 \pm 0.5)x10^10 (mu/15)^{-1}(X_CO/0.8) M_sun and (1.7 \pm 0.3)x10^10 (mu/20)^{-1}(X_CO/0.8) M_sun, and gas depletion timescales of 4.9x10^7 (X_CO/0.8) yr and 2.6x10^7 (X_CO/0.8) yr, respectively, where mu corresponds to the lens magnification and X_CO is the CO luminosity to gas mass conversion factor. In the case of SPT-S 053816-5030.8, we also obtained significant detections of the rest-frame 115.7 and 132.4 GHz radio continuum. Based on the radio to infrared spectral energy distribution and an assumed synchrotron spectral index, we find that 42 \pm 10 % and 55 \pm 13 % of the flux at rest-frame 115.7 and 132.4 GHz arise from free-free emission. We find a radio-derived intrinsic star formation rate (SFR) of 470 \pm 170 M_sun/yr, consistent within the uncertainties with the infrared estimate. Based on the morphology of this object in the source plane, the derived gas mass and the possible flattening of the radio spectral index towards low frequencies, we argue that SPT-S 053816-5030.8 exhibits properties compatible with a scaled-up local ultra-luminous infrared galaxy.

Star-forming galaxies at redshifts z>6 are likely responsible for the reionization of the universe, and it is important to study the nature of these galaxies. We present three candidates for z~7 Lyman-break galaxies (LBGs) from a 155 arcmin^2 area in the CANDELS/COSMOS field imaged by the deep FourStar Galaxy Evolution (zFourGE) survey. The FourStar medium-band filters provide the equivalent of R~10 spectroscopy, which cleanly distinguishes between z~7 LBGs and brown dwarf stars. The distinction between stars and galaxies based on an object's angular size can become unreliable even when using HST imaging; there exists at least one very compact z~7 candidate (FWHM~0.5-1 kpc) that is indistinguishable from a point source. The medium-band filters provide narrower redshift distributions compared with broad-band-derived redshifts. The UV luminosity function derived using the three z~7 candidates is consistent with previous studies, suggesting an evolution at the bright end (MUV -21.6 mag) from z~7 to z~5. Fitting the galaxies' spectral energy distributions, we predict Lyman-alpha equivalent widths for the two brightest LBGs, and find that the presence of a Lyman-alpha line affects the medium-band flux thereby changing the constraints on stellar masses and UV spectral slopes. This illustrates the limitations of deriving LBG properties using only broad-band photometry. The derived specific star-formation rates for the bright LBGs are ~13 per Gyr, slightly higher than the lower-luminosity LBGs, implying that the star-formation rate increases with stellar mass for these galaxies.

In the past decade, our understanding of galaxy evolution has been revolutionized by the discovery that luminous, dusty, starburst galaxies were 1,000 times more abundant in the early Universe than at present. It has, however, been difficult to measure the complete redshift 2 distribution of these objects, especially at the highest redshifts (z > 4). Here we report a redshift survey at a wavelength of three millimeters, targeting carbon monoxide line emission from the star-forming molecular gas in the direction of extraordinarily bright millimetrewave-selected sources. High-resolution imaging demonstrates that these sources are strongly gravitationally lensed by foreground galaxies. We detect spectral lines in 23 out of 26 sources and multiple lines in 12 of those 23 sources, from which we obtain robust, unambiguous redshifts. At least 10 of the sources are found to lie at z > 4, indicating that the fraction of dusty starburst galaxies at high redshifts is greater than previously thought. Models of lens geometries in the sample indicate that the background objects are ultra-luminous infrared galaxies, powered by extreme bursts of star formation.

We present Atacama Large Millimeter/submillimeter Array (ALMA) 860 micrometer imaging of four high-redshift (z=2.8-5.7) dusty sources that were detected using the South Pole Telescope (SPT) at 1.4 mm and are not seen in existing radio to far-infrared catalogs. At 1.5 arcsec resolution, the ALMA data reveal multiple images of each submillimeter source, separated by 1-3 arcsec, consistent with strong lensing by intervening galaxies visible in near-IR imaging of these sources. We describe a gravitational lens modeling procedure that operates on the measured visibilities and incorporates self-calibration-like antenna phase corrections as part of the model optimization, which we use to interpret the source structure. Lens models indicate that SPT0346-52, located at z=5.7, is one of the most luminous and intensely star-forming sources in the universe with a lensing corrected FIR luminosity of 3.7 X 10^13 L_sun and star formation surface density of 4200 M_sun yr^-1 kpc^-2. We find magnification factors of 5 to 22, with lens Einstein radii of 1.1-2.0 arcsec and Einstein enclosed masses of 1.6-7.2x10^11 M_sun. These observations confirm the lensing origin of these objects, allow us to measure the their intrinsic sizes and luminosities, and demonstrate the important role that ALMA will play in the interpretation of lensed submillimeter sources.

Using the Atacama Large Millimeter/submillimeter Array (ALMA), we have conducted a blind redshift survey in the 3 mm atmospheric transmission window for 26 strongly lensd dusty star-forming galaxies (DSFGs) selected with the South Pole Telescope (SPT). The sources were selected to have S_1.4mm>20 mJy and a dust-like spectrum and, to remove low-z sources, not have bright radio (S_843MHz<6mJy) or far-infrared counterparts (S_100um<1 Jy, S_60um<200mJy). We robustly detect 44 line features in our survey, which we identify as redshifted emission lines of 12CO, 13CO, [CI], H2O, and H2O+. We find one or more spectral features in 23 sources yielding a ~90% detection rate for this survey; in 12 of these sources we detect multiple lines, while in 11 sources we detect only a single line. For the sources with only one detected line, we break the redshift degeneracy with additional spectroscopic observations if available, or infer the most likely line identification based on photometric data. This yields secure redshifts for ~70% of the sample. The three sources with no lines detected are tentatively placed in the redshift desert between 1.7<z<2.0. The resulting mean redshift of our sample is <z>=3.5. This finding is in contrast to the redshift distribution of radio-identified DSFGs, which have a significantly lower mean redshift of <z>=2.3 and for which only 10-15% of the population is expected to be at z>3. We discuss the effect of gravitational lensing on the redshift distribution and compare our measured redshift distribution to that of models in the literature.

We measured and modeled the spectral energy distributions (SEDs) in 28 bands from the ultraviolet to the far-infrared (FIR) for 31 interacting galaxies in 14 systems. The sample is drawn from the Spitzer Interacting Galaxy Survey, which probes a range of galaxy interaction parameters at multiple wavelengths with an emphasis on the infrared bands. The subset presented in this paper consists of all galaxies for which FIR Herschel SPIRE observations are publicly available. Our SEDs combine the Herschel photometry with multi-wavelength data from Spitzer, GALEX, Swift UVOT, and 2MASS. While the shapes of the SEDs are broadly similar across our sample, strongly interacting galaxies typically have more mid-infrared emission relative to their near-infrared and FIR emission than weakly or moderately interacting galaxies. We modeled the full SEDs to derive host galaxy star formation rates (SFR), specific star formation rates (sSFR), stellar masses, dust temperatures, dust luminosities, and dust masses. We find increases in the dust luminosity and mass, SFR, and cold (15-25 K) dust temperature as the interaction progresses from moderately to strongly interacting and between non- interacting and strongly interacting galaxies. We also find increases in the SFR between weakly and strongly interacting galaxies. In contrast, the sSFR remains unchanged across all the interaction stages. The SFR derived from the SED modeling agrees well with rates estimated by proportionality relations that depend on infrared emission. (abridged)

We report the first investigation of cool-core properties of galaxy clusters selected via their Sunyaev--Zel'dovich (SZ) effect. We use 13 galaxy clusters uniformly selected from 178 deg^2 observed with the South Pole Telescope (SPT) and followed up by the Chandra X-ray Observatory. They form an approximately mass-limited sample (> 3 x 10^14 M_sun h^-1_70) spanning redshifts 0.3 < z < 1.1. Using previously published X-ray-selected cluster samples, we compare two proxies of cool-core strength: surface brightness concentration (cSB) and cuspiness ({\alpha}). We find that cSB is better constrained. We measure cSB for the SPT sample and find several new z > 0.5 cool-core clusters, including two strong cool cores. This rules out the hypothesis that there are no z > 0.5 clusters that qualify as strong cool cores at the 5.4{\sigma} level. The fraction of strong cool-core clusters in the SPT sample in this redshift regime is between 7% and 56% (95% confidence). Although the SPT selection function is significantly different from the X-ray samples, the high-z cSB distribution for the SPT sample is statistically consistent with that of X-ray-selected samples at both low and high redshifts. The cool-core strength is inversely correlated with the offset between the brightest cluster galaxy and the X-ray centroid, providing evidence that the dynamical state affects the cool-core strength of the cluster. Larger SZ-selected samples will be crucial in understanding the evolution of cluster cool cores over cosmic time.

We present a multi-wavelength study of a 3.6 $\mu$m-selected galaxy sample in the Extended Groth strip. The sample is complete for galaxies with stellar mass $>10^{9.5}$ \Msun and redshift $0.4<z<1.2$. In this redshift range, the IRAC 3.6 $\mu$m band measures the rest-frame near-infrared band, permitting nearly unbiased selection with respect to both quiescent and star-forming galaxies. The numerous spectroscopic redshifts available in the EGS are used to train an Artificial Neural Network to estimate photometric redshifts. The distribution of photometric redshift errors is Gaussian with standard deviation ${\sim}0.025(1+z)$, and the fraction of redshift failures (${>}3\sigma$ errors) is about 3.5%. A new method of validation based on pair statistics confirms the estimate of standard deviation even for galaxies lacking spectroscopic redshifts. Basic galaxy properties measured include rest-frame $U-B$ colors, $B$- and $K$-band absolute magnitudes, and stellar masses. We divide the sample into quiescent and star-forming galaxies according to their rest-frame $U-B$ colors and 24 to 3.6 \micron\ flux density ratios and derive rest $K$-band luminosity functions and stellar mass functions for quiescent, star forming, and all galaxies. The results show that massive, quiescent galaxies were in place by $z\approx1$, but lower mass galaxies generally ceased their star formation at later epochs.

One key goal of the Hubble Space Telescope Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey is to track galaxy evolution back to z ~ 8. Its two-tiered "wide and deep" strategy bridges significant gaps in existing near-infrared surveys. Here we report on z ~ 8 galaxy candidates selected as F105W-band dropouts in one of its deep fields, which covers 50.1 square arcmin to 4 ks depth in each of three near-infrared bands in the Great Observatories Origins Deep Survey southern field. Two of our candidates have J<26.2 mag, and are > 1 mag brighter than any previously known F105W-dropouts. We derive constraints on the bright-end of the rest-frame ultraviolet luminosity function of galaxies at z ~ 8, and show that the number density of such very bright objects is higher than expected from the previous Schechter luminosity function estimates at this redshift. Another two candidates are securely detected in Spitzer Infrared Array Camera images, which are the first such individual detections at z ~ 8. Their derived stellar masses are on the order of a few x 10^9 M_sun, from which we obtain the first measurement of the high-mass end of the galaxy stellar mass function at z ~ 8. The high number density of very luminous and very massive galaxies at z ~ 8, if real, could imply a large stellar-to-halo mass ratio and an efficient conversion of baryons to stars at such an early time.

We measure the clustering of Extremely Red Objects (EROs) in ~8 deg^2 of the NOAO Deep Wide Field Survey Bo\"otes field in order to establish robust links between ERO z~1.2 and local galaxy z<0.1 populations. Three different color selection criteria from the literature are analyzed to assess the consequences of using different criteria for selecting EROs. Specifically, our samples are (R-K_s)>5.0 (28,724 galaxies), (I-K_s)>4.0 (22,451 galaxies) and (I-[3.6])>5.0 (64,370 galaxies). Magnitude-limited samples show the correlation length (r_0) to increase for more luminous EROs, implying a correlation with stellar mass. We can separate star-forming and passive ERO populations using the (K_s-[24]) and ([3.6]-[24]) colors to K_s=18.4 and [3.6]=17.5, respectively. Star-forming and passive EROs in magnitude limited samples have different clustering properties and host dark halo masses, and cannot be simply understood as a single population. Based on the clustering, we find that bright passive EROs are the likely progenitors of >4L^* elliptical galaxies. Bright EROs with ongoing star formation were found to occupy denser environments than star-forming galaxies in the local Universe, making these the likely progenitors of >L^* local ellipticals. This suggests that the progenitors of massive >4L^* local ellipticals had stopped forming stars by z>1.2, but that the progenitors of less massive ellipticals (down to L^*) can still show significant star formation at this epoch.

We present the results of the ground- and space-based optical and near-infrared (NIR) follow-up of 224 galaxy cluster candidates detected with the Sunyaev-Zel'dovich (SZ) effect in the 720 deg^2 of the South Pole Telescope (SPT) survey completed in the 2008 and 2009 observing seasons. We use the optical/NIR data to establish whether each candidate is associated with an overdensity of galaxies and to estimate the cluster redshift. Most photometric redshifts are derived through a combination of three different cluster redshift estimators using red-sequence galaxies, resulting in an accuracy of \Delta z/(1+z)=0.017, determined through comparison with a subsample of 57 clusters for which we have spectroscopic redshifts. We successfully measure redshifts for 158 systems and present redshift lower limits for the remaining candidates. The redshift distribution of the confirmed clusters extends to z=1.35 with a median of z_{med}=0.57. Approximately 18% of the sample with measured redshifts lies at z>0.8. We estimate a lower limit to the purity of this SPT SZ-selected sample by assuming that all unconfirmed clusters are noise fluctuations in the SPT data. We show that the cumulative purity at detection significance \xi>5 (\xi>4.5) is >= 95 (>= 70%). We present the red brightest cluster galaxy (rBCG) positions for the sample and examine the offsets between the SPT candidate position and the rBCG. The radial distribution of offsets is similar to that seen in X-ray-selected cluster samples, providing no evidence that SZ-selected cluster samples include a different fraction of recent mergers than X-ray-selected cluster samples.

We compare cosmic microwave background lensing convergence maps derived from South Pole Telescope (SPT) data with galaxy survey data from the Blanco Cosmology Survey, the Wide-field Infrared Survey Explorer, and a new large Spitzer/IRAC field designed to overlap with the SPT survey. Using optical and infrared catalogs covering between 17 and 68 square degrees of sky, we detect correlation between the SPT convergence maps and each of the galaxy density maps at >4 sigma, with zero cross-correlation robustly ruled out in all cases. The amplitude and shape of the cross-power spectra are in good agreement with theoretical expectations and the measured galaxy bias is consistent with previous work. The detections reported here utilize a small fraction of the full 2500 square degree SPT survey data and serve as both a proof of principle of the technique and an illustration of the potential of this emerging cosmological probe.

We measure the faint end slope of the galaxy luminosity function (LF) for cluster galaxies at 1<z<1.5 using Spitzer IRAC data. We investigate whether this slope, alpha, differs from that of the field LF at these redshifts, and with the cluster LF at low redshifts. The latter is of particular interest as low-luminosity galaxies are expected to undergo significant evolution. We use seven high-redshift spectroscopically confirmed galaxy clusters drawn from the IRAC Shallow Cluster Survey to measure the cluster galaxy LF down to depths of M* + 3 (3.6 microns) and M* + 2.5 (4.5 microns). The summed LF at our median cluster redshift (z=1.35) is well fit by a Schechter distribution with alpha[3.6] = -0.97 +/- 0.14 and alpha[4.5] = -0.91 +/- 0.28, consistent with a flat faint end slope and is in agreement with measurements of the field LF in similar bands at these redshifts. A comparison to alpha in low-redshift clusters finds no statistically significant evidence of evolution. Combined with past studies which show that M* is passively evolving out to z~1.3, this means that the shape of the cluster LF is largely in place by z~1.3. This suggests that the processes that govern the build up of the mass of low-mass cluster galaxies have no net effect on the faint end slope of the cluster LF at z<1.3.

Unresolved near-infrared background anisotropies are expected to have contributions from the earliest galaxies during reionization and faint, dwarf galaxies at intermediate redshifts. Previous measurements were unable to conclusively pinpoint the dominant origin because they did not sample spatial scales that were sufficiently large to distinguish between these two possibilities. Here we report a measurement of the anisotropy power spectrum from sub-arcminute to one degree angular scales and find the clustering amplitude to be larger than the model predictions involving the two existing explanations. As the shot-noise level of the power spectrum is consistent with that expected from faint galaxies, a new source population on the sky is not necessary to explain the observations. A physical mechanism that increases the clustering amplitude, however, is needed. Motivated by recent results related to the extended stellar light profile in dark matter halos, we consider the possibility that the fluctuations originate from diffuse intrahalo stars of all galaxies. We find that the measured power spectrum can be explained by an intrahalo light fraction of 0.07 to 0.2 % relative to the total luminosity in dark matter halos of masses log(M/M_Sun) ~ 9 to 12 at redshifts of ~ 1 to 4.

The galaxy cluster SPT-CL J0205-5829 currently has the highest spectroscopically-confirmed redshift, z=1.322, in the South Pole Telescope Sunyaev-Zel'dovich (SPT-SZ) survey. XMM-Newton observations measure a core-excluded temperature of Tx=8.7keV producing a mass estimate that is consistent with the Sunyaev-Zel'dovich derived mass. The combined SZ and X-ray mass estimate of M500=(4.9+/-0.8)e14 h_{70}^{-1} Msun makes it the most massive known SZ-selected galaxy cluster at z>1.2 and the second most massive at z>1. Using optical and infrared observations, we find that the brightest galaxies in SPT-CL J0205-5829 are already well evolved by the time the universe was <5 Gyr old, with stellar population ages >3 Gyr, and low rates of star formation (<0.5Msun/yr). We find that, despite the high redshift and mass, the existence of SPT-CL J0205-5829 is not surprising given a flat LambdaCDM cosmology with Gaussian initial perturbations. The a priori chance of finding a cluster of similar rarity (or rarer) in a survey the size of the 2500 deg^2 SPT-SZ survey is 69%.

We use weak gravitational lensing to measure the masses of five galaxy clusters selected from the South Pole Telescope (SPT) survey, with the primary goal of comparing these with the SPT Sunyaev--Zel'dovich (SZ) and X-ray based mass estimates. The clusters span redshifts 0.28 < z < 0.43 and have masses M_500 > 2 x 10^14 h^-1 M_sun, and three of the five clusters were discovered by the SPT survey. We observed the clusters in the g'r'i' passbands with the Megacam imager on the Magellan Clay 6.5m telescope. We measure a mean ratio of weak lensing (WL) aperture masses to inferred aperture masses from the SZ data, both within an aperture of R_500,SZ derived from the SZ mass, of 1.04 +/- 0.18. We measure a mean ratio of spherical WL masses evaluated at R_500,SZ to spherical SZ masses of 1.07 +/- 0.18, and a mean ratio of spherical WL masses evaluated at R_500,WL to spherical SZ masses of 1.10 +/- 0.24. We explore potential sources of systematic error in the mass comparisons and conclude that all are subdominant to the statistical uncertainty, with dominant terms being cluster concentration uncertainty and N-body simulation calibration bias. Expanding the sample of SPT clusters with WL observations has the potential to significantly improve the SPT cluster mass calibration and the resulting cosmological constraints from the SPT cluster survey. These are the first WL detections using Megacam on the Magellan Clay telescope.

In the cores of some galaxy clusters the hot intracluster plasma is dense enough that it should cool radiatively in the cluster's lifetime, leading to continuous "cooling flows" of gas sinking towards the cluster center, yet no such cooling flow has been observed. The low observed star formation rates and cool gas masses for these "cool core" clusters suggest that much of the cooling must be offset by astrophysical feedback to prevent the formation of a runaway cooling flow. Here we report X-ray, optical, and infrared observations of the galaxy cluster SPT-CLJ2344-4243 at z = 0.596. These observations reveal an exceptionally luminous (L_2-10 keV = 8.2 x 10^45 erg/s) galaxy cluster which hosts an extremely strong cooling flow (dM/dt = 3820 +/- 530 Msun/yr). Further, the central galaxy in this cluster appears to be experiencing a massive starburst (740 +/- 160 Msun/yr), which suggests that the feedback source responsible for preventing runaway cooling in nearby cool core clusters may not yet be fully established in SPT-CLJ2344-4243. This large star formation rate implies that a significant fraction of the stars in the central galaxy of this cluster may form via accretion of the intracluster medium, rather than the current picture of central galaxies assembling entirely via mergers.

We report the results of a study exploring the stellar populations of 13 luminous (L>L*), spectroscopically confirmed, galaxies in the redshift interval 5.5<z<6.5, all with WFC3/IR and IRAC imaging from the HST/CANDELS and Spitzer/SEDS surveys. Based on fitting the observed photometry with SED templates covering a wide range of different star-formation histories, and a self consistent treatment of Lyman-alpha emission, we find that the derived stellar masses lie within the range 10^9 Msun < M*< 10^10 Msun and are robust to within a factor of two. In contrast, we confirm previous reports that the ages of the stellar populations are poorly constrained. Although the best-fitting models for three objects have ages >= 300 Myr, the degeneracies introduced by dust extinction mean that only two of these objects actually require a >300 Myr old stellar population to reproduce the observed photometry. Moreover, when considering only smoothly-varying star-formation histories, we observe a clear tension between the data and models such that a galaxy SED template with an old age is often chosen in order to try and fit objects with blue UV-slopes but red UV-to-optical colours. To break this tension we explore SED fitting with two-component models (burst plus on-going star-formation) and allow for nebular emission. On average, the inclusion of nebular emission leads to lower stellar-mass estimates (median offset 0.18 dex), moderately higher specific star-formation rates, and allows for a wider range of plausible stellar ages. However, based on our SED modelling, we find no strong evidence for extremely young ages in our sample (<50 Myr). Finally, considering all of the different star-formation histories explored, we find that the median best-fitting ages are of the order 200-300 Myr and that the objects with the tightest constraints indicate ages in the range 50-200 Myr (Abridged).

We combine high-resolution HST/WFC3 images with multi-wavelength photometry to track the evolution of structure and activity of massive (log(M*) > 10) galaxies at redshifts z = 1.4 - 3 in two fields of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). We detect compact, star-forming galaxies (cSFGs) whose number densities, masses, sizes, and star formation rates qualify them as likely progenitors of compact, quiescent, massive galaxies (cQGs) at z = 1.5 - 3. At z > 2 most cSFGs have specific star-formation rates (sSFR = 10^-9 yr^-1) half that of typical, massive SFGs at the same epoch, and host X-ray luminous AGN 30 times (~30%) more frequently. These properties suggest that cSFGs are formed by gas-rich processes (mergers or disk-instabilities) that induce a compact starburst and feed an AGN, which, in turn, quench the star formation on dynamical timescales (few 10^8 yr). The cSFGs are continuously being formed at z = 2 - 3 and fade to cQGs by z = 1.5. After this epoch, cSFGs are rare, thereby truncating the formation of new cQGs. Meanwhile, down to z = 1, existing cQGs continue to enlarge to match local QGs in size, while less-gas-rich mergers and other secular mechanisms shepherd (larger) SFGs as later arrivals to the red sequence. In summary, we propose two evolutionary scenarios of QG formation: an early (z > 2), fast-formation path of rapidly-quenched cSFGs that evolve into cQGs that later enlarge within the quiescent phase, and a slow, late-arrival (z < 2) path for SFGs to form QGs without passing through a compact state.

We present APEX SABOCA 350micron and LABOCA 870micron observations of 11 representative examples of the rare, extremely bright (S_1.4mm > 15mJy), dust-dominated millimeter-selected galaxies recently discovered by the South Pole Telescope (SPT). All 11 sources are robustly detected with LABOCA with 40 < S_870micron < 130mJy, approximately an order of magnitude higher than the canonical submillimeter galaxy (SMG) population. Six of the sources are also detected by SABOCA at >3sigma, with the detections or upper limits providing a key constraint on the shape of the spectral energy distribution (SED) near its peak. We model the SEDs of these galaxies using a simple modified blackbody and perform the same analysis on samples of SMGs of known redshift from the literature. These calibration samples inform the distribution of dust temperature for similar SMG populations, and this dust temperature prior allows us to derive photometric redshift estimates and far infrared luminosities for the sources. We find a median redshift of <z> = 3.0, higher than the <z> = 2.2 inferred for the normal SMG population. We also derive the apparent size of the sources from the temperature and apparent luminosity, finding them to appear larger than our unlensed calibration sample, which supports the idea that these sources are gravitationally magnified by massive structures along the line of sight.

Infrared 3.6 to 8 micron images of the Extended Groth Strip yield plausible counterpart identifications for all but one of 510 radio sources in the AEGIS20 S(1.4 GHz) > 50 micro-Jy sample. This is the first such deep sample that has been effectively 100% identified. Achieving the same identification rate at R-band would require observations reaching R_AB > 27. Spectroscopic redshifts are available for 46% of the sample and photometric redshifts for an additional 47%. Almost all of the sources with 3.6 micron AB magnitudes brighter than 19 have spectroscopic redshifts z < 1.1, while fainter objects predominantly have photometric redshifts with 1 \lapprox z \lapprox 3. Unlike more powerful radio sources that are hosted by galaxies having large stellar masses within a relatively narrow range, the AEGIS20 counterparts have stellar masses spanning more than a factor of 10 at z \sim 1. The sources are roughly 10--15% starbursts at z \lapprox 0.5 and 20--25% AGNs mostly at z > 1 with the remainder of uncertain nature.

The transneptunian objects (TNOs) trapped in mean-motion resonances with Neptune were likely emplaced there during planet migration late in the giant-planet formation process. We perform detailed modelling of the resonant objects detected in the Canada-France Ecliptic Plane Survey (CFEPS) in order to provide population estimates and, for some resonances, constrain the complex internal orbital element distribution. Detection biases play a critical role because phase relationships with Neptune make object discovery more likely at certain longitudes. This paper discusses the 3:2, 5:2, 2:1, 3:1, 5:1, 4:3, 5:3, 7:3, 5:4, and 7:4 mean-motion resonances, all of which had CFEPS detections, along with our upper limit on 1:1 Neptune Trojans (which is consistent with their small population estimated elsewhere). For the plutinos (TNOs in the 3:2 resonance) we refine the orbital element distribution given in Kavelaars et al. (2009) and show that steep H-magnitude distributions (N(H) proportional to 10aH, with a=0.8-0.9) are favoured in the range Hg=8-9, and confirm that this resonance does not share the inclination distribution of the classical Kuiper Belt. We give the first population estimate for the 5:2 resonance and find that, to within the uncertainties, the population is equal to that of the 3:2 (13,000 TNOs with Hg < 9.16), whereas the 2:1 population is smaller by a factor of 3-4 compared to the other two resonances. We also measure significant populations inhabiting the 4:3, 5:3, 7:3, 5:4, 7:4, 3:1, and 5:1 resonances, with Hg < 9.16 (D >100 km) populations in the thousands. We compare our intrinsic population and orbital-element distributions with several published models of resonant-TNO production; the most striking discrepancy is that resonances beyond the 2:1 are in reality more heavily populated than in published models.

We extend previous measurements of cosmic infrared background (CIB) fluctuations to ~ 1 deg using new data from the Spitzer Extended Deep Survey. Two fields, with depths of ~12 hr/pixel over 3 epochs, are analyzed at 3.6 and 4.5 mic. Maps of the fields were assembled using a self-calibration method uniquely suitable for probing faint diffuse backgrounds. Resolved sources were removed from the maps to a magnitude limit of AB mag ~ 25, as indicated by the level of the remaining shot noise. The maps were then Fourier-transformed and their power spectra were evaluated. Instrumental noise was estimated from the time-differenced data, and subtracting this isolates the spatial fluctuations of the actual sky. The power spectra of the source-subtracted fields remain identical (within the observational uncertainties) for the three epochs indicating that zodiacal light contributes negligibly to the fluctuations. Comparing to 8 mic power spectra shows that Galactic cirrus cannot account for the fluctuations. The signal appears isotropically distributed on the sky as required for an extragalactic origin. The CIB fluctuations continue to diverge to > 10 times those of known galaxy populations on angular scales out to < 1 deg. The low shot noise levels remaining in the diffuse maps indicate that the large scale fluctuations arise from the spatial clustering of faint sources well below the confusion noise. The spatial spectrum of these fluctuations is in reasonable agreement with an origin in populations clustered according to the standard cosmological model (LCDM) at epochs coinciding with the first stars era.

We have analysed a sample of 25 extremely red H-[4.5]>4 galaxies, selected using 4.5 micron data from the Spitzer SEDS survey and deep H-band data from the Hubble Space Telescope (HST) CANDELS survey, over ~180 square arcmin of the UKIDSS Ultra Deep Survey (UDS) field. Our aim is to investigate the nature of this rare population of mid-infrared (mid-IR) sources that display such extreme near-to-mid-IR colours. Using up to 17-band photometry (U through 8.0 microns), we have studied in detail their spectral energy distributions, including possible degeneracies in the photometric redshift/internal extinction (zphot-Av) plane. Our sample appears to include sources of very different nature. Between 45% and 75% of them are dust-obscured, massive galaxies at 3<zphot<5. All of the 24 micron-detected sources in our sample are in this category. Two of these have S(24 micron)>300 microJy, which at 3<zphot<5 suggests that they probably host a dust-obscured active galactic nucleus (AGN). Our sample also contains four highly obscured (Av>5) sources at zphot<1. Finally, we analyse in detail two zphot~6 galaxy candidates, and discuss their plausibility and implications. Overall, our red galaxy sample contains the tip of the iceberg of a larger population of z>3 galaxies to be discovered with the future James Webb Space Telescope.

We present a catalog of 224 galaxy cluster candidates, selected through their Sunyaev-Zel'dovich (SZ) effect signature in the first 720 deg2 of the South Pole Telescope (SPT) survey. This area was mapped with the SPT in the 2008 and 2009 austral winters to a depth of 18 uK-arcmin at 150 GHz; 550 deg2 of it was also mapped to 44 uK-arcmin at 95 GHz. Based on optical imaging of all candidates and near-infrared imaging of the majority of candidates, we have found optical and/or infrared counterparts for 158 clusters. Of these, 135 were first identified as clusters in SPT data, including 117 new discoveries reported in this work. This catalog triples the number of confirmed galaxy clusters discovered through the SZ effect. We report photometrically derived (and in some cases spectroscopic) redshifts for confirmed clusters and redshift lower limits for the remaining candidates. The catalog extends to high redshift with a median redshift of z = 0.55 and maximum redshift of z = 1.37. Based on simulations, we expect the catalog to be nearly 100% complete above M500 ~ 5e14 Msun h_{70}^-1 at z > 0.6. There are 121 candidates detected at signal-to-noise greater than five, at which the catalog purity is measured to be 95%. From this high-purity subsample, we exclude the z < 0.3 clusters and use the remaining 100 candidates to improve cosmological constraints following the method presented by Benson et al., 2011. Adding the cluster data to CMB+BAO+H0 data leads to a preference for non-zero neutrino masses while only slightly reducing the upper limit on the sum of neutrino masses to sum mnu < 0.38 eV (95% CL). For a spatially flat wCDM cosmological model, the addition of this catalog to the CMB+BAO+H0+SNe results yields sigma8=0.807+-0.027 and w = -1.010+-0.058, improving the constraints on these parameters by a factor of 1.4 and 1.3, respectively. [abbrev]

The South Pole Telescope (SPT) is a 10 meter telescope operating at mm wavelengths. It has recently completed a three-band survey covering 2500 sq. degrees. One of the survey's main goals is to detect galaxy clusters using Sunyaev-Zeldovich effect and use these clusters for a variety of cosmological and astrophysical studies such as the dark energy equation of state, the primordial non-gaussianity and the evolution of galaxy populations. Since 2005, we have been engaged in a comprehensive optical and near-infrared followup program (at wavelengths between 0.4 and 5 {\mu}m) to image high-significance SPT clusters, to measure their photometric redshifts, and to estimate the contamination rate of the candidate lists. These clusters are then used for various cosmological and astrophysical studies.

The evolution of galaxies is greatly influenced by their interactions. As part of a program to study interacting galaxies, we have measured and modeled the spectral energy distri- butions (SEDs) from the ultraviolet (UV) to the far-infrared (FIR). We describe the constraints imposed on star formation histories by these SEDs, and the variations therein seen across the interaction sequence, and we compare the results of different star formation rate prescriptions applied to the data. The sample itself is based on the Spitzer Interacting Galaxy Survey (SIGS) of 111 galaxies in 50 systems, a project designed to probe a range of galaxy interaction parameters in the infrared. Our SEDs combine the Spitzer results with multiwavelength data from other missions, in particular GALEX and Herschel. The subset presented here is the sample for which FIR Herschel observations are currently publicly available.

This paper describes the Hubble Space Telescope imaging data products and data reduction procedures for the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). This survey is designed to document the evolution of galaxies and black holes at $z\sim1.5-8$, and to study Type Ia SNe beyond $z>1.5$. Five premier multi-wavelength sky regions are selected, each with extensive multiwavelength observations. The primary CANDELS data consist of imaging obtained in the Wide Field Camera 3 / infrared channel (WFC3/IR) and UVIS channel, along with the Advanced Camera for Surveys (ACS). The CANDELS/Deep survey covers \sim125 square arcminutes within GOODS-N and GOODS-S, while the remainder consists of the CANDELS/Wide survey, achieving a total of \sim800 square arcminutes across GOODS and three additional fields (EGS, COSMOS, and UDS). We summarize the observational aspects of the survey as motivated by the scientific goals and present a detailed description of the data reduction procedures and products from the survey. Our data reduction methods utilize the most up to date calibration files and image combination procedures. We have paid special attention to correcting a range of instrumental effects, including CTE degradation for ACS, removal of electronic bias-striping present in ACS data after SM4, and persistence effects and other artifacts in WFC3/IR. For each field, we release mosaics for individual epochs and eventual mosaics containing data from all epochs combined, to facilitate photometric variability studies and the deepest possible photometry. A more detailed overview of the science goals and observational design of the survey are presented in a companion paper.

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) is designed to document the first third of galactic evolution, over the approximate redshift (z) range 8--1.5. It will image >250,000 distant galaxies using three separate cameras on the Hubble Space Telescope, from the mid-ultraviolet to the near-infrared, and will find and measure Type Ia supernovae at z>1.5 to test their accuracy as standardizable candles for cosmology. Five premier multi-wavelength sky regions are selected, each with extensive ancillary data. The use of five widely separated fields mitigates cosmic variance and yields statistically robust and complete samples of galaxies down to a stellar mass of 10^9 M_\odot to z \approx 2, reaching the knee of the ultraviolet luminosity function (UVLF) of galaxies to z \approx 8. The survey covers approximately 800 arcmin^2 and is divided into two parts. The CANDELS/Deep survey (5\sigma\ point-source limit H=27.7 mag) covers \sim 125 arcmin^2 within GOODS-N and GOODS-S. The CANDELS/Wide survey includes GOODS and three additional fields (EGS, COSMOS, and UDS) and covers the full area to a 5\sigma\ point-source limit of H \gtrsim 27.0 mag. Together with the Hubble Ultra Deep Fields, the strategy creates a three-tiered "wedding cake" approach that has proven efficient for extragalactic surveys. Data from the survey are nonproprietary and are useful for a wide variety of science investigations. In this paper, we describe the basic motivations for the survey, the CANDELS team science goals and the resulting observational requirements, the field selection and geometry, and the observing design. The Hubble data processing and products are described in a companion paper.

We identify an abundant population of extreme emission line galaxies (EELGs) at redshift z~1.7 in the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS) imaging from Hubble Space Telescope/Wide Field Camera 3 (HST/WFC3). 69 EELG candidates are selected by the large contribution of exceptionally bright emission lines to their near-infrared broad-band magnitudes. Supported by spectroscopic confirmation of strong [OIII] emission lines -- with rest-frame equivalent widths ~1000\AA -- in the four candidates that have HST/WFC3 grism observations, we conclude that these objects are galaxies with 10^8 Msol in stellar mass, undergoing an enormous starburst phase with M_*/(dM_*/dt) of only ~15 Myr. These bursts may cause outflows that are strong enough to produce cored dark matter profiles in low-mass galaxies. The individual star formation rates and the co-moving number density (3.7x10^-4 Mpc^-3) can produce in ~4 Gyr much of the stellar mass density that is presently contained in 10^8-10^9 Msol dwarf galaxies. Therefore, our observations provide a strong indication that many or even most of the stars in present-day dwarf galaxies formed in strong, short-lived bursts, mostly at z>1.

Chandra X-ray observations of the high redshift (z =1.532) radio-loud quasar 3C270.1 in 2008 February show the nucleus to have a power-law spectrum, Gamma = 1.66 +/- 0.08, typical of a radio-loud quasar, and a marginally-detected Fe Kalpha emission line. The data also reveal extended X-ray emission, about half of which is associated with the radio emission from this source. The southern emission is co-spatial with the radio lobe and peaks at the position of the double radio hotspot. Modeling this hotspot including Spitzer upper limits rules out synchrotron emission from a single power-law population of electrons, favoring inverse-Compton emission with a field of ~11nT, roughly a third of the equipartition value. The northern emission is concentrated close to the location of a 40 deg. bend where the radio jet is presumed to encounter external material. It can be explained by inverse Compton emission involving Cosmic Microwave Background photons with a field of ~3nT, roughly a factor of nine below the equipartition value. The remaining, more diffuse X-ray emission is harder (HR=-0.09 +/- 0.22). With only 22.8+/-5.6 counts, the spectral form cannot be constrained. Assuming thermal emission with a temperature of 4 keV yields an estimate for the luminosity of 1.8E44 erg/s, consistent with the luminosity-temperature relation of lower-redshift clusters. However deeper Chandra X-ray observations are required to delineate the spatial distribution, and better constrain the spectrum of the diffuse emission to verify that we have detected X-ray emission from a high-redshift cluster.

Star formation is arguably the most important physical process in the cosmos. It is a fundamental driver of galaxy evolution and the ultimate source of most of the energy emitted by galaxies. A correct interpretation of star formation rate (SFR) measures is therefore essential to our understanding of galaxy formation and evolution. Unfortunately, however, no single SFR estimator is universally available or even applicable in all circumstances: the numerous galaxies found in deep surveys are often too faint (or too distant) to yield significant detections with most standard SFR measures, and until now there have been no global, multi-band observations of nearby galaxies that span all the conditions under which star-formation is taking place. To address this need in a systematic way, we have undertaken a multi-band survey of all types of star-forming galaxies in the local Universe. This project, the Star Formation Reference Survey (SFRS), is based on a statistically valid sample of 369 nearby galaxies that span all existing combinations of dust temperature, SFR, and specific SFR. Furthermore, because the SFRS is blind with respect to AGN fraction and environment it serves as a means to assess the influence of these factors on SFR. Our panchromatic global flux measurements (including GALEX FUV+NUV, SDSS ugriz, 2MASS JHKs, Spitzer 3-8{\mu}m, and others) furnish uniform SFR measures and the context in which their reliability can be assessed. This paper describes the SFRS survey strategy, defines the sample, and presents the multi-band photometry collected to date.

We report the orbital distribution of the trans-neptunian objects (TNOs) discovered during the Canada-France Ecliptic Plane Survey, whose discovery phase ran from early 2003 until early 2007. The follow-up observations started just after the first discoveries and extended until late 2009. We obtained characterized observations of 321 sq.deg. of sky to depths in the range g ~ 23.5--24.4 AB mag. We provide a database of 169 TNOs with high-precision dynamical classification and known discovery efficiency. Using this database, we find that the classical belt is a complex region with sub-structures that go beyond the usual splitting of inner (interior to 3:2 mean-motion resonance [MMR]), outer (exterior to 2:1 MMR), and main (in between). The main classical belt (a=40--47 AU) needs to be modeled with at least three components: the `hot' component with a wide inclination distribution and two `cold' components (stirred and kernel) with much narrower inclination distributions. The hot component must have a significantly shallower absolute magnitude (Hg) distribution than the other two components. With 95% confidence, there are 8000+1800-1600 objects in the main belt with Hg <= 8.0, of which 50% are from the hot component, 40% from the stirred component and 10% from the kernel; the hot component's fraction drops rapidly with increasing Hg. Because of this, the apparent population fractions depend on the depth and ecliptic latitude of a trans-neptunian survey. The stirred and kernel components are limited to only a portion of the main belt, while we find that the hot component is consistent with a smooth extension throughout the inner, main and outer regions of the classical belt; the inner and outer belts are consistent with containing only hot-component objects. The Hg <= 8.0 TNO population estimates are 400 for the inner belt and 10,000 for the outer belt within a factor of two.

The South Pole Telescope (SPT) is currently surveying 2500 deg^2 of the southern sky to detect massive galaxy clusters out to the epoch of their formation using the Sunyaev-Zel'dovich (SZ) effect. This paper presents a catalog of the 26 most significant SZ cluster detections in the full survey region. The catalog includes 14 clusters which have been previously identified and 12 that are new discoveries. These clusters were identified in fields observed to two differing noise depths: 1500 deg^2 at the final SPT survey depth of 18 uK-arcmin at 150 GHz, and 1000 deg^2 at a depth of 54 uK-arcmin. Clusters were selected on the basis of their SZ signal-to-noise ratio (S/N) in SPT maps, a quantity which has been demonstrated to correlate tightly with cluster mass. The S/N thresholds were chosen to achieve a comparable mass selection across survey fields of both depths. Cluster redshifts were obtained with optical and infrared imaging and spectroscopy from a variety of ground- and space-based facilities. The redshifts range from 0.098 \leq z \leq 1.132 with a median of z_med = 0.40. The measured SZ S/N and redshifts lead to unbiased mass estimates ranging from 9.8 \times 10^14 M_sun/h_70 \leq M_200(rho_mean) \leq 3.1 \times 10^15 M_sun/h_70. Based on the SZ mass estimates, we find that none of the clusters are individually in significant tension with the LambdaCDM cosmological model. We also test for evidence of non-Gaussianity based on the cluster sample and find the data show no preference for non-Gaussian perturbations.

We report the results of an 87 square-degree point-source survey centered at R.A. 5h30m, decl. -55 deg. taken with the South Pole Telescope (SPT) at 1.4 and 2.0 mm wavelengths with arc-minute resolution and milli-Jansky depth. Based on the ratio of flux in the two bands, we separate the detected sources into two populations, one consistent with synchrotron emission from active galactic nuclei (AGN) and one consistent with thermal emission from dust. We present source counts for each population from 11 to 640 mJy at 1.4 mm and from 4.4 to 800 mJy at 2.0 mm. The 2.0 mm counts are dominated by synchrotron-dominated sources across our reported flux range; the 1.4 mm counts are dominated by synchroton-dominated sources above ~15 mJy and by dust-dominated sources below that flux level. We detect 141 synchrotron-dominated sources and 47 dust-dominated sources at S/N > 4.5 in at least one band. All of the most significantly detected members of the synchrotron-dominated population are associated with sources in previously published radio catalogs. Some of the dust-dominated sources are associated with nearby (z << 1) galaxies whose dust emission is also detected by the Infrared Astronomy Satellite (IRAS). However, most of the bright, dust-dominated sources have no counterparts in any existing catalogs. We argue that these sources represent the rarest and brightest members of the population commonly referred to as sub-millimeter galaxies (SMGs). Because these sources are selected at longer wavelengths than in typical SMG surveys, they are expected to have a higher mean redshift distribution and may provide a new window on galaxy formation in the early universe.

Based on the ultraviolet to far-infrared photometry already compiled and presented in a companion paper (Barro et al. 2011a, Paper I), we present a detailed SED analysis of nearly 80,000 IRAC 3.6+4.5 micron selected galaxies in the Extended Groth Strip. We estimate photometric redshifts, stellar masses, and star formation rates separately for each galaxy in this large sample. The catalog includes 76,936 sources with [3.6] < 23.75 (85% completeness level of the IRAC survey) over 0.48 square degrees. The typical photometric redshift accuracy is Delta z/(1+z)=0.034, with a catastrophic outlier fraction of just 2%. We quantify the systematics introduced by the use of different stellar population synthesis libraries and IMFs in the calculation of stellar masses. We find systematic offsets ranging from 0.1 to 0.4 dex, with a typical scatter of 0.3 dex. We also provide UV- and IR-based SFRs for all sample galaxies, based on several sets of dust emission templates and SFR indicators. We evaluate the systematic differences and goodness of the different SFR estimations using the deep FIDEL 70 micron data available in the EGS. Typical random uncertainties of the IR-bases SFRs are a factor of two, with non-negligible systematic effects at z$\gtrsim$1.5 observed when only MIPS 24 micron data is available. All data products (SEDs, postage stamps from imaging data, and different estimations of the photometric redshifts, stellar masses, and SFRs of each galaxy) described in this and the companion paper are publicly available, and they can be accessed through our the web-interface utility Rainbow-navigator

We present an IRAC 3.6+4.5 microns selected catalog in the Extended Groth Strip (EGS) containing photometry from the ultraviolet to the far-infrared and stellar parameters derived from the analysis of the multi-wavelength data. In this paper, we describe the method used to build coherent spectral energy distributions (SEDs) for all the sources. In a companion paper, we analyze those SEDs to obtain robust estimations of stellar parameters such as photometric redshifts, stellar masses, and star formation rates. The catalog comprises 76,936 sources with [3.6]<23.75 mag (85% completeness level of the IRAC survey in the EGS) over 0.48 square degrees. For approximately 16% of this sample, we are able to deconvolve the IRAC data to obtain robust fluxes for the multiple counterparts found in ground-based optical images. Typically, the SEDs of the IRAC sources in our catalog count with more than 15 photometric data points, spanning from the UV to the FIR. Approximately 95% and 90% of all IRAC sources are detected in the deepest optical and near-infrared bands. Only 10% of the sources have optical spectroscopy and redshift estimations. Almost 20% and 2% of the sources are detected by MIPS at 24 and 70 microns, respectively. We also cross-correlate our catalog with public X-ray and radio catalogs. Finally, we present the Rainbow Navigator public web-interface utility designed to browse all the data products resulting from this work, including images, spectra, photometry, and stellar parameters.

We present the J-band luminosity function of 1838 mid-infrared and X-ray selected AGNs in the redshift range 0<z<5.85. These luminosity functions are constructed by combining the deep multi-wavelength broad-band observations from the UV to the mid-IR of the NDWFS Bootes field with the X-ray observations of the XBootes survey and the spectroscopic observations of the same field by AGES. Our sample is primarily composed of IRAC-selected AGNs, targeted using modifications of the Stern et al.(2005) criteria, complemented by MIPS 24 microns and X-ray selected AGNs to alleviate the biases of IRAC mid-IR selection against z~4.5 quasars and AGNs faint with respect to their hosts. This sample provides an accurate link between low and high redshift AGN luminosity functions and does not suffer from the usual incompleteness of optical samples at z~3. We find that the space density of the brightest quasars strongly decreases from z=3 to z=0, while the space density of faint quasars is at least flat, and possibly increasing, over the same redshift range. At z>3 we observe a decrease in the space density of quasars of all brightnesses. We model the luminosity function by a double power-law and find that its evolution cannot be described by either pure luminosity or pure density evolution, but must be a combination of both. Our best-fit model has bright and faint power-law indices consistent with the low redshift measurements based on the 2QZ and 2SLAQ surveys and it generally agrees with the number of bright quasars predicted by other LFs at all redshifts. If we construct the QSO luminosity function using only the IRAC-selected AGNs, we find that the biases inherent to this selection method significantly modify the behavior of phi*(z) only for z<1 and have no significant impact upon the characteristic magnitude M*_J(z).

Using the South Pole Telescope (SPT), we have discovered the most massive known galaxy cluster at z > 1, SPT-CL J2106-5844. In addition to producing a strong Sunyaev-Zel'dovich effect signal, this system is a luminous X-ray source and its numerous constituent galaxies display spatial and color clustering, all indicating the presence of a massive galaxy cluster. VLT and Magellan spectroscopy of 18 member galaxies shows that the cluster is at z = 1.132^+0.002_-0.003. Chandra observations obtained through a combined HRC-ACIS GTO program reveal an X-ray spectrum with an Fe K line redshifted by z = 1.18 +/- 0.03. These redshifts are consistent with galaxy colors in extensive optical, near-infrared, and mid-infrared imaging. SPT-CL J2106-5844 displays extreme X-ray properties for a cluster, having a core-excluded temperature of kT = 11.0^+2.6_-1.9 keV and a luminosity (within r_500) of L_X (0.5 - 2.0 keV) = (13.9 +/- 1.0) x 10^44 erg/s. The combined mass estimate from measurements of the Sunyaev-Zel'dovich effect and X-ray data is M_200 = (1.27 +/- 0.21) x 10^15 M_sun. The discovery of such a massive gravitationally collapsed system at high redshift provides an interesting laboratory for galaxy formation and evolution, and is a powerful probe of extreme perturbations of the primordial matter density field. We discuss the latter, determining that, under the assumption of LambdaCDM cosmology with only Gaussian perturbations, there is only a 7% chance of finding a galaxy cluster similar to SPT-CL J2106-5844 in the 2500 deg^2 SPT survey region, and that only one such galaxy cluster is expected in the entire sky.

We report the X-ray detection of two z>1.4 infrared-selected galaxy clusters from the IRAC Shallow Cluster Survey (ISCS). We present new data from the Hubble Space Telescope and the W. M. Keck Observatory that spectroscopically confirm cluster ISCS J1432.4+3250 at z=1.49, the most distant of 18 confirmed z>1 clusters in the ISCS to date. We also present new spectroscopy for ISCS J1438.1+3414, previously reported at z = 1.41, and measure its dynamical mass. Clusters ISCS J1432.4+3250 and ISCS J1438.1+3414 are detected in 36ks and 143ks Chandra exposures at significances of 5.2 sigma and 9.7 sigma, from which we measure total masses of log(M_{200,Lx}/Msun) = 14.4 +/- 0.2 and 14.35^{+0.14}_{-0.11}, respectively. The consistency of the X-ray and dynamical properties of these high redshift clusters further demonstrates that the ISCS is robustly detecting massive clusters to at least z = 1.5.

We report the discovery of a six-month-long mid-infrared transient, SDWFS-MT-1 (aka SN 2007va), in the Spitzer Deep, Wide-Field Survey of the NOAO Deep Wide-Field Survey Bootes field. The transient, located in a z=0.19 low luminosity (M_[4.5]~-18.6 mag, L/L_MilkyWay~0.01) metal-poor (12+log(O/H)~7.8) irregular galaxy, peaked at a mid-infrared absolute magnitude of M_[4.5]~-24.2 in the 4.5 micron Spitzer/IRAC band and emitted a total energy of at least 10^51 ergs. The optical emission was likely fainter than the mid-infrared, although our constraints on the optical emission are poor because the transient peaked when the source was "behind" the Sun. The Spitzer data are consistent with emission by a modified black body with a temperature of ~1350 K. We rule out a number of scenarios for the origin of the transient such as a Galactic star, AGN activity, GRB, tidal disruption of a star by a black hole and gravitational lensing. The most plausible scenario is a supernova exploding inside a massive, optically thick circumstellar medium, composed of multiple shells of previously ejected material. If the proposed scenario is correct, then a significant fraction (~10%) of the most luminous supernova may be self-enshrouded by dust not only before but also after the supernova occurs. The spectral energy distribution of the progenitor of such a supernova would be a slightly cooler version of eta Carina, peaking at 20-30 microns.

Combined Spitzer, Chandra, XMM-Newton, and VLA observations of the giant radio galaxy NGC 1316 (Fornax A) show a radio jet and X-ray cavities from AGN outbursts likely triggered by a merger with a late-type galaxy at least 0.4 Gyr ago. We detect a weak nucleus with an SED typical of a low-luminosity AGN with a bolometric luminosity of 2.4x10^42 erg/s. The Spitzer IRAC and MIPS images show dust emission strongest in regions with little or no radio emission. The large infrared luminosity relative to the galaxy's K-band luminosity implies an external origin for the dust. The dust mass implies that the merger spiral galaxy had a stellar mass of 1-6x10^10 M_sun and a gas mass of 2-4x10^9 M_sun. Chandra images show a small ~15"(1.6 kpc) cavity coincident with the radio jet, while the XMM-Newton image shows two large X-ray cavities lying 320"(34.8 kpc) east and west of the nucleus, each approximately 230"(25 kpc) in radius. The radio lobes lie at radii of 14.3'(93 kpc) and 15.6'(101 kcp), more distant from the nucleus than the detected X-ray cavities. The relative morphology of the large scale 1.4 GHz and X-ray emission suggests they were products of two distinct outbursts, an earlier one creating the radio lobes and a later one producing the X-ray cavities. Alternatively, if a single outburst created both the X-ray cavities and the radio lobes, this would require that the morphology is not fully defined by the 1.4 GHz emission. For the more likely two outburst scenario, we use the buoyancy rise times to estimate an age for the more recent outburst of 0.1 Gyr and use the associated PV work done by the expanding plasma to create the X-ray cavities to estimate the outburst's energy of 10^58 ergs. The present size and location of the large radio lobes implies that the AGN outburst that created them happened ~0.4 Gyr ago and released ~5x10^58 ergs. (abridged)

We report the spectroscopic confirmation of SPT-CL J0546-5345 at <z> = 1.067. To date this is the most distant cluster to be spectroscopically confirmed from the 2008 South Pole Telescope (SPT) catalog, and indeed the first z > 1 cluster discovered by the Sunyaev-Zel'dovich Effect (SZE). We identify 21 secure spectroscopic members within 0.9 Mpc of the SPT cluster position, 18 of which are quiescent, early-type galaxies. From these quiescent galaxies we obtain a velocity dispersion of 1179^{+232}_{-167} km/s, ranking SPT-CL J0546-5345 as the most dynamically massive cluster yet discovered at z > 1. Assuming that SPT-CL J0546-5345 is virialized, this implies a dynamical mass of M_200 = 1.0^{+0.6}_{-0.4} x 10^{15} Msun, in agreement with the X-ray and SZE mass measurements. Combining masses from several independent measures leads to a best-estimate mass of M_200 = (7.95 +/- 0.92) x 10^{14} Msun. The spectroscopic confirmation of SPT-CL J0546-5345, discovered in the wide-angle, mass-selected SPT cluster survey, marks the onset of the high redshift SZE-selected galaxy cluster era.

We have observed warm molecular hydrogen in two nearby edge-on disk galaxies, NGC 4565 and NGC 5907, using the Spitzer high-resolution infrared spectrograph. The 0-0 S(0) 28.2 micron and 0-0 S(1) 17.0 micron pure rotational lines were detected out to 10 kpc from the center of each galaxy on both sides of the major axis, and in NGC 4565 the S(0) line was detected at r = 15 kpc on one side. This location lies beyond a steep drop in the radio continuum emission from cosmic rays in the disk. Despite indications that star formation activity decreases with radius, the H2 excitation temperature and the ratio of the H2 line and the far-IR luminosity surface densities, Sigma_L(H2}/Sigma_L(TIR}, change very little as a function of radius, even into the diffuse outer region of the disk of NGC 4565. This suggests that the source of excitation of the H2 operates over a large range of radii, and is broadly independent of the strength and relative location of UV emission from young stars. Although excitation in photodissociation regions is the most common explanation for the widespread H2 emission, cosmic ray heating or shocks cannot be ruled out. The inferred mass surface densities of warm molecular hydrogen in both edge-on galaxies differ substantially, being 4(-60) M_solar/pc^2 and 3(-50) M_solar/pc^2 at r = 10 kpc for NGC 4565 and NGC 5907, respectively. The higher values represent very unlikely point-source upper limits. The point source case is not supported by the observed emission distribution in the spectral slits. These mass surface densities cannot support the observed rotation velocities in excess of 200 km/s. Therefore, warm molecular hydrogen cannot account for dark matter in these disk galaxies, contrary to what was implied by a previous ISO study of the nearby edge-on galaxy NGC 891.

Using the Spitzer Space Telescope, we have obtained rest frame 9-16mu spectra of 11 quasars and 9 radio galaxies from the 3CRR catalog at redshifts 1.0<z<1.4. This complete flux-limited 178MHz-selected sample is unbiased with respect to orientation and therefore suited to study orientation-dependent effects in the most powerful active galactic nuclei (AGN). The mean radio galaxy spectrum shows a clear silicate absorption feature (tau_9.7mu = 1.1) whereas the mean quasar spectrum shows silicates in emission. The mean radio galaxy spectrum matches a dust-absorbed mean quasar spectrum in both shape and overall flux level. The data for individual objects conform to these results. The trend of the silicate depth to increase with decreasing core fraction of the radio source further supports that for this sample, orientation is the main driver for the difference between radio galaxies and quasars, as predicted by AGN unification. However, comparing our high-z sample with lower redshift 3CRR objects reveals that the absorption of the high-z radio galaxy MIR continuum is lower than expected from a scaled up version of lower luminosity sources, and we discuss some effects that may explain these trends.

We present a study of a large, statistically complete sample of star-forming dwarf galaxies using mid-infrared observations from the {\it Spitzer Space Telescope}. The relationships between metallicity, star formation rate (SFR) and mid-infrared color in these systems show that the galaxies span a wide range of properties. However, the galaxies do show a deficit of 8.0 \um\ polycyclic aromatic hydrocarbon emission as is apparent from the median 8.0 \um\ luminosity which is only 0.004 \lstarf\ while the median $B$-band luminosity is 0.05 \lstarb. Despite many of the galaxies being 8.0 \um\ deficient, there is about a factor of 4 more extremely red galaxies in the [3.6] $-$ [8.0] color than for a sample of normal galaxies with similar optical colors. We show correlations between the [3.6] $-$ [8.0] color and luminosity, metallicity, and to a lesser extent SFRs that were not evident in the original, smaller sample studied previously. The luminosity--metallicity relation has a flatter slope for dwarf galaxies as has been indicated by previous work. We also show a relationship between the 8.0 \um\ luminosity and the metallicity of the galaxy which is not expected given the competing effects (stellar mass, stellar population age, and the hardness of the radiation field) that influence the 8.0 \um\ emission. This larger sample plus a well-defined selection function also allows us to compute the 8.0 \um\ luminosity function and compare it with the one for the local galaxy population. Our results show that below 10$^{9}$ $L$\solar, nearly all the 8.0 \um\ luminosity density of the local universe arises from dwarf galaxies that exhibit strong \ha\ emission -- i.e., 8.0 \um\ and \ha\ selection identify similar galaxy populations despite the deficit of 8.0 \um\ emission observed in these dwarfs.

This paper analyzes a sample of 489 Spitzer/IRAC sources in the Extended Groth Strip whose spectral energy distributions fit a red power law from 3.6 to 8.0 \micron. The median for sources with known redshift is <z>=1.6. Though all or nearly all of the sample are likely to be active galactic nuclei, only 33% were detected in the EGS X-ray survey (AEGIS-X) using 200 ks Chandra observations. The detected sources are X-ray luminous with L_X > 10^43 erg/s and moderately to heavily obscured with N_H > 10^22 cm^-2. Stacking the X-ray-undetected sample members yields a statistically significant X-ray signal, suggesting that they are on average more distant or more obscured than sources with X-ray detections. The ratio of X-ray to mid-infrared fluxes suggests that a substantial fraction of the sources undetected in X-rays are obscured at the Compton-thick level, in contrast to the X-ray-detected sources, all of which appear to be Compton-thin. For the X-ray-detected power-law sources with redshifts, an X-ray luminosity L_X ~ 10^44 erg/s marks a transition between low-luminosity, blue sources dominated by the host galaxy to high-luminosity, red power-law sources dominated by nuclear activity. X-ray-to-optical ratios, infrared variability, and 24 micron properties of the sample are consistent with the identification of infrared power-law sources as active nuclei, but a rough estimate is that only 22% of AGNs are selected by the power law criteria. Comparison of the power-law selection technique and various IRAC color criteria for identifying AGNs confirms that high-redshift samples selected via simple IRAC colors may be heavily contaminated by starlight-dominated objects.

We use the multi-epoch, mid-infrared Spitzer Deep, Wide-Field Survey to investigate the variability of 474,179 objects in 8.1 deg^2 of the NDWFS Bootes field. We perform a Difference Image Analysis of the four available epochs between 2004 and 2008, focusing on the deeper 3.6 and 4.5 micron bands. We find that 1.1% of the studied sample meet our standard selection criteria for being classed as a variable source. We require that the 3.6 and 4.5 micron light-curves are strongly correlated (r>0.8) and that their joint variance exceeds that for all sources with the same magnitude by 2 sigma. We then examine the mid-IR colors of the variable sources and match them with X-ray sources from the XBootes survey, radio catalogs, 24 micron-selected AGN candidates, and spectroscopically identified AGNs from the AGN and Galaxy Evolution Survey (AGES). Based on their mid-IR colors, most of the variable sources are AGNs (76%), with smaller contributions from stars (11%), galaxies (6%), and unclassified objects. Most of the stellar, galaxy and unclassified sources are false positives. For our standard selection criteria, 11-12% of the mid-IR counterparts to X-ray sources, 24 micron-selected AGN candidates and spectroscopically identified AGNs show variability. Mid-IR AGN variability can be well described by a single power-law structure function with a power-law index of 0.5 at both 3.6 and 4.5 microns, and an amplitude of 0.1 mag on rest-frame time scales of 2 years. The variability amplitude is higher for shorter rest-frame wavelengths and lower luminosities. (Abridged)

We report the detection and successful modeling of the unusual 9.7\mum Si--O stretching silicate emission feature in the type 1 (i.e. face-on) LINER nucleus of M81. Using the Infrared Spectrograph (IRS) instrument on Spitzer, we determine the feature in the central 230 pc of M81 to be in strong emission, with a peak at ~10.5\mum. This feature is strikingly different in character from the absorption feature of the galactic interstellar medium, and from the silicate absorption or weak emission features typical of galaxies with active star formation. We successfully model the high signal-to-noise ratio IRS spectra with porous silicate dust using laboratory-acquired mineral spectra. We find that the most probable fit uses micron-sized, porous grains of amorphous silicate and graphite. In addition to silicate dust, there is weak PAH emission present (particularly at 11.3\mum, arising from the C--H out-of-plane bending vibration of relatively large PAHs of ~500--1000 C atoms) whose character reflects the low-excitation AGN environment, with some evidence that small PAHs of ~100--200 C atoms (responsible for the 7.7\mum C--C stretching band) in the immediate vicinity of the nucleus have been preferentially destroyed. (abstract continues)

We have identified a sample of cool field brown dwarf candidates using IRAC data from the Spitzer Deep, Wide-Field Survey (SDWFS). The candidates were selected from 400,000 SDWFS sources with [4.5] <= 18.5 mag and required to have [3.6]-[4.5] >= 1.5 and [4.5] - [8.0] <= 2.0 on the Vega system. The first color requirement selects objects redder than all but a handful of presently known brown dwarfs with spectral classes later than T7, while the second eliminates 14 probable reddened AGN. Optical detection of 4 of the remaining 18 sources implies they are likely also AGN, leaving 14 brown dwarf candidates. For two of the brightest candidates (SDWFS J143524.44+335334.6 and SDWFS J143222.82+323746.5), the spectral energy distributions including near-infrared detections suggest a spectral class of ~ T8. The proper motion is < 0.25 "/yr, consistent with expectations for a luminosity inferred distance of >70 pc. The reddest brown dwarf candidate (SDWFS J143356.62+351849.2) has [3.6] - [4.5]=2.24 and H - [4.5] > 5.7, redder than any published brown dwarf in these colors, and may be the first example of the elusive Y-dwarf spectral class. Models from Burrows et al. (2003) predict larger numbers of cool brown dwarfs should be found for a Chabrier (2003) mass function. Suppressing the model [4.5] flux by a factor of two, as indicated by previous work, brings the Burrows models and observations into reasonable agreement. The recently launched Wide-field Infrared Survey Explorer (WISE) will probe a volume ~40x larger and should find hundreds of brown dwarfs cooler than T7.

We present high-resolution 345 GHz interferometric observations of two extreme luminous (L_{IR}>10^{13} L_sun), submillimetre-selected galaxies (SMGs) in the COSMOS field with the Submillimeter Array (SMA). Both targets were previously detected as unresolved point-sources by the SMA in its compact configuration, also at 345 GHz. These new data, which provide a factor of ~3 improvement in resolution, allow us to measure the physical scale of the far-infrared in the submillimetre directly. The visibility functions of both targets show significant evidence for structure on 0.5-1 arcsec scales, which at z=1.5 translates into a physical scale of 5-8 kpc. Our results are consistent with the angular and physical scales of two comparably luminous objects with high-resolution SMA followup, as well as radio continuum and CO sizes. These relatively compact sizes (<5-10 kpc) argue strongly for merger-driven starbursts, rather than extended gas-rich disks, as the preferred channel for forming SMGs. For the most luminous objects, the derived sizes may also have important physical consequences; under a series of simplifying assumptions, we find that these two objects in particular are forming stars close to or at the Eddington limit for a starburst.

We present new observational determinations of the evolution of the 2-10keV X-ray luminosity function (XLF) of AGN. We utilise data from a number of surveys including both the 2Ms Chandra Deep Fields and the AEGIS-X 200ks survey, enabling accurate measurements of the evolution of the faint end of the XLF. We combine direct, hard X-ray selection and spectroscopic follow-up or photometric redshift estimates at z<1.2 with a rest-frame UV colour pre-selection approach at higher redshifts to avoid biases associated with catastrophic failure of the photometric redshifts. Only robust optical counterparts to X-ray sources are considered using a likelihood ratio matching technique. A Bayesian methodology is developed that considers redshift probability distributions, incorporates selection functions for our high redshift samples, and allows robust comparison of different evolutionary models. We find that the XLF retains the same shape at all redshifts, but undergoes strong luminosity evolution out to z~1, and an overall negative density evolution with increasing redshift, which thus dominates the evolution at earlier times. We do not find evidence that a Luminosity-Dependent Density Evolution, and the associated flattening of the faint-end slope, is required to describe the evolution of the XLF. We find significantly higher space densities of low-luminosity, high-redshift AGN than in prior studies, and a smaller shift in the peak of the number density to lower redshifts with decreasing luminosity. The total luminosity density of AGN peaks at z=1.2+/-0.1, but there is a mild decline to higher redshifts. We find >50% of black hole growth takes place at z>1, with around half in Lx<10^44 erg/s AGN.

[abridged] We present constraints on the nature of the first galaxies selected at 350 microns. The sample includes galaxies discovered in the deepest blank-field survey at 350 microns (in the Bootes Deep Field) and also later serendipitous detections in the Lockman Hole. Spectral energy distribution templates are fit to identified counterparts, and the sample is found to comprise IR-luminous galaxies at 1<z<3 predominantly powered by star formation. The first spectrum of a 350-micron-selected galaxy provides an additional confirmation, showing prominent dust grain features typically associated with star-forming galaxies. Compared to submillimeter galaxies selected at 850 and 1100 microns, galaxies selected at 350 microns have a similar range of far-infrared color temperatures. However, no 350-micron-selected sources are reliably detected at 850 or 1100 microns. Galaxies in our sample with redshifts 1<z<2 show a tight correlation between the far- and mid-infrared flux densities, but galaxies at higher redshifts show a large dispersion in their mid- to far-infrared colors. The 350-micron flux densities (15<S(350)<40 mJy) place these objects near the Herschel/SPIRE 350-micron confusion threshold, with the lower limit on the star formation rate density suggesting the bulk of the 350-micron contribution will come from less luminous infrared sources and normal galaxies.

We present 20 band photometry from the far-IR to X-ray in the Spitzer IRAC dark field. The bias for the near-IR camera on Spitzer is calibrated by observing a ~20 arcminute diameter "dark" field near the north ecliptic pole roughly every two-to-three weeks throughout the mission duration of Spitzer. The field is unique for its extreme depth, low background, high quality imaging, time-series information, and accompanying photometry including data taken with Akari, Palomar, MMT, KPNO, Hubble, and Chandra. This serendipitous survey contains the deepest mid-IR data taken to date. This dataset is well suited for studies of intermediate redshift galaxy clusters, high redshift galaxies, the first generation of stars, and the lowest mass brown dwarfs, among others. This paper provides a summary of the data characteristics and catalog generation from all bands collected to date as well as a discussion of photometric redshifts and initial and expected science results and goals. To illustrate the scientific potential of this unique dataset, we also present here IRAC color color diagrams.

We present results from a continuing interferometric survey of high-redshift submillimeter galaxies with the Submillimeter Array, including high-resolution (beam size ~2 arcsec) imaging of eight additional AzTEC 1.1mm selected sources in the COSMOS Field, for which we obtain six reliable (peak S/N>5 or peak S/N>4 with multiwavelength counterparts within the beam) and two moderate significance (peak S/N>4) detections. When combined with previous detections, this yields an unbiased sample of millimeter-selected SMGs with complete interferometric followup. With this sample in hand, we (1) empirically confirm the radio-submillimeter association, (2) examine the submillimeter morphology - including the nature of submillimeter galaxies with multiple radio counterparts and constraints on the physical scale of the far infrared - of the sample, and (3) find additional evidence for a population of extremely luminous, radio-dim submillimeter galaxies that peaks at higher redshift than previous, radio-selected samples. In particular, the presence of such a population of high-redshift sources has important consequences for models of galaxy formation - which struggle to account for such objects even under liberal assumptions - and dust production models given the limited time since the Big Bang.

Using the Infrared Spectrograph on board the Spitzer Space Telescope, we present low-resolution (64 < lambda / dlambda < 124), mid-infrared (20-38 micron) spectra of 23 high-redshift ULIRGs detected in the Bootes field of the NOAO Deep Wide-Field Survey. All of the sources were selected to have 1) fnu(24 micron) > 0.5 mJy; 2) R-[24] > 14 Vega mag; and 3) a prominent rest-frame 1.6 micron stellar photospheric feature redshifted into Spitzer's 3-8 micron IRAC bands. Of these, 20 show emission from polycyclic aromatic hydrocarbons (PAHs), usually interpreted as signatures of star formation. The PAH features indicate redshifts in the range 1.5 < z < 3.0, with a mean of <z>=1.96 and a dispersion of 0.30. Based on local templates, these sources have extremely large infrared luminosities, comparable to that of submillimeter galaxies. Our results confirm previous indications that the rest-frame 1.6 micron stellar bump can be efficiently used to select highly obscured starforming galaxies at z~2, and that the fraction of starburst-dominated ULIRGs increases to faint 24 micron flux densities. Using local templates, we find that the observed narrow redshift distribution is due to the fact that the 24 micron detectability of PAH-rich sources peaks sharply at z = 1.9. We can analogously explain the broader redshift distribution of Spitzer-detected AGN-dominated ULIRGs based on the shapes of their SEDs. Finally, we conclude that z~2 sources with a detectable 1.6 micron stellar opacity feature lack sufficient AGN emission to veil the 7.7 micron PAH band.

The Spitzer Deep, Wide-Field Survey (SDWFS) is a four-epoch infrared survey of ten square degrees in the Bootes field of the NOAO Deep Wide-Field Survey using the IRAC instrument on the Spitzer Space Telescope. SDWFS, a Cycle four Spitzer Legacy project, occupies a unique position in the area-depth survey space defined by other Spitzer surveys. The four epochs that make up SDWFS permit -- for the first time -- the selection of infrared-variable and high proper motion objects over a wide field on timescales of years. Because of its large survey volume, SDWFS is sensitive to galaxies out to z~3 with relatively little impact from cosmic variance for all but the richest systems. The SDWFS datasets will thus be especially useful for characterizing galaxy evolution beyond z~1.5. This paper explains the SDWFS observing strategy and data processing, presents the SDWFS mosaics and source catalogs, and discusses some early scientific findings. The publicly-released, full-depth catalogs contain 6.78, 5.23, 1.20, and 0.96 x 10e5 distinct sources detected to the average 5-sigma, 4" diameter, aperture-corrected limits of 19.77, 18.83, 16.50, and 15.82 Vega mag at 3.6, 4.5, 5.8, and 8.0 micron, respectively. The SDWFS number counts and color-color distribution are consistent with other, earlier Spitzer surveys. At the 6 min integration time of the SDWFS IRAC imaging, more than 50% of isolated FIRST radio sources and more than 80% of on-axis XBootes sources are detected out to 8.0 micron. Finally, we present the four highest proper motion IRAC-selected sources identified from the multi-epoch imaging, two of which are likely field brown dwarfs of mid-T spectral class.

We present $u'g'R$ optical images taken with the MMT/Megacam and the Subaru/Suprime of the Extended Groth Strip survey. The total survey covers an area of about $\sim 1$ degree$^2$, including four sub-fields and is optimized for the study of galaxies at $z\sim3$. Our methods for photometric calibration in AB magnitudes, the limiting magnitude and the galaxy number count are described. A sample of 1642 photometrically selected candidate LBGs to an apparent $R_{AB}$ magnitude limit of 25.0 is present. The average sky surface density of our LBGs sample is $\sim$ 1.0 arcmin$^{-2}$, slightly higher than the previous finding.

Differential K-band luminosity functions (LFs) are presented for a complete sample of 1613 nearby bright galaxies segregated by visible morphology. The LF for late-type spirals follows a power law that rises towards low luminosities whereas the LFs for ellipticals, lenticulars and bulge-dominated spirals are peaked and decline toward both higher and lower luminosities. Each morphological type (E, S0, S0/a-Sab, Sb-Sbc, Sc-Scd) contributes approximately equally to the overall K-band luminosity density of galaxies in the local universe. Type averaged bulge/disk ratios are used to subtract the disk component leading to the prediction that the K-band LF for bulges is bimodal with ellipticals dominating the high luminosity peak, comprising 60% of the bulge luminosity density in the local universe with the remaining 40% contributed by lenticulars and the bulges of spirals. Overall, bulges contribute 30% of the galaxy luminosity density at K in the local universe with spiral disks making up the remainder. If bulge luminosities indicate central black hole masses, then our results predict that the black hole mass function is also bimodal.

UV continuum and mid-IR emission constitute two widely used star formation indicators at intermediate and high redshifts. We study 2430 galaxies with z<1.4 in the Extended Groth Strip with MIPS 24 mic observations from FIDEL, spectroscopy from DEEP2, and UV, optical, and near-IR photometry from AEGIS. The data are coupled with stellar population models and Bayesian SED fitting to estimate dust-corrected SFRs. In order to probe the dust heating from stellar populations of various ages, the derived SFRs were averaged over various timescales--from 100 Myr for "current" SFR to 1--3 Gyr for long-timescale SFRs. These SED-based UV/optical SFRs are compared to total infrared luminosities extrapolated from 24 mic observations. We find that for the blue, actively star forming galaxies the correlation between the IR luminosity and the UV/optical SFR shows a decrease in scatter when going from shorter to longer SFR-averaging timescales. We interpret this as the greater role of intermediate age stellar populations in heating the dust than what is typically assumed. This holds over the entire redshift range. Many so-called green valley galaxies are simply dust-obscured actively star-forming galaxies. However, there exist 24 mic-detected galaxies, some with L>10^11 L_sun, yet with little current star formation. For them a reasonable amount of dust absorption of stellar light is sufficient to produce the observed levels of IR. In our sample optical and X-ray AGNs do not contribute on average more than ~50% to the mid-IR luminosity, and we see no evidence for a large population of "IR excess" galaxies (Abridged).

We present infrared luminosities, star formation rates, colors, morphologies, locations, and AGN properties of 24 micron-detected sources in photometrically detected high-redshift clusters in order to understand the impact of environment on star formation and AGN evolution in cluster galaxies. We use three newly-identified z=1 clusters selected from the IRAC dark field; the deepest ever mid-IR survey with accompanying, 14 band multiwavelength data including deep HST imaging and deep wide-area Spitzer MIPS 24 micron imaging. We find 90 cluster members with MIPS detections within two virial radii of the cluster centers, of which 17 appear to have spectral energy distributions dominated by AGN and the rest dominated by star formation. We find that 43 of the star forming are luminous infrared galaxies (LIRGs). The majority of sources (81%) are spirals or irregulars. A large fraction (at least 25%) show obvious signs of interactions. The MIPS -detected member galaxies have varied spatial distributions as compared to the MIPS-undetected members with one of the three clusters showing SF galaxies being preferentially located on the cluster outskirts, while the other 2 clusters show no such trend. Both the AGN fraction and the summed SFR of cluster galaxies increases from z=0 to 1, at a rate that is a few times faster in clusters than over the same redshift range in the field. Cluster environment does have an effect on the evolution of both AGN fraction and SFR from redshift one to the present, but does not effect the infrared luminosities or morphologies of the MIPS sample. Star formation happens in the same way regardless of environment making MIPS sources look the same in the cluster and field, however the cluster environment does encourage a more rapid evolution with time as compared to the field.

We have carried out a search for galaxies at z ~ 7-10 in ~14.4 sq. arcmin of new NICMOS parallel imaging taken in the Great Observatories Origins Deep Survey (GOODS, 5.9 sq. arcmin), the Cosmic Origins Survey (COSMOS, 7.2 sq. arcmin), and SSA22 (1.3 sq. arcmin). These images reach 5 sigma sensitivities of J110 = 26.0-27.5 (AB), and combined they increase the amount of deep near-infrared data by more than 60% in fields where the investment in deep optical data has already been made. We find no z>7 candidates in our survey area, consistent with the Bouwens et al. (2008) measurements at z~7 and 9 (over 23 sq. arcmin), which predict 0.7 galaxies at z~7 and <0.03 galaxies at z~9. We estimate that 10-20% of z>7 galaxies are missed by this survey, due to incompleteness from foreground contamination by faint sources. For the case of luminosity evolution, assuming a Schecter parameterization with a typical phi* = 10^-3 Mpc^-3, we find M* > -20.0 for z~7 and M* > -20.7 for z~9 (68% confidence). This suggests that the downward luminosity evolution of LBGs continues to z~7, although our result is marginally consistent with the z~6 LF of Bouwens et al.(2006, 2007). In addition we present newly-acquired deep MMT/Megacam imaging of the z~9 candidate JD2325+1433, first presented in Henry et al. (2008). The resulting weak but significant detection at i' indicates that this galaxy is most likely an interloper at z~2.7.

Differential 2.2um (K-band) luminosity functions are presented for a complete sample of 1570 nearby Vgsr < 3000 km/s, where Vgsr is the velocity measured with respect to the Galactic standard of rest), bright (K < 10 mag), galaxies segregated by visible morphology. The K-band luminosity function for late-type spirals follows a power law that rises towards low luminosities whereas the K-band luminosity functions for ellipticals, lenticulars and bulge-dominated spirals are peaked with a fall off at both high and low luminosities. However, each morphological type (E, S0, S0/a-Sab, Sb-Sbc, Sc-Scd) contributes approximately equally to the overall K-band luminosity density in the local universe, and by inference, the stellar mass density as well.

In the paradigm of hierarchical galaxy formation, luminous radio galaxies mark mass assembly peaks that should contain clusters of galaxies. Observations of the z=1.53 quasar 3C270.1 with the Spitzer Space Telescope at 3.6-24 micron and with the 6.5-m MMT in the z'- and Y-bands allow detection of potential cluster members via photometric redshifts. Compared with nearby control fields, there is an excess of 11 extremely red objects (EROs) at 1.33 < z_phot < 1.73, consistent with a proto-cluster around the quasar. The spectral energy distributions (SEDs) of 3/4 of the EROs are better fitted with passive elliptical galaxies than withdust-reddened starbursts, and of four sources well-detected on an archival HST snapshot image, all have undisturbed morphologies. However, one ERO, not covered by the HST image, is a double source with 0.8" separation on the z' image and a marginal (2sigma) 24 micron detection indicating a dust-enshrouded starburst. The EROs are more luminous than L* (H = -23.6 AB mag at z=1.5).

We present the results of a Spitzer IRAC and MIPS 24 micron study of extended Lyman-alpha clouds (or Lyman-alpha Blobs, LABs) within the SSA22 filamentary structure at z = 3.09. We detect 6/26 LABs in all IRAC filters, four of which are also detected at 24 micron, and find good correspondence with the 850 micron measurements of Geach et al. 2005. An analysis of the rest-frame ultraviolet, optical, near- and mid-infrared colors reveals that these six systems exhibit signs of nuclear activity (AGN)and/or extreme star formation. Notably, they have properties that bridge galaxies dominated by star formation (Lyman-break galaxies; LBGs) and those with AGNs (LBGs classified as QSOs). The LAB systems not detected in all four IRAC bands, on the other hand, are, as a group, consistent with pure star forming systems, similar to the majority of the LBGs within the filament. These results indicate that the galaxies within LABs do not comprise a homogeneous population, though they are also consistent with scenarios in which the gas halos are ionized through a common mechanism such as galaxy-scale winds driven by the galaxies within them, or gravitational heating of the collapsing cloud itself.

We combine Spitzer IRAC mid-infrared (MIR) and Chandra X-ray observations of the dominant galaxies NGC6872 and NGC6876 in the Pavo group with archival optical and HI data to study interaction-induced star formation. In spiral galaxy NGC6872, 8.0 and 5.8 micron nonstellar emission having colors consistent with polycyclic aromatic hydrocarbons (PAHs) is concentrated in clumps in three regions: in a 5 kpc radius outer ring about the center of the spiral galaxy, in a bridge of emission connecting NGC6872's northern spiral arm to IC4970, and along the full extent of NGC6872's tidal arms. PAH emission is correlated with young star clusters and dense HI regions. We find no strong differences in the MIR colors of star-forming regions in the spiral galaxy NGC6872 as a function of position relative to the tidally interacting companion galaxy IC4970. We find 11 very luminous X-ray sources (>~ (0.5 - 5) x 10^{39} ergs/s) clustered to the southwest in NGC6872, near bright star-forming regions. In NGC6872's tidal features, young star clusters form at the boundaries of diffuse X-ray gas, suggesting that stars form as gas stripped by the interactions cools. The nucleus of NGC6872 is a weak X-ray point source (0.5-8 keV luminosity of 8.5 x 10^{39} ergs/s), but there is little evidence in the inner 1 kpc of NGC6872 for PAH emission from recent star formation or nuclear activity. However, a 4 kpc `stream', leading from the outer ring of NGC6872 to the nucleus, may signal transport of interstellar matter into NGC6872's nuclear region. Nonstellar emission, consistent with PAH emission, is also found in the central region of elliptical galaxy NGC6877, companion to dominant Pavo group elliptical NGC6876. However, in the central region of NGC6876, the dust emission is more likely due to silicate emission from old AGB stars.

We present high resolution submillimeter interferometric imaging of two of the brightest high-redshift submillimeter galaxies known: GN20 and AzTEC1 at 0.8 and 0.3 arcsec resolution respectively. Our data - the highest resolution submillimeter imaging of high redshift sources accomplished to date - was collected in three different array configurations: compact, extended, and very extended. We derive angular sizes of 0.6 and 1.0 arcsec for GN20 and 0.3 and 0.4 arcsec for AzTEC1 from modeling their visibility functions as a Gaussian and elliptical disk respectively. Because both sources are B-band dropouts, they likely lie within a relatively narrow redshift window around z~4, which indicates their angular extent corresponds to physical scales of 4-8 and 1.5-3 kpc respectively for the starburst region. By way of a series of simple assumptions, we find preliminary evidence that these hyperluminous starbursts - with star formation rates >1000 $M_\odot$ yr$^{-1}$ - are radiating at or close to their Eddington limit. Should future high resolution observations indicate that these two objects are typical of a population of high redshift Eddington-limited starbursts, this could have important consequences for models of star formation and feedback in extreme environments.

High-resolution submm imaging of the HzRG, 4C60.07, at z=3.8, has revealed two dusty components. Spitzer imaging shows that one of these components (B) is coincident with an extremely red AGN, offset by ~4" (~30 kpc) from the HzRG core. The other submm component (A) - resolved by our beam and devoid of emission at 3.6-8.0um - lies between B and the HzRG core. Since the HzRG was discovered via its young, steep-spectrum lobes and their creation was likely triggered by the interaction, we argue that we are witnessing an early-stage merger, prior to its eventual equilibrium state. The interaction is between the host galaxy of an actively-fueled BH, and a gas-rich starburst/AGN (B) marked by the compact submm component and coincident with broad CO emission. `A' is a plume of cold, dusty gas, associated with a narrow (~150 km/s) CO feature, and may represent a short-lived tidal structure. It has been claimed that HzRGs and SMGs differ only in the activity of their AGNs, but such complex submm morphologies are seen only rarely amongst SMGs. Our study has important implications: where a galaxy's gas is not aligned with its central BH, CO may be an unreliable probe of dynamical mass, affecting work on the co-assembly of BHs and spheroids. Our data support the picture wherein close binary AGN are induced by mergers. They also raise the possibility that some supposedly jet-induced starbursts may have formed co-evally with (yet independently of) the radio jets, both triggered by the same interaction. We note that the HzRG host would have gone unnoticed without its jets/companion, so there may be many other unseen BHs at high redshift, lost in the sea of ~5 x 10^8 similarly bright IRAC sources - sufficiently massive to drive a >10^27-W radio source, yet practically invisible unless actively fueled (abridged).

Using the Spitzer Space Telescope, we have obtained 3.6--24 micron photometry of 38 radio galaxies and 24 quasars from the 3CR catalog at redshift 1<z<2.5. This 178 MHz-selected sample is unbiased with respect to orientation and therefore suited to study orientation-dependent effects in the most powerful active galactic nuclei (AGN). Quasar and radio galaxy subsamples matched in isotropic radio luminosity are compared. The quasars all have similar spectral energy distributions (SEDs), nearly constant in nu F_nu through the rest 1.6-10 micron range, consistent with a centrally heated dust distribution which outshines the host galaxy contribution. The radio galaxy SEDs show larger dispersion, but the mean radio galaxy SED declines from rest 1.6 to 3 micron and then rises from 3 to 8 micron. The radio galaxies are on average a factor 3-10 less luminous in this spectral range than the quasars. These characteristics are consistent with composite emission from a heavily reddened AGN plus starlight from the host galaxy. The mid-infrared colors and radio to mid-infrared spectral slopes of individual galaxies are also consistent with this picture. Individual galaxies show different amounts of extinction and host galaxy starlight, consistent with the orientation-dependent unified scheme.

Using three newly identified galaxy clusters at z~1 (photometric redshift) we measure the evolution of the galaxies within clusters from high redshift to the present day by studying the growth of the red cluster sequence. The clusters are located in the Spitzer Infrared Array Camera (IRAC) Dark Field, an extremely deep mid-infrared survey near the north ecliptic pole with photometry in 18 total bands from X-ray through far-IR. Two of the candidate clusters are additionally detected as extended emission in matching Chandra data in the survey area allowing us to measure their masses to be M_{500}= 6.2 \pm 1.0 \times 10^{13} and 3.6 \pm 1.1 \times 10^{13} solar masses. For all three clusters we create a composite color magnitude diagram in rest-frame B-K using our deep HST and Spitzer imaging. By comparing the fraction of low luminosity member galaxies on the composite red sequence with the corresponding population in local clusters at z=0.1 taken from the COSMOS survey, we examine the effect of a galaxy's mass on its evolution. We find a deficit of faint galaxies on the red sequence in our z~1 clusters which implies that more massive galaxies have evolved in clusters faster than less massive galaxies, and that the less massive galaxies are still forming stars in clusters such that they have not yet settled onto the red sequence.

We explore the stellar population of M31 in a Spitzer Space Telescope survey utilizing IRAC and MIPS observations. Red supergiants are the brightest objects seen in the infrared; they are a prominent evolutionary phase. Due to their circumstellar envelopes, many of these radiate the bulk of their luminosity at IRAC wavelengths and do not stand out in the near infrared or optically. Going fainter, we see large numbers of luminous asymptotic giant branch (AGB) stars, many of which are known long period variables. Relative to M33, the AGB carbon star population of M31 appears sparse, but this needs to be spectroscopically confirmed.

High-redshift submillimetre-bright galaxies identified by blank field surveys at millimetre and submillimetre wavelengths appear in the region of the IRAC colour-colour diagrams previously identified as the domain of luminous active galactic nuclei (AGNs). Our analysis using a set of empirical and theoretical dusty starburst spectral energy distribution (SED) models shows that power-law continuum sources associated with hot dust heated by young (<100 Myr old), extreme starbursts at z>2 also occupy the same general area as AGNs in the IRAC colour-colour plots. A detailed comparison of the IRAC colours and SEDs demonstrates that the two populations are distinct from each other, with submillimetre-bright galaxies having a systematically flatter IRAC spectrum (>1 mag bluer in the observed [4.5]-[8.0] colour). Only about 20% of the objects overlap in the colour-colour plots, and this low fraction suggests that submillimetre galaxies powered by a dust-obscured AGN are not common. The red IR colours of the submillimetre galaxies are distinct from those of the ubiquitous foreground IRAC sources, and we propose a set of IR colour selection criteria for identifying SMG counterparts that can be used even in the absence of radio or Spitzer MIPS 24 micron data.

The Extended Groth Strip (EGS) is one of the premier fields for extragalactic deep surveys. Deep observations of the EGS with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope cover an area of 0.38 square degrees to a 50% completeness limit of 1.5 uJy at 3.6 um. The catalog comprises 57434 objects detected at 3.6 um, with 84%, 28%, and 24% also detected at 4.5, 5.8, and 8.0 um. Number counts are consistent with results from other Spitzer surveys. Color distributions show that the EGS IRAC sources comprise a mixture of populations: low-redshift star-forming galaxies, quiescent galaxies dominated by stellar emission at a range of redshifts, and high redshift galaxies and AGN.

We combine the results of the Spitzer IRAC Shallow Survey and the Chandra XBootes Survey of the 8.5 square degrees Bootes field of the NOAO Deep Wide- Field Survey to produce the largest comparison of mid-IR and X-ray sources to date. The comparison is limited to sources with X-ray fluxes >8x10-15 erg cm-2s-1 in the 0.5-7.0 keV range and mid-IR sources with 3.6 um fluxes brighter than 18.4 mag (12.3 uJy). In this most sensitive IRAC band, 85% of the 3086 X-ray sources have mid-IR counterparts at an 80% confidence level based on a Bayesian matching technique. Only 2.5% of the sample have no IRAC counterpart at all based on visual inspection. Even for a smaller but a significantly deeper Chandra survey in the same field, the IRAC Shallow Survey recovers most of the X-ray sources. A majority (65%) of the Chandra sources detected in all four IRAC bands occupy a well-defined region of IRAC [3.6] - [4.5] vs [5.8] - [8.0] color-color space. These X-ray sources are likely infrared luminous, unobscured type I AGN with little mid-infrared flux contributed by the AGN host galaxy. Of the remaining Chandra sources, most are lower luminosity type I and type II AGN whose mid-IR emission is dominated by the host galaxy, while approximately 5% are either Galactic stars or very local galaxies.

The application of multi-wavelength selection techniques is crucial for discovering a complete and unbiased set of Active Galactic Nuclei (AGNs). Here, we select a sample of 72 AGN candidates in the Extended Groth Strip (EGS) using deep radio and mid-infrared data from the Very Large Array (VLA) and the Spitzer Space Telescope, and analyze their properties across other wavelengths. Only 30% of these sources are detected in deep 200 ks Chandra X-ray Observatory pointings. The X-ray detected sources demonstrate moderate obscuration with column densities of N_H > 10^22 cm^-2. A stacked image of sources undetected by Chandra shows low levels of X-ray activity, suggesting they may be faint or obscured AGNs. Less than 40% of our sample are selected as AGNs with optical broad lines, mid-infrared power laws, or X-ray detections. Thus, if our candidates are indeed AGNs, the radio/mid-infrared selection criteria we use provide a powerful tool for identifying sources missed by other surveys.

We present photometric analysis of deep mid-infrared observations obtained by Spitzer/IRAC covering the fields Q1422+2309, Q2233+1341, DSF2237a,b, HDFN, SSA22a,b and B20902+34, giving the number counts and the depths for each field. In a sample of 751 LBGs lying in those fields, 443, 448, 137 and 152 are identified at 3.6microns, 4.5microns, 5.8microns, 8.0microns IRAC bands respectively, expanding their spectral energy distribution to rest-near-infrared and revealing that LBGs display a variety of colours. Their rest-near-infrared properties are rather inhomogeneous, ranging from those that are bright in IRAC bands and exhibit [R]-[3.6] > 1.5 colours to those that are faint or not detected at all in IRAC bands with [R]-[3.6] < 1.5 colours and these two groups of LBGs are investigated. We compare the mid-IR colours of the LBGs with the colours of star-forming galaxies and we find that LBGs have colours consistent with star-foming galaxies at z~3. The properties of the LBGs detected in the 8microns IRAC band (rest frame K-band) are examined separately, showing that they exhibit redder [R]-[3.6] colours than the rest of the population and that IRAC 8microns band can be used as a diagnostic tool, to separate AGN dominated objects from normal star-forming galaxies at z~3

The ISO/LWS far-infrared spectrum of the ultraluminous galaxy Mkn 231 shows OH and H_2O lines in absorption from energy levels up to 300 K above the ground state, and emission in the [O I] 63 micron and [C II] 158 micron lines. Our analysis shows that OH and H_2O are radiatively pumped by the far-infrared continuum emission of the galaxy. The absorptions in the high-excitation lines require high far-infrared radiation densities, allowing us to constrain the properties of the underlying continuum source. The bulk of the far-infrared continuum arises from a warm (T_dust=70-100 K), optically thick (tau_100micron=1-2) medium of effective diameter 200-400 pc. In our best-fit model of total luminosity L_IR, the observed OH and H2O high-lying lines arise from a luminous (L/L_IR~0.56) region with radius ~100 pc. The high surface brightness of this component suggests that its infrared emission is dominated by the AGN. The derived column densities N(OH)>~10^{17} cm^{-2} and N(H_2O)>~6x10^{16} cm^{-2} may indicate XDR chemistry, although significant starburst chemistry cannot be ruled out. The lower-lying OH, [C II] 158 micron, and [O I] 63 micron lines arise from a more extended (~350 pc) starburst region. We show that the [C II] deficit in Mkn 231 is compatible with a high average abundance of C+ because of an extreme overall luminosity to gas mass ratio. Therefore, a [C II] deficit may indicate a significant contribution to the luminosity by an AGN, and/or by extremely efficient star formation.

We present results from Chandra X-ray and Spitzer mid-infrared observations of the interacting galaxy pair NGC6872/IC4970 in the Pavo galaxy group and show that the smaller companion galaxy IC4970 hosts a highly obscured active galactic nucleus (AGN). The 0.5-10 keV X-ray luminosity of the nucleus is variable, increasing by a factor 2.9 to 1.7 x 10^{42} erg/s (bright state) on ~100 ks timescales. The X-ray spectrum of the is heavily absorbed (N_H = 3 x 10^{23} cm^{-2}) for power law models with Gamma = 1.5-2.0 and shows a clear 6.4 keV Fe Kalpha line with equivalent width of 144-195 eV. Limits on the diffuse emission in IC4970 from Chandra X-ray data suggest that the available power from Bondi accretion of hot interstellar gas may be an order of magnitude too small to power the AGN. Spitzer images show that 8 micron nonstellar emission is concentrated in the central 1 kpc of IC4970, consistent with high obscuration in this region. The mid-infrared colors of the nucleus are consistent with those expected for a highly obscured AGN. Taken together these data suggest that the nucleus of IC4970 is a Seyfert 2, triggered and fueled by cold material supplied to the central supermassive black hole as a result of the off-axis collision of IC4970 with the cold-gas rich spiral galaxy NGC6872.

We have used the Spitzer Space Telescope to study the dust properties of a sample of star-forming dwarf galaxies. The differences in the mid-infrared spectral energy distributions for these galaxies which, in general, are low metallicity systems, indicate differences in the physical properties, heating, and/or distribution of the dust. Specifically, these galaxies have more hot dust and/or very small grains and less PAH emission than either spiral or higher luminosity starburst galaxies. As has been shown in previous studies, there is a gradual decrease in PAH emission as a function of metallicity. Because much of the energy from star formation in galaxies is re-radiated in the mid-infrared, star-formation rate indicators based on both line and continuum measurements in this wavelength range are coming into more common usage. We show that the variations in the interstellar medium properties of galaxies in our sample, as measured in the mid-infrared, result in over an order of magnitude spread in the computed star-formation rates.

We have used the Submillimeter Array to image a flux limited sample of seven submillimeter galaxies, selected by the AzTEC camera on the JCMT at 1.1 mm, in the COSMOS field at 890um with 2" resolution. All of the sources - two radio-bright and five radio-dim - are detected as single point-sources at high significance (> 6\sigma), with positions accurate to 0.2" that enable counterpart identification at other wavelengths observed with similarly high angular resolution. All seven have IRAC counterparts, but only two have secure counterparts in deep HST/ACS imaging. As compared to the two radio-bright sources in the sample, and those in previous studies, the five radio-dim sources in the sample (1) have systematically higher submillimeter-to-radio flux ratios, (2) have lower IRAC 3.6-8.0um fluxes, and (3) are not detected at 24um. These properties, combined with size constraints at 890um (\theta < 1.2"), suggest that the radio-dim submillimeter galaxies represent a population of very dusty starbursts, with physical scales similar to local ultraluminous infrared galaxies, and an average redshift higher than radio-bright sources.

A SST survey in the NOAO Deep-Wide Field in Bo\"otes provides a complete, 8-micron-selected sample of galaxies to a limiting (Vega) magnitude of 13.5. In the 6.88 deg$^2$ field sampled, 79% of the 4867 galaxies have spectroscopic redshifts, allowing an accurate determination of the local (z<0.3) galaxy luminosity function. Stellar and dust emission can be separated on the basis of observed galaxy colors. Dust emission (mostly PAH) accounts for 80% of the 8 micron luminosity, stellar photospheres account for 19%, and AGN emission accounts for roughly 1 %. A sub-sample of the 8 micron-selected galaxies have blue, early-type colors, but even most of these have significant PAH emission. The luminosity functions for the total 8 micron luminosity and for the dust emission alone are both well fit by Schechter functions. For the 8 micron luminosity function, the characteristic luminosity is \nu L_{\nu}^*(8.0 \micron) = 1.8 \times 10^{10}$ \Lsun while for the dust emission alone it is 1.6 x 10^{10}$ \Lsun\null. The average 8 \micron luminosity density at z<0.3 is 3.1 x 10^7 \Lsun Mpc^{-3}, and the average luminosity density from dust alone is 2.5 x 10^7 \Lsun Mpc^{-3}. This luminos ity arises predominantly from galaxies with 8 \micron luminosities ($\nu L_{\nu}$) between $2\times 10^9$ and $2 x 10^{10}$ \Lsun, i.e., normal galaxies, not LIRGs or ULIRGs.

We analyze five epochs of Spitzer Space Telescope/Infrared Array Camera (IRAC) observations of the nearby spiral galaxy M33. Each epoch covered nearly a square degree at 3.6, 4.5, and 8.0 microns. The point source catalog from the full dataset contains 37,650 stars. The stars have luminosities characteristic of the asymptotic giant branch and can be separated into oxygen-rich and carbon-rich populations by their [3.6] - [4.5] colors. The [3.6] - [8.0] colors indicate that over 80% of the stars detected at 8.0 microns have dust shells. Photometric comparison of epochs using conservative criteria yields a catalog of 2,923 variable stars. These variables are most likely long-period variables amidst an evolved stellar population. At least one-third of the identified carbon stars are variable.

We have imaged a $\sim$6 arcminute$^2$ region in the Bo\"otes Deep Field using the 350 $\mu$m-optimised second generation Submillimeter High Angular Resolution Camera (SHARC II), achieving a peak 1$\sigma$ sensitivity of $\sim$5 mJy. We detect three sources above 3$\sigma$, and determine a spurious source detection rate of 1.09 in our maps. In the absence of $5\sigma$ detections, we rely on deep 24 $\mu$m and 20 cm imaging to deduce which sources are most likely to be genuine, giving two real sources. From this we derive an integral source count of 0.84$^{+1.39}_{-0.61}$ sources arcmin$^{-2}$ at $S>13$ mJy, which is consistent with 350 $\mu$m source count models that have an IR-luminous galaxy population evolving with redshift. We use these constraints to consider the future for ground-based short-submillimetre surveys.

Mid-infrared observations of the Andromeda galaxy, M31, obtained with the Infrared Array Camera on board the Spitzer Space Telescope, are presented. The image mosaics cover areas of approximate 3.7deg x 1.6deg and include the satellite galaxies M32 and NGC 205. The appearance of M31 varies dramatically in the different mid-infrared bands, from the smooth bulge and disk of the old stellar population seen at 3.6um to the well-known '10 kpc ring' dominating the 8um image. The similarity of the 3.6um and optical isophotes and nearly constant optical-mid-infrared color over the inner 400 arcsec confirms that there is no significant extinction at optical wavelengths in M31's bulge. The nuclear colors indicate the presence of dust but not an infrared-bright nucleus. The integrated 8um non-stellar luminosity implies a star formation rate of 0.4 Msun/yr, consistent with other indicators that show M31 to be a quiescent galaxy.

We have measured the contribution of submillimeter and mid-infrared sources to the extragalactic background radiation at 70 and 160um. Specifically, we have stacked flux in 70 and 160um Spitzer Space Telescope (Spitzer) observations of the Canada-UK Deep Sub-millimeter Survey 14h field at the positions of 850um sources detected by SCUBA and also 8 and 24um sources detected by Spitzer. We find that per source, the SCUBA galaxies are the strongest and the 8um sources the weakest contributors to the background flux at both 70 and 160um. Our estimate of the contribution of the SCUBA sources is higher than previous estimates. However, expressed as a total contribution, the full 8um source catalogue accounts for twice the total 24um source contribution and ~10 times the total SCUBA source contribution. The 8um sources account for the majority of the background radiation at 160um with a flux of 0.87+/-0.16 MJy/sr and at least a third at 70um with a flux of 0.103+/-0.019 MJy/sr. These measurements are consistent with current lower limits on the background at 70 and 160um. Finally, we have investigated the 70 and 160um emission from the 8 and 24um sources as a function of redshift. We find that the average 70um flux per 24um source and the average 160um flux per 8 and 24um source is constant over all redshifts, up to z~4. In contrast, the low-redshift half (z<1) of the of 8um sample contributes approximately four times the total 70um flux of the high-redshift half. These trends can be explained by a single non-evolving SED.

We present 0.5 -160 micron Spectral Energy Distributions (SEDs) of galaxies, detected at 70microns with the Multiband Imaging Photometer for Spitzer (MIPS), using broadband imaging data from Spitzer and ground-based telescopes. Spectroscopic redshifts, in the range 0.2<z<1.5, have been measured as part of the Deep Extragalactic Evolutionary Probe2 (DEEP2) project. Based on the SEDs we explore the nature and physical properties of the sources. Using the optical spectra we derive Hbeta and [OII]-based Star Formation Rates (SFR) which are 10-100 times lower than SFR estimates based on IR and radio. The median offset in SFR between optical and IR is reduced by a factor of ~3 when we apply a typical extinction corrections. We investigate mid-to-far infrared correlations for low redshift (>0.5) and high redshift (0.5<z<1.2) bins. Using this unique ``far-infrared'' selected sample we derive an empirical mid to far-infrared relationship that can be used to estimate the infrared energy budget of galaxies in the high-redshift universe. Our sample can be used as a template to translate far-infrared luminosities into bolometric luminosities for high redshift objects.

The unusual morphology of the Andromeda Spiral (Messier 31, the closest spiral galaxy to the Milky Way) has long been an enigma. Although regarded for decades as showing little evidence of a violent history, M~31 has a well-known outer ring of star formation at a radius of 10 kpc whose center is offset from the galaxy nucleus. In addition, the outer galaxy disk is warped as seen at both optical and radio wavelengths. The halo contains numerous loops and ripples. Here we report the discovery, based on analysis of previously-obtained data, of a second, inner dust ring with projected dimensions 1.5 by 1 kpc and offset by ~0.5kpc from the center of the galaxy. The two rings appear to be density waves propagating in the disk. Numerical simulations offer a completely new interpretation for the morphology of M31: both rings result from a companion galaxy plunging head-on through the center of the disk of M31. The most likely interloper is M32. Head-on collisions between galaxies are rare, but it appears nonetheless that one took place 210 million years ago in our Local Group of galaxies.

As part of the "All Wavelength Extended Groth Strip International Survey" (AEGIS), we describe the panchromatic characterization of an X-ray luminous active galactic nucleus (AGN) in a merging galaxy at z=1.15. This object is detected at infrared (8mic, 24mic, 70mic, 160mic), submillimeter (850mic) and radio wavelengths, from which we derive a bolometric luminosity L_bol ~ 9x10^12 Lsol. We find that the AGN clearly dominates the hot dust emission below 40mic but its total energetic power inferred from the hard X-rays is substantially less than the bolometric output of the system. About 50% of the infrared luminosity is indeed produced by a cold dust component that probably originates from enshrouded star formation in the host galaxy. In the context of a coeval growth of stellar bulges and massive black holes, this source might represent a ``transition'' object sharing properties with both quasars and luminous starbursts. Study of such composite galaxies will help address how the star formation and disk-accretion phenomena may have regulated each other at high redshift and how this coordination may have participated to the build-up of the relationship observed locally between the masses of black holes and stellar spheroids.

We study 87 Extremely Red Objects (EROs), selected both to have color redder than R-[3.6] = 4.0, and to have confirmed spectroscopic redshifts. Together, these two constraints result in this sample populating a fairly narrow redshift range at 0.76 < z < 1.42. The key new ingredient included here is deep Spitzer Space Telescope InfraRed Array Camera (IRAC) data. Based on [3.6]-[8.0] color, we demonstrate that it is possible to classify EROs into early-type, dusty starburst, or power-law (AGN) types. We present ultraviolet to mid-infrared spectral energy distributions (SEDs) and Advanced Camera for Surveys (ACS) images, both of which support our simple IRAC color classification.

We report the detection of rest--frame 6.2 and 7.7 \micron emission features arising from Polycyclic Aromatic Hydrocarbons (PAH) in the Spitzer/IRS spectrum of an infrared-luminous Lyman break galaxy at z=3.01. This is currently the highest redshift galaxy where these PAH emission features have been detected. The total infrared luminosity inferred from the MIPS 24 \micron and radio flux density is 2$\times10^{13}$ L$_{\odot}$, which qualifies this object as a so--called hyperluminous infrared galaxy (HyLIRG). However, unlike local HyLIRGs which are generally associated with QSO/AGNs and have weak or absent PAH emission features, this HyLIRG has very strong 6.2 and 7.7 \micron PAH emission. We argue that intense star formation dominates the infrared emission of this source, although we cannot rule out the presence of a deeply obscured AGN. This LBG appears to be a distorted system in the HST ACS F606W and F814W images, possibly indicating that a significant merger or interaction is driving the large IR luminosity.

A new 6-cm survey of almost 0.6 square degrees to a limit of 0.55-mJy/beam (10-sigma) finds 37 isolated radio sources and 7 radio source pairs (not necessarily physical companions). IRAC counterparts are identified for at least 92% of the radio sources within the area of deep IRAC coverage, which includes 31 isolated sources and 6 pairs. This contrasts with an identification rate of <74% to R<23.95 in visible light. Eight of the IRAC galaxies have power law spectral energy distributions, implying that the mid-infrared emission comes from a powerful AGN. The remaining 26 IRAC galaxies show stellar emission in the mid-infrared, probably in most of these galaxies because the stellar emission is bright enough to outshine an underlying AGN. The infrared colors suggest that the majority of these galaxies are bulge-dominated and have redshifts between approximately 0.5 and 1.0. Visible spectra from the DEEP2 redshift survey, available for 11 galaxies, are consistent with this suggestion. The IRAC galaxies fall into two distinct groups in a color-magnitude diagram, one group (the "stripe") includes all the AGN. The other group (the "blue clump") has blue 3.6 to 8 micron colors and a small range of 8 micron magnitudes. This separation should be useful in classifying galaxies found in other radio surveys.

Accurate photometric redshifts are calculated for nearly 200,000 galaxies to a 4.5 micron flux limit of ~13 uJy in the 8.5 deg^2 Spitzer/IRAC Shallow survey. Using a hybrid photometric redshift algorithm incorporating both neural-net and template-fitting techniques, calibrated with over 15,000 spectroscopic redshifts, a redshift accuracy of \sigma = 0.06(1+z) is achieved for 95% of galaxies at 0<z<1.5. The accuracy is \sigma = 0.12(1+z) for 95% of AGN at 0<z<3. Redshift probability functions, central to several ongoing studies of the galaxy population, are computed for the full sample. We demonstrate that these functions accurately represent the true redshift probability density, allowing the calculation of valid confidence intervals for all objects. These probability functions have already been used to successfully identify a population of Spitzer-selected high redshift (z>1) galaxy clusters. We present one such spectroscopically confirmed cluster at <z>=1.24, ISCS J1434.5+3427. Finally, we present a measurement of the 4.5 micron-selected galaxy redshift distribution.