We present Spitzer/IRAC observations of nine $z'$-band dropouts highly magnified (2<mu<12) by the Bullet Cluster. We combine archival imaging with our Exploratory program (SURFS UP), which results in a total integration time of ~30 hr per IRAC band. We detect (>3sigma) in both IRAC bands the brightest of these high-redshift galaxies, with [3.6]=23.80+-0.28 mag, [4.5]=23.78+-0.25 mag, and (H-[3.6])=1.17+-0.32 mag. The remaining eight galaxies are undetected to [3.6]~26.4 mag and [4.5]~26.0 mag with stellar masses of ~5x10^7 M_sol. The detected galaxy has an estimated magnification of mu=12+-4, which implies this galaxy has an ultraviolet luminosity of L_1500~0.3 L*_{z=7} --- the lowest luminosity individual source detected in IRAC at z>7. By modeling the broadband photometry, we estimate the galaxy has an intrinsic star-formation rate of SFR~1.3 M_sol/yr and stellar mass of M~2x10^9 M_sol, which gives a specific star-formation rate of sSFR~0.7 Gyr^-1. If this galaxy had sustained this star-formation rate since z~20, it could have formed the observed stellar mass (to within a factor of ~2), we also discuss alternate star-formation histories and argue the exponentially-increasing model is unlikely. Finally, based on the intrinsic star-formation rate, we estimate this galaxy has a likely [C II] flux of <f_[C II]> = 10^{-17} erg/s/cm2.

SURFSUP is a joint Spitzer and HST Exploration Science program using 10 galaxy clusters as cosmic telescopes to study z >~ 7 galaxies at intrinsically lower luminosities, enabled by gravitational lensing, than blank field surveys of the same exposure time. Our main goal is to measure stellar masses and ages of these galaxies, which are the most likely sources of the ionizing photons that drive reionization. Accurate knowledge of the star formation density and star formation history at this epoch is necessary to determine whether these galaxies indeed reionized the universe. Determination of the stellar masses and ages requires measuring rest frame optical light, which only Spitzer can probe for sources at z >~ 7, for a large enough sample of typical galaxies. Our program consists of 550 hours of Spitzer/IRAC imaging covering 10 galaxy clusters with very well-known mass distributions, making them extremely precise cosmic telescopes. We combine our data with archival observations to obtain mosaics with ~30 hours exposure time in both 3.6$\mu$m and 4.5$\mu$m in the central 4 arcmin x 4 arcmin field and ~15 hours in the flanking fields. This results in 3-$\sigma$ sensitivity limits of ~26.6 and ~26.2AB magnitudes for the central field in the IRAC 3.6 and 4.5$\mu$m bands, respectively. To illustrate the survey strategy and characteristics we introduce the sample, present the details of the data reduction and demonstrate that these data are sufficient for in-depth studies of z >~ 7 sources (using a z=9.5 galaxy behind MACSJ1149.5+2223 as an example). For the first cluster of the survey (the Bullet Cluster) we have released all high-level data mosaics and IRAC empirical PSF models. In the future we plan to release these data products for the entire survey.

We present the first results of our spectroscopic follow-up of 6.5 < z < 10 candidate galaxies behind clusters of galaxies. We report the spectroscopic confirmation of an intrinsically faint Lyman break galaxy (LBG) identified as a z 850LP-band dropout behind the Bullet Cluster. We detect an emission line at {\lambda} = 9412 {\AA} at >5{\sigma} significance using a 16 hr long exposure with FORS2 VLT. Based on the absence of flux in bluer broadband filters, the blue color of the source, and the absence of additional lines, we identify the line as Ly{\alpha} at z = 6.740 \pm 0.003. The integrated line flux is f = (0.7 \pm 0.1 \pm 0.3) \times 10^{-17} erg^{-1} s^{-1} cm^{-2} (the uncertainties are due to random and flux calibration errors, respectively) making it the faintest Ly{\alpha} flux detected at these redshifts. Given the magnification of {\mu} = 3.0 \pm 0.2 the intrinsic (corrected for lensing) flux is f^int = (0.23 \pm 0.03 \pm 0.10 \pm 0.02) \times 10^{-17} erg^{-1} s^{-1} cm^{-2} (additional uncertainty due to magnification), which is ~2-3 times fainter than other such measurements in z ~ 7 galaxies. The intrinsic H 160W-band magnitude of the object is m^int(H_160W)=27.57 \pm 0.17, corresponding to 0.5 L* for LBGs at these redshifts. The galaxy is one of the two sub-L* LBG galaxies spectroscopically confirmed at these high redshifts (the other is also a lensed z = 7.045 galaxy), making it a valuable probe for the neutral hydrogen fraction in the early universe.

We use imaging obtained with the Hubble Space Telescope Wide Field Camera 3 to search for z_850 dropouts at z~7 and J_110 dropouts at z~9 lensed by the Bullet Cluster. In total we find 10 z_850 dropouts in our 8.27 arcmin^2 field. Using magnification maps from a combined weak and strong lensing mass reconstruction of the Bullet Cluster and correcting for estimated completeness levels, we calculate the surface density and luminosity function of our z_850 dropouts as a function of intrinsic (accounting for magnification) magnitude. We find results consistent with published blank field surveys, despite using much shallower data, and demonstrate the effectiveness of cluster surveys in the search for z~7 galaxies.

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.

We present a detection-significance-limited catalog of 21 Sunyaev-Zel'dovich selected galaxy clusters. These clusters, along with 1 unconfirmed candidate, were identified in 178 deg^2 of sky surveyed in 2008 by the South Pole Telescope to a depth of 18 uK-arcmin at 150 GHz. Optical imaging from the Blanco Cosmology Survey (BCS) and Magellan telescopes provided photometric (and in some cases spectroscopic) redshift estimates, with catalog redshifts ranging from z=0.15 to z>1, with a median z = 0.74. Of the 21 confirmed galaxy clusters, three were previously identified as Abell clusters, three were presented as SPT discoveries in Staniszewski et al, 2009, and three were first identified in a recent analysis of BCS data by Menanteau et al, 2010; the remaining 12 clusters are presented for the first time in this work. Simulated observations of the SPT fields predict the sample to be nearly 100% complete above a mass threshold of M_200 ~ 5x10^14 M_sun/h at z = 0.6. This completeness threshold pushes to lower mass with increasing redshift, dropping to ~4x10^14 M_sun/h at z=1. The size and redshift distribution of this catalog are in good agreement with expectations based on our current understanding of galaxy clusters and cosmology. In combination with other cosmological probes, we use the cluster catalog to improve estimates of cosmological parameters. Assuming a standard spatially flat wCDM cosmological model, the addition of our catalog to the WMAP 7-year analysis yields sigma_8 = 0.81 +- 0.09 and w = -1.07 +- 0.29, a ~50% improvement in precision on both parameters over WMAP7 alone.

We use data from the first 100 square-degree field observed by the South Pole Telescope (SPT) in 2008 to measure the angular power spectrum of temperature anisotropies contributed by the background of dusty star-forming galaxies (DSFGs) at millimeter wavelengths. From the auto and cross-correlation of 150 and 220 GHz SPT maps, we significantly detect both Poisson distributed and, for the first time at millimeter wavelengths, clustered components of power from a background of DSFGs. The spectral indices between 150 and 220 GHz of the Poisson and clustered components are found to be 3.86 +- 0.23 and 3.8 +- 1.3 respectively, implying a steep scaling of the dust emissivity index beta ~ 2. The Poisson and clustered power detected in SPT, BLAST (at 600, 860, and 1200 GHz), and Spitzer (1900 GHz) data can be understood in the context of a simple model in which all galaxies have the same graybody spectrum with dust emissivity index of beta = 2 and dust temperature T_d = 34 K. In this model, half of the 150 GHz background light comes from redshifts greater than 3.2. We also use the SPT data to place an upper limit on the amplitude of the kinetic Sunyaev-Zel'dovich power spectrum at l = 3000 of 13 uK^2 at 95% confidence.

We present Sunyaev-Zel'dovich measurements of 15 massive X-ray selected galaxy clusters obtained with the South Pole Telescope. The Sunyaev-Zel'dovich (SZ) cluster signals are measured at 150 GHz, and concurrent 220 GHz data are used to reduce astrophysical contamination. Radial profiles are computed using a technique that takes into account the effects of the beams and filtering. In several clusters, significant SZ decrements are detected out to a substantial fraction of the virial radius. The profiles are fit to the beta model and to a generalized NFW pressure profile, and are scaled and stacked to probe their average behavior. We find model parameters that are consistent with previous studies: beta=0.86 and r_core/r_500 = 0.20 for the beta model, and (alpha, beta, gamma, c_500)=(1.0,5.5,0.5,1.0) for the generalized NFW model. Both models fit the SPT data comparably well, and both are consistent with the average SZ profile out to the virial radius. The integrated Compton-y parameter Y_SZ is computed for each cluster using both model-dependent and model-independent techniques, and the results are compared to X-ray estimates of cluster parameters. We find that Y_SZ scales with Y_X and gas mass with low scatter. Since these observables have been found to scale with total mass, our results point to a tight mass-observable relation for the SPT cluster survey.

We report cosmic microwave background (CMB) power spectrum measurements from the first 100 sq. deg. field observed by the South Pole Telescope (SPT) at 150 and 220 GHz. On angular scales where the primary CMB anisotropy is dominant, ell ~< 3000, the SPT power spectrum is consistent with the standard LambdaCDM cosmology. On smaller scales, we see strong evidence for a point source contribution, consistent with a population of dusty, star-forming galaxies. After we mask bright point sources, anisotropy power on angular scales of 3000 < ell < 9500 is detected with a signal-to-noise > 50 at both frequencies. We combine the 150 and 220 GHz data to remove the majority of the point source power, and use the point source subtracted spectrum to detect Sunyaev-Zel'dovich (SZ) power at 2.6 sigma. At ell=3000, the SZ power in the subtracted bandpowers is 4.2 +/- 1.5 uK^2, which is significantly lower than the power predicted by a fiducial model using WMAP5 cosmological parameters. This discrepancy may suggest that contemporary galaxy cluster models overestimate the thermal pressure of intracluster gas. Alternatively, this result can be interpreted as evidence for lower values of sigma8. When combined with an estimate of the kinetic SZ contribution, the measured SZ amplitude shifts sigma8 from the primary CMB anisotropy derived constraint of 0.794 +/- 0.028 down to 0.773 +/- 0.025. The uncertainty in the constraint on sigma8 from this analysis is dominated by uncertainties in the theoretical modeling required to predict the amplitude of the SZ power spectrum for a given set of cosmological parameters.

The use of Type Ia supernovae (SNe Ia) as cosmological standard candles is a key to solving the mystery of dark energy. Improving the calibration of SNe Ia increases their power as cosmological standard candles. We find tentative evidence for a correlation between the late-time lightcurve slope and the peak luminosity of SNe Ia in the B band; brighter SNe Ia seem to have shallower lightcurve slopes between 100 and 150 days from maximum light. Using a Markov Chain Monte Carlo analysis in calibrating SNe Ia, we are able to simultaneously take into consideration the effect of dust extinction, the luminosity and lightcurve width correlation (parametrized by \Delta m_{15}), and the luminosity and late-time lightcurve slope correlation. For the available sample of 11 SNe Ia with well-measured late-time lightcurves, we find that correcting for the correlation between luminosity and late-time lightcurve slope of the SNe Ia leads to an intrinsic dispersion of 0.12 mag in the Hubble diagram. Our results have significant implications for future supernova surveys aimed to illuminate the nature of dark energy.

A comparative study of spectra of 21 Type Ia supernovae (SNe Ia) obtained about one week before maximum light, and 8 spectra obtained 11 or more days before maximum, is presented. To a large extent the premaximum spectra exhibit the defining characteristics of the four groups defined in Paper II (core-normal, broad-line, cool, and shallow-silicon). Comparisons with SYNOW synthetic spectra show that all strong features and most weak ones can be accounted for in a plausible way. The issues of detached high-velocity features, the possible ubiquity of carbon clumps, the maximum detectable ejecta velocities, and the possibility of blueshifted emission-line peaks are discussed.

Turbulence in an incompressible fluid with and without a magnetic field as well as moderately compressible MHD turbulence are compared. For the magnetohydrodynamic (MHD) models the probability distribution functions of the velocity components perpendicular to the external magnetic field are like the incompressible hydrodynamic (HD) model while those parallel to the field have a smaller range of velocities. The probability distribution functions of the transverse velocity increments for the MHD models decline slower than the incompressible HD model. The similarity of incompressible HD and both incompressible and compressible MHD turbulence persists over high order longitudinal structure function scaling exponents measured in the global reference frame as well as for motions perpendicular to the local mean field. In these two frames the longitudinal scaling exponents of both MHD models seem to follow theoretical incompressible HD dissipation structure predictions while the transverse scaling exponents of the incompressible MHD model seem to follow the predictions for incompressible MHD. In the local magnetic system the motions parallel to the local mean field for both MHD models are different from incompressible HD motions. The fields of the MHD simulations are decomposed into Alfvenic, fast, and slow modes and stuided separately. (abridged)

A comparative study of near-maximum-light optical spectra of 24 Type Ia supernovae (SNe Ia) is presented. The spectra are quantified in two ways, and assigned to four groups. Seven "core-normal" SNe Ia have very similar spectra, except for strong high-velocity CaII absorption in SN 2001el. Seven SNe Ia are assigned to a "broad-line" group, the most extreme of which is SN 1984A. Five SNe Ia, including SN 1991bg, are assigned to a "cool" group. Five SNe Ia, including SN 1991T, are assigned to a "shallow-silicon" group. Comparisons with SYNOW synthetic spectra provide a basis for discussion of line identifications, and an internally consistent quantification of the maximum-light spectroscopic diversity among SNe Ia. The extent to which SN Ia maximum-light spectra appear to have a continuous distribution of properties, rather than consisting of discrete subtypes, is discussed.

As the first step in a comprehensive, comparative, direct analysis of the spectra of Type Ia supernovae (SNe Ia), we use the parameterized supernova synthetic-spectrum code, SYNOW, to interpret 26 spectra of the well-observed SN 1994D. Our results are consistent with the traditional view that the composition structure (element abundance fractions versus ejection velocity) is radially stratified. We find that resonance-scattering features due to permitted lines of Ca II, Na I, and Fe II persist to more than 100 days after explosion. The fitting parameters for SN 1994D, together with those to be determined for other SNe Ia, will provide an internally consistent quantification of the spectroscopic diversity among SNe Ia, and shed light on how the various manifestations of observational diversity are related to their physical causes.