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.

The Probe of Inflation and Cosmic Origins (PICO) is a proposed probe-scale space mission consisting of an imaging polarimeter operating in frequency bands between 20 and 800 GHz. We describe the science achievable by PICO, which has sensitivity equivalent to more than 3300 Planck missions, the technical implementation, the schedule and cost.

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 compare the magnetic field orientation for the young giant molecular cloud Vela C inferred from 500-$\mu$m polarization maps made with the BLASTPol balloon-borne polarimeter to the orientation of structures in the integrated line emission maps from Mopra observations. Averaging over the entire cloud we find that elongated structures in integrated line-intensity, or zeroth-moment maps, for low density tracers such as $^{12}$CO and $^{13}$CO $J$ $\rightarrow$ 1 - 0 are statistically more likely to align parallel to the magnetic field, while intermediate or high density tracers show (on average) a tendency for alignment perpendicular to the magnetic field. This observation agrees with previous studies of the change in relative orientation with column density in Vela C, and supports a model where the magnetic field is strong enough to have influenced the formation of dense gas structures within Vela C. The transition from parallel to no preferred/perpendicular orientation appears to happen between the densities traced by $^{13}$CO and by C$^{18}$O $J$ $\rightarrow$ 1 - 0. Using RADEX radiative transfer models to estimate the characteristic number density traced by each molecular line we find that the transition occurs at a molecular hydrogen number density of approximately $10^3$ cm$^{-3}$. We also see that the Centre-Ridge (the highest column density and most active star-forming region within Vela C) appears to have a transition at a lower number density, suggesting that this may depend on the evolutionary state of the cloud.

(Abbreviated) We investigate the interplay between magnetic (B) field, gravity, and turbulence in the fragmentation process of cores within the filamentary infrared dark cloud G34.43+00.24. We observe the magnetic field (B) morphology across G34.43 and compare with the kinematics obtained from N2H+ across the filament. We derive local velocity gradients from N2H+, tracing motion in the plane of sky, and compare with the observed local B field orientations in the plane of sky. Besides a large-scale east-west velocity gradient, we find a close alignment between local B field orientations and local velocity gradients toward the MM1/MM2 ridge. This local correlation in alignment suggests that gas motions are influenced by the magnetic field morphology or vice versa. Additionally, this alignment seems to be getting even closer with increasing integrated emission in N2H+, possibly indicating that a growing gravitational pull is more and more aligning B field and gas motion. We analyze and quantify B field, gravity, turbulence, and their relative importance toward the MM1, MM2 and MM3 regions with various techniques over two scales, a larger clump area at 2 pc scale and the smaller core area at 0.6 pc scale. While gravitational energy, B field, and turbulent pressure all grow systematically from large to small scale, the ratios among the three constituents develop clearly differently over scale. We propose that this varying relative importance between B field, gravity, and turbulence over scale drives and explains the different fragmentation types seen at sub-pc scale (no fragmentation in MM1; aligned fragmentation in MM2; clustered fragmentation in MM3). We discuss uncertainties, subtleties, and the robustness of our conclusion, and we stress the need of a multi-scale joint analysis to understand the dynamics in these systems.

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.

The Probe of Inflation and Cosmic Origins (PICO) is an imaging polarimeter that will scan the sky for 5 years in 21 frequency bands spread between 21 and 799 GHz. It will produce full-sky surveys of intensity and polarization with a final combined-map noise level of 0.87 $\mu$K arcmin for the required specifications, equivalent to 3300 Planck missions, and with our current best-estimate would have a noise level of 0.61 $\mu$K arcmin (6400 Planck missions). PICO will either determine the energy scale of inflation by detecting the tensor to scalar ratio at a level $r=5\times 10^{-4}~(5\sigma)$, or will rule out with more than $5\sigma$ all inflation models for which the characteristic scale in the potential is the Planck scale. With LSST's data it could rule out all models of slow-roll inflation. PICO will detect the sum of neutrino masses at $>4\sigma$, constrain the effective number of light particle species with $\Delta N_{\rm eff}<0.06~(2\sigma)$, and elucidate processes affecting the evolution of cosmic structures by measuring the optical depth to reionization with errors limited by cosmic variance and by constraining the evolution of the amplitude of linear fluctuations $\sigma_{8}(z)$ with sub-percent accuracy. Cross-correlating PICO's map of the thermal Sunyaev-Zeldovich effect with LSST's gold sample of galaxies will precisely trace the evolution of thermal pressure with $z$. PICO's maps of the Milky Way will be used to determine the make up of galactic dust and the role of magnetic fields in star formation efficiency. With 21 full sky legacy maps in intensity and polarization, which cannot be obtained in any other way, the mission will enrich many areas of astrophysics. PICO is the only single-platform instrument with the combination of sensitivity, angular resolution, frequency bands, and control of systematic effects that can deliver this compelling, timely, and broad science.

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.

Linear polarization maps of the Carina Nebula were obtained at 250, 350, and 500 $\mu$m during the 2012 flight of the BLASTPol balloon-borne telescope. These measurements are combined with Planck 850 $\mu$m data in order to produce a submillimeter spectrum of the polarization fraction of the dust emission, averaged over the cloud. This spectrum is flat to within $\pm$15% (relative to the 350 $\mu$m polarization fraction). In particular, there is no evidence for a pronounced minimum of the spectrum near 350 $\mu$m, as suggested by previous ground-based measurements of other molecular clouds. This result of a flat polarization spectrum in Carina is consistent with recently-published BLASTPol measurements of the Vela C molecular cloud, and also agrees with a published model for an externally-illuminated, dense molecular cloud by Bethell and collaborators. The shape of the spectrum in Carina does not show any dependence on the radiative environment of the dust, as quantified by the Planck-derived dust temperature or dust optical depth at 353 GHz.

The Next Generation Balloon-borne Large Aperture Submillimeter Telescope (BLAST-TNG) is a submillimeter mapping experiment planned for a 28 day long-duration balloon (LDB) flight from McMurdo Station, Antarctica during the 2018-2019 season. BLAST-TNG will detect submillimeter polarized interstellar dust emission, tracing magnetic fields in galactic molecular clouds. BLAST-TNG will be the first polarimeter with the sensitivity and resolution to probe the $\sim$0.1 parsec-scale features that are critical to understanding the origin of structures in the interstellar medium. With three detector arrays operating at 250, 350, and 500 $\mu$m (1200, 857, and 600 GHz), BLAST-TNG will obtain diffraction-limited resolution at each waveband of 30, 41, and 59 arcseconds respectively. To achieve the submillimeter resolution necessary for its science goals, the BLAST-TNG telescope features a 2.5 m aperture carbon fiber composite primary mirror, one of the largest mirrors flown on a balloon platform. Successful performance of such a large telescope on a balloon-borne platform requires stiff, lightweight optical components and mounting structures. Through a combination of optical metrology and finite element modeling of thermal and mechanical stresses on both the telescope optics and mounting structures, we expect diffraction-limited resolution at all our wavebands. We expect pointing errors due to deformation of the telescope mount to be negligible. We have developed a detailed thermal model of the sun shielding, gondola, and optical components to optimize our observing strategy and increase the stability of the telescope over the flight. We present preflight characterization of the telescope and its platform.

The Next Generation Balloon-borne Large Aperture Submillimeter Telescope (BLAST-TNG) is a submillimeter mapping experiment planned for a 28 day long-duration balloon (LDB) flight from McMurdo Station, Antarctica during the 2018-2019 season. BLAST-TNG will detect submillimeter polarized interstellar dust emission, tracing magnetic fields in galactic molecular clouds. BLAST-TNG will be the first polarimeter with the sensitivity and resolution to probe the $\sim$0.1 parsec-scale features that are critical to understanding the origin of structures in the interstellar medium. BLAST-TNG features three detector arrays operating at wavelengths of 250, 350, and 500 $\mu$m (1200, 857, and 600 GHz) comprised of 918, 469, and 272 dual-polarization pixels, respectively. Each pixel is made up of two crossed microwave kinetic inductance detectors (MKIDs). These arrays are cooled to 275 mK in a cryogenic receiver. Each MKID has a different resonant frequency, allowing hundreds of resonators to be read out on a single transmission line. This inherent ability to be frequency-domain multiplexed simplifies the cryogenic readout hardware, but requires careful optical testing to map out the physical location of each resonator on the focal plane. Receiver-level optical testing was carried out using both a cryogenic source mounted to a movable xy-stage with a shutter, and a beam-filling, heated blackbody source able to provide a 10-50 $^\circ$C temperature chop. The focal plane array noise properties, responsivity, polarization efficiency, instrumental polarization were measured. We present the preflight characterization of the BLAST-TNG cryogenic system and array-level optical testing of the MKID detector arrays in the flight receiver.

BFORE is a high-altitude ultra-long-duration balloon mission to map the cosmic microwave background (CMB). During a 28-day mid-latitude flight launched from Wanaka, New Zealand, the instrument will map half the sky to improve measurements of the optical depth to reionization tau. This will break parameter degeneracies needed to detect neutrino mass. BFORE will also hunt for the gravitational wave B-mode signal, and map Galactic dust foregrounds. The mission will be the first near-space use of TES/mSQUID multichroic detectors (150/217 GHz and 280/353 GHz bands) with low-power readout electronics.

TolTEC is a new camera being built for the 50-meter Large Millimeter-wave Telescope (LMT) in Puebla, Mexico to survey distant galaxies and star-forming regions in the Milky Way. The optical design simultaneously couples the field of view onto focal planes at 150, 220, and 280 GHz. The optical design and detector properties, as well as a data-driven model of the atmospheric emission of the LMT site, inform the sensitivity model of the integrated instrument. This model is used to optimize the instrument design, and to calculate the mapping speed as an early forecast of the science reach of the instrument.

How important is the magnetic (B-) field when compared to gravity and turbulence in the star-formation process? Does its importance depend on scale and location? We summarize submm dust polarization observations towards the large filamentary infrared dark cloud G34 and towards a dense core in the high-mass star-forming region W51. We detect B-field orientations that are either perpendicular or parallel to the G34 filament axis. These B-field orientations further correlate with local velocity gradients. Towards three cores in G34 we find a varying importance between B-field, gravity, and turbulence that seems to dictate varying types of fragmentation. At highest resolution towards the gravity-dominated collapsing core W51 e2 we resolve new B-field features, such as converging B-field lines and possibly magnetic channels.

The relative orientation between filamentary structures in molecular clouds and the ambient magnetic field provides insight into filament formation and stability. To calculate the relative orientation, a measurement of filament orientation is first required. We propose a new method to calculate the orientation of the one pixel wide filament skeleton that is output by filament identification algorithms such as \textsc{filfinder}. We derive the local filament orientation from the direction of the intensity gradient in the skeleton image using the Sobel filter and a few simple post-processing steps. We call this the `Sobel-gradient method'. The resulting filament orientation map can be compared quantitatively on a local scale with the magnetic field orientation map to then find the relative orientation of the filament with respect to the magnetic field at each point along the filament. It can also be used in constructing radial profiles for filament width fitting. The proposed method facilitates automation in analysis of filament skeletons, which is imperative in this era of `big data'.

Polarized emission from aligned dust is a crucial tool for studies of magnetism in the ISM and a troublesome contaminant for studies of CMB polarization. In each case, an understanding of the significance of the polarization signal requires well-calibrated physical models of dust grains. Despite decades of progress in theory and observation, polarized dust models remain largely underconstrained. During its 2012 flight, the balloon-borne telescope BLASTPol obtained simultaneous broad-band polarimetric maps of a translucent molecular cloud at 250, 350, and 500 microns. Combining these data with polarimetry from the Planck 850 micron band, we have produced a submillimeter polarization spectrum for a cloud of this type for the first time. We find the polarization degree to be largely constant across the four bands. This result introduces a new observable with the potential to place strong empirical constraints on ISM dust polarization models in a previously inaccessible density regime. Comparing with models by Draine and Fraisse (2009), our result disfavors two of their models for which all polarization arises due only to aligned silicate grains. By creating simple models for polarized emission in a translucent cloud, we verify that extinction within the cloud should have only a small effect on the polarization spectrum shape compared to the diffuse ISM. Thus we expect the measured polarization spectrum to be a valid check on diffuse ISM dust models. The general flatness of the observed polarization spectrum suggests a challenge to models where temperature and alignment degree are strongly correlated across major dust components.

BFORE is a NASA high-altitude ultra-long-duration balloon mission proposed to measure the cosmic microwave background (CMB) across half the sky during a 28-day mid-latitude flight launched from Wanaka, New Zealand. With the unique access to large angular scales and high frequencies provided by the balloon platform, BFORE will significantly improve measurements of the optical depth to reionization tau, breaking parameter degeneracies needed for a measurement of neutrino mass with the CMB. The large angular scale data will enable BFORE to hunt for the large-scale gravitational wave B-mode signal, as well as the degree-scale signal, each at the r~0.01 level. The balloon platform allows BFORE to map Galactic dust foregrounds at frequencies where they dominate, in order to robustly separate them from CMB signals measured by BFORE, in addition to complementing data from ground-based telescopes. The combination of frequencies will also lead to velocity measurements for thousands of galaxy clusters, as well as probing how star-forming galaxies populate dark matter halos. The mission will be the first near-space use of TES multichroic detectors (150/217 GHz and 280/353 GHz bands) using highly-multiplexed mSQUID microwave readout, raising the technical readiness level of both technologies.

Different combinations of input parameters to filament identification algorithms, such as Disperse and FilFinder, produce numerous different output skeletons. The skeletons are a one pixel wide representation of the filamentary structure in the original input image. However, these output skeletons may not necessarily be a good representation of that structure. Furthermore, a given skeleton may not be as good a representation as another. Previously there has been no mathematical `goodness-of-fit' measure to compare output skeletons to the input image. Thus far this has been assessed visually, introducing visual bias. We propose the application of the mean structural similarity index (MSSIM) as a mathematical goodness-of-fit measure. We describe the use of the MSSIM to find the output skeletons most mathematically similar to the original input image (the optimum, or `best', skeletons) for a given algorithm, and independently of the algorithm. This measure makes possible systematic parameter studies, aimed at finding the subset of input parameter values returning optimum skeletons. It can also be applied to the output of non-skeleton based filament identification algorithms, such as the Hessian matrix method. The MSSIM removes the need to visually examine thousands of output skeletons, and eliminates the visual bias, subjectivity, and limited reproducibility inherent in that process, representing a major improvement on existing techniques. Importantly, it also allows further automation in the post-processing of output skeletons, which is crucial in this era of `big data'.

We present a large-scale combination of near-infrared (near-IR) interstellar polarization data from background starlight with polarized emission data at submillimeter (sub-mm) wavelengths for the Vela C molecular cloud. The near-IR data consist of more than 6700 detections probing a range of visual extinctions between $2$ and $20\,$mag in and around the cloud. The sub-mm data was collected in Antartica by the Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol). This is the first direct combination of near-IR and sub-mm polarization data for a molecular cloud aimed at measuring the "polarization efficiency ratio" ($R_{\mathrm{eff}}$), a quantity that is expected to depend only on grain intrinsic physical properties. It is defined as $p_{500}/(p_{I}/\tau_{V})$, where $p_{500}$ and $p_{I}$ are polarization fractions at $500\,\mu$m and $I$-band, respectively, and $\tau_{V}$ is the optical depth. To ensure that the same column density of material is producing both polarization from emission and from extinction, we conducted a careful selection of near-background stars using 2MASS, $Herschel$ and $Planck$ data. This selection excludes objects contaminated by the Galactic diffuse background material as well as objects located in the foreground. Accounting for statistical and systematic uncertainties, we estimate an average $R_{\mathrm{eff}}$ value of $2.4\pm0.8$, which can be used to test the predictions of dust grain models designed for molecular clouds when such predictions become available. $R_{\mathrm{eff}}$ appears to be relatively flat as a function of the cloud depth for the range of visual extinctions probed.

We statistically evaluate the relative orientation between gas column density structures, inferred from Herschel submillimetre observations, and the magnetic field projected on the plane of sky, inferred from polarized thermal emission of Galactic dust observed by BLASTPol at 250, 350, and 500 micron, towards the Vela C molecular complex. First, we find very good agreement between the polarization orientations in the three wavelength-bands, suggesting that, at the considered common angular resolution of 3.0 arcminutes that corresponds to a physical scale of approximately 0.61 pc, the inferred magnetic field orientation is not significantly affected by temperature or dust grain alignment effects. Second, we find that the relative orientation between gas column density structures and the magnetic field changes progressively with increasing gas column density, from mostly parallel or having no preferred orientation at low column densities to mostly perpendicular at the highest column densities. This observation is in agreement with previous studies by the Planck collaboration towards more nearby molecular clouds. Finally, we find a correspondence between the trends in relative orientation and the shape of the column density probability distribution functions. In the sub-regions of Vela C dominated by one clear filamentary structure, or "ridges", we find a sharp transition from preferentially parallel or having no preferred relative orientation at low column densities to preferentially perpendicular at highest column densities. In the sub-regions of Vela C dominated by several filamentary structures with multiple orientations, or "nests", such a transition is also present, but it is clearly less sharp than in the ridge-like sub-regions. Both of these results suggest that the magnetic field is dynamically important for the formation of density structures in this region.

Polarized thermal emission from interstellar dust grains can be used to map magnetic fields in star forming molecular clouds and the diffuse interstellar medium (ISM). The Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol) flew from Antarctica in 2010 and 2012 and produced degree-scale polarization maps of several nearby molecular clouds with arcminute resolution. The success of BLASTPol has motivated a next-generation instrument, BLAST-TNG, which will use more than 3000 linear polarization sensitive microwave kinetic inductance detectors (MKIDs) combined with a 2.5m diameter carbon fiber primary mirror to make diffraction-limited observations at 250, 350, and 500 $\mu$m. With 16 times the mapping speed of BLASTPol, sub-arcminute resolution, and a longer flight time, BLAST-TNG will be able to examine nearby molecular clouds and the diffuse galactic dust polarization spectrum in unprecedented detail. The 250 $\mu$m detector array has been integrated into the new cryogenic receiver, and is undergoing testing to establish the optical and polarization characteristics of the instrument. BLAST-TNG will demonstrate the effectiveness of kilo-pixel MKID arrays for applications in submillimeter astronomy. BLAST-TNG is scheduled to fly from Antarctica in December 2017 for 28 days and will be the first balloon-borne telescope to offer a quarter of the flight for "shared risk" observing by the community.

Polarization maps of the Vela C molecular cloud were obtained at 250, 350, and 500um during the 2012 flight of the balloon-borne telescope BLASTPol. These measurements are used in conjunction with 850um data from Planck to study the submillimeter spectrum of the polarization fraction for this cloud. The spectrum is relatively flat and does not exhibit a pronounced minimum at \lambda ~350um as suggested by previous measurements of other molecular clouds. The shape of the spectrum does not depend strongly on the radiative environment of the dust, as quantified by the column density or the dust temperature obtained from Herschel data. The polarization ratios observed in Vela C are consistent with a model of a porous clumpy molecular cloud being uniformly heated by the interstellar radiation field.

We present results for Vela C obtained during the 2012 flight of the Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol). We mapped polarized intensity across almost the entire extent of this giant molecular cloud, in bands centered at 250, 350, and 500 {\mu}m. In this initial paper, we show our 500 {\mu}m data smoothed to a resolution of 2.5 arcminutes (approximately 0.5 pc). We show that the mean level of the fractional polarization p and most of its spatial variations can be accounted for using an empirical three-parameter power-law fit, p = p_0 N^(-0.4) S^(-0.6), where N is the hydrogen column density and S is the polarization-angle dispersion on 0.5 pc scales. The decrease of p with increasing S is expected because changes in the magnetic field direction within the cloud volume sampled by each measurement will lead to cancellation of polarization signals. The decrease of p with increasing N might be caused by the same effect, if magnetic field disorder increases for high column density sightlines. Alternatively, the intrinsic polarization efficiency of the dust grain population might be lower for material along higher density sightlines. We find no significant correlation between N and S. Comparison of observed submillimeter polarization maps with synthetic polarization maps derived from numerical simulations provides a promising method for testing star formation theories. Realistic simulations should allow for the possibility of variable intrinsic polarization efficiency. The measured levels of correlation among p, N, and S provide points of comparison between observations and simulations.

The B-mode Foreground Experiment (BFORE) is a proposed NASA balloon project designed to make optimal use of the sub-orbital platform by concentrating on three dust foreground bands (270, 350, and 600 GHz) that complement ground-based cosmic microwave background (CMB) programs. BFORE will survey ~1/4 of the sky with 1.7 - 3.7 arcminute resolution, enabling precise characterization of the Galactic dust that now limits constraints on inflation from CMB B-mode polarization measurements. In addition, BFORE's combination of frequency coverage, large survey area, and angular resolution enables science far beyond the critical goal of measuring foregrounds. BFORE will constrain the velocities of thousands of galaxy clusters, provide a new window on the cosmic infrared background, and probe magnetic fields in the interstellar medium. We review the BFORE science case, timeline, and instrument design, which is based on a compact off-axis telescope coupled to >10,000 superconducting detectors.

AGN feedback from supermassive black holes (SMBHs) at the center of early type galaxies is commonly invoked as the explanation for the quenching of star formation in these systems. The situation is complicated by the significant amount of mass injected in the galaxy by the evolving stellar population over cosmological times. In absence of feedback, this mass would lead to unobserved galactic cooling flows, and to SMBHs two orders of magnitude more massive than observed. By using high-resolution 2D hydrodynamical simulations with radiative transport and star formation in state-of-the-art galaxy models, we show how the intermittent AGN feedback is highly structured on spatial and temporal scales, and how its effects are not only negative (shutting down the recurrent cooling episodes of the ISM), but also positive, inducing star formation in the inner regions of the host galaxy.

We report the first detection of linearly polarized emission at an observing wavelength of 350 $\mu$m from the radio-loud Active Galactic Nucleus 3C 279. We conducted polarization observations for 3C 279 using the SHARP polarimeter in the Caltech Submillimeter Observatory on 2014 March 13 and 14. For the first time, we detected the linear polarization with the degree of polarization of 13.3\%$\pm$3.4\% ($3.9\sigma$) and the Electric Vector Position Angle (EVPA) of 34.7$^\circ\pm5.6^\circ$. We also observed 3C~279 simultaneously at 13, 7, and 3.5 mm in dual polarization with the Korean very long baseline interferometry (VLBI) Network on 2014 March 6 (single dish) and imaged in milliarcsecond (mas) scales at 13, 7, 3.5, and 2.3 mm on March 22 (VLBI). We found that the degree of linear polarization increases from 10\% to 13\% at 13 mm to 350 $\mu$m and the EVPAs at all observing frequencies are parallel within $<10^\circ$ to the direction of the jet at mas scale, implying that the integrated magnetic fields are perpendicular to the jet in the innermost regions. We also found that the Faraday rotation measures RM are in a range of $-6.5\times10^2 \sim -2.7\times10^3$ rad m$^{-2}$ between 13 and 3.5 mm, and are scaled as a function of wavelength: $|{\rm RM}|\propto\lambda^{-2.2}$. These results indicate that the mm and sub-mm polarization emission are generated in the compact jet within 1~mas scale and affected by a Faraday screen in or in the close proximity of the jet.

We present first results of neutral carbon ([CI], 3P1 - 3P0 at 492 GHz) and carbon monoxide (13CO, J = 1 - 0) mapping in the Vela Molecular Ridge cloud C (VMR-C) and G333 giant molecular cloud complexes with the NANTEN2 and Mopra telescopes. For the four regions mapped in this work, we find that [CI] has very similar spectral emission profiles to 13CO, with comparable line widths. We find that [CI] has opacity of 0.1 - 1.3 across the mapped region while the [CI]/13CO peak brightness temperature ratio is between 0.2 to 0.8. The [CI] column density is an order of magnitude lower than that of 13CO. The H2 column density derived from [CI] is comparable to values obtained from 12CO. Our maps show CI is preferentially detected in gas with low temperatures (below 20 K), which possibly explains the comparable H2 column density calculated from both tracers (both CI and 12CO underestimate column density), as a significant amount of the CI in the warmer gas is likely in the higher energy state transition ([CI], 3P2 - 3P1 at 810 GHz), and thus it is likely that observations of both the above [CI] transitions are needed in order to recover the total H2 column density.

For the Class 0 protostar, L1527, we compare 131 polarization vectors from SCUPOL/JCMT, SHARP/CSO and TADPOL/CARMA observations with the corresponding model polarization vectors of four ideal-MHD, non-turbulent, cloud core collapse models. These four models differ by their initial magnetic fields before collapse; two initially have aligned fields (strong and weak) and two initially have orthogonal fields (strong and weak) with respect to the rotation axis of the L1527 core. Only the initial weak orthogonal field model produces the observed circumstellar disk within L1527. This is a characteristic of nearly all ideal-MHD, non-turbulent, core collapse models. In this paper we test whether this weak orthogonal model also has the best agreement between its magnetic field structure and that inferred from the polarimetry observations of L1527. We found that this is not the case; based on the polarimetry observations the most favored model of the four is the weak aligned model. However, this model does not produce a circumstellar disk, so our result implies that a non-turbulent, ideal-MHD global collapse model probably does not represent the core collapse that has occurred in L1527. Our study also illustrates the importance of using polarization vectors covering a large area of a cloud core to determine the initial magnetic field orientation before collapse; the inner core magnetic field structure can be highly altered by a collapse and so measurements from this region alone can give unreliable estimates of the initial field configuration before collapse.

The Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol) was a suborbital experiment designed to map magnetic fields in order to study their role in star formation processes. BLASTPol made detailed polarization maps of a number of molecular clouds during its successful flights from Antarctica in 2010 and 2012. We present the next-generation BLASTPol instrument (BLAST-TNG) that will build off the success of the previous experiment and continue its role as a unique instrument and a test bed for new technologies. With a 16-fold increase in mapping speed, BLAST-TNG will make larger and deeper maps. Major improvements include a 2.5 m carbon fiber mirror that is 40% wider than the BLASTPol mirror and ~3000 polarization sensitive detectors. BLAST-TNG will observe in three bands at 250, 350, and 500 microns. The telescope will serve as a pathfinder project for microwave kinetic inductance detector (MKID) technology, as applied to feedhorn coupled submillimeter detector arrays. The liquid helium cooled cryostat will have a 28-day hold time and will utilize a closed-cycle $^3$He refrigerator to cool the detector arrays to 270 mK. This will enable a detailed mapping of more targets with higher polarization resolution than any other submillimeter experiment to date. BLAST-TNG will also be the first balloon-borne telescope to offer shared risk observing time to the community. This paper outlines the motivation for the project and the instrumental design.

The Cosmology Large Angular Scale Surveyor (CLASS) is an experiment to measure the signature of a gravita-tional-wave background from inflation in the polarization of the cosmic microwave background (CMB). CLASS is a multi-frequency array of four telescopes operating from a high-altitude site in the Atacama Desert in Chile. CLASS will survey 70\% of the sky in four frequency bands centered at 38, 93, 148, and 217 GHz, which are chosen to straddle the Galactic-foreground minimum while avoiding strong atmospheric emission lines. This broad frequency coverage ensures that CLASS can distinguish Galactic emission from the CMB. The sky fraction of the CLASS survey will allow the full shape of the primordial B-mode power spectrum to be characterized, including the signal from reionization at low $\ell$. Its unique combination of large sky coverage, control of systematic errors, and high sensitivity will allow CLASS to measure or place upper limits on the tensor-to-scalar ratio at a level of $r=0.01$ and make a cosmic-variance-limited measurement of the optical depth to the surface of last scattering, $\tau$.

Two-dimensional integral field surveys such as ATLAS^3D are producing rich observational data sets yielding insights into galaxy formation. These new kinematic observations have highlighted the need to understand the evolutionary mechanisms leading to a spectrum of fast-rotators and slow-rotators in early-type galaxies. We address the formation of slow and fast rotators through a series of controlled, comprehensive hydrodynamical simulations sampling idealized galaxy merger scenarios constructed from model spiral galaxies. Idealized and controlled simulations of this sort complement the more 'realistic' cosmological simulations by isolating and analyzing the effects of specific parameters, as we do in this paper. We recreate minor and major binary mergers, binary merger trees with multiple progenitors, and multiple sequential mergers. Within each of these categories of formation history, we correlate progenitor gas fraction, mass ratio, orbital pericenter, orbital ellipticity, and spin with remnant kinematic properties. We create kinematic profiles of these 95 simulations comparable to ATLAS^3D data. By constructing remnant profiles of the projected specific angular momentum (lambda_R = <R|V|> / <sqrt(V^2+sigma^2)>, triaxiality, and measuring the incidences of kinematic twists and kinematically decoupled cores, we distinguish between varying formation scenarios. We find that binary mergers nearly always form fast rotators. Slow rotators can be formed from zero initial angular momentum configurations and gas-poor mergers, but are not as round as the ATLAS^3D galaxies. Remnants of binary merger trees are triaxial slow rotators. Sequential mergers form round slow rotators that most resemble the ATLAS^3D rotators.

An attitude determination system for balloon-borne experiments is presented. The system provides pointing information in azimuth and elevation for instruments flying on stratospheric balloons over Antarctica. In-flight attitude is given by the real-time combination of readings from star cameras, a magnetometer, sun sensors, GPS, gyroscopes, tilt sensors and an elevation encoder. Post-flight attitude reconstruction is determined from star camera solutions, interpolated by the gyroscopes using an extended Kalman Filter. The multi-sensor system was employed by the Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol), an experiment that measures polarized thermal emission from interstellar dust clouds. A similar system was designed for the upcoming flight of SPIDER, a Cosmic Microwave Background polarization experiment. The pointing requirements for these experiments are discussed, as well as the challenges in designing attitude reconstruction systems for high altitude balloon flights. In the 2010 and 2012 BLASTPol flights from McMurdo Station, Antarctica, the system demonstrated an accuracy of <5' rms in-flight, and <5" rms post-flight.

The Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol) is a suborbital mapping experiment, designed to study the role played by magnetic fields in the star formation process. BLASTPol observes polarized light using a total power instrument, photolithographic polarizing grids, and an achromatic half-wave plate to modulate the polarization signal. During its second flight from Antarctica in December 2012, BLASTPol made degree scale maps of linearly polarized dust emission from molecular clouds in three wavebands, centered at 250, 350, and 500 microns. The instrumental performance was an improvement over the 2010 BLASTPol flight, with decreased systematics resulting in a higher number of confirmed polarization vectors. The resultant dataset allows BLASTPol to trace magnetic fields in star-forming regions at scales ranging from cores to entire molecular cloud complexes.

The Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol) is a suborbital mapping experiment designed to study the role magnetic fields play in star formation. BLASTPol has had two science flights from McMurdo Station, Antarctica in 2010 and 2012. These flights have produced thousands of polarization vectors at 250, 350 and 500 microns in several molecular cloud targets. We present the design, specifications, and progress towards the next-generation BLASTPol experiment (BLAST-TNG). BLAST-TNG will fly a 40% larger diameter primary mirror, with almost 8 times the number of polarization-sensitive detectors resulting in a factor of 16 increase in mapping speed. With a spatial resolution of 22 arcseconds and four times the field of view of BLASTPol, BLAST-TNG will bridge the angular scales between Planck's low resolution all-sky maps and ALMA's ultra-high resolution narrow fields. The new receiver has a larger cryogenics volume, allowing for a 28 day hold time. BLAST-TNG employs three arrays of Microwave Kinetic Inductance Detectors (MKIDs) with 30% fractional bandwidth at 250, 350 and 500 microns. In this paper, we will present the new BLAST-TNG instrument and science objectives.

We present the thermal model of the Balloon-borne Large-Aperture Submillimeter Telescope for Polarimetry (BLASTPol). This instrument was successfully flown in two circumpolar flights from McMurdo, Antarctica in 2010 and 2012. During these two flights, BLASTPol obtained unprecedented information about the magnetic field in molecular clouds through the measurement of the polarized thermal emission of interstellar dust grains. The thermal design of the experiment addresses the stability and control of the payload necessary for this kind of measurement. We describe the thermal modeling of the payload including the sun-shielding strategy. We present the in-flight thermal performance of the instrument and compare the predictions of the model with the temperatures registered during the flight. We describe the difficulties of modeling the thermal behavior of the balloon-borne platform and establish landmarks that can be used in the design of future balloon-borne instruments.

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.

Based on two-dimensional high resolution hydrodynamic numerical simulation, we study the mechanical and radiative feedback effects from the central AGN on the cosmological evolution of an isolated elliptical galaxy. Physical processes such as star formation and supernovae are considered. The inner boundary of the simulation domain is carefully chosen so that the fiducial Bondi radius is resolved and the accretion rate of the black hole is determined self-consistently. In analogy to previous works, we assume that the specific angular momentum of the galaxy is low. It is well-known that when the accretion rates are high and low, the central AGNs will be in cold and hot accretion modes, which correspond to the radiative and kinetic feedback modes, respectively. The emitted spectrum from the hot accretion flows is harder than that from the cold accretion flows, which could result in a higher Compton temperature accompanied by a more efficient radiative heating, according to previous theoretical works. Such a difference of the Compton temperature between the two feedback modes, the focus of this study, has been neglected in previous works. Significant differences in the kinetic feedback mode are found as a result of the stronger Compton heating and accretion becomes more chaotic. More importantly, if we constrain models to correctly predict black hole growth and AGN duty cycle after cosmological evolution, we find that the favored model parameters are constrained: mechanical feedback efficiency diminishes with decreasing luminosity (the maximum efficiency being $\simeq 10^{-3.5}$) and X-ray Compton temperature increases with decreasing luminosity, although models with fixed mechanical efficiency and Compton temperature can be found that are satisfactory as well. We conclude that radiative feedback in the kinetic mode is much more important than previously thought.

We show that recently documented trends in galaxy sizes with mass and redshift can be understood in terms of the influence of underlying cosmic evolution; a holistic view which is complimentary to interpretations involving the accumulation of discreet evolutionary processes acting on individual objects. Using standard cosmology theory, supported with results from the Millennium simulations, we derive expected size trends for collapsed cosmic structures, emphasising the important distinction between these trends and the assembly paths of individual regions. We then argue that the observed variation in the stellar mass content of these structures can be understood to first order in terms of natural limitations of cooling and feedback. But whilst these relative masses vary by orders of magnitude, galaxy and host radii have been found to correlate linearly. We explain how these two aspects will lead to galaxy sizes that closely follow observed trends and their evolution, comparing directly with the COSMOS and SDSS surveys. Thus we conclude that the observed minimum radius for galaxies, the evolving trend in size as a function of mass for intermediate systems, and the observed increase in the sizes of massive galaxies, may all be considered an emergent consequence of the cosmic expansion.

Turbulence and magnetic fields are expected to be important for regulating molecular cloud formation and evolution. However, their effects on subparsec to 100 parsec scales, leading to the formation of starless cores, is not well understood. We investigate the prestellar core structure morphologies obtained from analysis of the Herschel-SPIRE 350 $\mu$m maps of the Lupus I cloud. This distribution is first compared on a statistical basis to the large scale shape of the main filament. We find the distribution of the elongation position angle of the cores to be consistent with a random distribution, which means no specific orientation of the morphology of the cores is observed with respect to a large-scale filament shape model for Lupus I, or relative to a large-scale bent filament model. This distribution is also compared to the mean orientation of the large-scale magnetic fields probed at 350 $\mu$m with the Balloon-borne Large Aperture Telescope for Polarimetry (BLASTPol) during its 2010 campaign. Here again we do not find any correlation between the core morphology distribution and the average orientation of the magnetic fields on parsec scales. Our main conclusion is that the local filament dynamics - including secondary filaments that often run orthogonally to the primary filament - and possibly small-scale variations in the local magnetic field direction, could be the dominant factors for explaining the final orientation of each core.

The question whether magnetic fields play an important role in the processes of molecular cloud and star formation has been debated for decades. Recent observations have revealed a simple picture that may help illuminate these questions: magnetic fields have a tendency to preserve their orientation at all scales that have been probed - from 100-pc scale inter-cloud media down to sub-pc scale cloud cores. This ordered morphology has implications for the way in which self-gravity and turbulence interact with magnetic fields: both gravitational contraction and turbulent velocities should be anisotropic, due to the influence of dynamically important magnetic fields. Such anisotropy is now observed. Here we review these recent observations and discuss how they can improve our understanding of cloud/star formation.

We present a comparison of the molecular gas properties in the outflow vs. in the ambient medium of the local prototype radio-loud and ultraluminous-infrared galaxy 4C12.50 (IRAS13451+1232), using new data from the IRAM Plateau de Bure interferometer and 30m telescope, and the Herschel space telescope. Previous H_2 (0-0) S(1) and S(2) observations with the Spitzer space telescope had indicated that the warm (~400K) molecular gas in 4C12.50 is made up of a 1.4(+-0.2)x10^8 M_sun ambient reservoir and a 5.2(+-1.7)x10^7 M_sun outflow. The new CO(1-0) data cube indicates that the corresponding cold (25K) H_2 gas mass is 1.0(+-0.1)x10^10 M_sun for the ambient medium and <1.3x10^8 M_sun for the outflow, when using a CO-intensity-to-H_2-mass conversion factor alpha of 0.8 M_sun /(K km/s pc^2). The combined mass outflow rate is high, 230-800 M_sun/yr, but the amount of gas that could escape the galaxy is low. A potential inflow of gas from a 3.3(+-0.3)x10^8 M_sun tidal tail could moderate any mass loss. The mass ratio of warm-to-cold molecular gas is >= 30 times higher in the outflow than in the ambient medium, indicating that a non-negligible fraction of the accelerated gas is heated to temperatures at which star formation is inefficient. This conclusion is robust against the use of different alpha factor values, and/or different warm gas tracers (H_2 vs. H_2 plus CO): with the CO-probed gas mass being at least 40 times lower at 400K than at 25K, the total warm-to-cold mass ratio is always lower in the ambient gas than in the entrained gas. Heating of the molecular gas could facilitate the detection of new outflows in distant galaxies by enhancing their emission in intermediate rotational number CO lines.

A cryogenic achromatic half-wave plate (HWP) for submillimetre astronomical polarimetry has been designed, manufactured, tested, and deployed in the Balloon-borne Large-Aperture Submillimeter Telescope for Polarimetry (BLASTPol). The design is based on the five-slab Pancharatnam recipe and it works in the wavelength range 200-600 micron, making it the broadest-band HWP built to date at (sub)millimetre wavelengths. The frequency behaviour of the HWP has been fully characterised at room and cryogenic temperatures with incoherent radiation from a polarising Fourier transform spectrometer. We develop a novel empirical model, complementary to the physical and analytical ones available in the literature, that allows us to recover the HWP Mueller matrix and phase shift as a function of frequency and extrapolated to 4K. We show that most of the HWP non-idealities can be modelled by quantifying one wavelength-dependent parameter, the position of the HWP equivalent axes, which is then readily implemented in a map-making algorithm. We derive this parameter for a range of spectral signatures of input astronomical sources relevant to BLASTPol, and provide a benchmark example of how our method can yield improved accuracy on measurements of the polarisation angle on the sky at submillimetre wavelengths.

A recent analysis of black hole scaling relations, used to estimate the local mass density in black holes, has indicated that the normalization of the scaling relations should be increased by approximately a factor of five. The local black hole mass density is connected to the mean radiative efficiency of accretion through the time integral of the quasar volume density. The correspondence between this estimate of the radiative efficiency and that expected theoretically from thin-disk accretion has long been used as an argument that most of the growth in black holes occurs via luminous accretion. The increase of the mass density in black holes pushes the mean observed radiative efficiency to values below that expected for thin-disk accretion for any value of the black hole spin, including retrograde accretion disks. This can be accommodated via black hole growth channels that are intrinsically radiatively inefficient, such as super-Eddington accretion, or via growth channels that are intrinsically radiatively efficient but for which few of the photons are observed, such as Compton thick accretion. Measurements of the 30 keV peak in the X-ray background indicate a significant population of Compton thick sources which can explain some, but not all, of the change in the local black hole mass density. If this result is taken as evidence that super-Eddington accretion is common, it greatly reduces the tension associated with the growth timescale of observed z=7 quasars compared to the age of the universe at that redshift.

The Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLASTPol) was created by adding polarimetric capability to the BLAST experiment that was flown in 2003, 2005, and 2006. BLASTPol inherited BLAST's 1.8 m primary and its Herschel/SPIRE heritage focal plane that allows simultaneous observation at 250, 350, and 500 {\mu}m. We flew BLASTPol in 2010 and again in 2012. Both were long duration Antarctic flights. Here we present polarimetry of the nearby filamentary dark cloud Lupus I obtained during the 2010 flight. Despite limitations imposed by the effects of a damaged optical component, we were able to clearly detect submillimeter polarization on degree scales. We compare the resulting BLASTPol magnetic field map with a similar map made via optical polarimetry (The optical data were published in 1998 by J. Rizzo and collaborators.). The two maps partially overlap and are reasonably consistent with one another. We compare these magnetic field maps to the orientations of filaments in Lupus I, and we find that the dominant filament in the cloud is approximately perpendicular to the large-scale field, while secondary filaments appear to run parallel to the magnetic fields in their vicinities. This is similar to what is observed in Serpens South via near-IR polarimetry, and consistent with what is seen in MHD simulations by F. Nakamura and Z. Li.

We present 450 {\mu}m polarimetric observations of the M17 molecular cloud obtained with the SHARP polarimeter at the Caltech Submillimeter Observatory. Across the observed region, the magnetic field orientation is consistent with previous submillimeter and far-infrared polarization measurements. Our observations are centered on a region of the molecular cloud that has been compressed by stellar winds from a cluster of OB stars. We have compared these new data with previous 350 {\mu}m polarimetry and find an anti-correlation between the 450 to 350 {\mu}m polarization magnitude ratio and the ratio of 21 cm to 450 {\mu}m intensity. The polarization ratio is lower near the east end of the studied region where the cloud is exposed to stellar winds and radiation. At the west end of the region, the polarization ratio is higher. We interpret the varying polarization spectrum as evidence supporting the radiative alignment torque (RAT) model for grain alignment, implying higher alignment efficiency in the region that is exposed to a higher anisotropic radiation field.

We present 350 micron polarization observations of four low-mass cores containing Class 0 protostars: L483, L1157, L1448-IRS2, and Serp-FIR1. This is the second paper in a larger survey aimed at testing magnetically regulated models for core-collapse. One key prediction of these models is that the mean magnetic field in a core should be aligned with the symmetry axis (minor axis) of the flattened YSO inner envelope (aka pseudodisk). Furthermore, the field should exhibit a pinched or hour-glass shaped morphology as gravity drags the field inward towards the central protostar. We combine our results for the four cores with results for three similar cores that were published in the first paper from our survey. An analysis of the 350 micron polarization data for the seven cores yields evidence of a positive correlation between mean field direction and pseudodisk symmetry axis. Our rough estimate for the probability of obtaining by pure chance a correlation as strong as the one we found is about 5%. In addition, we combine together data for multiple cores to create a source-averaged magnetic field map having improved signal-to-noise ratio, and this map shows good agreement between mean field direction and pseudodisk axis (they are within 15 degrees). We also see hints of a magnetic pinch in the source-averaged map. We conclude that core-scale magnetic fields appear to be strong enough to guide gas infall, as predicted by the magnetically regulated models. Finally, we find evidence of a positive correlation between core magnetic field direction and bipolar outflow axis.

We present the first detection of polarization around the Class 0 low-mass protostar L1157-mm at two different wavelengths. We show polarimetric maps at large scales (10" resolution at 350 um) from the SHARC-II Polarimeter and at smaller scales (1.2"-4.5" at 1.3 mm) from the Combined Array for Research in Millimeter-wave Astronomy (CARMA). The observations are consistent with each other and show inferred magnetic field lines aligned with the outflow. The CARMA observations suggest a full hourglass magnetic field morphology centered about the core; this is only the second well-defined hourglass detected around a low-mass protostar to date. We apply two different methods to CARMA polarimetric observations to estimate the plane-of-sky magnetic field magnitude, finding values of 1.4 and 3.4 mG.

The definition of homology for single-component galaxies is clear, but for multi-component (luminous and dark matter) galaxies there is some ambiguity. We attempt to clarify the situation by carefully separating the different concepts of homology that have been used to date. We argue that the most useful definition is that a set of galaxies is homologous if they are the same in all respects up to a set of three dimensional scaling constants which may differ from one galaxy to the next. Noting that we are free to choose the dimensional constants, we find that a set of hydrodynamic simulated galaxy merger remnants is significantly closer to homologous when the dimensional length constant is taken to be the radius containing equal amounts of dark and baryonic matter rather than the usual observationally motivated choice of the baryonic half-mass radius. Once the correct dimensional scaling constants are used, the stellar velocity dispersion anisotropy is essentially the sole source of the variation in the kinematic structure of these simulated merger remnants. In order to facilitate the use of these scaling constants to analyse observed galaxies, we calculated the relationship between our preferred dimensional scaling constants and the typical observationally accessible quantities.

To facilitate the study of black hole fueling, star formation, and feedback in galaxies, we outline a method for treating the radial forces on interstellar gas due to absorption of photons by dust grains. The method gives the correct behavior in all of the relevant limits (dominated by the central point source; dominated by the distributed isotropic source; optically thin; optically thick to UV/optical; optically thick to IR) and reasonably interpolates between the limits when necessary. The method is explicitly energy conserving so that UV/optical photons that are absorbed are not lost, but are rather redistributed to the IR where they may scatter out of the galaxy. We implement the radiative transfer algorithm in a two-dimensional hydrodynamical code designed to study feedback processes in the context of early-type galaxies. We find that the dynamics and final state of simulations are measurably but only moderately affected by radiative forces on dust, even when assumptions about the dust-to-gas ratio are varied from zero to a value appropriate for the Milky Way. In simulations with high gas densities designed to mimic ULIRGs with a star formation rate of several hundred solar masses per year, dust makes a more substantial contribution to the dynamics and outcome of the simulation. We find that, despite the large opacity of dust to UV radiation, the momentum input to the flow from radiation very rarely exceeds L/c due to two factors: the low opacity of dust to the re-radiated IR and the tendency for dust to be destroyed by sputtering in hot gas environments. We also develop a simplification of our radiative transfer algorithm that respects the essential physics but is much easier to implement and requires a fraction of the computational cost.

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.

We extend the black hole (BH) feedback models of Ciotti, Ostriker, and Proga to two dimensions. In this paper, we focus on identifying the differences between the one-dimensional and two-dimensional hydrodynamical simulations. We examine a normal, isolated $L_*$ galaxy subject to the cooling flow instability of gas in the inner regions. Allowance is made for subsequent star formation, Type Ia and Type II supernovae, radiation pressure, and inflow to the central BH from mildly rotating galactic gas which is being replenished as a normal consequence of stellar evolution. The central BH accretes some of the infalling gas and expels a conical wind with mass, momentum, and energy flux derived from both observational and theoretical studies. The galaxy is assumed to have low specific angular momentum in analogy with the existing one-dimensional case in order to isolate the effect of dimensionality. The code then tracks the interaction of the outflowing radiation and winds with the galactic gas and their effects on regulating the accretion. After matching physical modeling to the extent possible between the one-dimensional and two-dimensional treatments, we find essentially similar results in terms of BH growth and duty cycle (fraction of the time above a given fraction of the Eddington luminosity). In the two-dimensional calculations, the cool shells forming at 0.1--1 kpc from the center are Rayleigh--Taylor unstable to fragmentation, leading to a somewhat higher accretion rate, less effective feedback, and a more irregular pattern of bursting compared to the one-dimensional case.

We investigate the polarization modulation properties of a variable-delay polarization modulator (VPM). The VPM modulates polarization via a variable separation between a polarizing grid and a parallel mirror. We find that in the limit where the wavelength is much larger than the diameter of the metal wires that comprise the grid, the phase delay derived from the geometric separation between the mirror and the grid is sufficient to characterize the device. However, outside of this range, additional parameters describing the polarizing grid geometry must be included to fully characterize the modulator response. In this paper, we report test results of a VPM at wavelengths of 350 microns and 3 mm. Electromagnetic simulations of wire grid polarizers were performed and are summarized using a simple circuit model that incorporates the loss and polarization properties of the device.

We report new 350 micron polarization observations of the thermal dust emission from the cores surrounding the low-mass, Class 0 YSOs L1527, IC348-SMM2 and B335. We have inferred magnetic field directions from these observations, and have used them together with results in the literature to determine whether magnetically regulated core-collapse and star-formation models are consistent with the observations. These models predict a pseudo-disk with its symmetry axis aligned with the core magnetic field. The models also predict a magnetic field pinch structure on a scale less than or comparable to the infall radii for these sources. In addition, if the core magnetic field aligns (or nearly aligns) the core rotation axis with the magnetic field before core collapse, then the models predict the alignment (or near alignment) of the overall pinch field structure with the bipolar outflows in these sources. We show that if one includes the distorting effects of bipolar outflows on magnetic fields, then in general the observational results for L1527 and IC348-SMM2 are consistent with these magnetically regulated models. We can say the same for B335 only if we assume the distorting effects of the bipolar outflow on the magnetic fields within the B335 core are much greater than for L1527 and IC348-SMM2. We show that the energy densities of the outflows in all three sources are large enough to distort the magnetic fields predicted by magnetically regulated models.

We present new 6 and 20 cm Very Large Array (VLA) observations of polarized continuum emission of roughly 0.5 square degrees of the Galactic center (GC) region. The 6 cm observations detect diffuse linearly-polarized emission throughout the region with a brightness of roughly 1 mJy per 15"x10" beam. The Faraday rotation measure (RM) toward this polarized emission has structure on degree size scales and ranges from roughly +330 rad/m2 east of the dynamical center (Sgr A) to -880 rad/m2 west of the dynamical center. This RM structure is also seen toward several nonthermal radio filaments, which implies that they have a similar magnetic field orientation and constrains models for their origin. Modeling shows that the RM and its change with Galactic longitude are best explained by the high electron density and strong magnetic field of the GC region. Considering the emissivity of the GC plasma shows that while the absolute RM values are indirect measures of the GC magnetic field, the RM longitude structure directly traces the magnetic field in the central kiloparsec of the Galaxy. Combining this result with previous work reveals a larger RM structure covering the central ~2 degrees of the Galaxy. This RM structure is similar to that proposed by Novak and coworkers, but is shifted roughly 50 pc west of the dynamical center of the Galaxy. If this RM structure originates in the GC region, it shows that the GC magnetic field is organized on ~300 pc size scales. The pattern is consistent with a predominantly poloidal field geometry, pointing from south to north, that is perturbed by the motion of gas in the Galactic disk.

The deposition of mechanical feedback from a supermassive black hole (SMBH) in an active galactic nucleus (AGN) into the surrounding galaxy occurs via broad-line winds which must carry mass and radial momentum as well as energy. The effect can be summarized by the dimensionless parameter $\eta=dot{M_outflow}/dot{M_accretion}= (2 \epsilon_w c^2)/v_w^2$ where ($\epslion_w \equiv dot{E}_w/(dot{M_accretion} c^2)$) is the efficiency by which accreted matter is turned into wind energy in the disc surrounding the central SMBH. The outflowing mass and omentum are proportional to $\eta$, and many prior treatments have essentially assumed that $\eta=0$. We perform one- and two-dimensional simulations and find that the growth of the central SMBH is very sensitive to the inclusion of the mass and momentum driving but is insensitive to the assumed mechanical efficiency. For example in representative calculations, the omission of momentum and mass feedback leads to an hundred fold increase in the mass of the SMBH to over $10^{10} \Msun$. When allowance is made for momentum driving, the final SMBH mass is much lower and the wind efficiencies which lead to the most observationally acceptable results are relatively low with $\epsilon_w \lesssim 10^{-4}$.

The Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry (BLAST-Pol) is a suborbital mapping experiment designed to study the role played by magnetic fields in the star formation process. BLAST-Pol is the reconstructed BLAST telescope, with the addition of linear polarization capability. Using a 1.8 m Cassegrain telescope, BLAST-Pol images the sky onto a focal plane that consists of 280 bolometric detectors in three arrays, observing simultaneously at 250, 350, and 500 um. The diffraction-limited optical system provides a resolution of 30'' at 250 um. The polarimeter consists of photolithographic polarizing grids mounted in front of each bolometer/detector array. A rotating 4 K achromatic half-wave plate provides additional polarization modulation. With its unprecedented mapping speed and resolution, BLAST-Pol will produce three-color polarization maps for a large number of molecular clouds. The instrument provides a much needed bridge in spatial coverage between larger-scale, coarse resolution surveys and narrow field of view, and high resolution observations of substructure within molecular cloud cores. The first science flight will be from McMurdo Station, Antarctica in December 2010.

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

We present the first 350 micron polarization measurement for the disk of the T Tauri star (TTS) DG Tau. The data were obtained using the SHARP polarimeter at the Caltech Submillimeter Observatory. We did not detect polarization but we provide constraints on the normalized Stokes parameters q and u. We derive information about the polarization spectrum by comparing our 350 micron result with an 850 micron polarization detection previously published for this source. This empirical information concerning the polarization spectrum disagrees with the predictions of a recent model for TTS disk polarization. We conclude, however, that adding more mass to the model disk would probably result in model polarization spectra that agree much better with the 350 and 850 micron polarimetric observations. We suggest that multiwavelength polarimetry of TTS disk emission may provide a promising method for probing the opacity of TTS disks.

One can estimate the characteristic magnetic field strength in GMCs by comparing submillimeter polarimetric observations of these sources with simulated polarization maps developed using a range of different values for the assumed field strength. The point of comparison is the degree of order in the distribution of polarization position angles. In a recent paper by H. Li and collaborators, such a comparison was carried out using SPARO observations of two GMCs, and employing simulations by E. Ostriker and collaborators. Here we reexamine this same question, using the same data set and the same simulations, but using an approach that differs in several respects. The most important difference is that we incorporate new, higher angular resolution observations for one of the clouds, obtained using the Hertz polarimeter. We conclude that the agreement between observations and simulations is best when the total magnetic energy (including both uniform and fluctuating field components) is at least as large as the turbulent kinetic energy.

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) is a sub-orbital experiment designed to study the process of star formation in local galaxies (including the Milky Way) and in galaxies at cosmological distances. Using a 2-m Cassegrain telescope, BLAST images the sky onto a focal plane, which consists of 270 bolometric detectors split between three arrays, observing simultaneously in 30% wide bands, centered at 250, 350, and 500 microns. The diffraction-limited optical system provides a resolution of 30" at 250 microns. The pointing system enables raster-like scans with a positional accuracy of ~30", reconstructed to better than 5" rms in post-flight analysis. BLAST had two successful flights, from the Arctic in 2005, and from Antarctica in 2006, which provided the first high-resolution and large-area (~0.8-200 deg^2) submillimeter surveys at these wavelengths. As a pathfinder for the SPIRE instrument on Herschel, BLAST shares with the ESA satellite similar focal plane technology and scientific motivation. A third flight in 2009 will see the instrument modified to be polarization-sensitive (BLAST-Pol). With its unprecedented mapping speed and resolution, BLAST-Pol will provide insights into Galactic star-forming nurseries, and give the necessary link between the larger, coarse resolution surveys and the narrow, resolved observations of star-forming structures from space and ground based instruments being commissioned in the next 5 years.

We characterize the problem of artificial polarization for the Submillimeter High Angular Resolution Polarimeter (SHARP) through the use of simulated data and observations made at the Caltech Submillimeter Observatory (CSO). These erroneous, artificial polarization signals are introduced into the data through misalignments in the bolometer sub-arrays plus pointing drifts present during the data-taking procedure. An algorithm is outlined here to address this problem and correct for it, provided that one can measure the degree of the sub-array misalignments and telescope pointing drifts. Tests involving simulated sources of Gaussian intensity profile indicate that the level of introduced artificial polarization is highly dependent upon the angular size of the source. Despite this, the correction algorithm is effective at removing up to 60% of the artificial polarization during these tests. The analysis of Jupiter data taken in January 2006 and February 2007 indicates a mean polarization of 1.44%+/-0.04% and 0.95%+/-0.09%, respectively. The application of the correction algorithm yields mean reductions in the polarization of approximately 0.15% and 0.03% for the 2006 and 2007 data sets, respectively.

We study the redshift evolution of galaxy pair fractions and merger rates for different types of galaxies using kinematic pairs selected from the DEEP2 Redshift Survey. By parameterizing the evolution of the pair fraction as (1+z)^{m}, we find that the companion rate increases mildly with redshift with m = 0.41+-0.20 for all galaxies with -21 < M_B^{e} < -19. Blue galaxies show slightly faster evolution in the blue companion rate with m = 1.27+-0.35 while red galaxies have had fewer red companions in the past as evidenced by the negative slope m = -0.92+-0.59. We find that at low redshift the pair fraction within the red sequence exceeds that of the blue cloud, indicating a higher merger probability among red galaxies compared to that among the blue galaxies. With further assumptions on the merger timescale and the fraction of pairs that will merge, the galaxy major merger rates for 0.1 < z <1.2 are estimated to be ~10^{-3}h^{3}Mpc^{-3}Gyr^{-1} with a factor of 2 uncertainty. At z ~ 1.1, 68% of mergers are wet, 8% of mergers are dry, and 24% of mergers are mixed, compared to 31% wet mergers, 25% dry mergers, and 44% mixed mergers at z ~ 0.1. The growth of dry merger rates with decreasing redshift is mainly due to the increase in the co-moving number density of red galaxies over time. About 22% to 54% of present-day L^{*} galaxies have experienced major mergers since z ~ 1.2, depending on the definition of major mergers. Moreover, 24% of the red galaxies at the present epoch have had dry mergers with luminosity ratios between 1:4 and 4:1 since z ~ 1. Our results also suggest that the wet mergers and/or mixed mergers may be partially responsible for producing red galaxies with intermediate masses while a significant portion of massive red galaxies is assembled through dry mergers at later times.

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.

We present first results of Hertz/VPM, the first submillimeter polarimeter employing the dual Variable-delay Polarization Modulator (dual-VPM). This device differs from previously used polarization modulators in that it operates in translation rather than mechanical rotation. We discuss the basic theory behind this device, and its potential advantages over the commonly used half wave plate (HWP). The dual-VPM was tested both at the Submillimeter Telescope Observatory (SMTO) and in the lab. In each case we present a detailed description of the setup. We discovered nonideal behavior in the system. This is at least in part due to properties of the VPM wire grids (diameter, spacing) employed in our experiment. Despite this, we found that the dual-VPM system is robust, operating with high efficiency and low instrumental polarization. This device is well suited for air and space-borne applications.

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 examine the radius-luminosity (R-L) relation for blue galaxies in the Team Keck Redshift Survey (TKRS) of GOODS-N. We compare with a volume-limited, Sloan Digital Sky Survey sample and find that the R-L relation has evolved to lower surface brightness since z=1. Based on the detection limits of GOODS this can not be explained by incompleteness in low surface-brightness galaxies. Number density arguments rule out a pure radius evolution. It can be explained by a radius dependent decline in B-band luminosity with time. Assuming a linear shift in M_B with z, we use a maximum likelihood method to quantify the evolution. Under these assumptions, large (R_{1/2} > 5 kpc), and intermediate sized (3 < R_{1/2} < 5 kpc) galaxies, have experienced Delta M_B =1.53 (-0.10,+0.13) and 1.65 (-0.18, +0.08) magnitudes of dimming since z=1. A simple exponential decline in star formation with an e-folding time of 3 Gyr can result in this amount of dimming. Meanwhile, small galaxies, or some subset thereof, have experienced more evolution, 2.55 (+/- 0.38) magnitudes. This factor of ten decline in luminosity can be explained by sub-samples of starbursting dwarf systems that fade rapidly, coupled with a decline in burst strength or frequency. Samples of bursting, luminous, blue, compact galaxies at intermediate redshifts have been identified by various previous studies. If there has been some growth in galaxy size with time, these measurements are upper limits on luminosity fading.

Using a sample of 89 snapshots from 58 hydrodynamic binary galaxy major merger simulations, we find that stellar remnants are mostly oblate while dark matter halos are mostly prolate or triaxial. The stellar minor axis and the halo major axis are almost always nearly perpendicular. This can be understood by considering the influence of angular momentum and dissipation during the merger. If binary mergers of spiral galaxies are responsible for the formation of elliptical galaxies or some subpopulation thereof, these galaxies can be expected to be oblate and inhabit their halos with the predicted shapes and orientations. These predictions are potentially relevant to observational studies of weak gravitational lensing, where one must stack many optically aligned galaxies in order to determine the shape of the resulting stacked mass distribution. The simple relationship between the dark and luminous matter presented here can be used to guide the stacking of galaxies to minimize the information lost.

We present results from the Austral Winter 2003 observing campaign of SPARO, a 450 micron polarimeter used with a two-meter telescope at South Pole. We mapped large-scale magnetic fields in four Giant Molecular Clouds (GMCs) in the Galactic disk: NGC 6334, the Carina Nebula, G333.6-0.2 and G331.5-0.1. We find a statistically significant correlation of the inferred field directions with the orientation of the Galactic plane. Specifically, three of the four GMCs (NGC 6334 is the exception) have mean field directions that are within 15 degrees of the plane. The simplest interpretation is that the field direction tends to be preserved during the process of GMC formation. We have also carried out an analysis of published optical polarimetry data. For the closest of the SPARO GMCs, NGC 6334, we can compare the field direction in the cloud as measured by SPARO with the field direction in a larger region surrounding the cloud, as determined from optical polarimetry. For purposes of comparison, we also use optical polarimetry to determine field directions for other regions of similar size and distance. Overall, the results from this optical polarimetry analysis are consistent with our suggestion that field direction tends to be preserved during GMC formation. Finally, we compare the disorder in our magnetic field maps with the disorder seen in magnetic field maps derived from MHD turbulence simulations. We conclude from these comparisons that the magnetic energy density in our clouds is comparable to the turbulent energy density.

We address the correlations of black hole (BH) mass with four different host-galaxy properties from eleven existing data sets. To guide theoretical understanding, we first try to quantify the tightness of the intrinsic correlations. Given the estimated measurement errors, we evaluate the probability distribution of the residual variance in excess of that expected from the measurement errors. Our central result is that the current data sets do not allow definite conclusions regarding the quality of the true correlations because the obtained probability distributions for the residual variance overlap for most quantities. We then consider which of the relations offer the best inferences of BH mass when there is no direct measurement available. As with the residual variances, we find that the probability distribution of expected uncertainty in inferred BH masses overlaps significantly for most of the relations. Photometric methods would then be preferred because the data are easier to obtain, as long as bulge-disk decomposition or detailed modeling of the photometric profile (as in \citet{graham:01}) do not present problems. Determining which correlation offers the best inferences requires reducing the uncertainty in the expected error in the inferred BH masses (the ``error on the error''). This uncertainty is currently limited by uncertainty in the residual variance for all of the relations.

There is strong evidence that the mass in the Universe is dominated by dark matter, which exerts gravitational attraction but whose exact nature is unknown. In particular, all galaxies are believed to be embedded in massive haloes of dark matter. This view has recently been challenged by surprisingly low random stellar velocities in the outskirts of ordinary elliptical galaxies, which were interpreted as indicating a lack of dark matter (Mendez et al. 2001; Romanowsky et al. 2003). Here we show that the low velocities are in fact compatible with galaxy formation in dark-matter haloes. Using numerical simulations of disc-galaxy mergers, we find that the stellar orbits in the outer regions of the resulting ellipticals are very elongated. These stars were torn by tidal forces from their original galaxies during the first close passage and put on outgoing trajectories. The elongated orbits, combined with the steeply falling density profile of the observed tracers, explain the observed low velocities even in the presence of large amounts of dark matter. Projection effects when viewing a triaxial elliptical can lead to even lower observed velocities along certain lines of sight.

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.

Recent studies of the 55 Cancri system suggest the existence of three planets with periods of ~15 days, ~45 days, and ~5500 days (Marcy et al. 2002). The inner two planets are near the 3:1 mean motion commensurability and it is likely that these two planets became trapped in the resonance while farther from the star and migrated together. As the innermost planet begins to dissipate energy through tides the planets break out of the resonance. The final state of the system gives important information about its past history, such as the migration timescale that led to capture.

We have observed the linear polarization of 450 micron continuum emission from the Galactic center, using a new polarimetric detector system that is operated on a 2 m telescope at the South Pole. The resulting polarization map extends ~ 170 pc along the Galactic plane and ~ 30 pc in Galactic latitude, and thus covers a significant fraction of the central molecular zone. Our map shows that this region is permeated by large-scale toroidal magnetic fields. We consider our results together with radio observations that show evidence for poloidal fields in the Galactic center, and with Faraday rotation observations. We compare all of these observations with the predictions of a magnetodynamic model for the Galactic center that was proposed in order to explain the Galactic Center Radio Lobe as a magnetically driven gas outflow. We conclude that the observations are basically consistent with the model.

Observations of the microwave sky using the Python telescope in its fifth season of operation at the Amundsen-Scott South Pole Station in Antarctica are presented. The system consists of a 0.75 m off-axis telescope instrumented with a HEMT amplifier-based radiometer having continuum sensitivity from 37-45 GHz in two frequency bands. With a 0.91 deg x 1.02 deg beam the instrument fully sampled 598 deg^2 of sky, including fields measured during the previous four seasons of Python observations. Interpreting the observed fluctuations as anisotropy in the cosmic microwave background, we place constraints on the angular power spectrum of fluctuations in eight multipole bands up to l ~ 260. The observed spectrum is consistent with both the COBE experiment and previous Python results. There is no significant contamination from known foregrounds. The results show a discernible rise in the angular power spectrum from large (l ~ 40) to small (l ~ 200) angular scales. The shape of the observed power spectrum is not a simple linear rise but has a sharply increasing slope starting at l ~ 150.

We report on additional observations of degree scale anisotropy at 90~GHz from the Amundsen-Scott South Pole Station in Antarctica. Observations during the first season with the Python instrument yielded a statistically significant sky signal; in this paper we report the confirmation of that signal with data taken in the second year, and on results from an interleaving set of fields.