Far-infrared (FIR) dust polarimetry enables the study of interstellar magnetic fields via tracing of the polarized emission from dust grains that are partially aligned with the direction of the field. The advent of high quality polarimetric data has permitted the use of statistical methods to extract both the direction and magnitude of the magnetic field. This work presents an effort to produce maps (integrated along the line of sight) of the magnetic field strength in both the line-of-sight (LOS) direction and in the plane-of-sky (POS) in the Orion Molecular Cloud (OMC-1) using all four HAWC+ polarimetry bands (53, 89, 154 and 214 $\micron$). Using the statistics of the polarization data (angles and fractions) along with previously-reported Zeeman measurements, information about the three-dimensional magnetic field configuration was inferred over the OMC-1 region. POS magnetic field strengths of up to 2 mG were determined in OMC-1. The strongest fields were observed near the BN/KL object, while the OMC-1 bar shows strengths of up to few hundreds of $\mu$G. The LOS magnetic field has been determined to have comparable strengths to the POS components but it is differently distributed across the cloud. These estimates of the magnetic field components were used to produce more reliable maps of the mass-to-magnetic flux ratio ($M/\Phi$) -- a proxy for probing the conditions for star formation in molecular clouds -- and determine regions of sub- and super-criticality in OMC-1. Such maps can provide invaluable input and comparison to MHD simulations of star formation processes in filamentary structures of molecular clouds.

We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~ 1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multi-scale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~ 1 pc and remains continuous from the scales of filaments (~ 0.1 pc) to that of protostellar envelopes (~ 0.005 pc or ~ 1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data towards the filament network.

We report the first detection of galactic spiral structure by means of thermal emission from magnetically aligned dust grains. Our 89 $\mu$m polarimetric imaging of NGC 1068 with the High-resolution Airborne Wideband Camera/Polarimeter (HAWC+) on NASA's Stratospheric Observatory for Infrared Astronomy (SOFIA) also sheds light on magnetic field structure in the vicinity of the galaxy's inner-bar and active galactic nucleus (AGN). We find correlations between the 89 $\mu$m magnetic field vectors and other tracers of spiral arms, and a symmetric polarization pattern as a function of the azimuthal angle arising from the projection and inclination of the disk field component in the plane of the sky. The observations can be fit with a logarithmic spiral model with pitch angle of $16.9^{+2.7}_{-2.8}$$^{\circ}$ and a disk inclination of $48\pm2^{\circ}$. We infer that the bulk of the interstellar medium from which the polarized dust emission originates is threaded by a magnetic field that closely follows the spiral arms. Inside the central starburst disk ($<1.6$ kpc), the degree of polarization is found to be lower than for far-infrared sources in the Milky Way, and has minima at the locations of most intense star formation near the outer ends of the inner-bar. Inside the starburst ring, the field direction deviates from the model, becoming more radial along the leading edges of the inner-bar. The polarized flux and dust temperature peak $\sim 3-6$" NE of the AGN at the location of a bow shock between the AGN outflow and the surrounding interstellar medium, but the AGN itself is weakly polarized ($< 1$%) at both 53 and 89 \um.

We report on polarimetric maps made with HAWC+/SOFIA toward Rho Oph A, the densest portion of the Rho Ophiuchi molecular complex. We employed HAWC+ bands C (89 $\mu$m) and D (154 $\mu$m). The slope of the polarization spectrum was investigated by defining the quantity R_DC = p_D/p_C, where p_C and p_D represent polarization degrees in bands C and D, respectively. We find a clear correlation between R_DC and the molecular hydrogen column density across the cloud. A positive slope (R_DC > 1) dominates the lower density and well illuminated portions of the cloud, that are heated by the high mass star Oph S1, whereas a transition to a negative slope (R_DC < 1) is observed toward the denser and less evenly illuminated cloud core. We interpret the trends as due to a combination of: (1) Warm grains at the cloud outskirts, which are efficiently aligned by the abundant exposure to radiation from Oph S1, as proposed in the radiative torques theory; and (2) Cold grains deep in the cloud core, which are poorly aligned due to shielding from external radiation. To assess this interpretation, we developed a very simple toy model using a spherically symmetric cloud core based on Herschel data, and verified that the predicted variation of R_DC is consistent with the observations. This result introduces a new method that can be used to probe the grain alignment efficiency in molecular clouds, based on the analysis of trends in the far-infrared polarization spectrum.

We present the POL-2 850 $\mu$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We then use the polarization data at 850 $\mu$m to infer the plane-of-sky orientation of the large-scale magnetic field in Barnard 1. This magnetic field runs North-South across most of the cloud, with the exception of B1-c where it turns more East-West. From the dispersion of polarization angles, we calculate a turbulence correlation length of $5.0 \pm 2.5$ arcsec ($1500$ au), and a turbulent-to-total magnetic energy ratio of $0.5 \pm 0.3$ inside the cloud. We combine this turbulent-to-total magnetic energy ratio with observations of NH$_3$ molecular lines from the Green Bank Ammonia Survey (GAS) to estimate the strength of the plane-of-sky component of the magnetic field through the Davis-Chandrasekhar-Fermi method. With a plane-of-sky amplitude of $120 \pm 60$ $\mu$G and a criticality criterion $\lambda_c = 3.0 \pm 1.5$, we find that Barnard 1 is a supercritical molecular cloud with a magnetic field nearly dominated by its turbulent component.

We present the 850 $\mu$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scale ($\lesssim 1.0$ pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether or not the observed polarization traces the magnetic field in the HFS, we analyze the dependence between the observed polarization fraction and total intensity using a Bayesian approach with the polarization fraction described by the Rice likelihood function, which can correctly describe the probability density function (PDF) of the observed polarization fraction for low signal-to-noise ratio (SNR) data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in $A_V\sim$ 20--300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 regions. Our polarization maps reveal a curved magnetic field, possibly dragged by the contraction along the parsec-scale filament. We further obtain a magnetic field strength of 0.5$\pm$0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of $\sim$ $1.3\pm0.4$ and an Alfv\'{e}nic Mach number of $<$0.6. These results suggest that gravity and magnetic field is currently of comparable importance in the HFS, and turbulence is less important.

We present far-IR photometry and infrared spectrum of the z=3.9114 quasar/starburst composite system APM 08279+5255 obtained using the Stratospheric Observatory for Infrared Astronomy (SOFIA)/HAWC+ and the Spitzer Space Telescope Infrared Spectrograph (IRS). We decompose the IR-to-radio spectral energy distribution (SED), sampled in 51 bands, using (i) a model comprised of two-temperature modified blackbodies (MBB) and radio power-laws and (ii) a semi-analytic model, which also accounts for emission from a clumpy torus. The latter is more realistic but requires a well-sampled SED, as possible here. In the former model, we find temperatures of T_warm = 296^17_15 K and T_cold = 110^3_3 K for the warm and cold dust components, respectively. This model suggests that the cold dust component dominates the FIR energy budget (66%) but contributes only 17% to the total IR luminosity. Based on the torus models, we infer an inclination angle of i=15^8_8 degree and the presence of silicate emission, in accordance with the Type-1 active galactic nucleus nature of APM 08279+5255. Accounting for the torus' contribution to the FIR luminosity, we find a lensing-corrected star formation rate of SFR=3075x(4/mu_L) Msun yr^-1. We find that the central quasar contributes 30% to the FIR luminosity but dominates the total IR luminosity (93%). The 30% correction is in contrast to the 90% reported in previous work. In addition, the IR luminosity inferred from the torus model is a factor of two higher. These differences highlight the importance of adopting physically motivated models to properly account for IR emission in high-z quasars, which is now possible with SOFIA/HAWC+.

Understanding the physics of how stars form is a highly-prioritized goal of modern Astrophysics, in part because star formation is linked to both galactic dynamics on large scales and to the formation of planets on small scales. It is well-known that stars form from the gravitational collapse of molecular clouds, which are in turn formed out of the turbulent interstellar medium. Star formation is highly inefficient, with one of the likely culprits being the regulation against gravitational collapse provided by magnetic fields. Measurement of the polarized emission from interstellar dust grains, which are partially aligned with the magnetic field, provides a key tool for understanding the role these fields play in the star formation process. Over the past decade, much progress has been made by the most recent generation of polarimeters operating over a range of wavelengths (from the far-infrared through the millimeter part of the spectrum) and over a range of angular resolutions (from less than an arcsecond through fractions of a degree). Future developments in instrument sensitivity for ground-based, airborne, and space-borne polarimeters operating over range of spatial scales are critical for enabling revolutionary steps forward in our understanding of the magnetized turbulence from which stars are formed.

We report 850~$\mu$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $\rho$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $\mu$G, 136 $\pm$ 69 $\mu$G, and 213 $\pm$ 115 $\mu$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties.

We report new polarimetric and photometric maps of the massive star-forming region OMC-1 using the HAWC+ instrument on the Stratospheric Observatory for Infrared Astronomy (SOFIA). We present continuum polarimetric and photometric measurements of this region at 53, 89, 154, and 214 microns at angular resolutions of 5.1, 7.9, 14.0, and 18.7 arcseconds for the four bands, respectively. The photometric maps enable the computation of improved SEDs for the region. We find that at the longer wavelengths, the inferred magnetic field configuration matches the `hourglass' configuration seen in previous studies, indicating magnetically-regulated star formation. The field morphology differs at the shorter wavelengths. The magnetic field inferred at these wavelengths traces the bipolar structure of the explosive Becklin-Neugebauer (BN)/Kleinman-Low (KL) outflow emerging from OMC-1 behind the Orion Nebula. Using statistical methods to estimate the field strength in the region, we find that the explosion dominates the magnetic field near the center of the feature. Farther out, the magnetic field is close to energetic equilibrium with the ejecta and may be providing confinement to the explosion. The correlation between polarization fraction and the local polarization angle dispersion indicates that the depolarization as a function of unpolarized intensity is a result of intrinsic field geometry as opposed to decreases in grain alignment efficiency in denser regions.

We present Far-Infrared polarimetry observations of M82 at 53 and $154~\mu \rm{m}$ and NGC 253 at $89~\mu \rm{m}$, which were taken with HAWC+ in polarimetry mode on the Stratospheric Observatory for Infrared Astronomy (SOFIA). The polarization of M82 at $53~\mu \rm{m}$ clearly shows a magnetic field geometry perpendicular to the disk in the hot dust emission. For M82 the polarization at $154~\mu \rm{m}$ shows a combination of field geometry perpendicular to the disk in the nuclear region, but closer to parallel to the disk away from the nucleus. The fractional polarization at $53~\mu \rm{m}$ $(154~\mu \rm{m})$ ranges from 7% (3%) off nucleus to 0.5% (0.3%) near the nucleus. A simple interpretation of the observations of M82 invokes a massive polar outflow, dragging the field along, from a region $\sim 700$~pc in diameter that has entrained some of the gas and dust, creating a vertical field geometry seen mostly in the hotter $(53~\mu \rm{m})$ dust emission. This outflow sits within a larger disk with a more typical planar geometry that more strongly contributes to the cooler $(154~\mu \rm{m})$ dust emission. For NGC 253, the polarization at $89~\mu \rm{m}$ is dominated by a planar geometry in the tilted disk, with weak indication of a vertical geometry above and below the plane from the nucleus. The polarization observations of NGC 253 at $53~\mu \rm{m}$ were of insufficient S/N for detailed analysis.

We present the detection at 89 $\mu$m (observed frame) of the {\it Herschel}-selected gravitationally lensed starburst galaxy HATLASJ1429-0028 (also known as G15v2.19) in 15 minutes with the High-resolution Airborne Wideband Camera-plus (HAWC+) onboard the Stratospheric Observatory for Infrared Astronomy (SOFIA). The spectacular lensing system consists of an edge-on foreground disk galaxy at $z$ = 0.22 and a nearly complete Einstein ring of an intrinsic ultra-luminous infrared galaxy at $z$ = 1.03. Is this high IR luminosity powered by pure star formation (SF) or also an active galactic nucleus (AGN)? Previous nebular line diagnostics indicate that it is star-formation dominated. We perform a 27-band multi-wavelength spectral energy distribution modeling (SED) including the new SOFIA/HAWC+ data to constrain the fractional AGN contribution to the total IR luminosity. The AGN fraction in the IR turns out to be negligible. In addition, J1429-0028 serves as a testbed for comparing SED results from different models/templates and SED codes (MAGPHYS, SED3FIT, and CIGALE). We stress that star formation history is the dominant source of uncertainty in the derived stellar mass (as high as a factor of $\sim$ 10) even in the case of extensive photometric coverage. Furthermore, the detection of a source at $z$ $\sim$ 1 with SOFIA/HAWC+ demonstrates the potential of utilizing this facility for distant galaxy studies including the decomposition of SF/AGN components, which cannot be accomplished with other current facilities.

We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fields in star formation for nearby star-forming molecular clouds. We present a first look at the geometry and strength of magnetic fields in Oph-B. The field geometry is traced over $\sim$0.2 pc, with clear detection of both of the sub-clumps of Oph-B. The field pattern appears significantly disordered in sub-clump Oph-B1. The field geometry in Oph-B2 is more ordered, with a tendency to be along the major axis of the clump, parallel to the filamentary structure within which it lies. The degree of polarization decreases systematically towards the dense core material in the two sub-clumps. The field lines in the lower density material along the periphery are smoothly joined to the large scale magnetic fields probed by NIR polarization observations. We estimated a magnetic field strength of 630$\pm$410 $\mu$G in the Oph-B2 sub-clump using a Davis-Chandeasekhar-Fermi analysis. With this magnetic field strength, we find a mass-to-flux ratio $\lambda$= 1.6$\pm$1.1, which suggests that the Oph-B2 clump is slightly magnetically supercritical.

We present 850 $\mu$m imaging polarimetry data of the $\rho$ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (B-fields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We identify 10 subregions with distinct polarization fractions and angles in the 0.2 pc $\rho$ Oph A core; some of them can be part of a coherent magnetic field structure in the $\rho$ Oph region. The results are consistent with previous observations of the brightest regions of $\rho$ Oph-A, where the degrees of polarization are at a level of a few percents, but our data reveal for the first time the magnetic field structures in the fainter regions surrounding the core where the degree of polarization is much higher ($> 5 \%$). A comparison with previous near-infrared polarimetric data shows that there are several magnetic field components which are consistent at near-infrared and submillimeter wavelengths. Using the Davis-Chandrasekhar-Fermi method, we also derive magnetic field strengths in several sub-core regions, which range from approximately 0.2 to 5 mG. We also find a correlation between the magnetic field orientations projected on the sky with the core centroid velocity components.

We present the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will aim to answer. The most important of these is the role of magnetic fields in the star formation process on the scale of individual filaments and cores in dense regions. We describe the data acquisition and reduction processes for POL-2, demonstrating both repeatability and consistency with previous data. We present a first-look analysis of the first results from the BISTRO survey in the OMC 1 region. We see that the magnetic field lies approximately perpendicular to the famous 'integral filament' in the densest regions of that filament. Furthermore, we see an 'hour-glass' magnetic field morphology extending beyond the densest region of the integral filament into the less-dense surrounding material, and discuss possible causes for this. We also discuss the more complex morphology seen along the Orion Bar region. We examine the morphology of the field along the lower-density north-eastern filament. We find consistency with previous theoretical models that predict magnetic fields lying parallel to low-density, non-self-gravitating filaments, and perpendicular to higher-density, self-gravitating filaments.

Variable emission from Sgr~A*, the luminous counterpart to the super-massive black hole at the center of our Galaxy, arises from the innermost portions of the accretion flow. Better characterization of the variability is important for constraining models of the low-luminosity accretion mode powering Sgr~A*, and could further our ability to use variable emission as a probe of the strong gravitational potential in the vicinity of the $4\times10^{6}\mathrm{M}_{\odot}$ black hole. We use the \textit{Herschel} Spectral and Photometric Imaging Receiver (SPIRE) to monitor Sgr~A* at wavelengths that are difficult or impossible to observe from the ground. We find highly significant variations at 0.25, 0.35, and 0.5 mm, with temporal structure that is highly correlated across these wavelengths. While the variations correspond to $<$1% changes in the total intensity in the \textit{Herschel} beam containing Sgr~A*, comparison to independent, simultaneous observations at 0.85 mm strongly supports the reality of the variations. The lowest point in the light curves, $\sim$0.5 Jy below the time-averaged flux density, places a lower bound on the emission of Sgr~A* at 0.25 mm, the first such constraint on the THz portion of the SED. The variability on few hour timescales in the SPIRE light curves is similar to that seen in historical 1.3 mm data, where the longest time series is available, but the distribution of variations in the sub-mm do not show a tail of large-amplitude variations seen at 1.3 mm. Simultaneous X-ray photometry from XMM-Newton shows no significant variation within our observing period, which may explain the lack of very large variations if X-ray and submillimeter flares are correlated.

In a companion paper we have reported a $>5\sigma$ detection of degree scale $B $-mode polarization at 150 GHz by the BICEP2 experiment. Here we provide a detailed study of potential instrumental systematic contamination to that measurement. We focus extensively on spurious polarization that can potentially arise from beam imperfections. We present a heuristic classification of beam imperfections according to their symmetries and uniformities, and discuss how resulting contamination adds or cancels in maps that combine observations made at multiple orientations of the telescope about its boresight axis. We introduce a technique, which we call "deprojection", for filtering the leading order beam-induced contamination from time ordered data, and show that it removes power from BICEP2's $BB$ spectrum consistent with predictions using high signal-to-noise beam shape measurements. We detail the simulation pipeline that we use to directly simulate instrumental systematics and the calibration data used as input to that pipeline. Finally, we present the constraints on $BB$ contamination from individual sources of potential systematics. We find that systematics contribute $BB$ power that is a factor $\sim10\times$ below BICEP2's 3-year statistical uncertainty, and negligible compared to the observed $BB$ signal. The contribution to the best-fit tensor/scalar ratio is at a level equivalent to $r=(3-6)\times10^{-3}$.

We used wide area surveys over 39 deg$^2$ by the HerMES collaboration, performed with the Herschel Observatory SPIRE multi-wavelength camera, to estimate the low-redshift, $0.02<z<0.5$, monochromatic luminosity functions (LFs) of galaxies at 250, 350 and 500$\,\mu$m. SPIRE flux densities were also combined with Spitzer photometry and multi-wavelength archival data to perform a complete SED fitting analysis of SPIRE detected sources to calculate precise k-corrections, as well as the bolometric infrared (8-1000$\,\mu$m) luminosity functions and their low-$z$ evolution from a combination of statistical estimators. Integration of the latter prompted us to also compute the local luminosity density (LLD) and the comoving star formation rate density (SFRD) for our sources, and to compare them with theoretical predictions of galaxy formation models. The luminosity functions show significant and rapid luminosity evolution already at low redshifts, $0.02<z<0.2$, with L$_{IR}^* \propto (1+z)^{6.0\pm0.4}$ and $\Phi_{IR}^* \propto (1+z)^{-2.1\pm0.4}$, L$_{250}^* \propto (1+z)^{5.3\pm0.2}$ and $\Phi_{250}^* \propto (1+z)^{-0.6\pm0.4}$ estimated using the IR bolometric and the 250$\,\mu$m LFs respectively. Converting our IR LD estimate into an SFRD assuming a standard Salpeter IMF and including the unobscured contribution based on the UV dust-uncorrected emission from local galaxies, we estimate a SFRD scaling of SFRD$_0+0.08 z$, where SFRD$_0\simeq (1.9\pm 0.03)\times 10^{-2} [\mathrm{M}_\odot\,\mathrm{Mpc}^{-3}]$ is our total SFRD estimate at $z\sim0.02$.

The Keck Array is a system of cosmic microwave background (CMB) polarimeters, each similar to the BICEP2 experiment. In this paper we report results from the 2012 and 2013 observing seasons, during which the Keck Array consisted of five receivers all operating in the same (150 GHz) frequency band and observing field as BICEP2. We again find an excess of B-mode power over the lensed-$\Lambda$CDM expectation of $> 5 \sigma$ in the range $30 < \ell < 150$ and confirm that this is not due to systematics using jackknife tests and simulations based on detailed calibration measurements. In map difference and spectral difference tests these new data are shown to be consistent with BICEP2. Finally, we combine the maps from the two experiments to produce final Q and U maps which have a depth of 57 nK deg (3.4 $\mu$K arcmin) over an effective area of 400 deg$^2$ for an equivalent survey weight of 250,000 $\mu$K$^{-2}$. The final BB band powers have noise uncertainty a factor of 2.3 times better than the previous results, and a significance of detection of excess power of $> 6\sigma$.

BICEP2 and the Keck Array are polarization-sensitive microwave telescopes that observe the cosmic microwave background (CMB) from the South Pole at degree angular scales in search of a signature of inflation imprinted as B-mode polarization in the CMB. BICEP2 was deployed in late 2009, observed for three years until the end of 2012 at 150 GHz with 512 antenna-coupled transition edge sensor bolometers, and has reported a detection of B-mode polarization on degree angular scales. The Keck Array was first deployed in late 2010 and will observe through 2016 with five receivers at several frequencies (95, 150, and 220 GHz). BICEP2 and the Keck Array share a common optical design and employ the field-proven BICEP1 strategy of using small-aperture, cold, on-axis refractive optics, providing excellent control of systematics while maintaining a large field of view. This design allows for full characterization of far-field optical performance using microwave sources on the ground. Here we describe the optical design of both instruments and report a full characterization of the optical performance and beams of BICEP2 and the Keck Array at 150 GHz.

BICEP1 is a millimeter-wavelength telescope designed specifically to measure the inflationary B-mode polarization of the Cosmic Microwave Background (CMB) at degree angular scales. We present results from an analysis of the data acquired during three seasons of observations at the South Pole (2006 to 2008). This work extends the two-year result published in Chiang et al. (2010), with additional data from the third season and relaxed detector-selection criteria. This analysis also introduces a more comprehensive estimation of band-power window functions, improved likelihood estimation methods and a new technique for deprojecting monopole temperature-to-polarization leakage which reduces this class of systematic uncertainty to a negligible level. We present maps of temperature, E- and B-mode polarization, and their associated angular power spectra. The improvement in the map noise level and polarization spectra error bars are consistent with the 52% increase in integration time relative to Chiang et al. (2010). We confirm both self-consistency of the polarization data and consistency with the two-year results. We measure the angular power spectra at 21 <= l <= 335 and find that the EE spectrum is consistent with Lambda Cold Dark Matter (LCDM) cosmology, with the first acoustic peak of the EE spectrum now detected at 15sigma. The BB spectrum remains consistent with zero. From B-modes only, we constrain the tensor-to-scalar ratio to r = 0.03+0.27-0.23, or r < 0.70 at 95% confidence level.

We report the results of a joint analysis of data from BICEP2/Keck Array and Planck. BICEP2 and Keck Array have observed the same approximately 400 deg$^2$ patch of sky centered on RA 0h, Dec. $-57.5\deg$. The combined maps reach a depth of 57 nK deg in Stokes $Q$ and $U$ in a band centered at 150 GHz. Planck has observed the full sky in polarization at seven frequencies from 30 to 353 GHz, but much less deeply in any given region (1.2 $\mu$K deg in $Q$ and $U$ at 143 GHz). We detect 150$\times$353 cross-correlation in $B$-modes at high significance. We fit the single- and cross-frequency power spectra at frequencies $\geq 150$ GHz to a lensed-$\Lambda$CDM model that includes dust and a possible contribution from inflationary gravitational waves (as parameterized by the tensor-to-scalar ratio $r$), using a prior on the frequency spectral behavior of polarized dust emission from previous \planck\ analysis of other regions of the sky. We find strong evidence for dust and no statistically significant evidence for tensor modes. We probe various model variations and extensions, including adding a synchrotron component in combination with lower frequency data, and find that these make little difference to the $r$ constraint. Finally we present an alternative analysis which is similar to a map-based cleaning of the dust contribution, and show that this gives similar constraints. The final result is expressed as a likelihood curve for $r$, and yields an upper limit $r_{0.05}<0.12$ at 95% confidence. Marginalizing over dust and $r$, lensing $B$-modes are detected at $7.0\,\sigma$ significance.

We have developed antenna-coupled transition-edge sensor (TES) bolometers for a wide range of cosmic microwave background (CMB) polarimetry experiments, including BICEP2, Keck Array, and the balloon borne SPIDER. These detectors have reached maturity and this paper reports on their design principles, overall performance, and key challenges associated with design and production. Our detector arrays repeatedly produce spectral bands with 20%-30% bandwidth at 95, 150, or 220~GHz. The integrated antenna arrays synthesize symmetric co-aligned beams with controlled side-lobe levels. Cross-polarized response on boresight is typically ~0.5%, consistent with cross-talk in our multiplexed readout system. End-to-end optical efficiencies in our cameras are routinely 35% or higher, with per detector sensitivities of NET~300 uKrts. Thanks to the scalability of this design, we have deployed 2560 detectors as 1280 matched pairs in Keck Array with a combined instantaneous sensitivity of ~9 uKrts, as measured directly from CMB maps in the 2013 season. Similar arrays have recently flown in the SPIDER instrument, and development of this technology is ongoing.

Photometric instruments operating at far infrared to millimetre wavelengths often have broad spectral passbands (central wavelength/bandwidth ~ 3 or less), especially those operating in space. A broad passband can result in significant variation of the beam profile and aperture efficiency across the passband, effects which thus far have not generally been taken into account in the flux calibration of such instruments. With absolute calibration uncertainties associated with the brightness of primary calibration standards now in the region of 5% or less, variation of the beam properties across the passband can be a significant contributor to the overall calibration accuracy for extended emission. We present a calibration framework which takes such variations into account for both antenna-coupled and absorber-coupled focal plane architectures. The scheme covers point source and extended source cases, and also the intermediate case of a semi-extended source profile. We apply the new method to the Herschel-SPIRE space-borne photometer.

We report on the design and performance of the BICEP2 instrument and on its three-year data set. BICEP2 was designed to measure the polarization of the cosmic microwave background (CMB) on angular scales of 1 to 5 degrees ($\ell$=40-200), near the expected peak of the B-mode polarization signature of primordial gravitational waves from cosmic inflation. Measuring B-modes requires dramatic improvements in sensitivity combined with exquisite control of systematics. The BICEP2 telescope observed from the South Pole with a 26~cm aperture and cold, on-axis, refractive optics. BICEP2 also adopted a new detector design in which beam-defining slot antenna arrays couple to transition-edge sensor (TES) bolometers, all fabricated on a common substrate. The antenna-coupled TES detectors supported scalable fabrication and multiplexed readout that allowed BICEP2 to achieve a high detector count of 500 bolometers at 150 GHz, giving unprecedented sensitivity to B-modes at degree angular scales. After optimization of detector and readout parameters, BICEP2 achieved an instrument noise-equivalent temperature of 15.8 $\mu$K sqrt(s). The full data set reached Stokes Q and U map depths of 87.2 nK in square-degree pixels (5.2 $\mu$K arcmin) over an effective area of 384 square degrees within a 1000 square degree field. These are the deepest CMB polarization maps at degree angular scales to date. The power spectrum analysis presented in a companion paper has resulted in a significant detection of B-mode polarization at degree scales.

We present the second generation BLASTbus electronics. The primary purposes of this system are detector readout, attitude control, and cryogenic housekeeping, for balloon-borne telescopes. Readout of neutron transmutation doped germanium (NTD-Ge) bolometers requires low noise and parallel acquisition of hundreds of analog signals. Controlling a telescope's attitude requires the capability to interface to a wide variety of sensors and motors, and to use them together in a fast, closed loop. To achieve these different goals, the BLASTbus system employs a flexible motherboard-daughterboard architecture. The programmable motherboard features a digital signal processor (DSP) and field-programmable gate array (FPGA), as well as slots for three daughterboards. The daughterboards provide the interface to the outside world, with versions for analog to digital conversion, and optoisolated digital input/output. With the versatility afforded by this design, the BLASTbus also finds uses in cryogenic, thermometry, and power systems. For accurate timing control to tie everything together, the system operates in a fully synchronous manner. BLASTbus electronics have been successfully deployed to the South Pole, and flown on stratospheric balloons.

(abridged for arXiv) We report results from the BICEP2 experiment, a cosmic microwave background (CMB) polarimeter specifically designed to search for the signal of inflationary gravitational waves in the B-mode power spectrum around $\ell\sim80$. The telescope comprised a 26 cm aperture all-cold refracting optical system equipped with a focal plane of 512 antenna coupled transition edge sensor 150 GHz bolometers each with temperature sensitivity of $\approx300\mu\mathrm{K}_\mathrm{CMB}\sqrt{s}$. BICEP2 observed from the South Pole for three seasons from 2010 to 2012. A low-foreground region of sky with an effective area of 380 square deg was observed to a depth of 87 nK deg in Stokes $Q$ and $U$. We find an excess of $B$-mode power over the base lensed-LCDM expectation in the range $30< \ell< 150$, inconsistent with the null hypothesis at a significance of $> 5\sigma$. Through jackknife tests and simulations we show that systematic contamination is much smaller than the observed excess. We also examine a number of available models of polarized dust emission and find that at their default parameter values they predict power $\sim(5-10)\times$ smaller than the observed excess signal. However, these models are not sufficiently constrained to exclude the possibility of dust emission bright enough to explain the entire excess signal. Cross correlating BICEP2 against 100 GHz maps from the BICEP1 experiment, the excess signal is confirmed and its spectral index is found to be consistent with that of the CMB, disfavoring dust at $1.7\sigma$. The observed $B$-mode power spectrum is well fit by a lensed-LCDM + tensor theoretical model with tensor-to-scalar ratio $r=0.20^{+0.07}_{-0.05}$, with $r=0$ disfavored at $7.0\sigma$. Accounting for the contribution of foreground dust will shift this value downward by an amount which will be better constrained with upcoming data sets.

Cosmic Microwave Background (CMB) polarimeters aspire to measure the faint $B$-mode signature predicted to arise from inflationary gravitational waves. They also have the potential to constrain cosmic birefringence which would produce non-zero expectation values for the CMB's $TB$ and $EB$ spectra. However, instrumental systematic effects can also cause these $TB$ and $EB$ correlations to be non-zero. In particular, an overall miscalibration of the polarization orientation of the detectors produces $TB$ and $EB$ spectra which are degenerate with isotropic cosmological birefringence, while also introducing a small but predictable bias on the $BB$ spectrum. The \bicep three-year spectra, which use our standard calibration of detector polarization angles from a dielectric sheet, are consistent with a polarization rotation of $\alpha = -2.77^\circ \pm 0.86^\circ \text{(statistical)} \pm 1.3^\circ \text{(systematic)}$. We revise the estimate of systematic error on the polarization rotation angle from the two-year analysis by comparing multiple calibration methods. We investigate the polarization rotation for the \bicep 100 GHz and 150 GHz bands separately to investigate theoretical models that produce frequency-dependent cosmic birefringence. We find no evidence in the data supporting either these models or Faraday rotation of the CMB polarization by the Milky Way galaxy's magnetic field. If we assume that there is no cosmic birefringence, we can use the $TB$ and $EB$ spectra to calibrate detector polarization orientations, thus reducing bias of the cosmological $B$-mode spectrum from leaked $E$-modes due to possible polarization orientation miscalibration. After applying this "self-calibration" process, we find that the upper limit on the tensor-to-scalar ratio decreases slightly, from $r<0.70$ to $r<0.65$ at $95\%$ confidence.

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

The old, red stars which constitute the bulges of galaxies, and the massive black holes at their centres, are the relics of a period in cosmic history when galaxies formed stars at remarkable rates and active galactic nuclei (AGN) shone brightly from accretion onto black holes. It is widely suspected, but unproven, that the tight correlation in mass of the black hole and stellar components results from the AGN quenching the surrounding star formation as it approaches its peak luminosity. X-rays trace emission from AGN unambiguously, while powerful star-forming galaxies are usually dust-obscured and are brightest at infrared to submillimetre wavelengths. Here we report observations in the submillimetre and X-ray which show that rapid star formation was common in the host galaxies of AGN when the Universe was 2-6 Gyrs old, but that the most vigorous star formation is not observed around black holes above an X-ray luminosity of 10^44 erg/s. This suppression of star formation in the host galaxies of powerful AGN is a key prediction of models in which the AGN drives a powerful outflow, expelling the interstellar medium of its host galaxy and transforming the galaxy's properties in a brief period of cosmic time.

We exploit the deep and extended far infrared data sets (at 70, 100 and 160 um) of the Herschel GTO PACS Evolutionary Probe (PEP) Survey, in combination with the HERschel Multi tiered Extragalactic Survey (HerMES) data at 250, 350 and 500 um, to derive the evolution of the restframe 35 um, 60 um, 90 um, and total infrared (IR) luminosity functions (LFs) up to z~4. We detect very strong luminosity evolution for the total IR LF combined with a density evolution. In agreement with previous findings, the IR luminosity density increases steeply to z~1, then flattens between z~1 and z~3 to decrease at z greater than 3. Galaxies with different SEDs, masses and sSFRs evolve in very different ways and this large and deep statistical sample is the first one allowing us to separately study the different evolutionary behaviours of the individual IR populations contributing to the IR luminosity density. Galaxies occupying the well established SFR/stellar mass main sequence (MS) are found to dominate both the total IR LF and luminosity density at all redshifts, with the contribution from off MS sources (0.6 dex above MS) being nearly constant (~20% of the total IR luminosity density) and showing no significant signs of increase with increasing z over the whole 0.8<z<2.2 range. Sources with mass in the 10< log(M/Msun) <11 range are found to dominate the total IR LF, with more massive galaxies prevailing at the bright end of the high-z LF. A two-fold evolutionary scheme for IR galaxies is envisaged: on the one hand, a starburst-dominated phase in which the SMBH grows and is obscured by dust, is followed by an AGN dominated phase, then evolving toward a local elliptical. On the other hand, moderately starforming galaxies containing a low-luminosity AGN have various properties suggesting they are good candidates for systems in a transition phase preceding the formation of steady spiral galaxies.

We have observed four massive galaxy clusters with the SPIRE instrument on the Herschel Space Observatory and measure a deficit of surface brightness within their central region after subtracting sources. We simulate the effects of instrumental sensitivity and resolution, the source population, and the lensing effect of the clusters to estimate the shape and amplitude of the deficit. The amplitude of the central deficit is a strong function of the surface density and flux distribution of the background sources. We find that for the current best fitting faint end number counts, and excellent lensing models, the most likely amplitude of the central deficit is the full intensity of the cosmic infrared background (CIB). Our measurement leads to a lower limit to the integrated total intensity of the CIB of I(250 microns) > 0.69_(-0.03)^(+0.03) (stat.)_(-0.06)^(+0.11) (sys.) MJy/sr, with more CIB possible from both low-redshift sources and from sources within the target clusters. It should be possible to observe this effect in existing high angular resolution data at other wavelengths where the CIB is bright, which would allow tests of models of the faint source component of the CIB.

A set of six debris disk candidates identified with IRAS or WISE excesses were observed at either 350 um or 450 um with the CSO. Five of the targets - HIP 51658, HIP 68160, HIP 73512, HIP 76375, and HIP 112460 - have among the largest measured excess emission from cold dust from IRAS in the 25-100 um bands. Single temperature blackbody fits to the excess dust emission of these sources predict 350-450 um fluxes above 240 mJy. The final target - HIP 73165 - exhibits weak excess emission above the stellar photosphere from WISE measurements at 22 um, indicative of a population of warm circumstellar dust. None of the six targets were detected, with 3 sigma upper limits ranging from 51-239 mJy. These limits are significantly below the expected fluxes from SED fitting. Two potential causes of the null detections were explored - companion stars and contamination. To investigate the possible influence of companion stars, imaging data were analyzed from new AO data from the MMT and archival HST, NIRI, and POSS/2MASS data. The images are sensitive to all stellar companions beyond a radius of 1-94 AU. One target is identified as a binary system, but with a separation too large to impact the disk. While the gravitational effects of a companion do not appear to provide an explanation for the submm upper limits, the majority of the IRAS excess targets show evidence for contaminating sources, based on investigation of higher resolution WISE and archival Spitzer and Herschel images. Finally, the exploratory submm measurements of the WISE excess source suggest that the hot dust present around these targets is not matched by a comparable population of colder, outer dust. More extensive and more sensitive Herschel observations of WISE excess sources will build upon this initial example to further define the characteristics of warm debris disks sources.

We describe the procedure used to flux calibrate the three-band submillimetre photometer in the Spectral and Photometric Imaging REceiver (SPIRE) instrument on the Herschel Space Observatory. This includes the equations describing the calibration scheme, a justification for using Neptune as the primary calibration source, a description of the observations and data processing procedures used to derive flux calibration parameters (for converting from voltage to flux density) for every bolometer in each array, an analysis of the error budget in the flux calibration for the individual bolometers, and tests of the flux calibration on observations of primary and secondary calibrators. The procedure for deriving the flux calibration parameters is divided into two parts. In the first part, we use observations of astronomical sources in conjunction with the operation of the photometer internal calibration source to derive the unscaled derivatives of the flux calibration curves. To scale the calibration curves in Jy/beam/V, we then use observations of Neptune in which the beam of each bolometer is mapped using Neptune observed in a very fine scan pattern. The total instrumental uncertainties in the flux calibration for the individual bolometers is ~0.5% for most bolometers, although a few bolometers have uncertainties of ~1-5% because of issues with the Neptune observations. Based on application of the flux calibration parameters to Neptune observations performed using typical scan map observing modes, we determined that measurements from each array as a whole have instrumental uncertainties of 1.5%. This is considerably less than the absolute calibration uncertainty associated with the model of Neptune, which is estimated at 4%.

We present measurements of the auto- and cross-frequency power spectra of the cosmic infrared background (CIB) at 250, 350, and 500um (1200, 860, and 600 GHz) from observations totaling ~ 70 deg^2 made with the SPIRE instrument aboard the Herschel Space Observatory. We measure a fractional anisotropy dI / I = 14 +- 4%, detecting signatures arising from the clustering of dusty star-forming galaxies in both the linear (2-halo) and non-linear (1-halo) regimes; and that the transition from the 2- to 1-halo terms, below which power originates predominantly from multiple galaxies within dark matter halos, occurs at k_theta ~ 0.1 - 0.12 arcmin^-1 (l ~ 2160 - 2380), from 250 to 500um. New to this paper is clear evidence of a dependence of the Poisson and 1-halo power on the flux-cut level of masked sources --- suggesting that some fraction of the more luminous sources occupy more massive halos as satellites, or are possibly close pairs. We measure the cross-correlation power spectra between bands, finding that bands which are farthest apart are the least correlated, as well as hints of a reduction in the correlation between bands when resolved sources are more aggressively masked. In the second part of the paper we attempt to interpret the measurements in the framework of the halo model. With the aim of fitting simultaneously with one model the power spectra, number counts, and absolute CIB level in all bands, we find that this is achievable by invoking a luminosity-mass relationship, such that the luminosity-to-mass ratio peaks at a particular halo mass scale and declines towards lower and higher mass halos. Our best-fit model finds that the halo mass which is most efficient at hosting star formation in the redshift range of peak star-forming activity, z ~ 1-3, is log(M_peak/M_sun) ~ 12.1 +- 0.5, and that the minimum halo mass to host infrared galaxies is log(M_min/M_sun) ~ 10.1 +- 0.6.

ABRIGED Herschel/SPIRE has provided confusion limited maps of deep fields at 250, 350, and 500um, as part of the HerMES survey. Due to confusion, only a small fraction of the Cosmic Infrared Background can be resolved into individually-detected sources. Our goal is to produce deep galaxy number counts and redshift distributions below the confusion limit, which we then use to place strong constraints on the origins of the cosmic infrared background and on models of galaxy evolution. We individually extracted the bright SPIRE with a method using the positions, the flux densities, and the redshifts of the 24um sources as a prior, and derived the number counts and redshift distributions of the bright SPIRE sources. For fainter SPIRE sources, we reconstructed the number counts and the redshift distribution below the confusion limit using the deep 24um catalogs associated with photometric redshift and information provided by the stacking of these sources into the deep SPIRE maps. Finally, by integrating all these counts, we studied the contribution of the galaxies to the CIB as a function of their flux density and redshift. Through stacking, we managed to reconstruct the source counts per redshift slice down to ~2 mJy in the three SPIRE bands, which lies about a factor 10 below the 5sigma confusion limit. None of the pre-existing population models are able to reproduce our results at better than 3sigma. Finally, we extrapolate our counts to zero flux density in order to derive an estimate of the total contribution of galaxies to the CIB, finding 10.1, 6.5, and 2.8 nW/m2/sr at 250, 350, and 500um, respectively. These values agree well with FIRAS absolute measurements, suggesting our number counts and their extrapolation are sufficient to explain the CIB. Finally, combining our results with other works, we estimate the energy budget contained in the CIB between 8 and 1000um: 26 nW/m2/sr.

Massive present-day early-type (elliptical and lenticular) galaxies probably gained the bulk of their stellar mass and heavy elements through intense, dust-enshrouded starbursts - that is, increased rates of star formation - in the most massive dark matter halos at early epochs. However, it remains unknown how soon after the Big Bang such massive starburst progenitors exist. The measured redshift distribution of dusty, massive starbursts has long been suspected to be biased low in redshift owing to selection effects, as confirmed by recent findings of systems out to redshift z~5. Here we report the identification of a massive starburst galaxy at redshift 6.34 through a submillimeter color-selection technique. We unambiguously determined the redshift from a suite of molecular and atomic fine structure cooling lines. These measurements reveal a hundred billion solar masses of highly excited, chemically evolved interstellar medium in this galaxy, which constitutes at least 40% of the baryonic mass. A "maximum starburst" converts the gas into stars at a rate more than 2,000 times that of the Milky Way, a rate among the highest observed at any epoch. Despite the overall downturn of cosmic star formation towards the highest redshifts, it seems that environments mature enough to form the most massive, intense starbursts existed at least as early as 880 million years after the Big Bang.

Using Herschel data from the deepest SPIRE and PACS surveys (HerMES and PEP) in COSMOS and GOODS (N+S), we examine the dust properties of IR-luminous (L_IR>10^10 L_sun) galaxies at 0.1<z<2 and determine how these evolve with cosmic time. The unique angle of this work is the rigorous analysis of survey selection effects, making this the first study of the star-formation-dominated, IR-luminous population within a framework almost entirely free of selection biases. We find that IR-luminous galaxies have SEDs with broad far-IR peaks characterised by cool/extended dust emission and average dust temperatures in the 25-45K range. Hot (T>45K) SEDs and cold (T<25K), cirrus-dominated SEDs are rare, with most sources being within the range occupied by warm starbursts such as M82 and cool spirals such as M51. We observe a luminosity-temperature (L-T) relation, where the average dust temperature of log[L_IR/L_sun]=12.5 galaxies is about 10K higher than that of their log[L_IR/L_sun]=10.5 counterparts. However, although the increased dust heating in more luminous systems is the driving factor behind the L-T relation, the increase in dust mass and/or starburst size with luminosity plays a dominant role in shaping it. Our results show that the dust conditions in IR-luminous sources evolve with cosmic time: at high redshift, dust temperatures are on average up to 10K lower than what is measured locally. This is manifested as a flattening of the L-T relation, suggesting that (U)LIRGs in the early Universe are typically characterised by a more extended dust distribution and/or higher dust masses than local equivalent sources. Interestingly, the evolution in dust temperature is luminosity dependent, with the fraction of LIRGs with T<35K showing a 2-fold increase from z~0 to z~2, whereas that of ULIRGs with T<35K shows a 6-fold increase.

The design and performance of a wide bandwidth linear polarization modulator based on the Faraday effect is described. Faraday Rotation Modulators (FRMs) are solid-state polarization switches that are capable of modulation up to ~10 kHz. Six FRMs were utilized during the 2006 observing season in the Background Imaging of Cosmic Extragalactic Polarization (BICEP) experiment; three FRMs were used at each of BICEP's 100 and 150 GHz frequency bands. The technology was verified through high signal-to-noise detection of Galactic polarization using two of the six FRMs during four observing runs in 2006. The features exhibit strong agreement with BICEP's measurements of the Galaxy using non-FRM pixels and with the Galactic polarization models. This marks the first detection of high signal-to-noise mm-wave celestial polarization using fast, active optical modulation. The performance of the FRMs during periods when they were not modulated was also analyzed and compared to results from BICEP's 43 pixels without FRMs.

We present a list of 13 candidate gravitationally lensed submillimeter galaxies (SMGs) from 95 square degrees of the Herschel Multi-tiered Extragalactic Survey, a surface density of 0.14\pm0.04deg^{-2}. The selected sources have 500um flux densities (S_500) greater than 100mJy. Gravitational lensing is confirmed by follow-up observations in 9 of the 13 systems (70%), and the lensing status of the four remaining sources is undetermined. We also present a supplementary sample of 29 (0.31\pm0.06deg^{-2}) gravitationally lensed SMG candidates with S_500=80--100mJy, which are expected to contain a higher fraction of interlopers than the primary candidates. The number counts of the candidate lensed galaxies are consistent with a simple statistical model of the lensing rate, which uses a foreground matter distribution, the intrinsic SMG number counts, and an assumed SMG redshift distribution. The model predicts that 32--74% of our S_500>100mJy candidates are strongly gravitationally lensed (mu>2), with the brightest sources being the most robust; this is consistent with the observational data. Our statistical model also predicts that, on average, lensed galaxies with S_500=100mJy are magnified by factors of ~9, with apparently brighter galaxies having progressively higher average magnification, due to the shape of the intrinsic number counts. 65% of the sources are expected to have intrinsic 500micron flux densities less than 30mJy. Thus, samples of strongly gravitationally lensed SMGs, such as those presented here, probe below the nominal Herschel detection limit at 500 micron. They are good targets for the detailed study of the physical conditions in distant dusty, star-forming galaxies, due to the lensing magnification, which can lead to spatial resolutions of ~0.01" in the source plane.

Submillimeter cameras now have up to $10^4$ pixels (SCUBA 2). The proposed CCAT 25-meter submillimeter telescope will feature a 1 degree field-of-view. Populating the focal plane at 350 microns would require more than $10^6$ photon-noise limited pixels. To ultimately achieve this scaling, simple detectors and high-density multiplexing are essential. We are addressing this long-term challenge through the development of frequency-multiplexed superconducting microresonator detector arrays. These arrays use lumped-element, direct-absorption resonators patterned from titanium nitride films. We will discuss our progress toward constructing a scalable 350 micron pathfinder instrument focusing on fabrication simplicity, multiplexing density, and ultimately a low per-pixel cost.

Between the BICEP2 and Keck Array experiments, we have deployed over 1500 dual polarized antenna coupled bolometers to map the Cosmic Microwave Background's polarization. We have been able to rapidly deploy these detectors because they are completely planar with an integrated phased-array antenna. Through our experience in these experiments, we have learned of several challenges with this technology- specifically the beam synthesis in the antenna- and in this paper we report on how we have modified our designs to mitigate these challenges. In particular, we discus differential steering errors between the polarization pairs' beam centroids due to microstrip cross talk and gradients of penetration depth in the niobium thin films of our millimeter wave circuits. We also discuss how we have suppressed side lobe response with a Gaussian taper of our antenna illumination pattern. These improvements will be used in Spider, Polar-1, and this season's retrofit of Keck Array.

The Keck Array (SPUD) began observing the cosmic microwave background's polarization in the winter of 2011 at the South Pole. The Keck Array follows the success of the predecessor experiments Bicep and Bicep2, using five on-axis refracting telescopes. These have a combined imaging array of 2500 antenna-coupled TES bolometers read with a SQUID-based time domain multiplexing system. We will discuss the detector noise and the optimization of the readout. The achieved sensitivity of the Keck Array is 11.5 {\mu}K_(CMB)*sqrt{s} in the 2012 configuration.

The Keck Array (SPUD) is a set of microwave polarimeters that observes from the South Pole at degree angular scales in search of a signature of Inflation imprinted as B-mode polarization in the Cosmic Microwave Background (CMB). The first three Keck Array receivers were deployed during the 2010-2011 Austral summer, followed by two new receivers in the 2011-2012 summer season, completing the full five-receiver array. All five receivers are currently observing at 150 GHz. The Keck Array employs the field-proven BICEP/BICEP2 strategy of using small, cold, on-axis refractive optics, providing excellent control of systematics while maintaining a large field of view. This design allows for full characterization of far-field optical performance using microwave sources on the ground. We describe our efforts to characterize the main beam shape and beam shape mismatch between co-located orthogonally-polarized detector pairs, and discuss the implications of measured differential beam parameters on temperature to polarization leakage in CMB analysis.

The BICEP2 and Keck Array experiments are designed to measure the polarization of the cosmic microwave background (CMB) on angular scales of 2-4 degrees (l=50-100). This is the region in which the B-mode signal, a signature prediction of cosmic inflation, is expected to peak. BICEP2 was deployed to the South Pole at the end of 2009 and is in the middle of its third year of observing with 500 polarization-sensitive detectors at 150 GHz. The Keck Array was deployed to the South Pole at the end of 2010, initially with three receivers--each similar to BICEP2. An additional two receivers have been added during the 2011-12 summer. We give an overview of the two experiments, report on substantial gains in the sensitivity of the two experiments after post-deployment optimization, and show preliminary maps of CMB polarization from BICEP2.

We measured polarized dust emission at 350um towards the high-mass star forming massive dense clump IRAS 20126+4104 using the SHARC II Polarimeter, SHARP, at the Caltech Submillimeter Observatory. Most of the observed magnetic field vectors agree well with magnetic field vectors obtained from a numerical simulation for the case when the global magnetic field lines are inclined with respect to the rotation axis of the dense clump. The results of the numerical simulation show that rotation plays an important role on the evolution of the massive dense clump and its magnetic field. The direction of the cold CO 1-0 bipolar outflow is parallel to the observed magnetic field within the dense clump as well as the global magnetic field, as inferred from optical polarimetry data, indicating that the magnetic field also plays a critical role in an early stage of massive star formation. The large-scale Keplerian disk of the massive (proto)star rotates in almost opposite sense to the clump's envelope. The observed magnetic field morphology and the counter-rotating feature of the massive dense clump system provide hints to constrain the role of magnetic fields in the process of high mass star formation.

The Herschel Multi-tiered Extragalactic Survey, HerMES, is a legacy program designed to map a set of nested fields totalling ~380 deg^2. Fields range in size from 0.01 to ~20 deg^2, using Herschel-SPIRE (at 250, 350 and 500 \mu m), and Herschel-PACS (at 100 and 160 \mu m), with an additional wider component of 270 deg^2 with SPIRE alone. These bands cover the peak of the redshifted thermal spectral energy distribution from interstellar dust and thus capture the re-processed optical and ultra-violet radiation from star formation that has been absorbed by dust, and are critical for forming a complete multi-wavelength understanding of galaxy formation and evolution. The survey will detect of order 100,000 galaxies at 5\sigma in some of the best studied fields in the sky. Additionally, HerMES is closely coordinated with the PACS Evolutionary Probe survey. Making maximum use of the full spectrum of ancillary data, from radio to X-ray wavelengths, it is designed to: facilitate redshift determination; rapidly identify unusual objects; and understand the relationships between thermal emission from dust and other processes. Scientific questions HerMES will be used to answer include: the total infrared emission of galaxies; the evolution of the luminosity function; the clustering properties of dusty galaxies; and the properties of populations of galaxies which lie below the confusion limit through lensing and statistical techniques. This paper defines the survey observations and data products, outlines the primary scientific goals of the HerMES team, and reviews some of the early results.

We investigate the potential of submm-mm and submm-mm-radio photometric redshifts using a sample of mm-selected sources as seen at 250, 350 and 500 {\mu}m by the SPIRE instrument on Herschel. From a sample of 63 previously identified mm-sources with reliable radio identifications in the GOODS-N and Lockman Hole North fields 46 (73 per cent) are found to have detections in at least one SPIRE band. We explore the observed submm/mm colour evolution with redshift, finding that the colours of mm-sources are adequately described by a modified blackbody with constant optical depth {\tau} = ({\nu}/{\nu}0)^{\beta} where {\beta} = +1.8 and {\nu}0 = c/100 {\mu}m. We find a tight correlation between dust temperature and IR luminosity. Using a single model of the dust temperature and IR luminosity relation we derive photometric redshift estimates for the 46 SPIRE detected mm-sources. Testing against the 22 sources with known spectroscopic, or good quality optical/near-IR photometric, redshifts we find submm/mm photometric redshifts offer a redshift accuracy of |z|/(1+z) = 0.16 (< |z| >= 0.51). Including constraints from the radio-far IR correlation the accuracy is improved to |z|/(1 + z) = 0.15 (< |z| >= 0.45). We estimate the redshift distribution of mm-selected sources finding a significant excess at z > 3 when compared to ~ 850 {\mu}m selected samples.

We describe the generation of single-band point source catalogues from submillimetre Herschel-SPIRE observations taken as part of the Science Demonstration Phase of the Herschel Multi-tiered Extragalactic Survey (HerMES). Flux densities are found by means of peak-finding and the fitting of a Gaussian point-response function. With highly-confused images, careful checks must be made on the completeness and flux density accuracy of the detected sources. This is done by injecting artificial sources into the images and analysing the resulting catalogues. Measured flux densities at which 50 per cent of injected sources result in good detections at (250, 350, 500) {\mu}m range from (11.6, 13.2, 13.1) mJy to (25.7, 27.1, 35.8) mJy, depending on the depth of the observation (where a `good' detection is taken to be one with positional offset less than one full-width half-maximum of the point-response function, and with the measured flux density within a factor of 2 of the flux density of the injected source). This paper acts as a reference for the 2010 July HerMES public data release.

Dynamical interactions between planets and debris disks may sculpt the disk structure and impact planetary orbits, but only a few systems with both imaged planets and spatially resolved debris disks are known. With the Caltech Submm Observatory (CSO), we have observed the HR 8799 debris disk at 350{\mu}m. The 350{\mu}m map is the first spatially resolved measurement of the debris disk encircling the HR 8799 planetary system at this wavelength. Both the flux and size of the emission are consistent with a Kuiper belt of dust extending from ~100-300 AU. Although the resolution of the current map is limited, the map shows an indication of offset asymmetric emission, and several scenarios for this possibility are explored with radiative transfer calculations of a star-disk system and N-body numerical simulations of planet-disk interactions with parameters representative of the HR 8799 system.

The Keck Array is a cosmic microwave background (CMB) polarimeter that will begin observing from the South Pole in late 2010. The initial deployment will consist of three telescopes similar to BICEP2 housed in ultra-compact, pulse tube cooled cryostats. Two more receivers will be added the following year. In these proceedings we report on the design and performance of the Keck cryostat. We also report some initial results on the performance of antenna-coupled TES detectors operating in the presence of a pulse tube. We find that the performance of the detectors is not seriously impacted by the replacement of BICEP2's liquid helium cryostat with a pulse tube cooled cryostat.

We report the detection of CO 5-4, 3-2, and 1-0 emission in the strongly lensed, Herschel/SPIRE-selected submillimeter galaxy (SMG) HLSW-01 at z=2.9574+/-0.0001, using the Plateau de Bure Interferometer, the Combined Array for Research in Millimeter-wave Astronomy, and the Green Bank Telescope. The observations spatially resolve the molecular gas into four lensed images with a maximum separation of ~9", and reveal the internal gas dynamics in this system. We derive lensing-corrected CO line luminosities of L'(CO 1-0) = (4.17+/-0.41), L'(CO 3-2) = (3.96+/-0.20) and L'(CO 5-4) = (3.45+/-0.20) x 10^10 (mu_L/10.9)^-1 Kkm/s pc^2, corresponding to luminosity ratios of r_31 = 0.95+/-0.10, r_53 = 0.87+/-0.06, and r_51 = 0.83+/-0.09. This suggests a total molecular gas mass of Mgas = 3.3 x 10^10 (alpha_CO/0.8) (mu_L/10.9)^-1 Msun. The gas mass, gas mass fraction, gas depletion timescale, star formation efficiency, and specific star formation rate are typical for an SMG. The velocity structure of the gas reservoir suggests that the brightest two lensed images are dynamically resolved projections of the same dust-obscured region in the galaxy that are kinematically offset from the unresolved fainter images. The resolved kinematics appear consistent with the complex velocity structure observed in major, `wet' (i.e., gas-rich) mergers. Major mergers are commonly observed in SMGs, and are likely to be responsible for fueling their intense starbursts at high gas consumption rates. This study demonstrates the level of detail to which galaxies in the early universe can be studied by utilizing the increase in effective spatial resolution and sensitivity provided by gravitational lensing.

We report the discovery of a bright ($f(250\mum) > 400$ mJy), multiply-lensed submillimeter galaxy \obj\ in {\it Herschel}/SPIRE Science Demonstration Phase data from the HerMES project. Interferometric 880\mum\ Submillimeter Array observations resolve at least four images with a large separation of $\sim 9\arcsec$. A high-resolution adaptive optics $K_p$ image with Keck/NIRC2 clearly shows strong lensing arcs. Follow-up spectroscopy gives a redshift of $z=2.9575$, and the lensing model gives a total magnification of $\mu \sim 11 \pm 1$. The large image separation allows us to study the multi-wavelength spectral energy distribution (SED) of the lensed source unobscured by the central lensing mass. The far-IR/millimeter-wave SED is well described by a modified blackbody fit with an unusually warm dust temperature, $88 \pm 3$ K. We derive a lensing-corrected total IR luminosity of $(1.43 \pm 0.09) \times 10^{13}\, \mathrm{L}_{\odot}$, implying a star formation rate of $\sim 2500\, \mathrm{M}_{\odot}\, \mathrm{yr}^{-1}$. However, models primarily developed from brighter galaxies selected at longer wavelengths are a poor fit to the full optical-to-millimeter SED. A number of other strongly lensed systems have already been discovered in early {\it Herschel} data, and many more are expected as additional data are collected.

We present the results of a gravitational lensing analysis of the bright $\zs=2.957$ sub-millimeter galaxy (SMG), HERMES J105751.1+573027 found in {\it Herschel}/SPIRE Science Demonstration Phase data from the Herschel Multi-tiered Extragalactic Survey (HerMES) project. The high resolution imaging available in optical and Near-IR channels, along with CO emission obtained with the Plateau de Bure Interferometer, allow us to precisely estimate the intrinsic source extension and hence estimate the total lensing magnification to be $\mu=10.9\pm 0.7$. We measure the half-light radius $R_{\rm eff}$ of the source in the rest-frame Near-UV and $V$ bands that characterize the unobscured light coming from stars and find $R_{\rm eff,*}= [2.0 \pm 0.1]$ kpc, in good agreement with recent studies on the Submillimeter Galaxy population. This lens model is also used to estimate the size of the gas distribution ($R_{\rm eff,gas}= [1.1\pm0.5]$) kpc by mapping back in the source plane the CO (J=5-4) transition line emission. The lens modeling yields a relatively large Einstein radius $R_{\rm Ein}= 4\farcs10 \pm 0\farcs02$, corresponding to a deflector velocity dispersion of [$483\pm 16] \,\kms$. This shows that HERMES J105751.1+573027 is lensed by a {\it galaxy group-size} dark matter halo at redshift $\zl\sim 0.6$. The projected dark matter contribution largely dominates the mass budget within the Einstein radius with $f_{\rm dm}(<R_{\rm Ein})\sim 80%$. This fraction reduces to $f_{\rm dm}(<R_{\rm eff,G1}\simeq 4.5\kpc)\sim 47%$ within the effective radius of the main deflecting galaxy of stellar mass $M_{\rm *,G1}=[8.5\pm 1.6] \times 10^{11}\msun$. At this smaller scale the dark matter fraction is consistent with results already found for massive lensing ellipticals at $z\sim0.2$ from the SLACS survey.

We report on the redshift measurement and CO line excitation of HERMES J105751.1+573027 (HLSW-01), a strongly lensed submillimeter galaxy discovered in Herschel/SPIRE observations as part of the Herschel Multi-tiered Extragalactic Survey (HerMES). HLSW-01 is an ultra-luminous galaxy with an intrinsic far-infrared luminosity of 1.4x10^(13) solar luminosities, and is lensed by a massive group of galaxies into at least four images with a total magnification of 10.9+/-0.7. With the 100 GHz instantaneous bandwidth of the Z-Spec instrument on the Caltech Submillimeter Observatory, we robustly identify a redshift of z=2.958+/-0.007 for this source, using the simultaneous detection of four CO emission lines (J = 7-6, J = 8-7, J = 9-8, and J = 10-9). Combining the measured line fluxes for these high-J transitions with the J = 1-0, J = 3-2 and J = 5-4 line fluxes measured with the Green Bank Telescope, the Combined Array for Research in Millimeter Astronomy, and the Plateau de Bure Interferometer, respectively, we model the physical properties of the molecular gas in this galaxy. We find that the full CO spectral line energy distribution is well described by warm, moderate-density gas with Tkin = 86-235 K and n(H2) = (1.1-3.5)x10^3 cm^(-3). However, it is possible that the highest-J transitions are tracing a small fraction of very dense gas in molecular cloud cores, and two-component models that include a warm/dense molecular gas phase with Tkin ~ 200 K, n(H2) ~ 10^5 cm^(-3) are also consistent with these data. Higher signal-to-noise measurements of the J(upper) > 7 transitions with high spectral resolution, combined with high spatial resolution CO maps, are needed to improve our understanding of the gas excitation, morphology, and dynamics of this interesting high-redshift galaxy.

In addition to its potential to probe the Inflationary cosmological paradigm, millimeter-wave polarimetry is a powerful tool for studying the Milky Way galaxy's composition and magnetic field structure. Towards this end, presented here are Stokes I, Q, and U maps of the Galactic plane from the millimeter-wave polarimeter BICEP covering the Galactic longitude range 260 - 340 degrees in three atmospheric transmission windows centered on 100, 150, and 220 GHz. The maps sample an optical depth 1 < AV < 30, and are consistent with previous characterizations of the Galactic millimeter-wave frequency spectrum and the large-scale magnetic field structure permeating the interstellar medium. Polarized emission is detected over the entire region within two degrees of the Galactic plane and indicates that the large-scale magnetic field is oriented parallel to the plane of the Galaxy. An observed trend of decreasing polarization fraction with increasing total intensity rules out the simplest model of a constant Galactic magnetic field throughout the Galaxy. Including WMAP data in the analysis, the degree-scale frequency spectrum of Galactic polarization fraction is plotted between 23 and 220 GHz for the first time. A generally increasing trend of polarization fraction with electromagnetic frequency is found, which varies from 0.5%-1.5%at frequencies below 50 GHz to 2.5%-3.5%above 90 GHz. The BICEP and WMAP data are fit to a two-component (synchrotron and dust) model showing that the higher frequency BICEP data are necessary to tightly constrain the amplitude and spectral index of Galactic dust. Furthermore, the dust amplitude predicted by this two-component fit is consistent with model predictions of dust emission in the BICEP bands.

The extragalactic background light at far-infrared wavelengths originates from optically-faint, dusty, star-forming galaxies in the universe with star-formation rates at the level of a few hundred solar masses per year. Due to the relatively poor spatial resolution of far-infrared telescopes, the faint sub-millimetre galaxies are challenging to study individually. Instead, their average properties can be studied using statistics such as the angular power spectrum of the background intensity variations. A previous attempt at measuring this power spectrum resulted in the suggestion that the clustering amplitude is below the level computed with a simple ansatz based on a halo model. Here we report a clear detection of the excess clustering over the linear prediction at arcminute angular scales in the power spectrum of brightness fluctuations at 250, 350, and 500 microns. From this excess, we find that sub-millimetre galaxies are located in dark matter halos with a minimum mass of log[M_min/M_sun ]= 11.5^+0.7_-0.2 at 350 microns. This minimum dark matter halo mass corresponds to the most efficient mass scale for star formation in the universe, and is lower than that predicted by semi-analytical models for galaxy formation.

We investigate the effect of the Milky Way's magnetic field in star forming regions using archived 350 micron polarization data on 52 Galactic star formation regions from the Hertz polarimeter module. The polarization angles and percentages for individual telescope beams were combined in order to produce a large-scale average for each source and for complexes of sources. In more than 80% of the sources, we find a meaningful mean magnetic field direction, implying the existence of an ordered magnetic field component at the scale of these sources. The average polarization angles were analyzed with respect to the Galactic coordinates in order to test for correlations between polarization percentage, polarization angle, intensity, and Galactic location. No correlation was found, which suggests that the magnetic field in dense molecular clouds is decoupled from the large-scale Galactic magnetic field. Finally, we show that the magnetic field directions in the complexes are consistent with a random distribution on the sky.

We have investigated the nature of flare emission from Sgr A* during multi-wavelength observations of this source that took place in 2004, 2005 and 2006. We present evidence for dimming of submm and radio flux during the peak of near-IR flares. This suggests that the variability of Sgr A* across its wavelength spectrum is phenomenologically related. The model explaining this new behavior of flare activity could be consistent with adiabatically cooling plasma blobs that are expanding but also partially eclipsing the background quiescent emission from Sgr A*. When a flare is launched, the plasma blob is most compact and is brightest in the optically thin regime whereas the emission in radio/submm wavelengths has a higher opacity. Absorption in the observed light curve of Sgr A* at radio/submm flux is due to the combined effects of lower brightness temperature of plasma blobs with respect to the quiescent brightness temperature and high opacity of plasma blobs. This implies that plasma blobs are mainly placed in the magnetosphere of a disk-like flow or further out in the flow. The depth of the absorption being larger in submm than in radio wavelengths implies that the intrinsic size of the quiescent emission increases with increasing wavelength which is consistent with previous size measurements of Sgr A*. Lastly, we believe that occultation of the quiescent emission of Sgr A* at radio/submm by IR flares can be used as a powerful tool to identify flare activity at its earliest phase of its evolution.

We describe the production and verification of sky maps of the five SPIRE fields observed as part of the Herschel Multi-tiered Extragalactic Survey (HerMES) during the Science Demonstration Phase (SDP) of the Herschel mission. We have implemented an iterative map-making algorithm (SHIM; The SPIRE-HerMES Iterative Mapper) to produce high fidelity maps that preserve extended diffuse emission on the sky while exploiting the repeated observations of the same region of the sky with many detectors in multiple scan directions to minimize residual instrument noise. We specify here the SHIM algorithm and outline the various tests that were performed to determine and characterize the quality of the maps and verify that the astrometry, point source flux and power on all relevant angular scales meets the needs of the HerMES science goals. These include multiple jackknife tests, determination of the map transfer function and detailed examination of the power spectra of both sky and jackknife maps. The map transfer function is approximately unity on scales from one arcminute to one degree. Final maps (v1.0), including multiple jackknives, as well as the SHIM pipeline, have been used by the HerMES team for the production of SDP papers.

Dusty, star forming galaxies contribute to a bright, currently unresolved cosmic far-infrared background. Deep Herschel-SPIRE images designed to detect and characterize the galaxies that comprise this background are highly confused, such that the bulk lies below the classical confusion limit. We analyze three fields from the HerMES programme in all three SPIRE bands (250, 350, and 500 microns); parameterized galaxy number count models are derived to a depth of ~2 mJy/beam, approximately 4 times the depth of previous analyses at these wavelengths, using a P(D) (probability of deflection) approach for comparison to theoretical number count models. Our fits account for 64, 60, and 43 per cent of the far-infrared background in the three bands. The number counts are consistent with those based on individually detected SPIRE sources, but generally inconsistent with most galaxy number counts models, which generically overpredict the number of bright galaxies and are not as steep as the P(D)-derived number counts. Clear evidence is found for a break in the slope of the differential number counts at low flux densities. Systematic effects in the P(D) analysis are explored. We find that the effects of clustering have a small impact on the data, and the largest identified systematic error arises from uncertainties in the SPIRE beam.

We present first results of a study of the submillimetre (rest frame far-infrared) properties of z~3 Lyman Break Galaxies (LBGs) and their lower-redshift counterparts BX/BM galaxies, based on Herschel-SPIRE observations of the Northern field of the Great Observatories Origins Deep Survey (GOODS-N). We use stacking analysis to determine the properties of LBGs well below the current limit of the survey. Although LBGs are not detected individually, stacking the infrared luminous LBGs (those detected with Spitzer at 24 microns yields a statistically significant submm detection with mean flux <S_{250}>= 5.9+/-1.4 mJy confirming the power of SPIRE in detecting UV-selected high-redshift galaxies at submillimetre wavelengths. In comparison, the Spitzer 24 microns detected BX/BM galaxies appear fainter with a stacked value of <S_{250}> = 2.7 +/-0.8 mJy. By fitting the Spectral Energy Distributions (SEDs) we derive median infrared luminosities, L_{IR}, of 2.8x10^{12} Lsun and 1.5x10^{11} Lsun for z~3 LBGs and BX/BMs, respectively. We find that $L_{IR} estimates derived from present measurements are in good agreement with those based on UV data for z~2 BX/BM galaxies, unlike the case for z~3 infrared luminous LBGs where the UV underestimates the true $L_{IR}. Although sample selection effects may influence this result we suggest that differences in physical properties (such as morphologies, dust distribution and extent of star-forming regions) between z ~3 LBGs and z~2 BX/BMs may also play a significant role.

We present the cross-identification and source photometry techniques used to process Herschel SPIRE imaging taken as part of the Herschel Multi-Tiered Extragalactic Survey (HerMES). Cross-identifications are performed in map-space so as to minimise source blending effects. We make use of a combination of linear inversion and model selection techniques to produce reliable cross-identification catalogues based on Spitzer MIPS 24 micron source positions. Testing on simulations and real Herschel observations show that this approach gives robust results for even the faintest sources S250~10 mJy. We apply our new technique to HerMES SPIRE observations taken as part of the science demostration phase of Herschel. For our real SPIRE observations we show that, for bright unconfused sources, our flux density estimates are in good agreement with those produced via more traditional point source detection methods (SussExtractor; Savage & Oliver et al. 2006) by Smith et al. 2010. When compared to the measured number density of sources in the SPIRE bands, we show that our method allows the recovery of a larger fraction of faint sources than these traditional methods. However this completeness is heavily dependant on the relative depth of the existing 24 micron catalogues and SPIRE imaging. Using our deepest multi-wavelength dataset in GOODS-N, we estimate that the use of shallow 24 micron in our other fields introduces an incompleteness at faint levels of between 20-40 per cent at 250 micron.

The Spectral and Photometric Imaging Receiver (SPIRE) on Herschel has been carrying out deep extragalactic surveys, one of whose aims is to establish spectral energy distributions (SED)s of individual galaxies spanning the infrared/submillimeter (IR/SMM) wavelength region. We report observations of the (IR/SMM) emission from the Lockman North field (LN) and Great Observatories Origins Deep Survey field North (GOODS-N). Because galaxy images in the wavelength range covered by Herschel generally represent a blend with contributions from neighboring galaxies, we present sets of galaxies in each field especially free of blending at 250, 350, and 500 microns. We identify the cumulative emission of these galaxies and the fraction of the far infrared cosmic background radiation they contribute. Our surveys reveal a number of highly luminous galaxies at redshift z ~< 3 and a novel relationship between infrared and visible emission that shows a dependence on luminosity and redshift.

We study the dust properties of galaxies in the redshift range 0.1<z<2.8 observed by the Herschel Space Observatory in the field of the Great Observatories Origins Deep Survey-North as part of PEP and HerMES key programmes. Infrared (IR) luminosity (L_IR) and dust temperature (T_dust) of galaxies are derived from the spectral energy distribution (SED) fit of the far-infrared (FIR) flux densities obtained with PACS and SPIRE instruments onboard Herschel. As a reference sample, we also obtain IR luminosities and dust temperatures of local galaxies at z<0.1 using AKARI and IRAS data in the field of the Sloan Digital Sky Survey. We compare the L_IR-T_dust relation between the two samples and find that: the median T_dust of Herschel-selected galaxies at z>0.5 with L_IR>5x10^{10} L_\odot, appears to be 2-5 K colder than that of AKARI-selected local galaxies with similar luminosities; and the dispersion in T_dust for high-z galaxies increases with L_IR due to the existence of cold galaxies that are not seen among local galaxies. We show that this large dispersion of the L_IR-T_dust relation can bridge the gap between local star-forming galaxies and high-z submillimeter galaxies (SMGs). We also find that three SMGs with very low T_dust (<20 K) covered in this study have close neighbouring sources with similar 24-\mum brightness, which could lead to an overestimation of FIR/(sub)millimeter fluxes of the SMGs.

We present a method of cross-calibrating the polarization angle of a polarimeter using BICEP Galactic observations. \bicep\ was a ground based experiment using an array of 49 pairs of polarization sensitive bolometers observing from the geographic South Pole at 100 and 150 GHz. The BICEP polarimeter is calibrated to +/-0.01 in cross-polarization and less than +/-0.7 degrees in absolute polarization orientation. BICEP observed the temperature and polarization of the Galactic plane (R.A= 100 degrees ~ 270 degrees and Dec. = -67 degrees ~ -48 degrees). We show that the statistical error in the 100 GHz BICEP Galaxy map can constrain the polarization angle offset of WMAP Wband to 0.6 degrees +\- 1.4 degrees. The expected 1 sigma errors on the polarization angle cross-calibration for Planck or EPIC are 1.3 degrees and 0.3 degrees at 100 and 150 GHz, respectively. We also discuss the expected improvement of the BICEP Galactic field observations with forthcoming BICEP2 and Keck observations.

The reliability of infrared (IR) and ultraviolet (UV) emissions to measure star formation rates in galaxies is investigated for a large sample of galaxies observed with the SPIRE and PACS instruments on Herschel as part of the HerMES project. We build flux-limited 250 micron samples of sources at redshift z<1, cross-matched with the Spitzer/MIPS and GALEX catalogues. About 60 % of the Herschel sources are detected in UV. The total IR luminosities, L_IR, of the sources are estimated using a SED-fitting code that fits to fluxes between 24 and 500 micron. Dust attenuation is discussed on the basis of commonly-used diagnostics: the L_IR/L_UV ratio and the slope, beta, of the UV continuum. A mean dust attenuation A_UV of ~ 3 mag is measured in the samples. L_IR/L_UV is found to correlate with L_IR. Galaxies with L_IR > 10 ^{11} L_sun and 0.5< z<1 exhibit a mean dust attenuation A_UV about 0.7 mag lower than that found for their local counterparts, although with a large dispersion. Our galaxy samples span a large range of beta and L_IR/L_UV values which, for the most part, are distributed between the ranges defined by the relations found locally for starburst and normal star-forming galaxies. As a consequence the recipe commonly applied to local starbursts is found to overestimate the dust attenuation correction in our galaxy sample by a factor ~2-3 .

We present SHARC-2 350 micron data on 20 luminous z ~ 2 starbursts with S(1.2mm) > 2 mJy from the Spitzer-selected samples of Lonsdale et al. and Fiolet et al. All the sources were detected, with S(350um) > 25 mJy for 18 of them. With the data, we determine precise dust temperatures and luminosities for these galaxies using both single-temperature fits and models with power-law mass--temperature distributions. We derive appropriate formulae to use when optical depths are non-negligible. Our models provide an excellent fit to the 6um--2mm measurements of local starbursts. We find characteristic single-component temperatures T1 ~ 35.5+-2.2 K and integrated infrared (IR) luminosities around 10^(12.9+-0.1) Lsun for the SWIRE-selected sources. Molecular gas masses are estimated at 4 x 10^(10) Msun, assuming kappa(850um)=0.15 m^2/kg and a submillimeter-selected galaxy (SMG)-like gas-to-dust mass ratio. The best-fit models imply >~2 kpc emission scales. We also note a tight correlation between rest-frame 1.4 GHz radio and IR luminosities confirming star formation as the predominant power source. The far-IR properties of our sample are indistinguishable from the purely submillimeter-selected populations from current surveys. We therefore conclude that our original selection criteria, based on mid-IR colors and 24 um flux densities, provides an effective means for the study of SMGs at z ~ 1.5--2.5.

The Spectral and Photometric Imaging Receiver (SPIRE), is the Herschel Space Observatory`s submillimetre camera and spectrometer. It contains a three-band imaging photometer operating at 250, 350 and 500 microns, and an imaging Fourier Transform Spectrometer (FTS) which covers simultaneously its whole operating range of 194-671 microns (447-1550 GHz). The SPIRE detectors are arrays of feedhorn-coupled bolometers cooled to 0.3 K. The photometer has a field of view of 4' x 8', observed simultaneously in the three spectral bands. Its main operating mode is scan-mapping, whereby the field of view is scanned across the sky to achieve full spatial sampling and to cover large areas if desired. The spectrometer has an approximately circular field of view with a diameter of 2.6'. The spectral resolution can be adjusted between 1.2 and 25 GHz by changing the stroke length of the FTS scan mirror. Its main operating mode involves a fixed telescope pointing with multiple scans of the FTS mirror to acquire spectral data. For extended source measurements, multiple position offsets are implemented by means of an internal beam steering mirror to achieve the desired spatial sampling and by rastering of the telescope pointing to map areas larger than the field of view. The SPIRE instrument consists of a cold focal plane unit located inside the Herschel cryostat and warm electronics units, located on the spacecraft Service Module, for instrument control and data handling. Science data are transmitted to Earth with no on-board data compression, and processed by automatic pipelines to produce calibrated science products. The in-flight performance of the instrument matches or exceeds predictions based on pre-launch testing and modelling: the photometer sensitivity is comparable to or slightly better than estimated pre-launch, and the spectrometer sensitivity is also better by a factor of 1.5-2.

We present preliminary results about the detection of high redshift (U)LIRGs in the Bullet cluster field by the PACS and SPIRE instruments within the Herschel Lensing Survey (HLS) Program. We describe in detail a photometric procedure designed to recover robust fluxes and deblend faint Herschel sources near the confusion noise. The method is based on the use of the positions of Spitzer/MIPS 24 um sources as priors. Our catalogs are able to reliably (5 sigma) recover galaxies with fluxes above 6 and 10 mJy in the PACS 100 and 160 um channels, respectively, and 12 to 18 mJy in the SPIRE bands. We also obtain spectral energy distributions covering the optical through the far-infrared/millimeter spectral ranges of all the Herschel detected sources, and analyze them to obtain independent estimations of the photometric redshift based on either stellar population or dust emission models. We exemplify the potential of the combined use of Spitzer position priors plus independent optical and IR photometric redshifts to robustly assign optical/NIR counterparts to the sources detected by Herschel and other (sub-)mm instruments.

The Sunyaev-Zel'dovich (SZ) effect is a spectral distortion of the cosmic microwave background as observed through the hot plasma in galaxy clusters. This distortion is a decrement in the CMB intensity for lambda > 1.3 mm, an increment at shorter wavelengths, and small again by lambda ~250 um. As part of the Herschel Lensing Survey (HLS) we have mapped 1E0657-56 (the Bullet cluster) with SPIRE with bands centered at 250, 350 and 500 um and have detected the SZ effect at the two longest wavelengths. The measured SZ effect increment central intensities are Delta I_{0} = 0.097 +- 0.019 MJy sr^{-1} at 350 um and Delta I_{0} = 0.268 +- 0.031 MJy sr^{-1} at 500 um, consistent with the SZ effect spectrum derived from previous measurements at 2 mm. No other diffuse emission is detected. The presence of the finite temperature SZ effect correction is preferred by the SPIRE data at a significance of 2.1 sigma, opening the possibility that the relativistic SZ effect correction can be constrained by SPIRE in a sample of clusters. The results presented here have important ramifications for both sub-mm measurements of galaxy clusters and blank field surveys with SPIRE.

We use deep, five band (100-500um) data from the Herschel Lensing Survey (HLS) to fully constrain the obscured star formation rate, SFR_FIR, of galaxies in the Bullet cluster (z=0.296), and a smaller background system (z=0.35) in the same field. Herschel detects 23 Bullet cluster members with a total SFR_FIR = 144 +/- 14 M_sun yr^-1. On average, the background system contains brighter far-infrared (FIR) galaxies, with ~50% higher SFR_FIR (21 galaxies; 207 +/- 9 M_sun yr^-1). SFRs extrapolated from 24um flux via recent templates (SFR_24) agree well with SFR_FIR for ~60% of the cluster galaxies. In the remaining ~40%, SFR_24 underestimates SFR_FIR due to a significant excess in observed S_100/S_24 (rest frame S_75/S_18) compared to templates of the same FIR luminosity.

The Herschel Lensing Survey (HLS) takes advantage of gravitational lensing by massive galaxy clusters to sample a population of high-redshift galaxies which are too faint to be detected above the confusion limit of current far-infrared/submillimeter telescopes. Measurements from 100-500 micron bracket the peaks of the far-infrared spectral energy distributions of these galaxies, characterizing their infrared luminosities and star formation rates. We introduce initial results from our science demonstration phase observations, directed toward the Bullet cluster (1E0657-56). By combining our observations with LABOCA 870 micron and AzTEC 1.1 mm data we fully constrain the spectral energy distributions of 19 MIPS 24 micron selected galaxies which are located behind the cluster. We find that their colors are best fit using templates based on local galaxies with systematically lower infrared luminosities.This suggests that our sources are not like local ultra-luminous infrared galaxies in which vigorous star formation is contained in a compact highly dust-obscured region. Instead, they appear to be scaled up versions of lower luminosity local galaxies with star formation occurring on larger physical scales.

The Herschel Lensing Survey (HLS) will conduct deep PACS and SPIRE imaging of ~40 massive clusters of galaxies. The strong gravitational lensing power of these clusters will enable us to penetrate through the confusion noise, which sets the ultimate limit on our ability to probe the Universe with Herschel. Here, we present an overview of our survey and a summary of the major results from our Science Demonstration Phase (SDP) observations of the Bullet Cluster (z=0.297). The SDP data are rich, allowing us to study not only the background high-redshift galaxies (e.g., strongly lensed and distorted galaxies at z=2.8 and 3.2) but also the properties of cluster-member galaxies. Our preliminary analysis shows a great diversity of far-infrared/submillimeter spectral energy distributions (SEDs), indicating that we have much to learn with Herschel about the properties of galaxy SEDs. We have also detected the Sunyaev-Zel'dovich (SZ) effect increment with the SPIRE data. The success of this SDP program demonstrates the great potential of the Herschel Lensing Survey to produce exciting results in a variety of science areas.

We measure the angular correlation function, w(theta), from 0.5 to 30 arcminutes of detected sources in two wide fields of the Herschel Multi-tiered Extragalactic Survey (HerMES). Our measurements are consistent with the expected clustering shape from a population of sources that trace the dark matter density field, including non-linear clustering at arcminute angular scales arising from multiple sources that occupy the same dark matter halos. By making use of the halo model to connect the spatial clustering of sources to the dark matter halo distribution, we estimate source bias and halo occupation number for dusty sub-mm galaxies at z ~ 2. We find that sub-mm galaxies with 250 micron flux densities above 30 mJy reside in dark matter halos with mass above (5\pm4) x 10^12 M_sun, while (14\pm8)% of such sources appear as satellites in more massive halos.

The Herschel Space Observatory enables us to accurately measure the bolometric output of starburst galaxies and active galactic nuclei (AGN) by directly sampling the peak of their far-infrared (IR) emission. Here we examine whether the spectral energy distribution (SED) and dust temperature of galaxies have strongly evolved since z~2.5. We use Herschel deep extragalactic surveys from 100 to 500um to compute total IR luminosities in galaxies down to the faintest levels, using PACS and SPIRE in the GOODS-North field (PEP and HerMES key programs). We show that measurements in the SPIRE bands can be used below the statistical confusion limit if information at higher spatial resolution is used to identify isolated galaxies whose flux is not boosted by bright neighbors. Below z~1.5, mid-IR extrapolations are correct for star-forming galaxies with a dispersion of only 40% (0.15dex), therefore similar to z~0 galaxies. This narrow distribution is puzzling when considering the range of physical processes that could have affected the SED of these galaxies. Extrapolations from only one of the 160um, 250um or 350um bands alone tend to overestimate the total IR luminosity. This may be explained by the lack of far-IR constraints around and above ~150um (rest-frame) on local templates. We also note that the dust temperature of luminous IR galaxies around z~1 is mildly colder by 10-15% than their local analogs and up to 20% for ULIRGs at z~1.6. Above z=1.5, distant galaxies are found to exhibit a substantially larger mid- over far-IR ratio, which could either result from stronger broad emission lines or warm dust continuum heated by a hidden AGN. Two thirds of the AGNs identified in the field with a measured redshift exhibit the same behavior as purely star-forming galaxies. Hence a large fraction of AGNs harbor star formation at very high SFR and in conditions similar to purely star-forming galaxies.

Local luminosity functions are fundamental benchmarks for high-redshift galaxy formation and evolution studies as well as for models describing these processes. Determining the local luminosity function in the submillimeter range can help to better constrain in particular the bolometric luminosity density in the local Universe, and Herschel offers the first opportunity to do so in an unbiased way by imaging large sky areas at several submillimeter wavelengths. We present the first Herschel measurement of the submillimeter 0<z<0.2 local luminosity function and infrared bolometric (8-1000 $\mu$m) local luminosity density based on SPIRE data from the HerMES Herschel Key Program over 14.7 deg^2. Flux measurements in the three SPIRE channels at 250, 350 and 500 \mum are combined with Spitzer photometry and archival data. We fit the observed optical-to-submillimeter spectral energy distribution of SPIRE sources and use the 1/V_{max} estimator to provide the first constraints on the monochromatic 250, 350 and 500 \mum as well as on the infrared bolometric (8-1000 \mum) local luminosity function based on Herschel data. We compare our results with modeling predictions and find a slightly more abundant local submillimeter population than predicted by a number of models. Our measurement of the infrared bolometric (8-1000 \mum) local luminosity function suggests a flat slope at low luminosity, and the inferred local luminosity density, 1.31_-0.21^+0.24 x 10^8 Lsun Mpc^-3, is consistent with the range of values reported in recent literature.

We present colours of sources detected with the Herschel/SPIRE instrument in deep extragalactic surveys of the Lockman Hole, Spitzer-FLS, and GOODS-N fields in three photometric bands at 250, 350 and 500 micrometers. We compare these with expectations from the literature and discuss associated uncertainties and biases in the SPIRE data. We identify a 500 micrometer flux limited selection of sources from the HerMES point source catalogue that appears free from neighbouring/blended sources in all three SPIRE bands. We compare the colours with redshift tracks of various contemporary models. Based on these spectral templates we show that regions corresponding to specific population types and redshifts can be identified better in colour-flux space. The redshift tracks as well as the colour-flux plots imply a majority of detected objects with redshifts at 1<z<3.5, somewhat depending on the group of model SEDs used. We also find that a population of S_250/S_350<0.8 at fluxes above 50 mJy as observed by SPIRE is not well represented by contemporary models and could consist of a mix of cold and lensed galaxies.

We have carried out two extremely deep surveys with SPIRE, one of the two cameras on Herschel, at 250 microns, close to the peak of the far-infrared background. We have used the results to investigate the evolution of the rest-frame 250-micron luminosity function out to z=2. We find evidence for strong evolution out to a redshift of around 1 but evidence for at most weak evolution beyond this redshift. Our results suggest that a significant part of the stars and metals in the Universe today were formed at z<1.4 in spiral galaxies.

We report on the sensitivity of SPIRE photometers on the Herschel Space Observatory. Specifically, we measure the confusion noise from observations taken during the Science Demonstration Phase of the Herschel Multi-tiered Extragalactic Survey. Confusion noise is defined to be the spatial variation of the sky intensity in the limit of infinite integration time, and is found to be consistent among the different fields in our survey at the level of 5.8, 6.3 and 6.8 mJy/beam at 250, 350 and 500 microns, respectively. These results, together with the measured instrument noise, may be used to estimate the integration time required for confusion-limited maps, and provide a noise estimate for maps obtained by SPIRE.

Nuclear and starburst activity are known to often occur concomitantly. Herschel-SPIRE provides sampling of the FIR SEDs of type 1 and type 2 AGN, allowing for the separation between the hot dust (torus) and cold dust (starburst) emission. We study large samples of spectroscopically confirmed type 1 and type 2 AGN lying within the Herschel Multi-tiered Extragalactic Survey (HerMES) fields observed during the science demonstration phase, aiming to understand their FIR colour distributions and constrain their starburst contributions. We find that one third of the spectroscopically confirmed AGN in the HerMES fields have 5-sigma detections at 250um, in agreement with previous (sub)mm AGN studies. Their combined Spitzer-MIPS and Herschel-SPIRE colours - specifically S(250)/S(70) vs. S(70)/S(24) - quite clearly separate them from the non-AGN, star-forming galaxy population, as their 24-um flux is dominated by the hot torus emission. However, their SPIRE colours alone do not differ from those of non-AGN galaxies. SED fitting shows that all those AGN need a starburst component to fully account for their FIR emission. For objects at z > 2, we find a correlation between the infrared luminosity attributed to the starburst component, L(SB), and the AGN accretion luminosity, L(acc), with L(SB) propto L(acc)^0.35. Type 2 AGN detected at 250um show on average higher L(SB) than type 1 objects but their number is still too low to establish whether this trend indicates stronger star-formation activity.

Emission at far-infrared wavelengths makes up a significant fraction of the total light detected from galaxies over the age of Universe. Herschel provides an opportunity for studying galaxies at the peak wavelength of their emission. Our aim is to provide a benchmark for models of galaxy population evolution and to test pre-existing models of galaxies. With the Herschel Multi-tiered Extra-galactic survey, HerMES, we have observed a number of fields of different areas and sensitivity using the SPIRE instrument on Herschel. We have determined the number counts of galaxies down to ~20 mJy. Our constraints from directly counting galaxies are consistent with, though more precise than, estimates from the BLAST fluctuation analysis. We have found a steep rise in the Euclidean normalised counts at <100 mJy. We have directly resolved 15% of the infrared extra-galactic background at the wavelength near where it peaks.

We set out to determine the ratio, q(IR), of rest-frame 8-1000um flux, S(IR), to monochromatic radio flux, S(1.4GHz), for galaxies selected at far-IR and radio wavelengths, to search for signs that the ratio evolves with redshift, luminosity or dust temperature, and to identify any far-IR-bright outliers - useful laboratories for exploring why the far-IR/radio correlation is generally so tight when the prevailing theory suggests variations are almost inevitable. We use flux-limited 250-um and 1.4-GHz samples, obtained in GOODS-N using Herschel (HerMES; PEP) and the VLA. We determine bolometric IR output using ten bands spanning 24-1250um, exploiting data from PACS and SPIRE, as well as Spitzer, SCUBA, AzTEC and MAMBO. We also explore the properties of an L(IR)-matched sample, designed to reveal evolution of q(IR) with z, spanning log L(IR) = 11-12 L(sun) and z=0-2, by stacking into the radio and far-IR images. For 1.4-GHz-selected galaxies, we see tentative evidence of a break in the flux ratio, q(IR), at L(1.4GHz) ~ 10^22.7 W/Hz, where AGN are starting to dominate the radio power density, and of weaker correlations with z and T(d). From our 250-um-selected sample we identify a small number of far-IR-bright outliers, and see trends of q(IR) with L(1.4GHz), L(IR), T(d) and z, noting that some of these are inter-related. For our L(IR)-matched sample, there is no evidence that q(IR) changes significantly as we move back into the epoch of galaxy formation: we find q(IR) goes as (1+z)^gamma, where gamma = -0.04 +/- 0.03 at z=0-2; however, discounting the least reliable data at z<0.5 we find gamma = -0.26 +/- 0.07, modest evolution which may be related to the radio background seen by ARCADE2, perhaps driven by <10uJy radio activity amongst ordinary star-forming galaxies at z>1.

Background Imaging of Cosmic Extragalactic Polarization (BICEP) is a bolometric polarimeter designed to measure the inflationary B-mode polarization of the cosmic microwave background (CMB) at degree angular scales. During three seasons of observing at the South Pole (2006 through 2008), BICEP mapped ~2% of the sky chosen to be uniquely clean of polarized foreground emission. Here we present initial results derived from a subset of the data acquired during the first two years. We present maps of temperature, Stokes Q and U, E and B modes, and associated angular power spectra. We demonstrate that the polarization data are self-consistent by performing a series of jackknife tests. We study potential systematic errors in detail and show that they are sub-dominant to the statistical errors. We measure the E-mode angular power spectrum with high precision at 21 < ell < 335, detecting for the first time the peak expected at ell ~ 140. The measured E-mode spectrum is consistent with expectations from a LCDM model, and the B-mode spectrum is consistent with zero. The tensor-to-scalar ratio derived from the B-mode spectrum is r = 0.03+0.31-0.26, or r < 0.72 at 95% confidence, the first meaningful constraint on the inflationary gravitational wave background to come directly from CMB B-mode polarization.

The BICEP experiment was designed specifically to search for the signature of inflationary gravitational waves in the polarization of the cosmic microwave background (CMB). Using a novel small-aperture refractor and 49 pairs of polarization-sensitive bolometers, BICEP has completed 3 years of successful observations at the South Pole beginning in 2006 February. To constrain the amplitude of the inflationary B-mode polarization, which is expected to be at least 7 orders of magnitude fainter than the 3 K CMB intensity, precise control of systematic effects is essential. This paper describes the characterization of potential systematic errors for the BICEP experiment, supplementing a companion paper on the initial cosmological results. Using the analysis pipelines for the experiment, we have simulated the impact of systematic errors on the B-mode polarization measurement. Guided by these simulations, we have established benchmarks for the characterization of critical instrumental properties including bolometer relative gains, beam mismatch, polarization orientation, telescope pointing, sidelobes, thermal stability, and timestream noise model. A comparison of the benchmarks with the measured values shows that we have characterized the instrument adequately to ensure that systematic errors do not limit BICEP's 2-year results, and identifies which future refinements are likely necessary to probe inflationary B-mode polarization down to levels below a tensor-to-scalar ratio r = 0.1.

We present a summary of data obtained with the 350 micron polarimeter, Hertz, at the Caltech Submillimeter Observatory. We give tabulated results and maps showing polarization vectors and flux contours. The summary includes over 4300 individual measurements in 56 Galactic sources and 2 galaxies. Of these measurements, 2153 have P >= 3\sigma_p statistical significance. The median polarization of the entire data set is 1.46%.

We report the detection of variable emission from Sgr A* in almost all wavelength bands (i.e. centimeter, millimeter, submillimeter, near-IR and X-rays) during a multi-wavelength observing campaign. Three new moderate flares are detected simultaneously in both near-IR and X-ray bands. The ratio of X-ray to near-IR flux in the flares is consistent with inverse Compton scattering of near-IR photons by submillimeter emitting relativistic particles which follow scaling relations obtained from size measurements of Sgr A*. We also find that the flare statistics in near-IR wavelengths is consistent with the probability of flare emission being inversely proportional to the flux. At millimeter wavelengths, the presence of flare emission at 43 GHz (7mm) using VLBA with milli-arcsecond spatial resolution indicates the first direct evidence that hourly time scale flares are localized within the inner 30$\times$70 Schwarzschild radii of Sgr A*. We also show several cross correlation plots between near-IR, millimeter and submillimeter light curves that collectively demonstrate the presence of time delays between the peaks of emission up to three hours. The evidence for time delays at millimeter and submillimeter wavelengths are consistent with the source of emission being optically thick initially followed by a transition to an optically thin regime. In particular, there is an intriguing correlation between the optically thin near-IR and X-ray flare and optically thick radio flare at 43 GHz that occurred on 2007 April 4. This would be the first evidence of a radio flare emission at 43 GHz delayed with respect to the near-IR and X-ray flare emission.

One of the key problems in star formation research is to determine the role of magnetic fields. Starting from the atomic inter-cloud medium (ICM) which has density nH ~ 1 per cubic cm, gas must accumulate from a volume several hundred pc across in order to form a typical molecular cloud. Star formation usually occurs in cloud cores, which have linear sizes below 1 pc and densities nH2 > 10^5 per cubic cm. With current technologies, it is hard to probe magnetic fields at scales lying between the accumulation length and the size of cloud cores, a range corresponds to many levels of turbulent eddy cascade, and many orders of magnitude of density amplification. For field directions detected from the two extremes, however, we show here that a significant correlation is found. Comparing this result with molecular cloud simulations, only the sub-Alfvenic cases result in field orientations consistent with our observations.

We present a novel statistical analysis aimed at deriving the intrinsic shapes and magnetic field orientations of molecular clouds using dust emission and polarization observations by the Hertz polarimeter. Our observables are the aspect ratio of the projected plane-of-the-sky cloud image, and the angle between the mean direction of the plane-of-the-sky component of the magnetic field and the short axis of the cloud image. To overcome projection effects due to the unknown orientation of the line-of-sight, we combine observations from 24 clouds, assuming that line-of-sight orientations are random and all are equally probable. Through a weighted least-squares analysis, we find that the best-fit intrinsic cloud shape describing our sample is an oblate disk with only small degrees of triaxiality. The best-fit intrinsic magnetic field orientation is close to the direction of the shortest cloud axis, with small (~24 deg) deviations toward the long/middle cloud axes. However, due to the small number of observed clouds, the power of our analysis to reject alternative configurations is limited.

We present single-dish 350 micron dust continuum polarimetry as well as HCN and HCO+ J=4-3 rotational emission spectra obtained on NGC 1333 IRAS 4. The polarimetry indicates a uniform field morphology over a 20" radius from the peak continuum flux of IRAS 4A, in agreement with models of magnetically supported cloud collapse. The field morphology around IRAS 4B appears to be quite distinct however, with indications of depolarization observed towards the peak flux of this source. Inverse P-Cygni profiles are observed in the HCN J=4-3 line spectra towards IRAS 4A, providing a clear indication of infall gas motions. Taken together, the evidence gathered here appears to support the scenario that IRAS 4A is a cloud core in a critical state of support against gravitational collapse.

Context: Among late-type red giants, an interesting change occurs in the structure of the outer atmospheric layers as one moves to later spectral types in the Hertzsprung-Russell diagram: a chromosphere is always present, but the coronal emission diminishes and a cool massive wind steps in. Aims: Where most studies have focussed on short-wavelength observations, this article explores the influence of the chromosphere and the wind on long-wavelength photometric measurements. Methods: The observational spectral energy distributions are compared with the theoretical predictions of the MARCS atmosphere models for a sample of 9 K- and M-giants. The discrepancies found are explained using basic models for flux emission originating from a chromosphere or an ionized wind. Results: For 7 out of 9 sample stars, a clear flux excess is detected at (sub)millimeter and/or centimeter wavelengths. The precise start of the excess depends upon the star under consideration. The flux at wavelengths shorter than about 1 mm is most likely dominated by an optically thick chromosphere, where an optically thick ionized wind is the main flux contributor at longer wavelengths. Conclusions: Although the optical to mid-infrared spectrum of the studied K- and M-giants is well represented by a radiative equilibrium atmospheric model, the presence of a chromosphere and/or ionized stellar wind at higher altitudes dominates the spectrum in the (sub)millimeter and centimeter wavelength ranges. The presence of a flux excess also has implications on the role of these stars as fiducial spectrophotometric calibrators in the (sub)millimeter and centimeter wavelength range.

When and how did galaxies form and their metals accumulate? Over the last decade, this has moved from an archeological question to a live investigation: there is now a broad picture of the evolution of galaxies in dark matter halos: their masses, stars, metals and supermassive blackholes. Galaxies have been found and studied in which these formation processes are taking place most vigorously, all the way back in cosmic time to when the intergalactic medium (IGM) was still largely neutral. However, the details of how and why the interstellar medium (ISM) in distant galaxies cools, is processed, recycled and enriched in metals by stars, and fuels active galactic nuclei (AGNs) remain uncertain. In particular, the cooling of gas to fuel star formation, and the chemistry and physics of the most intensely active regions is hidden from view at optical wavelengths, but can be seen and diagnosed at mid- & far-infrared (IR) wavelengths. Rest-frame IR observations are important first to identify the most luminous, interesting and important galaxies, secondly to quantify accurately their total luminosity, and finally to use spectroscopy to trace the conditions in the molecular and atomic gas out of which stars form. In order to map out these processes over the full range of environments and large-scale structures found in the universe - from the densest clusters of galaxies to the emptiest voids - we require tools for deep, large area surveys, of millions of galaxies out to z~5, and for detailed follow-up spectroscopy. The necessary tools can be realized technically. Here, we outline the requirements for gathering the crucial information to build, validate and challenge models of galaxy evolution.

Spitzer and Caltech Submillimeter Observatory (CSO) images and spectrophotometry of epsilon Eridani at wavelengths from 3.5 to 350 um reveal new details of its bright debris disk. The 350 um map confirms the presence of a ring at r = 11-28 arcsec (35-90 AU) observed previously at longer sub-mm wavelengths. The Spitzer mid- and far-IR images do not show the ring, but rather a featureless disk extending from within a few arcsec of the star across the ring to r ~ 34 arcsec (110 AU). The spectral energy distribution (SED) of the debris system implies a complex structure. A model constrained by the surface brightness profiles and the SED indicates that the sub-mm ring emission is primarily from large (a ~ 135 um) grains, with smaller (a ~ 15 um) grains also present in and beyond the ring. The Spitzer IRS and MIPS SED-mode spectrophotometry data clearly show the presence of spatially compact excess emission at lambda > 15 um that requires the presence of two additional narrow belts of dust within the sub-mm ring's central void. The innermost belt at r ~ 3 AU is composed of silicate dust. A simple dynamical model suggests that dust produced collisionally by a population of about 11 M_Earth of planetesimals in the sub-mm ring could be the source of the emission from both in and beyond the sub-mm ring. Maintaining the inner belts and the inner edge to the sub-mm ring may require the presence of three planets in this system including the candidate radial velocity object.

BICEP is a ground-based millimeter-wave bolometric array designed to target the primordial gravity wave signature on the polarization of the cosmic microwave background (CMB) at degree angular scales. Currently in its third year of operation at the South Pole, BICEP is measuring the CMB polarization with unprecedented sensitivity at 100 and 150 GHz in the cleanest available 2% of the sky, as well as deriving independent constraints on the diffuse polarized foregrounds with select observations on and off the Galactic plane. Instrument calibrations are discussed in the context of rigorous control of systematic errors, and the performance during the first two years of the experiment is reviewed.

Sgr A*, the massive black hole at the center of the Galaxy, varies in radio through X-ray emission on hourly time scales. The flare activity is thought to arise from the innermost region of an accretion flow onto Sgr A*. We present simultaneous light curves of Sgr A* in radio, sub-mm and X-rays that show a possible time delay of 110$\pm17$ minutes between X-ray and 850 $\mu$m suggesting that the sub-mm flare emission is optically thick. At radio wavelengths, we detect time lags of of $20.4\pm6.8, 30\pm12$ and 20$\pm6$ minutes between the flare peaks observed at 13 and 7 mm (22 and 43 GHz) in three different epochs using the VLA. Linear polarization of 1$\pm0.2$% and 0.7$\pm0.1$% is detected at 7 and 13 mm, respectively, when averaged over the entire observation on 2006 July 17. A simple model of a bubble of synchrotron emitting electrons cooling via adiabatic expansion can explain the time delay between various wavelengths, the asymmetric shape of the light curves, and the observed polarization of the flare emission at 43 and 22 GHz. The derived physical quantities that characterize the emission give an expansion speed of v$_{exp} \sim 0.003-0.1$c, magnetic field of B$\sim$10-70 Gauss and particle spectral index p$\sim$1-2. These parameters suggest that the associated plasma cannot escape from Sgr A* unless it has a large bulk motion.

Energetic flares are observed in the Galactic supermassive black hole Sagittarius A* from radio to X-ray wavelengths. On a few occasions, simultaneous flares have been detected in IR and X-ray observations, but clear counterparts at longer wavelengths have not been seen. We present a flare observed over several hours on 2006 July 17 with the Chandra X-Ray Observatory, the Keck II telescope, the Caltech Submillimeter Observatory, and the Submillimeter Array. All telescopes observed strong flare events, but the submillimeter peak is found to occur nearly 100 minutes after the X-ray peak. Submillimeter polarization data show linear polarization in the excess flare emission, increasing from 9% to 17% as the flare passes through its peak, consistent with a transition from optically thick to thin synchrotron emission. The temporal and spectral behavior of the flare require that the energetic electrons responsible for the emission cool faster than expected from their radiative output. This is consistent with adiabatic cooling in an expanding emission region, with X-rays produced through self-Compton scattering, although not consistent with the simplest model of such expansion. We also present a submillimeter flare that followed a bright IR flare on 2005 July 31. Compared to 2006, this event had a larger peak IR flux and similar submillimeter flux, but it lacked measurable X-ray emission. It also showed a shorter delay between the IR and submillimeter peaks. Based on these events we propose a synchrotron and self-Compton model to relate the submillimeter lag and the variable IR/X-ray luminosity ratio.

We have used the SHARP polarimeter at the Caltech Submillimeter Observatory to map the polarization at wavelengths of 350 and 450 micron in a ~2 x 3 arcmin region of the Orion Molecular Cloud. The map covers the brightest region of the OMC-1 ridge including the Kleinmann-Low (KL) nebula and the submillimeter source Orion-south. The ratio of 450-to-350 micron polarization is ~ 1.3 +/- 0.3 in the outer parts of the cloud and drops by a factor of 2 towards KL. The outer cloud ratio is consistent with measurements in other clouds at similar wavelengths and confirms previous measurements placing the minimum of the polarization ratio in dusty molecular clouds at a wavelength ~ 350 micron.

Continuum observations at 350um are presented of seven nearby elliptical galaxies, for which CO-gas disks have recently been resolved with interferometry mapping. These SHARCII mapping results provide the first clearly resolved far-infrared(FIR) to submillimeter(submm) continuum emission from cold dust (with temperatures 32K > T > 22K) of any elliptical galaxy at a distance >40Mpc. The measured FIR excess shows that the most likely and dominant heating source of this dust is not dilute stellar radiation or cooling flows, but rather star-formation, that could have been triggered by an accretion or merger event and fueled by dust-rich material that has settled in a dense region co-spatial with the central CO-gas disks. The dust is detected even in two cluster ellipticals that are deficient in HI, showing that, unlike the HI, cold dust and CO in ellipticals can survive in the presence of hot X-ray gas, even in galaxy clusters. No dust cooler than 20K, either distributed outside the CO disks, or co-spatial with and heated by the entire dilute stellar optical galaxy (or very extended HI), is currently evident.

We present observations of four z>= SDSS quasars at 350 micron with the SHARC-II bolometer camera on the Caltech Submillimeter Observatory. These are among the deepest observations that have been made by SHARC-II at 350 micron, and three quasars are detected at >=3 sigma significance, greatly increasing the sample of 350 micron (corresponds to rest frame wavelengths of <60 micron at z>=5), detected high-redshift quasars. The derived rest frame far-infrared (FIR) emission in the three detected sources is about five to ten times stronger than that expected from the average SED of the local quasars given the same 1450A luminosity. Combining the previous submillimeter and millimeter observations at longer wavelengths, the temperatures of the FIR-emitting warm dust from the three quasar detections are estimated to be in the range of 39 to 52 K. Additionally, the FIR-to-radio SEDs of the three 350 micron detections are consistent with the emission from typical star forming galaxies. The FIR luminosities are ~10^{13} L_solar and the dust masses are >= 10^{8}M_solar. These results confirm that huge amounts of warm dust can exist in the host galaxies of optically bright quasars as early as z~6. The universe is so young at these epochs (~1 Gyr) that a rapid dust formation mechanism is required. We estimate the size of the FIR dust emission region to be about a few kpc, and further provide a comparison of the SEDs among different kinds of dust emitting sources to investigate the dominant dust heating mechanism.

Multi-wavelength imaging polarimetry at far-infrared wavelengths has proven to be an excellent tool for studying the physical properties of dust, molecular clouds, and magnetic fields in the interstellar medium. Although these wavelengths are only observable from airborne or space-based platforms, no first-generation instrument for the Stratospheric Observatory for Infrared Astronomy (SOFIA) is presently designed with polarimetric capabilities. We study several options for upgrading the High-resolution Airborne Wideband Camera (HAWC) to a sensitive FIR polarimeter. HAWC is a 12 x 32 pixel bolometer camera designed to cover the 53 - 215 micron spectral range in 4 colors, all at diffraction-limited resolution (5 - 21 arcsec). Upgrade options include: (1) an external set of optics which modulates the polarization state of the incoming radiation before entering the cryostat window; (2) internal polarizing optics; and (3) a replacement of the current detector array with two state-of-the-art superconducting bolometer arrays, an upgrade of the HAWC camera as well as polarimeter. We discuss a range of science studies which will be possible with these upgrades including magnetic fields in star-forming regions and galaxies and the wavelength-dependence of polarization.

We have developed a fore-optics module that converts the SHARC-II camera at the Caltech Submillimeter Observatory into a sensitive imaging polarimeter at wavelengths of 350 and 450 microns. We refer to this module as "SHARP". SHARP splits the incident radiation into two orthogonally polarized beams that are then re-imaged onto opposite ends of the 32 x 12 pixel detector array in SHARC-II. A rotating half-wave plate is used just upstream from the polarization-splitting optics. The effect of SHARP is to convert SHARC-II into a dual-beam 12 x 12 pixel polarimeter. A novel feature of SHARP's design is the use of a crossed grid in a submillimeter polarimeter. Here we describe the detailed optical design of SHARP and present results of tests carried out during our first few observing runs. At 350 microns, the beam size (9 arcseconds), throughput (75%), and instrumental polarization (< 1%) are all very close to our design goals.

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.

We present 350micron observations of 36 ultraluminous infrared galaxies (ULIRGs) at intermediate redshifts (0.089 <= z <= 0.926) using the Submillimeter High Angular Resolution Camera II (SHARC-II) on the Caltech Submillimeter Observatory (CSO). In total, 28 sources are detected at S/N >= 3, providing the first flux measurements longward of 100micron for a statistically significant sample of ULIRGs in the redshift range of 0.1 < z < 1.0. Combining our 350micron flux measurements with the existing IRAS 60 and 100micron data, we fit a single-temperature model to the spectral energy distribution (SED), and thereby estimate dust temperatures and far-IR luminosities. Assuming an emissivity index of beta = 1.5, we find a median dust temperature and far-IR luminosity of Td = 42.8+-7.1K and log(Lfir/Lsolar) = 12.2+-0.5, respectively. The far-IR/radio correlation observed in local star-forming galaxies is found to hold for ULIRGs in the redshift range 0.1 < z < 0.5, suggesting that the dust in these sources is predominantly heated by starbursts. We compare the far-IR luminosities and dust temperatures derived for dusty galaxy samples at low and high redshifts with our sample of ULIRGs at intermediate redshift. A general Lfir-Td relation is observed, albeit with significant scatter, due to differing selection effects and variations in dust mass and grain properties. The relatively high dust temperatures observed for our sample compared to that of high-z submillimeter-selected starbursts with similar far-IR luminosities suggest that the dominant star formation in ULIRGs at moderate redshifts takes place on smaller spatial scales than at higher redshifts.

Sgr A* is considered to be a massive black hole at the Galactic center and is known to be variable in radio, millimeter, near-IR and X-rays. Recent multi-wavelength observing campaigns show a simultaneous X-ray and near-IR flare, as well as sub-millimeter and near-IR flares from Sgr A*. The flare activity is thought to be arising from the innermost region of Sgr A*. We have recently argued that the duration of flares in near-IR and submillimeter wavelengths implies that the burst of emission expands and cools on a dynamical time scale before the flares leave Sgr A*. The detection of radio flares with a time delay in the range of 20 and 40 minutes between 7 and 12mm peak emission implies adiabatic expansion of a uniform, spherical hot blob due to flare activity. We suspect that this simple outflow picture shows some of the characteristics that are known to take place in microquasars, thus we may learn much from comparative study of Sgr A* and its environment vs. microquasars.

We present 350 micron observations of 15 Chapman et al. submillimeter galaxies (SMGs) with radio counterparts and optical redshifts. We detect 12 and obtain sensitive upper limits for three, providing direct, precise measurements of their far-infrared luminosities and characteristic dust temperatures. With these, we verify the linear radio--far-infrared correlation at redshifts of z ~ 1--3 and luminosities of 10^11--10^13 L_sun, with a power-law index of 1.02+-0.12 and rms scatter of 0.12 dex. However, either the correlation constant q or the dust emissivity index beta is lower than measured locally. The best fitting q ~ 2.14 is consistent with SMGs being predominantly starbust galaxies, without significant AGN contribution, at far-infrared wavelengths. Gas-to-dust mass ratios are estimated at 54+-13 (kappa_850 / 0.15 m^2 kg^-1), depending on the absoption efficiency kappa, with intrinsic dispersion ~40% around the mean value. Dust temperatures consistent with 34.6+-3 K (beta/1.5)^-0.71, at z ~ 1.5--3.5, suggest that far-infrared photometric redshifts may be viable, and perhaps accurate to 10% < dz/(1+z), for up to 80% of the SMG population in this range, if the above temperature characterizes the full range of SMGs. However, observed temperature evolution of T_d ~ (1+z) is also plausible and could result from selection effects. From the observed luminosity-temperature (L-T) relation, L ~ T_obs^(2.82+-0.29), we derive scaling relations for dust mass versus dust temperature, and we identify expressions to inter-relate the observed quantities. These suggest that measurements at a single wavelength, in the far-infrared, submillimeter or radio wave bands, might constrain dust temperatures and far-infrared luminosities for most SMGs with redshifts at z ~ 0.5--4.

We have used the SHARC II camera at Caltech Submillimeter Observatory to make 350 and 450 micron images of the Vega dust disk at spatial resolutions (FWHM) of 9.7" and 11.1", respectively. The images show a ring-like morphology (radius \~ 100 AU) with inhomogeneous structure that is qualitatively different from that previously reported at 850 microns and longer wavelengths. We attribute the 350/450 micron emission to a grain population whose characteristic size (~ 1 mm) is intermediate between that of the cm-sized grains responsible for emission longward of 850 microns and the much smaller grains (< 18 microns) in the extensive halo, visible at 70 microns, discussed by Su et al. (2005). We have combined our submillimeter images with Spitzer data at 70 microns to produce 2-d maps of line-of-sight optical depth (relative column density). These "tau maps" suggest that the mm-sized grains are located preferentially in three symmetrically-located concentrations. If so, then this structure could be understood in terms of the Wyatt (2003) model in which planetesimals are trapped in the mean motion resonances of a Neptune-mass planet at 65 AU, provided allowance is made for the spatial distribution of dust grains to differ from that of the parent planetesimals. The peaks of the tau maps are, in fact, located near the expected positions corresponding to the 4:3 resonance. If this identification is confirmed by future observations, it would resolve an ambiguity with regard to the location of the planet.

The Robinson Telescope (BICEP) is a ground-based millimeter-wave bolometric array designed to study the polarization of the cosmic microwave background radiation (CMB) and galactic foreground emission. Such measurements probe the energy scale of the inflationary epoch, tighten constraints on cosmological parameters, and verify our current understanding of CMB physics. Robinson consists of a 250-mm aperture refractive telescope that provides an instantaneous field-of-view of 17 degrees with angular resolution of 55 and 37 arcminutes at 100 GHz and 150 GHz, respectively. Forty-nine pair of polarization-sensitive bolometers are cooled to 250 mK using a 4He/3He/3He sorption fridge system, and coupled to incoming radiation via corrugated feed horns. The all-refractive optics is cooled to 4 K to minimize polarization systematics and instrument loading. The fully steerable 3-axis mount is capable of continuous boresight rotation or azimuth scanning at speeds up to 5 deg/s. Robinson has begun its first season of observation at the South Pole. Given the measured performance of the instrument along with the excellent observing environment, Robinson will measure the E-mode polarization with high sensitivity, and probe for the B-modes to unprecedented depths. In this paper we discuss aspects of the instrument design and their scientific motivations, scanning and operational strategies, and the results of initial testing and observations.

We report detections of six high-redshift (1.8 < z < 6.4), optically luminous, radio-quiet quasars at 350 micron, using the SHARC II bolometer camera at the Caltech Submillimeter Observatory. Our observations double the number of high-redshift quasars for which 350 micron photometry is available. By combining the 350 micron measurements with observations at other submillimeter/millimeter wavelengths, for each source we have determined the temperature of the emitting dust (ranging from 40 to 60 K) and the far-infrared luminosity (0.6 to 2.2 x 10(13) Lo). The combined mean spectral energy distribution (SED) of all high-redshift quasars with two or more rest frame far-infrared photometric measurements is best fit with a greybody with temperature of 47 +- 3 K and a dust emissivity power-law spectral index of beta = 1.6 +- 0.1. This warm dust component is a good tracer of the starburst activity of the quasar host galaxy. The ratio of the far-infrared to radio luminosities of infrared luminous, radio-quiet high-redshift quasars is consistent with that found for local star-forming galaxies.

(abridged) We describe highlights of the results of two observing campaigns in 2004 to investigate the correlation of flare activity in Sgr A* in different wavelength regimes, using a total of nine ground and space-based telescopes. We report the detection of several new near-IR flares during the campaign based on {\it HST} observations. The level of near-IR flare activity can be as low as $\sim0.15$ mJy at 1.6 $\mu$m and continuous up to about 40% of the total observing time. Using the NICMOS instrument on the {\it HST}, the {\it XMM-Newton} and CSO observatories, we also detect simultaneous bright X-ray and near-IR flare in which we observe for the first time correlated substructures as well as simultaneous submillimeter and near-IR flaring. X-ray emission is arising from the population of near-IR-synchrotron-emitting relativistic particles which scatter submillimeter seed photons within the inner 10 Schwarzschild radii of Sgr A* up to X-ray energies. In addition, using the inverse Compton scattering picture, we explain the high energy 20-120 keV emission from the direction toward Sgr A*, and the lack of one-to-one X-ray counterparts to near-IR flares, by the variation of the magnetic field and the spectral index distributions of this population of nonthermal particles. In this picture, the evidence for the variability of submillimeter emission during a near-IR flare is produced by the low-energy component of the population of particles emitting synchrotron near-IR emission. Based on the measurements of the duration of flares in near-IR and submillimeter wavelengths, we argue that the cooling could be due to adiabatic expansion with the implication that flare activity may drive an outflow.

SHARC II, 350-micron continuum and archival HST J-H band maps are presented of NGC 3656, the brightest of our sample of six elliptical galaxies for which resolved CO gas disks have recently been detected with 7''-spatial-resolution, interferometry mapping. These gas disks confirm the conclusions of earlier results showing optical dust lanes and unresolved CO that implied the common existence of molecular gas in ellipticals and the disk-like structure of this gas. The presented SHARC II mapping results provide the best to date resolved FIR-submm extent of NGC 3656 and of any elliptical galaxy > 40 Mpc, showing that dust of 29 K exists out to at least 1.8 kpc in this galaxy. These new data are used in conjunction with the archival HST maps and other published data to determine dust properties and associations with galactic structures, including dominant heating sources such as nuclear-activity, star-formation or diffuse-stellar radiation.

Using the SHARC-II camera at the Caltech Submillimeter Observatory to obtain 350 micron images of sources detected with the MIPS instrument on Spitzer, we have discovered a remarkable object at z=1.325+/-0.002 with an apparent Far-Infrared luminosity of 3.2(+/-0.7) x 10^13 Lsun. Unlike other z>1 sources of comparable luminosity selected from mid-IR surveys, MIPS J142824.0+352619 lacks any trace of AGN activity, and is likely a luminous analog of galaxies selected locally by IRAS, or at high redshift in the submillimeter. This source appears to be lensed by a foreground elliptical galaxy at z=1.034, although the amplification is likely modest (~10). We argue that the contribution to the observed optical/Near-IR emission from the foreground galaxy is small, and hence are able to present the rest-frame UV through radio Spectral Energy Distribution of this galaxy. Due to its unusually high luminosity, MIPS J142824.0+352619 presents a unique chance to study a high redshift dusty starburst galaxy in great detail.

We have imaged the circumstellar disk of Fomalhaut at 350 microns wavelength, using SHARC II at the Caltech Submillimeter Observatory. The spatial resolution of the raw images (9 arcsec) has been enhanced by a factor of three using the HiRes deconvolution procedure. We find that at this wavelength and signal to noise ratio (~ 12), the observed morphology is that of a simple inclined ring (i ~ 70 deg), with little or no other apparent structure; this is the first observation that shows clearly the ring morphology of the disk. We have combined our 350 micron data with Spitzer Space Telescope images at 24, 70, and 160 microns in order to estimate the 2-dimensional spatial variation of relative column density ("tau map") using our DISKFIT procedure. The tau map is based on the following physical assumptions: (1) the wavelength variation of opacity is the same throughout the disk, (2) the radial variation of dust temperature is dictated by the energy balance of individual grains in the stellar radiation field, and (3) the vertical scale height of the disk follows a power-law radial variation. The results confirm the ring-like morphology, but also show that the geometric center is displaced from the star by about 8 AU and that the ring has an apocentric enhancement of approximately 14% in integrated column density. If we interpret the displacement in terms of elliptical orbital motion due to gravitational perturbation by an unseen planet, then the implied forced eccentricity is ~ 0.06; dynamical modeling then predicts an apocentric density enhancement consistent with that inferred from the tau map.

We use extensive 350 um polarimetry and continuum maps obtained with Hertz and SHARC II along with HCN and HCO+ spectroscopic data to trace the orientation of the magnetic field in the Orion A star-forming region. Using the polarimetry data, we find that the direction of the projection of the magnetic field in the plane of the sky relative to the orientation of the integral-shaped filament varies considerably as one moves from north to south. While in IRAS 05327-0457 and OMC-3 MMS 1-6 the projection of the field is primarily perpendicular to the filament it becomes better aligned with it at OMC-3 MMS 8-9 and well aligned with it at OMC-2 FIR 6. The OMC-2 FIR 4 cloud, located between the last two, is a peculiar object where we find almost no polarization. The projected angle of the field is more complicated in OMC-1 where it exhibits smooth variations in its orientation across the face of this massive complex. By combining the polarimetry and spectroscopic data we were able to measure a set of average values for the inclination angle of the magnetic field relative to the line of sight. We find that the field is oriented quite close to the plane of the sky in most places. More precisely, the inclination of the magnetic field is ~73 deg. around OMC-3 MMS 6, ~74 deg. at OMC-3 MMS 8-9, ~80 deg. at OMC-2 FIR 4, ~65 deg. in the northeastern part of OMC-1, and ~49 deg. in the Bar. We also present polarimetry data for the OMC-4 region located some 13 arcminutes south of OMC-1.

We demonstrate that the combination of Zeeman, polarimetry and ion-to-neutral molecular line width ratio measurements permits the determination of the magnitude and orientation of the magnetic field in the weakly ionized parts of molecular clouds. Zeeman measurements provide the strength of the magnetic field along the line of sight, polarimetry measurements give the field orientation in the plane of the sky and the ion-to-neutral molecular line width ratio determines the angle between the magnetic field and the line of sight. We apply the technique to the M17 star-forming region using a HERTZ 350 um polarimetry map and HCO+-to-HCN molecular line width ratios to provide the first three-dimensional view of the magnetic field in M17.

We report the detection of polarized emission in the vicinity of the Galactic center for 158 positions within eight different pointings of the Hertz polarimeter operating on the Caltech Submillimeter Observatory. These pointings include positions 2 arcminutes offset to the E, NE, and NW of M-0.02-0.07; positions to the SE and NW of the 20 km/s cloud (M-0.13-0.08), CO+0.02-0.02, M+0.07-0.08, and M+0.11-0.08. We use these data in conjunction with previous far-infrared and submillimeter polarization results to find that the direction of the inferred magnetic field is related to the density of the molecular material in the following way: in denser regions, the field is generally parallel to the Galactic plane, whereas in regions with lower densities, the field is generally perpendicular to the plane. This finding is consistent with a model in which an initially poloidal field has been sheared into a toroidal configuration in regions that are dense enough such that the gravitational energy density is greater than the energy density of the magnetic field. Using this model, we estimate the characteristic strength of the magnetic field in the central 30 pc of our Galaxy to be a few mG.

Using the 20-pixel camera SHARC, we have surveyed the 350 micron emission from a 60 arcminute by 12 arcminute region in the Galactic center. A comparable region has been observed at 800 microns by Lis & Carlstrom (1994); the SHARC map has better spatial resolution and sensitivity to extended emission, however. This paper introduces several features not visible in prior maps of dust emission. We have detected a 10+ pc band of emission peaking at l = -0.14 degrees, b = +0.02 degrees, which is probably associated with the 'negative velocity arc' observed in 13CO (Bally et al. 1988). We suggest an association of a rounded 350 micron feature near the Radio Arc with an expanding shell observed in CS (Tsuboi, Ukita, & Handa 1997). Faint emission at the same Galactic longitude as the Dust Ridge (Lis & Carlstrom 1994) but at opposite latitude is visible as well as a compact source at the base of a thermal radio filament at negative Galactic latitude. Diffuse dust emission within a few arcminutes of the Sgr D core is detected for the first time. 350 micron sources are associated with both H2O masers in Sgr D, but only one out of four OH masers (Mehringer et al. 1998). The mean 350 micron / 800 micron flux ratio for the Galactic center is approximately 17 (beta = 2.0) over the map but is higher (beta = 2.5) in parts of the Dust Ridge and lower in Sgr B2 (N).

The advent of submillimeter wavelength array cameras operating on large ground-based telescopes is revolutionizing imaging at these wavelengths, enabling high-resolution submillimeter surveys of dust emission in star-forming regions and galaxies. Here we present a recent 350 micron image of the edge-on galaxy NGC 891, which was obtained with the Submillimeter High Angular Resolution Camera (SHARC) at the Caltech Submillimeter Observatory (CSO). We find that high resolution submillimeter data is a vital complement to shorter wavelength satellite data, which enables a reliable separation of the cold dust component seen at millimeter wavelengths from the warmer component which dominates the far-infrared (FIR) luminosity.