We present observations of the H-alpha, H-beta, [SII] 6716, 6731 and [NII] 6583 emission lines in the galactic HII region Sh2-235 with the Mapper of Narrow Galaxy Lines (MaNGaL), a tunable filter at the 1-m telescope of Special Astrophysical Observatory of the Russian Academy of Sciences. We show that the HII region is obscured by neutral material with AV = 2-4 mag. The area with the highest AV is situated to the south-west from the ionizing star and coincides with a maximum detected electron density of >=300 cm(-3). The combination of these results with archive AKARI far-infrared data allows us to estimate the contribution of the front and rear walls to the total column density of neutral material in S235, and explain the three-dimensional structure of the region. The HII region consist of a denser, more compact portion deeply embedded in the neutral medium and the less dense and obscured gas. The front and rear walls of the HII region are inhomogeneous, with the material in the rear wall having a higher column density. We find a two-sided photodissociation region in the dense clump S235 East 1, illuminated by a UV field with G0=50-70 and 200 Habing units in the western and eastern parts, respectively.

Future mm-wave and sub-mm space missions will employ large arrays of multiplexed Transition Edge Sensor (TES) bolometers. Such instruments must contend with the high flux of cosmic rays beyond our atmosphere that induce "glitches" in bolometer data, which posed a challenge to data analysis from the Planck bolometers. Future instruments will face the additional challenges of shared substrate wafers and multiplexed readout wiring. In this work we explore the susceptibility of modern TES arrays to the cosmic ray environment of space using two data sets: the 2015 long-duration balloon flight of the SPIDER cosmic microwave background polarimeter, and a laboratory exposure of SPIDER flight hardware to radioactive sources. We find manageable glitch rates and short glitch durations, leading to minimal effect on SPIDER analysis. We constrain energy propagation within the substrate through a study of multi-detector coincidences, and give a preliminary look at pulse shapes in laboratory data.

Families of cosmic inflation models predict a primordial gravitational-wave background that imprints B-mode polarization pattern in the Cosmic Microwave Background (CMB). High sensitivity instruments with wide frequency coverage and well-controlled systematic errors are needed to constrain the faint B-mode amplitude. We have developed antenna-coupled Transition Edge Sensor (TES) arrays for high-sensitivity polarized CMB observations over a wide range of millimeter-wave bands. BICEP Array, the latest phase of the BICEP/Keck experiment series, is a multi-receiver experiment designed to search for inflationary B-mode polarization to a precision $\sigma$(r) between 0.002 and 0.004 after 3 full years of observations, depending on foreground complexity and the degree of lensing removal. We describe the electromagnetic design and measured performance of BICEP Array low-frequency 40-GHz detector, their packaging in focal plane modules, and optical characterization including efficiency and beam matching between polarization pairs. We summarize the design and simulated optical performance, including an approach to improve the optical efficiency due to mismatch losses. We report the measured beam maps for a new broad-band corrugation design to minimize beam differential ellipticity between polarization pairs caused by interactions with the module housing frame, which helps minimize polarized beam mismatch that converts CMB temperature to polarization ($T \rightarrow P$) anisotropy in CMB maps.

Branches of cosmic inflationary models, such as slow-roll inflation, predict a background of primordial gravitational waves that imprints a unique odd-parity B-mode pattern in the Cosmic Microwave Background (CMB) at amplitudes that are within experimental reach. The BICEP/Keck (BK) experiment targets this primordial signature, the amplitude of which is parameterized by the tensor-to-scalar ratio r, by observing the polarized microwave sky through the exceptionally clean and stable atmosphere at the South Pole. B-mode measurements require an instrument with exquisite sensitivity, tight control of systematics, and wide frequency coverage to disentangle the primordial signal from the Galactic foregrounds. BICEP Array represents the most recent stage of the BK program, and comprises four BICEP3-class receivers observing at 30/40, 95, 150 and 220/270 GHz. The 30/40 GHz receiver will be deployed at the South Pole during the 2019/2020 austral summer. After 3 full years of observations with 30,000+ detectors, BICEP Array will measure primordial gravitational waves to a precision $\sigma(r)$ between 0.002 and 0.004, depending on foreground complexity and the degree of lensing removal. In this paper we give an overview of the instrument, highlighting the design features in terms of cryogenics, magnetic shielding, detectors and readout architecture as well as reporting on the integration and tests that are ongoing with the first receiver at 30/40 GHz.

Luminosity outbursts of the FUOri type stars, which have a magnitude of $\sim100 L_{\odot}$ and last for decades, may affect chemical composition of the surrounding protoplanetary disk. Using astrochemical modeling we analyze the changes induced by the outburst and search for species sensitive to the luminosity rise. Some changes in the disk molecular composition appear not only during the outburst itself but can also retain for decades after the end of the outburst. We analyze main chemical processes responsible for these effects and assess timescales at which chemically inert species return to the pre-outburst abundances.

The BICEP/Keck experiment (BK) is a series of small-aperture refracting telescopes observing degree-scale Cosmic Microwave Background (CMB) polarization from the South Pole in search of a primordial $B$-mode signature. This $B$-mode signal arises from primordial gravitational waves interacting with the CMB, and has amplitude parametrized by the tensor-to-scalar ratio $r$. Since 2016, BICEP3 and the Keck Array have been observing with 4800 total antenna-coupled transition-edge sensor detectors, with frequency bands spanning 95, 150, 220, and 270 GHz. Here we present the optical performance of these receivers from 2016 to 2019, including far-field beams measured in situ with an improved chopped thermal source and instrument spectral response measured with a field-deployable Fourier Transform Spectrometer. As a pair differencing experiment, an important systematic that must be controlled is the differential beam response between the co-located, orthogonally polarized detectors. We generate per-detector far-field beam maps and the corresponding differential beam mismatch that is used to estimate the temperature-to-polarization leakage in our CMB maps and to give feedback on detector and optics fabrication. The differential beam parameters presented here were estimated using improved low-level beam map analysis techniques, including efficient removal of non-Gaussian noise as well as improved spatial masking. These techniques help minimize systematic uncertainty in the beam analysis, with the goal of constraining the bias on $r$ induced by temperature-to-polarization leakage to be subdominant to the statistical uncertainty. This is essential as we progress to higher detector counts in the next generation of CMB experiments.

The BICEP/Keck (BK) experiment aims to detect the imprint of primordial gravitational waves in the Cosmic Microwave Background polarization, which would be direct evidence of the inflation theory. While the tensor-to-scalar ratio has been constrained to be r_0.05 < 0.06 at 95% c.l., further improvements on this upper limit are hindered by polarized Galactic foreground emissions and removal of gravitational lensing polarization. The 30/40 GHz receiver of the BICEP Array (BA) will deploy at the end of 2019 and will constrain the synchrotron foreground with unprecedented accuracy within the BK sky patch. We will show the design of the 30/40 GHz detectors and test results summarizing its performance. The low optical and atmospheric loading at these frequencies requires our TES detectors to have low saturation power in order to be photon-noise dominated. To realize the low thermal conductivity required from a 250 mK base temperature, we developed new bolometer leg designs. We will present the relevant measured detector parameters: G, Tc, Rn, Psat , and spectral bands, and noise spectra. We achieved a per bolometer NEP including all noise components of 2.07E-17 W/sqrt(Hz), including an anticipated photon noise level 1.54E-17 W/sqrt(Hz).

We describe an on-sky demonstration of a microwave-multiplexing readout system in one of the receivers of the Keck Array, a polarimetry experiment observing the cosmic microwave background at the South Pole. During the austral summer of 2018-2019, we replaced the time-division multiplexing readout system with microwave-multiplexing components including superconducting microwave resonators coupled to radio-frequency superconducting quantum interference devices at the sub-Kelvin focal plane, coaxial-cable plumbing and amplification between room temperature and the cold stages, and a SLAC Microresonator Radio Frequency system for the warm electronics. In the range 5-6 GHz, a single coaxial cable reads out 528 channels. The readout system is coupled to transition-edge sensors, which are in turn coupled to 150-GHz slot-dipole phased-array antennas. Observations began in April 2019, and we report here on an initial characterization of the system performance.

Precision measurements of cosmic microwave background (CMB) polarization require extreme control of instrumental systematics. In a companion paper we have presented cosmological constraints from observations with the BICEP2 and Keck Array experiments up to and including the 2015 observing season (BK15), resulting in the deepest CMB polarization maps to date and a statistical sensitivity to the tensor-to-scalar ratio of $\sigma(r) = 0.020$. In this work we characterize the beams and constrain potential systematic contamination from main beam shape mismatch at the three BK15 frequencies (95, 150, and 220 GHz). Far-field maps of 7,360 distinct beam patterns taken from 2010-2015 are used to measure differential beam parameters and predict the contribution of temperature-to-polarization leakage to the BK15 B-mode maps. In the multifrequency, multicomponent likelihood analysis that uses BK15, Planck, and WMAP maps to separate sky components, we find that adding this predicted leakage to simulations induces a bias of $\Delta r = 0.0027 \pm 0.0019$. Future results using higher-quality beam maps and improved techniques to detect such leakage in CMB data will substantially reduce this uncertainty, enabling the levels of systematics control needed for BICEP Array and other experiments that plan to definitively probe large-field inflation.

The H II region RCW120 is a well-known object, which is often considered as a target to verify theoretical models of gas and dust dynamics in the interstellar medium. However, the exact geometry of RCW120 is still a matter of debate. In this work, we analyse observational data on molecular emission in RCW120 and show that 13CO(2-1) and C18O(2-1) lines are fitted by a 2D model representing a ring-like face-on structure. The changing of the C18O(3-2) line profile from double-peaked to single-peaked from the dense molecular Condensation 1 might be a signature of stalled expansion in this direction. In order to explain a self-absorption dip of the 13CO(2-1) and 13CO(3-2) lines, we suggest that RCW120 is surrounded by a diffuse molecular cloud, and find confirmation of this cloud on a map of interstellar extinction. Optically thick 13CO(2-1) emission and the infrared 8 um PAH band form a neutral envelope of the H II region resembling a ring, while the envelope breaks into separate clumps on images made with optically thin C18O(2-1) line and far-infrared dust emission.

The work is devoted to the adaptation of the results of laboratory studies of the laser-induced dissociation of molecules of benzene adsorbed on a quartz substrate to the conditions of the interstellar medium. Adsorption was performed under conditions of low temperature and deep vacuum. The difference between the photolysis of adsorbed molecules and molecules in the gas phase is identified. Significance of process of photolytic desorption in the interstellar conditions is analyzed, in particular, in the conditions of photodissociation regions. It is shown that the efficiency and dissociation channels of photolysis of adsorbed and gas phase benzene differ substantially. It is concluded that the photolysis of aromatic hydrocarbons adsorbed on the interstellar dust grains contributes a negligible fraction to the abundance of small hydrocarbons in the interstellar medium.

We analyze parameters of the interstellar matter emission in star-forming complexes in the high metallicity galaxies NGC~628, NGC~2976, and NGC~3351, which have different morphological types. The relation between H$\alpha$ emission and emission in CO and HI lines is considered along with the relation between H$\alpha$ emission and dust emission in the infrared range (IR). The fluxes and surface brightnesses in the UV and IR correlate well with H$\alpha$ emission. The HI emission also correlates well with H$\alpha$, while the correlation between the CO and H$\alpha$ emission is much less prominent. The ratio of the fluxes at 8 and 24 $\mu$m decreases with increasing H$\alpha$ flux. This may be due to changes in the properties of the dust ensemble (a decrease in the mass fraction of polycyclic aromatic hydrocarbons) or to changes in excitation conditions. Analysis of the kinematics of the CO lines shows that the CO flux grows with increasing velocity scatter $\Delta V$ when $\Delta V\lesssim70$~km/s. Preliminary evidence for the existence of star-forming complexes with higher values of $\Delta V$ is presented, and the increase in the velocity scatter is accompanied by a decrease in the CO line luminosity of the complex.

We present a numerical tool Shiva designed to simulate the dust destruction in warm neutral, warm ionized, and hot ionized media under the influence of photo-processing, sputtering, and shattering. The tool is designed primarily to study the evolution of hydrogenated amorphous carbons (HACs), but options to simulate polycyclic aromatic hydrocarbons (PAHs), silicate and graphite grains are also implemented. HAC grain photo-processing includes both dehydrogenation and carbon atom loss. Dehydrogenation leads to material transformation from aliphatic to aromatic structure. Simultaneously, some other physical properties (band gap energy, optical properties, etc.) of the material change as well. The Shiva tool allows calculating the time-dependent evolution of the dust size distribution depending on hydrogen, helium, and carbon number densities and ionization state, gas temperature, radiation flux, relative gas-dust and grain-grain velocities. For HAC grains the evolution of band gap energy distribution is also computed. We describe a dust evolution model, on which the tool relies, and present evolutionary time-scales for dust grains of different sizes depending on external conditions. This allows a user to estimate quickly a lifetime of a specific dust grain under relevant conditions. As an example of the tool usage, we demonstrate how grain properties and corresponding infrared spectra evolve in photo-dissociation regions, H II regions, and supernova remnant shocks.

We present the strongest constraints to date on anisotropies of cosmic microwave background (CMB) polarization rotation derived from 150 GHz data taken by the BICEP2/Keck Array CMB experiments up to and including the 2014 observing season (BK14). The definition of the polarization angle in BK14 maps has gone through self-calibration in which the overall angle is adjusted to minimize the observed TB and EB power spectra. After this procedure, the QU maps lose sensitivity to a uniform polarization rotation but are still sensitive to anisotropies of polarization rotation. This analysis places constraints on the anisotropies of polarization rotation, which could be generated by CMB photons interacting with axionlike pseudoscalar fields or Faraday rotation induced by primordial magnetic fields. The sensitivity of BK14 maps ($\sim 3\mu$K-arcmin) makes it possible to reconstruct anisotropies of the polarization rotation angle and measure their angular power spectrum much more precisely than previous attempts. Our data are found to be consistent with no polarization rotation anisotropies, improving the upper bound on the amplitude of the rotation angle spectrum by roughly an order of magnitude compared to the previous best constraints. Our results lead to an order of magnitude better constraint on the coupling constant of the Chern-Simons electromagnetic term $g_{a\gamma}\leq 7.2\times 10^{-2}/H_I$ (95% confidence) than the constraint derived from the B-mode spectrum, where $H_I$ is the inflationary Hubble scale. This constraint leads to a limit on the decay constant of $10^{-6}\lesssim f_a/M_{\rm pl}$ at mass range of $10^{-33}< m_a< 10^{-28}$ eV for $r=0.01$, assuming $g_{a\gamma}\sim\alpha/(2\pi f_a)$ with $\alpha$ denoting the fine structure constant. The upper bound on the amplitude of the primordial magnetic fields is 30nG (95% confidence) from the polarization rotation anisotropies.

Dust is known to drift and grow in protoplanetary discs, which results in dust segregation over the disc extent. Maps of the spectral index $\alpha$ are a common tool for studying the dust content in protoplanetary discs. The analysis of observationally derived maps reveals significant gradients of the spectral index, confirming that dust evolves in the disc, but a more detailed information about the dust redistribution is required to make inferences about the early stages of dust growth. We calculated the spectral index maps based on the results of numerical hydrodynamical simulations using the FEOSAD code, which allows studying a long-term dynamics of a self-gravitating viscous disc populated with coagulating, drifting, and fragmenting dust. Here we demonstrate that values of the spectral index estimated for different wavelength intervals within the far-infrared and radio bands reveal the presence of dust grains of various sizes. Specifically, we show that the disc regions with the maximal spectral index in a specific wavelength interval are the regions with the prevalence of dust grains of a specific size. Thus, a set of spectral index maps derived using different wavelength intervals can be used to recover the dust size-distribution over the disc extent.

High-resolution ALMA images towards the Orion Bar show no discernible offset between the peak of H2 emission in the photodissociation region (PDR) and the CO(3-2) and HCO+(4-3) emission in the molecular region. This implies that positions of H2 and CO dissociation fronts are indistinguishable in the limit of ALMA resolution. We use the chemo-dynamical model MARION to show that the ALMA view of the Orion Bar, namely, no appreciable offset between the CO(3-2) and HCO+(4-3) peaks, merged H2 and CO dissociation fronts, and high intensity of HCO+(4-3) emission, can only be explained by the ongoing propagation of the dissociation fronts through the molecular cloud, coupled to the dust motion driven by the stellar radiation pressure, and are not reproduced in the model where the dissociation fronts are assumed to be stationary. Modelling line intensities, we demonstrate that after the fronts have merged, the angular separation of the CO(3-2) and HCO+(4-3) peaks is indeed unresolvable with the ALMA observations. Our model predictions are consistent with the results of the ALMA observations about the relation of the bright HCO+(4-3) emission to the compressed gas at the border of the PDR in the sense that the theoretical HCO+(4-3) peak does correspond to the gas density enhancement, which naturally appears in the dynamical simulation, and is located near the H2 dissociation front at the illuminated side of the CO dissociation front.

We present an observational study of small-scale feedback processes operating in the star-forming region located in the wall of the expanding supershell around the Cyg OB1 association. The interstellar gas and dust content and pre-stellar populations in the vicinity of the open star cluster vdB 130 are analysed based on new optical and IR observations performed with the 6-m (3D spectroscopic mapping in the [SII] doublet) and 2.5-m (optical and NIR images) Russian telescopes along with the archival data of Spitzer and Herschel space telescopes. Analysing ionized gas kinematics and emission spectra, we discovered a compact region with supersonic motions. These motions may be caused either by stellar wind, or a bipolar outflow from a protostellar disc. Young stellar objects were identified and classified in the area under study. Two star-forming regions were identified. One of them is a region of ongoing star formation in the head of the molecular cloud observed there and another one is a burst of star formation in the cloud tail.

The chemical influence of luminosity outbursts on the environments of young solar-type stars is explored. Species are categorised into several types according to their response to the outburst. The first and second types imply chemical changes only during the outburst (with slightly different behaviours). These response types are mostly observed close to the star and are caused by icy mantle evaporation. However, mantles recover after the outburst almost immediately. A notable exception is benzene ice, which is accumulated on dust surfaces during and after the outburst, so that its abundance exceeds the pre-outburst level by orders of magnitude. The third type of response is mostly seen at the disc periphery and implies alteration of abundances during the outburst and preservation of these `abnormal' abundances for centuries. This behaviour is typical of organic compounds, like HCOOCH$_3$, CH$_3$CN, CH$_2$CO. Their presence in the dark disc regions can be a manifestation of the past outburst. CO and CO$_2$ only trace past outbursts at the remote disc regions. The outburst changes the C/O ratio, but it quickly returns to the pre-outburst value almost everywhere in the disc. An important factor determining the sensitivity of molecular composition to the outburst is the dust size distribution. The duration of the pre-outburst stage and of the outburst itself influence the chemical effects, if the burst duration is shorter than 50 yr and the duration of the quiescent phase between the bursts is shorter than 100 kyr.

The BICEP and Keck Array experiments are a suite of small-aperture refracting telescopes observing the microwave sky from the South Pole. They target the degree-scale B-mode polarization signal imprinted in the Cosmic Microwave Background (CMB) by primordial gravitational waves. Such a measurement would shed light on the physics of the very early universe. While BICEP2 observed for the first time a B-mode signal at 150 GHz, higher frequencies from the Planck satellite showed that it could be entirely due to the polarized emission from Galactic dust, though uncertainty remained high. Keck Array has been observing the same region of the sky for several years, with an increased detector count, producing the deepest polarized CMB maps to date. New detectors at 95 GHz were installed in 2014, and at 220 GHz in 2015. These observations enable a better constraint of galactic foreground emissions, as presented here. In 2015, BICEP2 was replaced by BICEP3, a 10 times higher throughput telescope observing at 95 GHz, while Keck Array is now focusing on higher frequencies. In the near future, BICEP Array will replace Keck Array, and will allow unprecedented sensitivity to the gravitational wave signal. High resolution observations from the South Pole Telescope (SPT) will also be used to remove the lensing contribution to B-modes.

We present results from an analysis of all data taken by the BICEP2/Keck CMB polarization experiments up to and including the 2015 observing season. This includes the first Keck Array observations at 220 GHz and additional observations at 95 & 150 GHz. The $Q/U$ maps reach depths of 5.2, 2.9 and 26 $\mu$K$_{cmb}$ arcmin at 95, 150 and 220 GHz respectively over an effective area of $\approx 400$ square degrees. The 220 GHz maps achieve a signal-to-noise on polarized dust emission approximately equal to that of Planck at 353 GHz. We take auto- and cross-spectra between these maps and publicly available WMAP and Planck maps at frequencies from 23 to 353 GHz. We evaluate the joint likelihood of the spectra versus a multicomponent model of lensed-$\Lambda$CDM+$r$+dust+synchrotron+noise. The foreground model has seven parameters, and we impose priors on some of these using external information from Planck and WMAP derived from larger regions of sky. The model is shown to be an adequate description of the data at the current noise levels. The likelihood analysis yields the constraint $r_{0.05}<0.07$ at 95% confidence, which tightens to $r_{0.05}<0.06$ in conjunction with Planck temperature measurements and other data. The lensing signal is detected at $8.8 \sigma$ significance. Running maximum likelihood search on simulations we obtain unbiased results and find that $\sigma(r)=0.020$. These are the strongest constraints to date on primordial gravitational waves.

The FU Ori-type young stellar objects are characterized by a sudden increase in luminosity by 1$-$2 orders of magnitude, followed by slow return to the pre-outburst state on timescales of $\sim$10$-$100 yr. The outburst strongly affects the entire disk, changing its thermal structure and radiation field. In this paper, using a detailed physical-chemical model we study the impact of the FU Ori outburst on the disk chemical inventory. Our main goal is to identify gas-phase molecular tracers of the outburst activity that could be observed after the outburst with modern telescopes such as ALMA and NOEMA. We find that the majority of molecules experience a considerable increase in the total disk gas-phase abundances due to the outburst, mainly due to the sublimation of their ices. Their return to the pre-outburst chemical state takes different amounts of time, from nearly instantaneous to very long. Among the former ones we identify CO, NH$_3$, C$_2$H$_6$, C$_3$H$_4$, etc. Their abundance evolution tightly follows the luminosity curve. For CO the abundance increase does not exceed an order of magnitude, while for other tracers the abundances increase by 2$-$5 orders of magnitude. Other molecules like H$_2$CO and NH$_2$OH have longer retention timescales, remaining in the gas phase for $\sim10-10^3$ yr after the end of the outburst. Thus H$_2$CO could be used as an indicator of the previous outbursts in the post-outburst FU Ori systems. We investigate the corresponding time-dependent chemistry in detail and present the most favorable transitions and ALMA configurations for future observations.

Bicep Array is a cosmic microwave background (CMB) polarization experiment that will begin observing at the South Pole in early 2019. This experiment replaces the five Bicep2 style receivers that compose the Keck Array with four larger Bicep3 style receivers observing at six frequencies from 30 to 270GHz. The 95GHz and 150GHz receivers will continue to push the already deep Bicep/Keck CMB maps while the 30/40GHz and 220/270GHz receivers will constrain the synchrotron and galactic dust foregrounds respectively. Here we report on the design and performance of the Bicep Array instruments focusing on the mount and cryostat systems.

BICEP3 is a 520mm aperture on-axis refracting telescope observing the polarization of the cosmic microwave background (CMB) at 95GHz in search of the B-mode signal originating from inflationary gravitational waves. BICEP3's focal plane is populated with modularized tiles of antenna-coupled transition edge sensor (TES) bolometers. BICEP3 was deployed to the South Pole during 2014-15 austral summer and has been operational since. During the 2016-17 austral summer, we implemented changes to optical elements that lead to better noise performance. We discuss this upgrade and show the performance of BICEP3 at its full mapping speed from the 2017 and 2018 observing seasons. BICEP3 achieves an order-of-magnitude improvement in mapping speed compared to a Keck 95GHz receiver. We demonstrate $6.6\mu K\sqrt{s}$ noise performance of the BICEP3 receiver.

Targeting faint polarization patterns arising from Primordial Gravitational Waves in the Cosmic Microwave Background requires excellent observational sensitivity. Optical elements in small aperture experiments such as Bicep3 and Keck Array are designed to optimize throughput and minimize losses from transmission, reflection and scattering at millimeter wavelengths. As aperture size increases, cryostat vacuum windows must withstand larger forces from atmospheric pressure and the solution has often led to a thicker window at the expense of larger transmission loss. We have identified a new candidate material for the fabrication of vacuum windows: with a tensile strength two orders of magnitude larger than previously used materials, woven high-modulus polyethylene could allow for dramatically thinner windows, and therefore significantly reduced losses and higher sensitivity. In these proceedings we investigate the suitability of high-modulus polyethylene windows for ground-based CMB experiments, such as current and future receivers in the Bicep/Keck Array program. This includes characterizing their optical transmission as well as their mechanical behavior under atmospheric pressure. We find that such ultra-thin materials are promising candidates to improve the performance of large-aperture instruments at millimeter wavelengths, and outline a plan for further tests ahead of a possible upcoming field deployment of such a science-grade window.

BICEP Array is the newest multi-frequency instrument in the BICEP/Keck Array program. It is comprised of four 550 mm aperture refractive telescopes observing the polarization of the cosmic microwave background (CMB) at 30/40, 95, 150 and 220/270 GHz with over 30,000 detectors. We present an overview of the receiver, detailing the optics, thermal, mechanical, and magnetic shielding design. BICEP Array follows BICEP3's modular focal plane concept, and upgrades to 6" wafer to reduce fabrication with higher detector count per module. The first receiver at 30/40 GHz is expected to start observing at the South Pole during the 2019-20 season. By the end of the planned BICEP Array program, we project $\sigma(r) \sim 0.003$, assuming current modeling of polarized Galactic foreground and depending on the level of delensing that can be achieved with higher resolution maps from the South Pole Telescope.

BICEP Array is a degree-scale Cosmic Microwave Background (CMB) experiment that will search for primordial B-mode polarization while constraining Galactic foregrounds. BICEP Array will be comprised of four receivers to cover a broad frequency range with channels at 30/40, 95, 150 and 220/270 GHz. The first low-frequency receiver will map synchrotron emission at 30 and 40 GHz and will deploy to the South Pole at the end of 2019. In this paper, we give an overview of the BICEP Array science and instrument, with a focus on the detector module. We designed corrugations in the metal frame of the module to suppress unwanted interactions with the antenna-coupled detectors that would otherwise deform the beams of edge pixels. This design reduces the residual beam systematics and temperature-to-polarization leakage due to beam steering and shape mismatch between polarized beam pairs. We report on the simulated performance of single- and wide-band corrugations designed to minimize these effects. Our optimized design alleviates beam differential ellipticity caused by the metal frame to about 7% over 57% bandwidth (25 to 45 GHz), which is close to the level due the bare antenna itself without a metal frame. Initial laboratory measurements are also presented.

The Canadian Hydrogen Intensity Mapping Experiment (CHIME) is a novel transit radio telescope operating across the 400-800-MHz band. CHIME is comprised of four 20-m x 100-m semi-cylindrical paraboloid reflectors, each of which has 256 dual-polarization feeds suspended along its axis, giving it a >200 square degree field-of-view. This, combined with wide bandwidth, high sensitivity, and a powerful correlator makes CHIME an excellent instrument for the detection of Fast Radio Bursts (FRBs). The CHIME Fast Radio Burst Project (CHIME/FRB) will search beam-formed, high time-and frequency-resolution data in real time for FRBs in the CHIME field-of-view. Here we describe the CHIME/FRB backend, including the real-time FRB search and detection software pipeline as well as the planned offline analyses. We estimate a CHIME/FRB detection rate of 2-42 FRBs/sky/day normalizing to the rate estimated at 1.4-GHz by Vander Wiel et al. (2016). Likely science outcomes of CHIME/FRB are also discussed. CHIME/FRB is currently operational in a commissioning phase, with science operations expected to commence in the latter half of 2018.

A search for infrared ring nebulae associated with regions of ionized hydrogen has been carried out. The New GPS Very Large Array survey at 20 cm forms the basis of the search, together with observations obtained with the Spitzer Space Telescope at 8 and 24 $\mu$m and the Herschel Space Telescope at 70 $\mu$m. Objects having ring-like morphologies at 8 $\mu$m and displaying extended emission at 20 cm were selected visually. Emission at 24 $\mu$m having the form of an inner ring or central peak is also observed in the selected objects. A catalog of 99 ring nebulae whose shapes at 8 and 70 $\mu$m are well approximated by ellipses has been compiled. The catalog contains 32 objects whose shapes are close to circular (eccentricities of the fitted ellipses at 8 $\mu$m no greater than 0.6, angular radius exceeding 20). These objects are promising for comparisons with the results of one-dimensional hydrodynamical simulations of expanding regions of ionized hydrogen.

Total infrared fluxes are estimated for 99 HII regions around massive stars. The following wavebands have been used for the analysis: 8 and 24 $\mu$m, based on data from Spitzer space telescope (IRAC and MIPS, respectively); 70, 160, 250, 350, and 500 $\mu$m, based on data from Herschel Space Observatory (PACS and SPIRE). The estimated fluxes are used to evaluate the mass fraction of polycyclic aromatic hydrocarbons ($q_{\rm PAH}$) and the intensity of the ultraviolet emission in the studied objects. It is shown that the PAH mass fraction, $q_{\rm PAH}$, is much lower in these objects than the average Galactic value, implying effective destruction of aromatic particles in HII regions. Estimated radiation field intensities ($U$) are close to those derived for extragalactic HII complexes. Color indices [$F_{24}$/$F_{8}$], [$F_{70}$/$F_{24}$], [$F_{160}$/$F_{24}$], [$F_{160}$/$F_{70}$] are compared to criteria proposed to distinguish between regions of ionized hydrogen and planetary nebulae. Also, we relate our results to analogous color indices for extragalactic complexes of ionized hydrogen.

The relationships between atomic and molecular hydrogen and dust of various sizes in extragalactic star-forming regions are considered, based on observational data from the Spitzer and Herschel infrared space telescopes, the Very Large Array (atomic hydrogen emission) and IRAM (CO emission). The source sample consists of approximately 300 star-forming regions in 11 nearby galaxies. Aperture photometry has been applied to measure the fluxes in eight infrared bands (3.6, 4.5, 5.8, 8, 24, 70, 100, and 160$\mu$m), the atomic hydrogen (21cm) line and CO (2--1) lines. The parameters of the dust in the starforming regions were determined via synthetic-spectra fitting, such as the total dust mass, the fraction of polycyclic aromatic hydrocarbons (PAHs), etc. Comparison of the observed fluxes with the measured parameters shows that the relationships between atomic hydrogen, molecular hydrogen, and dust are different in low- and high-metallicity regions. Low-metallicity regions contain more atomic gas, but less molecular gas and dust, including PAHs. The mass of dust constitutes about $1\%$ of the mass of molecular gas in all regions considered. Fluxes produced by atomic and molecular gas do not correlate with the parameters of the stellar radiation, whereas the dust fluxes grow with increasing mean intensity of stellar radiation and the fraction of enhanced stellar radiation. The ratio of the fluxes at 8 and 24$\mu$m, which characterizes the PAH content, decreases with increasing intensity of the stellar radiation, possibly indicating evolutionary variations of the PAH content. The results confirm that the contribution of the 24$\mu$m emission to the total IR luminosity of extragalactic star-forming regions does not depend on the metallicity.

Galaxies with polar rings consist of two subsystems, a disk and a ring, which rotate almost in orthogonal planes. In this paper, we analyze the parameters characterizing the composition of the interstellar medium and star formation in star-forming complexes, belonging to a polar ring galaxy NGC660. We show that star-forming regions in the ring of the galaxy are distinctively different from those in the galaxy disk. They possess substantially lower infrared luminosities, indicative of less dust mass in these regions than in a typical disk star-forming region. UV and H$\alpha$ luminosities also appear to be lower in the ring, probably, being a consequence of its relatively recent formation.

SPIDER is a balloon-borne instrument designed to map the polarization of the millimeter-wave sky at large angular scales. SPIDER targets the B-mode signature of primordial gravitational waves in the cosmic microwave background (CMB), with a focus on mapping a large sky area with high fidelity at multiple frequencies. SPIDER's first longduration balloon (LDB) flight in January 2015 deployed a total of 2400 antenna-coupled Transition Edge Sensors (TESs) at 90 GHz and 150 GHz. In this work we review the design and in-flight performance of the SPIDER instrument, with a particular focus on the measured performance of the detectors and instrument in a space-like loading and radiation environment. SPIDER's second flight in December 2018 will incorporate payload upgrades and new receivers to map the sky at 285 GHz, providing valuable information for cleaning polarized dust emission from CMB maps.

We describe the construction and characterization of the 280 GHz bolometric focal plane units (FPUs) to be deployed on the second flight of the balloon-borne SPIDER instrument. These FPUs are vital to SPIDER's primary science goal of detecting or placing an upper limit on the amplitude of the primordial gravitational wave signature in the cosmic microwave background (CMB) by constraining the B-mode contamination in the CMB from Galactic dust emission. Each 280 GHz focal plane contains a 16 x 16 grid of corrugated silicon feedhorns coupled to an array of aluminum-manganese transition-edge sensor (TES) bolometers fabricated on 150 mm diameter substrates. In total, the three 280 GHz FPUs contain 1,530 polarization sensitive bolometers (765 spatial pixels) optimized for the low loading environment in flight and read out by time-division SQUID multiplexing. In this paper we describe the mechanical, thermal, and magnetic shielding architecture of the focal planes and present cryogenic measurements which characterize yield and the uniformity of several bolometer parameters. The assembled FPUs have high yields, with one array as high as 95% including defects from wiring and readout. We demonstrate high uniformity in device parameters, finding the median saturation power for each TES array to be ~3 pW at 300 mK with a less than 6% variation across each array at one standard deviation. These focal planes will be deployed alongside the 95 and 150 GHz telescopes in the SPIDER-2 instrument, slated to fly from McMurdo Station in Antarctica in December 2018.

Protoplanetary disk mass is a key parameter controlling the process of planetary system formation. CO molecular emission is often used as a tracer of gas mass in the disk. In this study we consider the ability of CO to trace the gas mass over a wide range of disk structural parameters and search for chemical species that could possibly be used as alternative mass tracers to CO. Specifically, we apply detailed astrochemical modeling to a large set of models of protoplanetary disks around low-mass stars, to select molecules with abundances correlated with the disk mass and being relatively insensitive to other disk properties. We do not consider sophisticated dust evolution models, restricting ourselves with the standard astrochemical assumption of $0.1~\mu $m dust. We find that CO is indeed the best molecular tracer for total gas mass, despite the fact that it is not the main carbon carrier, provided reasonable assumptions about CO abundance in the disk are used. Typically, chemical reprocessing lowers the abundance of CO by a factor of 3, compared to the case of photo-dissociation and freeze-out as the only ways of CO depletion. On average only 13% C-atoms reside in gas-phase CO, albeit with variations from 2 to 30%. CO$_2$, H$_2$O and H$_2$CO can potentially serve as alternative mass tracers, the latter two being only applicable if disk structural parameters are known.

We present a new upper limit on CMB circular polarization from the 2015 flight of SPIDER, a balloon-borne telescope designed to search for $B$-mode linear polarization from cosmic inflation. Although the level of circular polarization in the CMB is predicted to be very small, experimental limits provide a valuable test of the underlying models. By exploiting the non-zero circular-to-linear polarization coupling of the HWP polarization modulators, data from SPIDER's 2015 Antarctic flight provide a constraint on Stokes $V$ at 95 and 150 GHz from $33<\ell<307$. No other limits exist over this full range of angular scales, and SPIDER improves upon the previous limit by several orders of magnitude, providing 95% C.L. constraints on $\ell (\ell+1)C_{\ell}^{VV}/(2\pi)$ ranging from 141 $\mu K ^2$ to 255 $\mu K ^2$ at 150 GHz for a thermal CMB spectrum. As linear CMB polarization experiments become increasingly sensitive, the techniques described in this paper can be applied to obtain even stronger constraints on circular polarization.

In this paper, we extend the study initiated in PaperI by modelling grain ensemble evolution in a dynamical model of an expanding HII region and checking the effects of momentum transfer from dust to gas. The radiation pressure on the dust, the dust drift, and the lug on the gas by the dust are all important process that should be considered simultaneously to describe the dynamics of HII regions. With accounting for the momentum transfer from the dust to the gas, the expansion time of the HII region is notably reduced (for our model of RCW120, the time to reach the observed radius of the HII region is reduced by a factor of 1.5). Under the common approximation of frozen dust, where there is no relative drift between the dust and gas, the radiation pressure from the ionizing star drives the formation of the very deep gas cavity near the star. Such a cavity is much less pronounced when the dust drift is taken into account. The dust drift leads to the two-peak morphology of the dust density distribution and significantly reduces the dust-to-gas ratio in the ionized region (by a factor of 2 to 10). The dust-to-gas ratio is larger for higher temperatures of the ionizing star since the dust grains have a larger electric charge and are more strongly coupled to the gas.

We present results from a new incoherent-beam Fast Radio Burst (FRB) search on the Canadian Hydrogen Intensity Mapping Experiment (CHIME) Pathfinder. Its large instantaneous field of view (FoV) and relative thermal insensitivity allow us to probe the ultra-bright tail of the FRB distribution, and to test a recent claim that this distribution's slope, $\alpha\equiv-\frac{\partial \log N}{\partial \log S}$, is quite small. A 256-input incoherent beamformer was deployed on the CHIME Pathfinder for this purpose. If the FRB distribution were described by a single power-law with $\alpha=0.7$, we would expect an FRB detection every few days, making this the fastest survey on sky at present. We collected 1268 hours of data, amounting to one of the largest exposures of any FRB survey, with over 2.4\,$\times$\,10$^5$\,deg$^2$\,hrs. Having seen no bursts, we have constrained the rate of extremely bright events to $<\!13$\,sky$^{-1}$\,day$^{-1}$ above $\sim$\,220$\sqrt{(\tau/\rm ms)}$ Jy\,ms for $\tau$ between 1.3 and 100\,ms, at 400--800\,MHz. The non-detection also allows us to rule out $\alpha\lesssim0.9$ with 95$\%$ confidence, after marginalizing over uncertainties in the GBT rate at 700--900\,MHz, though we show that for a cosmological population and a large dynamic range in flux density, $\alpha$ is brightness-dependent. Since FRBs now extend to large enough distances that non-Euclidean effects are significant, there is still expected to be a dearth of faint events and relative excess of bright events. Nevertheless we have constrained the allowed number of ultra-intense FRBs. While this does not have significant implications for deeper, large-FoV surveys like full CHIME and APERTIF, it does have important consequences for other wide-field, small dish experiments.

We present evolutionary calculations for the size and aromatization degree distributions of interstellar dust grains, driven by their destruction by radiation, collisions with gas particles, and shattering due to grain-grain collisions. Based on these calculations we model dust emission spectra. The initial grain size distribution play the important role in the evolution of an ensemble of dust particles. Radiation in the considered intensity range mostly aromatizes grains. The smallest grains are mainly destructed via sputtering by collisions with gas particles. There are no grains smaller than 20~\AA\ in the medium at relative gas-dust velocities more than 50~km/s, which is typical for shocks in supernova remnants. The IR emission spectrum changes significantly due to the dust evolution depending on the adopted grain properties, in particular, on the energy of the C--C bonds ($E_0$). Aromatic bands in the near-IR (2--15~$\mu$m) are absent, if $E_0$ is low, even when the medium properties correspond to the average interstellar medium in our Galaxy. As in reality these bands are observed, high $E_0$ values are more preferable. We consider dependence of the emission intensity ratios for various photometric bands on the medium properties. The aromatization degree of small grains shows up most strongly in the $I_{3.4}/I_{11.3}$ intensity ratio while the fraction of aromatic grains in the total dust mass influences the $I_{3.3}/I_{70+160}$ ratio.

A model of key processes influencing the evolution of a hydrocarbon grain of an arbitrary size under astrophysical conditions corresponding to ionized hydrogen regions (HII regions) and supernova remnants is presented. The considered processes include aromatization and photodestruction, sputtering by electrons and ions, and shattering due to collisions between grains. The model can be used to simulate the grain size distribution and the aromatization degree during the evolution of HII regions and supernova remnants for a specified radiation field, relative velocity of gas and dust, etc. The contribution of various processes to the evolution of hydrocarbon dust grains for parameters typical for the interstellar medium of our Galaxy is presented. Small grains (less than 50 carbon atoms) should be fully aromatized in the general interstellar medium. If larger grains initially have an aliphatic structure, it is preserved to a substantial extent. Variations in the size distribution of the grains due to their mutual collisions depend appreciably on the adopted initial size distribution. For the MRN initial distribution a significant redistribution of grain sizes is obtained, which increases the mass fraction of smaller grains. Characteristic for an initial distribution from the work of Jones et al. (2013), with high initial fraction of small grains, is a general decrease in the number of grains of all sizes.

BICEP3, the latest telescope in the BICEP/Keck program, started science observations in March 2016. It is a 550mm aperture refractive telescope observing the polarization of the cosmic microwave background at 95 GHz. We show the focal plane design and detector performance, including spectral response, optical efficiency and preliminary sensitivity of the upgraded BICEP3. We demonstrate 9.72$\mu$K$\sqrt{\textrm{s}}$ noise performance of the BICEP3 receiver.

BICEP3 is a 520 mm aperture, compact two-lens refractor designed to observe the polarization of the cosmic microwave background (CMB) at 95 GHz. Its focal plane consists of modularized tiles of antenna-coupled transition edge sensors (TESs), similar to those used in BICEP2 and the Keck Array. The increased per-receiver optical throughput compared to BICEP2/Keck Array, due to both its faster f/1.7 optics and the larger aperture, more than doubles the combined mapping speed of the BICEP/Keck program. The BICEP3 receiver was recently upgraded to a full complement of 20 tiles of detectors (2560 TESs) and is now beginning its second year of observation (and first science season) at the South Pole. We report on its current performance and observing plans. Given its high per-receiver throughput while maintaining the advantages of a compact design, BICEP3-class receivers are ideally suited as building blocks for a 3rd-generation CMB experiment, consisting of multiple receivers spanning 35 GHz to 270 GHz with total detector count in the tens of thousands. We present plans for such an array, the new "BICEP Array" that will replace the Keck Array at the South Pole, including design optimization, frequency coverage, and deployment/observing strategies.

We describe 280 GHz bolometric detector arrays that instrument the balloon-borne polarimeter SPIDER. A primary science goal of SPIDER is to measure the large-scale B-mode polarization of the cosmic microwave background in search of the cosmic-inflation, gravitational-wave signature. 280 GHz channels aid this science goal by constraining the level of B-mode contamination from galactic dust emission. We present the focal plane unit design, which consists of a 16$\times$16 array of conical, corrugated feedhorns coupled to a monolithic detector array fabricated on a 150 mm diameter silicon wafer. Detector arrays are capable of polarimetric sensing via waveguide probe-coupling to a multiplexed array of transition-edge-sensor (TES) bolometers. The SPIDER receiver has three focal plane units at 280 GHz, which in total contains 765 spatial pixels and 1,530 polarization sensitive bolometers. By fabrication and measurement of single feedhorns, we demonstrate 14.7$^{\circ}$ FHWM Gaussian-shaped beams with $<$1% ellipticity in a 30% fractional bandwidth centered at 280 GHz. We present electromagnetic simulations of the detection circuit, which show 94% band-averaged, single-polarization coupling efficiency, 3% reflection and 3% radiative loss. Lastly, we demonstrate a low thermal conductance bolometer, which is well-described by a simple TES model and exhibits an electrical noise equivalent power (NEP) = 2.6 $\times$ 10$^{-17}$ W/$\sqrt{\mathrm{Hz}}$, consistent with the phonon noise prediction.

The Canadian Hydrogen Intensity Mapping Experiment (CHIME) Pathfinder radio telescope is currently surveying the northern hemisphere between 400 and 800 MHz. By mapping the large scale structure of neutral hydrogen through its redshifted 21 cm line emission between $z \sim 0.8-2.5$ CHIME will contribute to our understanding of Dark Energy. Bright astrophysical foregrounds must be separated from the neutral hydrogen signal, a task which requires precise characterization of the polarized telescope beams. Using the DRAO John A. Galt 26 m telescope, we have developed a holography instrument and technique for mapping the CHIME Pathfinder beams. We report the status of the instrument and initial results of this effort.

A linear polarization field on the sphere can be uniquely decomposed into an E-mode and a B-mode component. These two components are analytically defined in terms of spin-2 spherical harmonics. Maps that contain filtered modes on a partial sky can also be decomposed into E-mode and B-mode components. However, the lack of full sky information prevents orthogonally separating these components using spherical harmonics. In this paper, we present a technique for decomposing an incomplete map into E and B-mode components using E and B eigenmodes of the pixel covariance in the observed map. This method is found to orthogonally define E and B in the presence of both partial sky coverage and spatial filtering. This method has been applied to the BICEP2 and the Keck Array maps and results in reducing E to B leakage from LCDM E-modes to a level corresponding to a tensor-to-scalar ratio of $r<1\times10^{-4}$.

We present measurements of polarization lensing using the 150 GHz maps which include all data taken by the BICEP2 & Keck Array CMB polarization experiments up to and including the 2014 observing season (BK14). Despite their modest angular resolution ($\sim 0.5^\circ$), the excellent sensitivity ($\sim 3\mu$K-arcmin) of these maps makes it possible to directly reconstruct the lensing potential using only information at larger angular scales ($\ell\leq 700$). From the auto-spectrum of the reconstructed potential we measure an amplitude of the spectrum to be $A^{\phi\phi}_{\rm L}=1.15\pm 0.36$ (Planck $\Lambda$CDM prediction corresponds to $A^{\phi\phi}_{\rm L}=1$), and reject the no-lensing hypothesis at 5.8$\sigma$, which is the highest significance achieved to date using an EB lensing estimator. Taking the cross-spectrum of the reconstructed potential with the Planck 2015 lensing map yields $A^{\phi\phi}_{\rm L}=1.13\pm 0.20$. These direct measurements of $A^{\phi\phi}_{\rm L}$ are consistent with the $\Lambda$CDM cosmology, and with that derived from the previously reported BK14 B-mode auto-spectrum ($A^{\rm BB}_{\rm L}=1.20\pm 0.17$). We perform a series of null tests and consistency checks to show that these results are robust against systematics and are insensitive to analysis choices. These results unambiguously demonstrate that the B-modes previously reported by BICEP / Keck at intermediate angular scales ($150\lesssim\ell\lesssim 350$) are dominated by gravitational lensing. The good agreement between the lensing amplitudes obtained from the lensing reconstruction and B-mode spectrum starts to place constraints on any alternative cosmological sources of B-modes at these angular scales.

We present results from an analysis of all data taken by the BICEP2 & Keck Array CMB polarization experiments up to and including the 2014 observing season. This includes the first Keck Array observations at 95 GHz. The maps reach a depth of 50 nK deg in Stokes $Q$ and $U$ in the 150 GHz band and 127 nK deg in the 95 GHz band. We take auto- and cross-spectra between these maps and publicly available maps from WMAP and Planck at frequencies from 23 GHz to 353 GHz. An excess over lensed-LCDM is detected at modest significance in the 95x150 $BB$ spectrum, and is consistent with the dust contribution expected from our previous work. No significant evidence for synchrotron emission is found in spectra such as 23x95, or for correlation between the dust and synchrotron sky patterns in spectra such as 23x353. We take the likelihood of all the spectra for a multi-component model including lensed-LCDM, dust, synchrotron and a possible contribution from inflationary gravitational waves (as parametrized by the tensor-to-scalar ratio $r$), using priors on the frequency spectral behaviors of dust and synchrotron emission from previous analyses of WMAP and Planck data in other regions of the sky. This analysis yields an upper limit $r_{0.05}<0.09$ at 95% confidence, which is robust to variations explored in analysis and priors. Combining these $B$-mode results with the (more model-dependent) constraints from Planck analysis of CMB temperature plus BAO and other data, yields a combined limit $r_{0.05}<0.07$ at 95% confidence. These are the strongest constraints to date on inflationary gravitational waves.

We describe the cryogenic half-wave plate rotation mechanisms built for and used in Spider, a polarization-sensitive balloon-borne telescope array that observed the Cosmic Microwave Background at 95 GHz and 150 GHz during a stratospheric balloon flight from Antarctica in January 2015. The mechanisms operate at liquid helium temperature in flight. A three-point contact design keeps the mechanical bearings relatively small but allows for a large (305 mm) diameter clear aperture. A worm gear driven by a cryogenic stepper motor allows for precise positioning and prevents undesired rotation when the motors are depowered. A custom-built optical encoder system monitors the bearing angle to an absolute accuracy of +/- 0.1 degrees. The system performed well in Spider during its successful 16 day flight.

BICEP3 is a $550~mm$ aperture telescope with cold, on-axis, refractive optics designed to observe at the $95~GHz$ band from the South Pole. It is the newest member of the BICEP/Keck family of inflationary probes specifically designed to measure the polarization of the cosmic microwave background (CMB) at degree-angular scales. BICEP3 is designed to house 1280 dual-polarization pixels, which, when fully-populated, totals to $\sim$9$\times$ the number of pixels in a single Keck $95~GHz$ receiver, thus further advancing the BICEP/Keck program's $95~GHz$ mapping speed. BICEP3 was deployed during the austral summer of 2014-2015 with 9 detector tiles, to be increased to its full capacity of 20 in the second season. After instrument characterization measurements were taken, CMB observation commenced in April 2015. Together with multi-frequency observation data from Planck, BICEP2, and the Keck Array, BICEP3 is projected to set upper limits on the tensor-to-scalar ratio to $r$ $\lesssim 0.03$ at $95\%$ C.L..

We describe the SPIDER flight cryostat, which is designed to cool six millimeter-wavelength telescopes during an Antarctic long-duration balloon flight. The cryostat, one of the largest to have flown on a stratospheric payload, uses liquid helium-4 to deliver cooling power to stages at 4.2 and 1.6 K. Stainless steel capillaries facilitate a high flow impedance connection between the main liquid helium tank and a smaller superfluid tank, allowing the latter to operate at 1.6 K as long as there is liquid in the 4.2 K main tank. Each telescope houses a closed cycle helium-3 adsorption refrigerator that further cools the focal planes down to 300 mK. Liquid helium vapor from the main tank is routed through heat exchangers that cool radiation shields, providing negative thermal feedback. The system performed successfully during a 17 day flight in the 2014-2015 Antarctic summer. The cryostat had a total hold time of 16.8 days, with 15.9 days occurring during flight.

Stellar population and the interstellar gas-dust medium in the vicinity of the open star cluster vdB 130 are analysed using optical observations taken with the 6-m telescope of the SAO RAS and the 125-cm telescope of the SAI MSU along with the data of Spitzer and Herschel. Based on proper motions and BV and JHKs 2MASS photometric data, we select additional 36 stars as probable cluster members. Some stars in vdB 130 are classified as B stars. Our estimates of minimum colour excess, apparent distance modulus and the distance are consistent with young age (from 5 to 10 Myrs) of the cluster vdB 130. We suppose the large deviations from the conventional extinction law in the cluster direction, with $R_V$ ~ 4 - 5. The cluster vdB 130 appears to be physically related to the supershell around Cyg OB1, a cometary CO cloud, ionized gas, and regions of infrared emission. There are a few regions of bright mid-infrared emission in the vicinity of vdB 130. The largest of them is also visible on H-alpha and [SII] emission maps. We suggest that the infrared blobs that coincide in projection with the head of the molecular cloud are HII regions, excited by the cluster B-stars. Some signatures of a shock front are identified between these IR-bright regions.

We consider dust drift under the influence of stellar radiation pressure during the pressure-driven expansion of an HII region using the chemo-dynamical model MARION. Dust size distribution is represented by four dust types: conventional polycyclic aromatic hydrocarbons (PAHs), very small grains (VSGs), big grains (BGs) and also intermediate-sized grains (ISGs), which are larger than VSGs and smaller than BGs. The dust is assumed to move at terminal velocity determined locally from the balance between the radiation pressure and gas drag. As Coulomb drag is an important contribution to the overall gas drag, we evaluate a grain charge evolution within the HII region for each dust type. BGs are effectively swept out of the HII region. The spatial distribution of ISGs within the HII region has a double peak structure, with a smaller inner peak and a higher outer peak. PAHs and VSGs are mostly coupled to the gas. The mean charge of PAHs is close to zero, so they can become neutral from time to time because of charge fluctuations. These periods of neutrality occur often enough to cause the removal of PAHs from the very interior of the HII region. For VSGs, the effect of charge fluctuations is less pronounced but still significant. We conclude that accounting for charge dispersion is necessary to describe the dynamics of small grains.

This paper describes outstanding issues in astrophysics and cosmology that can be solved by astronomical observations in a broad spectral range from far infrared to millimeter wavelengths. The discussed problems related to the formation of stars and planets, galaxies and the interstellar medium, studies of black holes and the development of the cosmological model can be addressed by the planned space observatory Millimetron (the "Spectr-M" project) equipped with a cooled 10-m mirror. Millimetron can operate both as a single-dish telescope and as a part of a space-ground interferometer with very long baseline.

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

We present an analysis of optical spectroscopy and infrared aperture photometry of more than 100 H II complexes in nine galaxies. Spectra obtained with the 6-m telescope of SAO RAS are used along with archival data from Spitzer and several ground-based telescopes to infer a strength of polycyclic aromatic hydrocarbon (PAH) emission, age, properties of the UV radiation field, and metallicity of studied H II complexes. Physical properties (age, radiation field parameters, metallicity) are related to the $F_{8}/F_{24}$ ratio used as a proxy for the PAH abundance in order to reveal factors that may influence the PAH evolution in H II complexes. The well-known correlation between the $F_{8}/F_{24}$ ratio and metallicity is confirmed in the studied complexes. The infrared flux ratio also correlates with the [O III]$\lambda 5007/\mathrm{H\beta}$ ratio which is often considered as an indicator of the radiation field hardness, but this correlation seems to be a mere reflection of a correlation between [O III]$\lambda 5007/\mathrm{H\beta}$ and metallicity. In separate metallicity bins, the $F_{8}/F_{24}$ ratio is found to correlate with an age of an H II complex, which is estimated from the equivalent width of $\mathrm{H}\beta$ line. The correlation is positive for low metallicity complexes and negative for high metallicity complexes. Analysing various mechanisms of PAH formation and destruction in the context of found correlations, we suggest that PAH abundance is likely altered by the UV radiation within H II complexes, but this is not necessarily due to their destruction. If PAHs can also form in H II complexes due to some processes like aromatisation, photodestruction, shattering and sputtering of very small grains, the net $F_{8}/F_{24}$ ratio is determined by a balance between all these processes that can be different at different metallicities.

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.

A detailed photometric study of star-forming regions (SFRs) in the galaxy Holmberg II has been carried out using archival observational data from the far infrared to ultraviolet obtained with the GALEX, Spitzer, and Herschel telescopes. Spectroscopic observations with the 6-m telescope of Special Astrophysical Observatory of the Russian Academy of Sciences are used to estimate ages and metallicities of SFRs. For the first time, the ages of SFRs have been related to their emission parameters in a wide spectral range and with the physical parameters determined by fitting the observed spectra. It is shown that fluxes at 8 and 24 micron characterizing the emission of polycyclic aromatic hydrocarbons (PAHs) and hot dust grains decrease with age, but their ratio increases. This implies that the relative PAH contribution to the total infrared flux increases with age. It is suggested that the detected increase in the ratio of the fluxes at 8 and 24 micron is related to the growth in the PAH mass due to destruction of larger grains.

We present the results of integration and characterization of the SPIDER instrument after the 2013 pre-flight campaign. SPIDER is a balloon-borne polarimeter designed to probe the primordial gravitational wave signal in the degree-scale $B$-mode polarization of the cosmic microwave background. With six independent telescopes housing over 2000 detectors in the 94 GHz and 150 GHz frequency bands, SPIDER will map 7.5% of the sky with a depth of 11 to 14 $\mu$K$\cdot$arcmin at each frequency, which is a factor of $\sim$5 improvement over Planck. We discuss the integration of the pointing, cryogenic, electronics, and power sub-systems, as well as pre-flight characterization of the detectors and optical systems. SPIDER is well prepared for a December 2014 flight from Antarctica, and is expected to be limited by astrophysical foreground emission, and not instrumental sensitivity, over the survey region.

We introduce the light-weight carbon fiber and aluminum gondola designed for the SPIDER balloon-borne telescope. SPIDER is designed to measure the polarization of the Cosmic Microwave Background radiation with unprecedented sensitivity and control of systematics in search of the imprint of inflation: a period of exponential expansion in the early Universe. The requirements of this balloon-borne instrument put tight constrains on the mass budget of the payload. The SPIDER gondola is designed to house the experiment and guarantee its operational and structural integrity during its balloon-borne flight, while using less than 10% of the total mass of the payload. We present a construction method for the gondola based on carbon fiber reinforced polymer tubes with aluminum inserts and aluminum multi-tube joints. We describe the validation of the model through Finite Element Analysis and mechanical tests.

We present the technology and control methods developed for the pointing system of the SPIDER experiment. SPIDER is a balloon-borne polarimeter designed to detect the imprint of primordial gravitational waves in the polarization of the Cosmic Microwave Background radiation. We describe the two main components of the telescope's azimuth drive: the reaction wheel and the motorized pivot. A 13 kHz PI control loop runs on a digital signal processor, with feedback from fibre optic rate gyroscopes. This system can control azimuthal speed with < 0.02 deg/s RMS error. To control elevation, SPIDER uses stepper-motor-driven linear actuators to rotate the cryostat, which houses the optical instruments, relative to the outer frame. With the velocity in each axis controlled in this way, higher-level control loops on the onboard flight computers can implement the pointing and scanning observation modes required for the experiment. We have accomplished the non-trivial task of scanning a 5000 lb payload sinusoidally in azimuth at a peak acceleration of 0.8 deg/s$^2$, and a peak speed of 6 deg/s. We can do so while reliably achieving sub-arcminute pointing control accuracy.

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.

A pathfinder version of CHIME (the Canadian Hydrogen Intensity Mapping Experiment) is currently being commissioned at the Dominion Radio Astrophysical Observatory (DRAO) in Penticton, BC. The instrument is a hybrid cylindrical interferometer designed to measure the large scale neutral hydrogen power spectrum across the redshift range 0.8 to 2.5. The power spectrum will be used to measure the baryon acoustic oscillation (BAO) scale across this poorly probed redshift range where dark energy becomes a significant contributor to the evolution of the Universe. The instrument revives the cylinder design in radio astronomy with a wide field survey as a primary goal. Modern low-noise amplifiers and digital processing remove the necessity for the analog beamforming that characterized previous designs. The Pathfinder consists of two cylinders 37\,m long by 20\,m wide oriented north-south for a total collecting area of 1,500 square meters. The cylinders are stationary with no moving parts, and form a transit instrument with an instantaneous field of view of $\sim$100\,degrees by 1-2\,degrees. Each CHIME Pathfinder cylinder has a feedline with 64 dual polarization feeds placed every $\sim$30\,cm which Nyquist sample the north-south sky over much of the frequency band. The signals from each dual-polarization feed are independently amplified, filtered to 400-800\,MHz, and directly sampled at 800\,MSps using 8 bits. The correlator is an FX design, where the Fourier transform channelization is performed in FPGAs, which are interfaced to a set of GPUs that compute the correlation matrix. The CHIME Pathfinder is a 1/10th scale prototype version of CHIME and is designed to detect the BAO feature and constrain the distance-redshift relation.

The Canadian Hydrogen Intensity Mapping Experiment (CHIME) is a transit interferometer currently being built at the Dominion Radio Astrophysical Observatory (DRAO) in Penticton, BC, Canada. We will use CHIME to map neutral hydrogen in the frequency range 400 -- 800\,MHz over half of the sky, producing a measurement of baryon acoustic oscillations (BAO) at redshifts between 0.8 -- 2.5 to probe dark energy. We have deployed a pathfinder version of CHIME that will yield constraints on the BAO power spectrum and provide a test-bed for our calibration scheme. I will discuss the CHIME calibration requirements and describe instrumentation we are developing to meet these requirements.

We report first science results from our new spectrometer, the 2nd generation z(Redshift) and Early Universe Spectrometer (ZEUS-2), recently commissioned on the Atacama Pathfinder Experiment telescope (APEX). ZEUS-2 is a submillimeter grating spectrometer optimized for detecting the faint and broad lines from distant galaxies that are redshifted into the telluric windows from 200 to 850 microns. It utilizes a focal plane array of transition-edge sensed bolometers, the first use of these arrays for astrophysical spectroscopy. ZEUS-2 promises to be an important tool for studying galaxies in the years to come due to its synergy with ALMA and its capabilities in the short submillimeter windows that are unique in the post Herschel era. Here we report on our first detection of the [CII] 158 $\mu m$ line with ZEUS-2. We detect the line at z ~ 1.8 from H-ATLAS J091043.1-000322 with a line flux of $(6.44 \pm 0.42) \times 10^{-18} W m^{-2}$. Combined with its far-infrared luminosity and a new Herschel-PACS detection of the [OI] 63 $\mu m $ line we model the line emission as coming from a photo-dissociation region with far-ultraviolet radiation field, $G \approx 2 \times 10^{4} G_{0}$, gas density, $n \approx 1 \times 10^{3} cm^{-3}$ and size between ~ 0.4 and 1 kpc. Based on this model, we conclude that H-ATLAS J091043.1-000322 is a high redshift analogue of a local ultra-luminous infrared galaxy, i.e. it is likely the site of a compact starburst due to a major merger. Further identification of these merging systems is important for constraining galaxy formation and evolution models.

Gas density and temperature in star forming regions around Sh2-235 are derived from ammonia line observations. This information is used to evaluate formation scenarios and to determine evolutionary stages of the young embedded clusters S235 East1, S235 East2, and S235 Central. We also estimate the gas mass in the embedded clusters and its ratio to the stellar mass. S235 East1 appears to be less evolved than S235 East2 and S235 Central. In S235 East1 the molecular gas mass exceeds that in the other clusters. Also, this cluster is more embedded in the parent gas cloud than the other two. Comparison with a theoretical model shows that the formation of these three clusters could have been stimulated by the expansion of the Sh2-235 HII region (hereafter S235) via a collect-and-collapse process, provided the density in the surrounding gas exceeds $3\cdot10^3$ cm$^{-3}$, or via collapse of pre-existing clumps. The expansion of S235 cannot be responsible for star formation in the southern S235 A-B region. However, formation of the massive stars in this region might have been triggered by a large-scale supernova shock. Thus, triggered star formation in the studied region may come in three varieties, namely collect-and-collapse and collapse of pre-existing clumps, both initiated by expansion of the local HII regions, and triggering by an external large-scale shock. We argue that the C235 A HII region expands into a highly non-uniform medium with increasing density. It is too young to trigger star formation in its vicinity by a collect-and-collapse process. There is an age spread inside the S235 A-B region. Massive stars in the S235 A-B region are considerably younger than lower mass stars in the same area. This follows from the estimates of their ages and the ages of associated HII regions.

The generation of infrared (IR) radiation and the observed IR intensity distribution at wavelengths of 8, 24, and 100 micron in the ionized hydrogen region around a young, massive star is investigated. The evolution of the HII region is treated using a self-consistent chemical-dynamical model in which three dust populations are included -- large silicate grains, small graphite grains, and polycyclic, aromatic hydrocarbons (PAHs). A radiative transfer model taking into account stochastic heating of small grains and macromolecules is used to model the IR spectral energy distribution. The computational results are compared with Spitzer and Herschel observations of the RCW 120 nebula. The contributions of collisions with gas particles and the radiation field of the star to stochastic heating of small grains are investigated. It is shown that a model with a homogeneous PAH content cannot reproduce the ring-like IR-intensity distribution at 8 micron. A model in which PAHs are destroyed in the ionized region provides a means to explain this intensity distribution. This model is in agreement with observations for realistic characteristic destruction times for the PAHs.

We evaluate the ability of SPIDER, a balloon-borne polarimeter, to detect a divergence-free polarization pattern ("B-modes") in the Cosmic Microwave Background (CMB). In the inflationary scenario, the amplitude of this signal is proportional to that of the primordial scalar perturbations through the tensor-to-scalar ratio r. We show that the expected level of systematic error in the SPIDER instrument is significantly below the amplitude of an interesting cosmological signal with r=0.03. We present a scanning strategy that enables us to minimize uncertainty in the reconstruction of the Stokes parameters used to characterize the CMB, while accessing a relatively wide range of angular scales. Evaluating the amplitude of the polarized Galactic emission in the SPIDER field, we conclude that the polarized emission from interstellar dust is as bright or brighter than the cosmological signal at all SPIDER frequencies (90 GHz, 150 GHz, and 280 GHz), a situation similar to that found in the "Southern Hole." We show that two ~20-day flights of the SPIDER instrument can constrain the amplitude of the B-mode signal to r<0.03 (99% CL) even when foreground contamination is taken into account. In the absence of foregrounds, the same limit can be reached after one 20-day flight.

The abundance of polycyclic aromatic hydrocarbons (PAHs) in low- and high-metallicity galaxies has been widely discussed since the time when detailed infrared data for extragalactic objects were first obtained. On the scales of entire galaxies, a smaller PAH abundance in lower-metallicity galaxies is often observed. We study this relationship for star-forming regions in nearby galaxies, for a sample containing more than 200 HII complexes, using spatially-resolved observations from the Herschel Space Observatory and Spitzer Space Telescope. We use a model for the dust emission to estimate the physical parameters (PAH abundance, metallicity, ultraviolet radiation field, etc.) of these complexes. The same correlation of PAH abundance with metallicity, as seen for entire galaxies, is apparently preserved at smaller scales, at least when the Kobulnicky & Kewley metallicity calibration is used. We discuss possible reasons for this correlation, noting that traces of less-effective PAH formation in low-metallicity AGB stars should be smeared out by radial mixing in galactic disks. Effective destruction by the harder and more intensive ultraviolet field in low-metallicity environments is qualitatively consistent with our data, as the ultraviolet field intensity, derived from the infrared photometry, is indeed smaller in HII complexes with lower metallicity.

We present a self-consistent model of a protoplanetary disk: 'ANDES' ('AccretioN disk with Dust Evolution and Sedimentation'). ANDES is based on a flexible and extendable modular structure that includes 1) a 1+1D frequency-dependent continuum radiative transfer module, 2) a module to calculate the chemical evolution using an extended gas-grain network with UV/X-ray-driven processes surface reactions, 3) a module to calculate the gas thermal energy balance, and 4) a 1+1D module that simulates dust grain evolution. For the first time, grain evolution and time-dependent molecular chemistry are included in a protoplanetary disk model. We find that grain growth and sedimentation of large grains to the disk midplane lead to a dust-depleted atmosphere. Consequently, dust and gas temperatures become higher in the inner disk (R < 50 AU) and lower in the outer disk (R > 50 AU), in comparison with the disk model with pristine dust. The response of disk chemical structure to the dust growth and sedimentation is twofold. First, due to higher transparency a partly UV-shielded molecular layer is shifted closer to the dense midplane. Second, the presence of big grains in the disk midplane delays the freeze-out of volatile gas-phase species such as CO there, while in adjacent upper layers the depletion is still effective. Molecular concentrations and thus column densities of many species are enhanced in the disk model with dust evolution, e.g., CO2, NH2CN, HNO, H2O, HCOOH, HCN, CO. We also show that time-dependent chemistry is important for a proper description of gas thermal balance.

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

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

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

The submillimeter opacity of dust in the diffuse Galactic interstellar medium (ISM) has been quantified using a pixel-by-pixel correlation of images of continuum emission with a proxy for column density. We used three BLAST bands at 250, 350, and 500 \mu m and one IRAS at 100 \mu m. The proxy is the near-infrared color excess, E(J-Ks), obtained from 2MASS. Based on observations of stars, we show how well this color excess is correlated with the total hydrogen column density for regions of moderate extinction. The ratio of emission to column density, the emissivity, is then known from the correlations, as a function of frequency. The spectral distribution of this emissivity can be fit by a modified blackbody, whence the characteristic dust temperature T and the desired opacity \sigma_e(1200) at 1200 GHz can be obtained. We have analyzed 14 regions near the Galactic plane toward the Vela molecular cloud, mostly selected to avoid regions of high column density (N_H > 10^{22} cm^-2) and small enough to ensure a uniform T. We find \sigma_e(1200) is typically 2 to 4 x 10^{-25} cm^2/H and thus about 2 to 4 times larger than the average value in the local high Galactic latitude diffuse atomic ISM. This is strong evidence for grain evolution. There is a range in total power per H nucleon absorbed (re-radiated) by the dust, reflecting changes in the interstellar radiation field and/or the dust absorption opacity. These changes affect the equilibrium T, which is typically 15 K, colder than at high latitudes. Our analysis extends, to higher opacity and lower T, the trend of increasing opacity with decreasing T that was found at high latitudes. The recognition of changes in the emission opacity raises a cautionary flag because all column densities deduced from dust emission maps, and the masses of compact structures within them, depend inversely on the value adopted.

Stochastic heating of small grains is often mentioned as a primary cause of large infrared (IR) fluxes from star-forming galaxies, e.g. at 24\mu m. If the mechanism does work at a galaxy-wide scale, it should show up at smaller scales as well. We calculate temperature probability density distributions within a model protostellar core for four dust components: large silicate and graphite grains, small graphite grains, and polycyclic aromatic hydrocarbon particles. The corresponding spectral energy distributions are calculated and compared with observations of a representative infrared dark cloud core. We show that stochastic heating, induced by the standard interstellar radiation field, cannot explain high mid-IR emission toward the centre of the core. In order to reproduce the observed emission from the core projected centre, in particular, at 24\mu m, we need to increase the ambient radiation field by a factor of about 70. However, the model with enhanced radiation field predicts even higher intensities at the core periphery, giving it a ring-like appearance, that is not observed. We discuss possible implications of this finding and also discuss a role of other non-radiative dust heating processes.

We present measurements of the mean mid-infrared-to-submillimeter flux densities of massive (M\ast \approx 2 \times 10^11 Msun) galaxies at redshifts 1.7 < z < 2.9, obtained by stacking positions of known objects taken from the GOODS NICMOS Survey (GNS) catalog on maps: at 24 {\mu}m (Spitzer/MIPS); 70, 100, and 160{\mu}m (Herschel/PACS); 250, 350, 500{\mu}m (BLAST); and 870{\mu}m (LABOCA). A modified blackbody spectrum fit to the stacked flux densities indicates a median [interquartile] star-formation rate of SFR = 63 [48, 81] Msun yr^-1 . We note that not properly accounting for correlations between bands when fitting stacked data can significantly bias the result. The galaxies are divided into two groups, disk-like and spheroid-like, according to their Sersic indices, n. We find evidence that most of the star formation is occurring in n \leq 2 (disk-like) galaxies, with median [interquartile] SFR = 122 [100,150] Msun yr^-1, while there are indications that the n > 2 (spheroid-like) population may be forming stars at a median [interquartile] SFR = 14 [9,20] Msun yr^-1, if at all. Finally, we show that star formation is a plausible mechanism for size evolution in this population as a whole, but find only marginal evidence that it is what drives the expansion of the spheroid-like galaxies.

This is the second paper in a series where we study the influence of transport processes on the chemical evolution of protoplanetary disks. Our analysis is based on a flared alpha-model of the DM Tau system, coupled to a large gas-grain chemical network. To account for production of complex molecules, the chemical network is supplied with an extended set of surface reactions and photo-processes in ice mantles. Our disk model covers a wide range of radii, 10-800 AU (from a Jovian planet-forming zone to the outer disk edge). Turbulent transport of gases and ices is implicitly modeled in full 2D along with the time-dependent chemistry. Two regimes are considered, with high and low efficiency of turbulent mixing. The results of the chemical model with suppressed turbulent diffusion are close to those from the laminar model, but not completely. A simple analysis for the laminar chemical model to highlight potential sensitivity of a molecule to transport processes is performed. It is shown that the higher the ratio of the characteristic chemical timescale to the turbulent transport timescale for a given molecule, the higher the probability that its column density will be affected by diffusion. We find that turbulent transport enhances abundances and column densities of many gas-phase species and ices, particularly, complex ones. For such species a chemical steady-state is not reached due to long timescales associated with evaporation and surface photoprocessing and recombination. In contrast, simple radicals and molecular ions, which chemical evolution is fast and proceeds solely in the gas phase, are not much affected by dynamics. All molecules are divided into three groups according to the sensitivity of their column densities to the turbulent diffusion. [Abridged]

We constrain the evolution of the rest-frame far-infrared (FIR) luminosity function out to high redshift, by combining several pieces of complementary information provided by the deep Balloon-borne Large-Aperture Submillimeter Telescope surveys at 250, 350 and 500 micron, as well as other FIR and millimetre data. Unlike most other phenomenological models, we characterise the uncertainties in our fitted parameters using Monte Carlo Markov Chains. We use a bivariate local luminosity function that depends only on FIR luminosity and 60-to-100 micron colour, along with a single library of galaxy spectral energy distributions indexed by colour, and apply simple luminosity and density evolution. We use the surface density of sources, Cosmic Infrared Background (CIB) measurements and redshift distributions of bright sources, for which identifications have been made, to constrain this model. The precise evolution of the FIR luminosity function across this crucial range has eluded studies at longer wavelengths (e.g., using SCUBA and MAMBO) and at shorter wavelengths (e.g., Spitzer), and should provide a key piece of information required for the study of galaxy evolution. Our adoption of Monte Carlo methods enables us not only to find the best-fit evolution model, but also to explore correlations between the fitted parameters. Our model-fitting approach allows us to focus on sources of tension coming from the combination of data-sets. We specifically find that our choice of parameterisation has difficulty fitting the combination of CIB measurements and redshift distribution of sources near 1 mm. Existing and future data sets will be able to dramatically improve the fits, as well as break strong degeneracies among the models. [abridged]

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

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

We study details of the UV radiation transfer in a protoplanetary disk, paying attention to the influence of dust growth and sedimentation on the disk density and temperature. Also, we show how the dust evolution affects photoreaction rates of key molecules, like CN and CS.

We identify and investigate the nature of the 20 brightest 250-micron sources detected by the BLAST within the central 150 sq. arcmin of the GOODS-South field. Aided by the available deep VLA radio imaging, reaching S_1.4 = 30 micro-Jy, we have identified radio counterparts for 17/20 of the 250-micron sources. The resulting enhanced positional accuracy of ~1 arcsec has then allowed us to exploit the deep multi-frequency imaging of GOODS-South to establish secure galaxy counterparts for the 17 radio-identified sources, and plausible galaxy candidates for the 3 radio-unidentified sources. Confusion is a serious issue for this deep BLAST 250-micron survey, due to the large size of the beam. Nevertheless, we argue that our chosen counterparts are significant, and often dominant contributors to the measured BLAST flux densities. For all of these 20 galaxies we have been able to determine spectroscopic (8) or photometric (12) redshifts. The result is the first `complete' redshift distribution for a deep 250-micron selected galaxy sample. This reveals that 250-micron surveys reaching detection limits of ~30 mJy contain not only low-redshift spirals/LIRGs, but also the extreme z~2 dust-enshrouded starburst galaxies previously discovered at sub-millimetre wavelengths. Based on their IRAC colours, we find that virtually all of the BLAST galaxy identifications appear better described as analogues of the M82 starburst galaxy, or Sc star-forming discs rather than highly obscured ULIRGs. Inspection of the LABOCA 870-micron imaging of the GOODS-South field yields detections of 7/11 of the z>1 BLAST sources, and reveals 250/870 flux-density ratios consistent with a standard 40K modified black-body fit with a dust emissivity index beta=1.5.

Infrared observations from the Spitzer Space Telescope archive are used to study the dust component of the interstellar medium in the IC~10 irregular galaxy. Dust distribution in the galaxy is compared to the distributions of H$\alpha$ and [SII] emission, neutral hydrogen and CO clouds, and ionizing radiation sources. The distribution of polycyclic aromatic hydrocarbons (PAH) in the galaxy is shown to be highly non-uniform with the mass fraction of these particles in the total dust mass reaching 4%. PAHs tend to avoid bright HII regions and correlate well with atomic and molecular gas. This pattern suggests that PAHs form in the dense interstellar gas. We propose that the significant decrease of the PAH abundance at low metallicity is observed not only globally (at the level of entire galaxies), but also locally (at least, at the level of individual HII regions). We compare the distribution of the PAH mass fraction to the distribution of high-velocity features, that we have detected earlier in wings of H$\alpha$ and SII lines, over the entire available galaxy area. No conclusive evidence for shock destruction of PAHs in the IC~10 galaxy could be found.

Spider is a balloon-borne instrument designed to map the polarization of the cosmic microwave background (CMB) with degree-scale resolution over a large fraction of the sky. Spider's main goal is to measure the amplitude of primordial gravitational waves through their imprint on the polarization of the CMB if the tensor-to-scalar ratio, r, is greater than 0.03. To achieve this goal, instrumental systematic errors must be controlled with unprecedented accuracy. Here, we build on previous work to use simulations of Spider observations to examine the impact of several systematic effects that have been characterized through testing and modeling of various instrument components. In particular, we investigate the impact of the non-ideal spectral response of the half-wave plates, coupling between focal plane components and the Earth's magnetic field, and beam mismatches and asymmetries. We also present a model of diffuse polarized foreground emission based on a three-dimensional model of the Galactic magnetic field and dust, and study the interaction of this foreground emission with our observation strategy and instrumental effects. We find that the expected level of foreground and systematic contamination is sufficiently low for Spider to achieve its science goals.

We have carried out the first general submillimeter analysis of the field towards GRSMC 45.46+0.05, a massive star forming region in Aquila. The deconvolved 6 deg^2 (3\degree X 2\degree) maps provided by BLAST in 2005 at 250, 350, and 500 micron were used to perform a preliminary characterization of the clump population previously investigated in the infrared, radio, and molecular maps. Interferometric CORNISH data at 4.8 GHz have also been used to characterize the Ultracompact HII regions (UCHIIRs) within the main clumps. By means of the BLAST maps we have produced an initial census of the submillimeter structures that will be observed by Herschel, several of which are known Infrared Dark Clouds (IRDCs). Our spectral energy distributions of the main clumps in the field, located at ~7 kpc, reveal an active population with temperatures of T~35-40 K and masses of ~10^3 Msun for a dust emissivity index beta=1.5. The clump evolutionary stages range from evolved sources, with extended HII regions and prominent IR stellar population, to massive young stellar objects, prior to the formation of an UCHIIR.The CORNISH data have revealed the details of the stellar content and structure of the UCHIIRs. In most cases, the ionizing stars corresponding to the brightest radio detections are capable of accounting for the clump bolometric luminosity, in most cases powered by embedded OB stellar clusters.

We study the impact of dust evolution in a protoplanetary disk around a T Tauri star on the disk chemical composition. For the first time we utilize a comprehensive model of dust evolution which includes growth, fragmentation and sedimentation. Specific attention is paid to the influence of grain evolution on the penetration of the UV field in the disk. A chemical model that includes a comprehensive set of gas phase and grain surface chemical reactions is used to simulate the chemical structure of the disk. The main effect of the grain evolution on the disk chemical composition comes from sedimentation, and, to a lesser degree, from the reduction of the total grain surface area. The net effect of grain growth is suppressed by the fragmentation process which maintains a population of small grains, dominating the total grain surface area. We consider three models of dust properties. In model GS both growth and sedimentation are taken into account. In models A5 and A4 all grains are assumed to have the same size (10(-5) cm and 10(-4) cm, respectively) with constant gas-to-dust mass ratio of 100. Like in previous studies, the "three-layer" pattern (midplane, molecular layer, hot atmosphere) in the disk chemical structure is preserved in all models, but shifted closer to the midplane in models with increased grain size (GS and A4). Unlike other similar studies, we find that in models GS and A4 column densities of most gas-phase species are enhanced by 1-3 orders of magnitude relative to those in a model with pristine dust (A5), while column densities of their surface counterparts are decreased. We show that column densities of certain species, like C2H, HC(2n+1)N (n=0-3), H2O and some other molecules, as well as the C2H2/HCN abundance ratio which are accessible with Herschel and ALMA can be used as observational tracers of early stages of the grain evolution process in protoplanetary disks.

We present Cygnus X in a new multi-wavelength perspective based on an unbiased BLAST survey at 250, 350, and 500 micron, combined with rich datasets for this well-studied region. Our primary goal is to investigate the early stages of high mass star formation. We have detected 184 compact sources in various stages of evolution across all three BLAST bands. From their well-constrained spectral energy distributions, we obtain the physical properties mass, surface density, bolometric luminosity, and dust temperature. Some of the bright sources reaching 40 K contain well-known compact H II regions. We relate these to other sources at earlier stages of evolution via the energetics as deduced from their position in the luminosity-mass (L-M) diagram. The BLAST spectral coverage, near the peak of the spectral energy distribution of the dust, reveals fainter sources too cool (~ 10 K) to be seen by earlier shorter-wavelength surveys like IRAS. We detect thermal emission from infrared dark clouds and investigate the phenomenon of cold ``starless cores" more generally. Spitzer images of these cold sources often show stellar nurseries, but these potential sites for massive star formation are ``starless" in the sense that to date there is no massive protostar in a vigorous accretion phase. We discuss evolution in the context of the L-M diagram. Theory raises some interesting possibilities: some cold massive compact sources might never form a cluster containing massive stars; and clusters with massive stars might not have an identifiable compact cold massive precursor.

We carry out a multi-wavelength study of individual galaxies detected by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) and identified at other wavelengths, using data spanning the radio to the ultraviolet (UV). We develop a Monte Carlo method to account for flux boosting, source blending, and correlations among bands, which we use to derive deboosted far-infrared (FIR) luminosities for our sample. We estimate total star-formation rates for BLAST counterparts with z < 0.9 by combining their FIR and UV luminosities. Star formation is heavily obscured at L_FIR > 10^11 L_sun, z > 0.5, but the contribution from unobscured starlight cannot be neglected at L_FIR < 10^11 L_sun, z < 0.25. We assess that about 20% of the galaxies in our sample show indication of a type-1 active galactic nucleus (AGN), but their submillimeter emission is mainly due to star formation in the host galaxy. We compute stellar masses for a subset of 92 BLAST counterparts; these are relatively massive objects, with a median mass of ~10^11 M_sun, which seem to link the 24um and SCUBA populations, in terms of both stellar mass and star-formation activity. The bulk of the BLAST counterparts at z<1 appear to be run-of-the-mill star-forming galaxies, typically spiral in shape, with intermediate stellar masses and practically constant specific star-formation rates. On the other hand, the high-z tail of the BLAST counterparts significantly overlaps with the SCUBA population, in terms of both star-formation rates and stellar masses, with observed trends of specific star-formation rate that support strong evolution and downsizing.

A one-dimensional method for reconstructing the structure of prestellar and protostellar clouds is presented. The method is based on radiative transfer computations and a comparison of theoretical and observed intensity distributions at both millimeter and infrared wavelengths. The radiative transfer of dust emission is modeled for specified parameters of the density distribution, central star, and external background, and the theoretical distribution of the dust temperature inside the cloud is determined. The intensity distributions at millimeter and IR wavelengths are computed and quantitatively compared with observational data. The best-fit model parameters are determined using a genetic minimization algorithm, which makes it possible to reveal the ranges of parameter degeneracy as well. The method is illustrated by modeling the structure of the two infrared dark clouds IRDC-320.27+029 (P2) and IRDC-321.73+005 (P2). The derived density and temperature distributions can be used to model the chemical structure and spectral maps in molecular lines.

We present an implementation of the iterative flux-conserving Lucy-Richardson (L-R) deconvolution method of image restoration for maps produced by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST). We have analyzed its performance and convergence extensively through simulations and cross-correlations of the deconvolved images with available highresolution maps. We present new science results from two BLAST surveys, in the Galactic regions K3-50 and IC 5146, further demonstrating the benefits of performing this deconvolution. We have resolved three clumps within a radius of 4.'5 inside the star-forming molecular cloud containing K3-50. Combining the well-resolved dust emission map with available multi-wavelength data, we have constrained the Spectral Energy Distributions (SEDs) of five clumps to obtain masses (M), bolometric luminosities (L), and dust temperatures (T). The L-M diagram has been used as a diagnostic tool to estimate the evolutionary stages of the clumps. There are close relationships between dust continuum emission and both 21-cm radio continuum and 12CO molecular line emission. The restored extended large scale structures in the Northern Streamer of IC 5146 have a strong spatial correlation with both SCUBA and high resolution extinction images. A dust temperature of 12 K has been obtained for the central filament. We report physical properties of ten compact sources, including six associated protostars, by fitting SEDs to multi-wavelength data. All of these compact sources are still quite cold (typical temperature below ~ 16 K) and are above the critical Bonner-Ebert mass. They have associated low-power Young Stellar Objects (YSOs). Further evidence for starless clumps has also been found in the IC 5146 region.

We present results from a survey carried out by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) on a 9 deg^2 field near the South Ecliptic Pole at 250, 350 and 500 {\mu}m. The median 1{\sigma} depths of the maps are 36.0, 26.4 and 18.4 mJy, respectively. We apply a statistical method to estimate submillimeter galaxy number counts and find that they are in agreement with other measurements made with the same instrument and with the more recent results from Herschel/SPIRE. Thanks to the large field observed, the new measurements give additional constraints on the bright end of the counts. We identify 132, 89 and 61 sources with S/N>4 at 250, 350, 500 {\mu}m, respectively and provide a multi-wavelength combined catalog of 232 sources with a significance >4{\sigma} in at least one BLAST band. The new BLAST maps and catalogs are available publicly at http://blastexperiment.info.

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

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.

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

We use new large area far infrared maps ranging from 65 - 500 microns obtained with the AKARI and the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) missions to characterize the dust emission toward the Cassiopeia A supernova remnant (SNR). Using the AKARI high resolution data we find a new "tepid" dust grain population at a temperature of ~35K and with an estimated mass of 0.06 solar masses. This component is confined to the central area of the SNR and may represent newly-formed dust in the unshocked supernova ejecta. While the mass of tepid dust that we measure is insufficient by itself to account for the dust observed at high redshift, it does constitute an additional dust population to contribute to those previously reported. We fit our maps at 65, 90, 140, 250, 350, and 500 microns to obtain maps of the column density and temperature of "cold" dust (near 16 K) distributed throughout the region. The large column density of cold dust associated with clouds seen in molecular emission extends continuously from the surrounding interstellar medium to project on the SNR, where the foreground component of the clouds is also detectable through optical, X-ray, and molecular extinction. At the resolution available here, there is no morphological signature to isolate any cold dust associated only with the SNR from this confusing interstellar emission. Our fit also recovers the previously detected "hot" dust in the remnant, with characteristic temperature 100 K.

We present observations at 250, 350, and 500 um of the nearby galaxy cluster Abell 3112 (z=0.075) carried out with BLAST, the Balloon-borne Large Aperture Submillimeter Telescope. Five cluster members are individually detected as bright submillimetre sources. Their far-infrared SEDs and optical colours identify them as normal star-forming galaxies of high mass, with globally evolved stellar populations. They all have B-R colours of 1.38+/-0.08, transitional between the blue, active population and the red, evolved galaxies that dominate the cluster core. We stack to determine the mean submillimetre emission from all cluster members, which is determined to be 16.6+/-2.5, 6.1+/-1.9, and 1.5+/-1.3 mJy at 250, 350, and 500 um, respectively. Stacking analyses of the submillimetre emission of cluster members reveal trends in the mean far-infrared luminosity with respect to cluster-centric radius and Ks-band magnitude. We find that a large fraction of submillimetre emission comes from the boundary of the inner, virialized region of the cluster, at cluster-centric distances around R_500. Stacking also shows that the bulk of the submillimetre emission arises in intermediate-mass galaxies (L<L*), with Ks magnitude ~1 mag fainter than the giant ellipticals. The results and constraints obtained in this work will provide a useful reference for the forthcoming surveys to be conducted on galaxy clusters by Herschel.

The Balloon-borne Large-Aperture Submillimeter Telescope (BLAST) carried out a 250, 350 and 500 micron survey of the galactic plane encompassing the Vela Molecular Ridge, with the primary goal of identifying the coldest dense cores possibly associated with the earliest stages of star formation. Here we present the results from observations of the Vela-D region, covering about 4 square degrees, in which we find 141 BLAST cores. We exploit existing data taken with the Spitzer MIPS, IRAC and SEST-SIMBA instruments to constrain their (single-temperature) spectral energy distributions, assuming a dust emissivity index beta = 2.0. This combination of data allows us to determine the temperature, luminosity and mass of each BLAST core, and also enables us to separate starless from proto-stellar sources. We also analyze the effects that the uncertainties on the derived physical parameters of the individual sources have on the overall physical properties of starless and proto-stellar cores, and we find that there appear to be a smooth transition from the pre- to the proto-stellar phase. In particular, for proto-stellar cores we find a correlation between the MIPS24 flux, associated with the central protostar, and the temperature of the dust envelope. We also find that the core mass function of the Vela-D cores has a slope consistent with other similar (sub)millimeter surveys.

We investigate the correlation between far-infrared (FIR) and radio luminosities in distant galaxies, a lynchpin of modern astronomy. We use data from the Balloon-borne Large Aperture Submillimetre Telescope (BLAST), Spitzer, the Large Apex BOlometer CamerA (LABOCA), the Very Large Array (VLA) and the Giant Metre-wave Radio Telescope (GMRT) in the Extended Chandra Deep Field South (ECDFS). For a catalogue of BLAST 250-micron-selected galaxies, we re-measure the 70--870-micron flux densities at the positions of their most likely 24-micron counterparts, which have a median [interquartile] redshift of 0.74 [0.25, 1.57]. From these, we determine the monochromatic flux density ratio, q_250 = log_10 (S_250micron / S_1400MHz), and the bolometric equivalent, q_IR. At z ~= 0.6, where our 250-micron filter probes rest-frame 160-micron emission, we find no evolution relative to q_160 for local galaxies. We also stack the FIR and submm images at the positions of 24-micron- and radio-selected galaxies. The difference between q_IR seen for 250-micron- and radio-selected galaxies suggests star formation provides most of the IR luminosity in ~< 100-uJy radio galaxies, but rather less for those in the mJy regime. For the 24-micron sample, the radio spectral index is constant across 0 < z < 3, but q_IR exhibits tentative evidence of a steady decline such that q_IR is proportional to (1+z)^(-0.15 +/- 0.03) - significant evolution, spanning the epoch of galaxy formation, with major implications for techniques that rely on the FIR/radio correlation. We compare with model predictions and speculate that we may be seeing the increase in radio activity that gives rise to the radio background.

We describe the application of a statistical method to estimate submillimeter galaxy number counts from confusion limited observations by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST). Our method is based on a maximum likelihood fit to the pixel histogram, sometimes called 'P(D)', an approach which has been used before to probe faint counts, the difference being that here we advocate its use even for sources with relatively high signal-to-noise ratios. This method has an advantage over standard techniques of source extraction in providing an unbiased estimate of the counts from the bright end down to flux densities well below the confusion limit. We specifically analyse BLAST observations of a roughly 10 sq. deg. map centered on the Great Observatories Origins Deep Survey South (GOODS-S) field. We provide estimates of number counts at the three BLAST wavelengths, 250, 350, and 500 microns; instead of counting sources in flux bins we estimate the counts at several flux density nodes connected with power-laws. We observe a generally very steep slope for the counts of about -3.7 at 250 microns and -4.5 at 350 and 500 microns, over the range ~0.02-0.5 Jy, breaking to a shallower slope below about 0.015 Jy at all three wavelengths. We also describe how to estimate the uncertainties and correlations in this method so that the results can be used for model-fitting. This method should be well-suited for analysis of data from the Herschel satellite.

We detect correlations in the cosmic far-infrared background due to the clustering of star-forming galaxies in observations made with the Balloon-borne Large Aperture Submillimeter Telescope, BLAST, at 250, 350, and 500 microns. We perform jackknife and other tests to confirm the reality of the signal. The measured correlations are well fit by a power law over scales of 5-25 arcminutes, with Delta I/I = 15.1 +/- 1.7%. We adopt a specific model for submillimeter sources in which the contribution to clustering comes from sources in the redshift ranges 1.3 <= z <= 2.2, 1.5 <= z <= 2.7, and 1.7 <= z <= 3.2, at 250, 350, and 500 microns, respectively. With these distributions, our measurement of the power spectrum, P(k_theta), corresponds to linear bias parameters, b = 3.8 +/- 0.6, 3.9 +/- 0.6 and 4.4 +/- 0.7, respectively. We further interpret the results in terms of the halo model, and find that at the smaller scales, the simplest halo model fails to fit our results. One way to improve the fit is to increase the radius at which dark matter halos are artificially truncated in the model, which is equivalent to having some star-forming galaxies at z >= 1 located in the outskirts of groups and clusters. In the context of this model we find a minimum halo mass required to host a galaxy is log (M_min / M_sun) = 11.5 (+0.4/-0.1), and we derive effective biases $b_eff = 2.2 +/- 0.2, 2.4 +/- 0.2, and 2.6 +/- 0.2, and effective masses log (M_eff / M_sun) = 12.9 +/- 0.3, 12.8 +/- 0.2, and 12.7 +/- 0.2, at 250, 350, and 500 microns, corresponding to spatial correlation lengths of r_0 = 4.9, 5.0, and 5.2 +/- 0.7 h^-1 Mpc, respectively. Finally, we discuss implications for clustering measurement strategies with Herschel and Planck.

Over the course of two flights, the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) made resolved maps of seven nearby (<25 Mpc) galaxies at 250, 350, and 500 microns. During its June 2005 flight from Sweden, BLAST observed a single nearby galaxy, NGC 4565. During the December 2006 flight from Antarctica, BLAST observed the nearby galaxies NGC 1097, NGC 1291, NGC 1365, NGC 1512, NGC 1566, and NGC 1808. We fit physical dust models to a combination of BLAST observations and other available data for the galaxies observed by Spitzer. We fit a modified blackbody to the remaining galaxies to obtain total dust mass and mean dust temperature. For the four galaxies with Spitzer data, we also produce maps and radial profiles of dust column density and temperature. We measure the fraction of BLAST detected flux originating from the central cores of these galaxies and use this to calculate a "core fraction," an upper limit on the "AGN fraction" of these galaxies. We also find our resolved observations of these galaxies give a dust mass estimate 5-19 times larger than an unresolved observations would predict. Finally, we are able to use these data to derive a value for the dust mass absorption co-efficient of kappa = 0.29 +/-0.03 m^2 kg^-1 at 250 microns. This study is an introduction to future higher-resolution and higher-sensitivity studies to be conducted by Herschel and SCUBA-II.

We report multi-wavelength power spectra of diffuse Galactic dust emission from BLAST observations at 250, 350, and 500 microns in Galactic Plane fields in Cygnus X and Aquila. These submillimeter power spectra statistically quantify the self-similar structure observable over a broad range of scales and can be used to assess the cirrus noise which limits the detection of faint point sources. The advent of submillimeter surveys with the Herschel Space Observatory makes the wavelength dependence a matter of interest. We show that the observed relative amplitudes of the power spectra can be related through a spectral energy distribution (SED). Fitting a simple modified black body to this SED, we find the dust temperature in Cygnus X to be 19.9 +/- 1.3 K and in the Aquila region 16.9 +/- 0.7 K. Our empirical estimates provide important new insight into the substantial cirrus noise that will be encountered in forthcoming observations.

We directly measure redshift evolution in the mean physical properties (far-infrared luminosity, temperature, and mass) of the galaxies that produce the cosmic infrared background (CIB), using measurements from the Balloon-borne Large Aperture Sub-millimeter Telescope (BLAST), and Spitzer which constrain the CIB emission peak. This sample is known to produce a surface brightness in the BLAST bands consistent with the full CIB, and photometric redshifts are identified for all of the objects. We find that most of the 70 micron background is generated at z <~ 1 and the 500 micron background generated at z >~ 1. A significant growth is observed in the mean luminosity from ~ 10^9 - 10^12 L_sun, and in the mean temperature by 10 K, from redshifts 0< z < 3. However, there is only weak positive evolution in the comoving dust mass in these galaxies across the same redshift range. We also measure the evolution of the far-infrared luminosity density, and the star-formation rate history for these objects, finding good agreement with other infrared studies up to z ~1, exceeding the contribution attributed to optically-selected galaxies.

We present the 250, 350, and 500 micron detection of bright submillimeter emission in the direction of the Bullet Cluster measured by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST). The 500 micron centroid is coincident with an AzTEC 1.1 mm point-source detection at a position close to the peak lensing magnification produced by the cluster. However, the 250 micron and 350 micron centroids are elongated and shifted toward the south with a differential shift between bands that cannot be explained by pointing uncertainties. We therefore conclude that the BLAST detection is likely contaminated by emission from foreground galaxies associated with the Bullet Cluster. The submillimeter redshift estimate based on 250-1100 micron photometry at the position of the AzTEC source is z_phot = 2.9 (+0.6 -0.3), consistent with the infrared color redshift estimation of the most likely IRAC counterpart. These flux densities indicate an apparent far-infrared luminosity of L_FIR = 2E13 Lsun. When the amplification due to the gravitational lensing of the cluster is removed, the intrinsic far-infrared luminosity of the source is found to be L_FIR <= 10^12 Lsun, consistent with typical luminous infrared galaxies.

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) has made one square degree, deep, confusion limited maps at three different bands, centered on the Great Observatories Origins Deep Survey South field. By calculating the covariance of these maps with catalogs of 24 micron sources from the Far-Infrared Deep Extragalactic Legacy Survey (FIDEL), we have determined that the total submillimeter intensities are 8.60 +/- 0.59, 4.93 +/- 0.34, and 2.27 +/- 0.20 nW m^2 sr^(-1) at 250, 350, and 500 micron, respectively. These numbers are more precise than previous estimates of the cosmic infrared background (CIB) and are consistent with 24 micron-selected galaxies generating the full intensity of the CIB. We find that the fraction of the CIB that originates from sources at z >= 1.2 increases with wavelength, with 60% from high redshift sources at 500 micron. At all BLAST wavelengths, the relative intensity of high-z sources is higher for 24 micron-faint sources than it is for 24 micron-bright sources. Galaxies identified as active galactic nuclei (AGN) by their Infrared Array Camera (IRAC) colors are 1.6-2.6 times brighter than the average population at 250-500 micron, consistent with what is found for X-ray-selected AGN. BzK-selected galaxies are found to be moderately brighter than typical 24 micron-selected galaxies in the BLAST bands. These data provide high precision constraints for models of the evolution of the number density and intensity of star forming galaxies at high redshift.

We present first results from an unbiased 50 deg^2 submillimeter Galactic survey at 250, 350, and 500 micron from the 2006 flight of the Balloon-borne Large Aperture Submillimeter Telescope (BLAST). The map has resolution ranging from 36 arcsec to 60 arcsec in the three submillimeter bands spanning the thermal emission peak of cold starless cores. We determine the temperature, luminosity, and mass of more than one thousand compact sources in a range of evolutionary stages and an unbiased statistical characterization of the population. From comparison with C^(18)O data, we find the dust opacity per gas mass, kappa r = 0.16 cm^2 g^(-1) at 250 micron, for cold clumps. We find that 2% of the mass of the molecular gas over this diverse region is in cores colder than 14 K, and that the mass function for these cold cores is consistent with a power law with index alpha = -3.22 +/- 0.14 over the mass range 14 M_sun < M < 80 M_sun. Additionally, we infer a mass-dependent cold core lifetime of t_c(M) = 4E6 (M/20 M_sun)^(-0.9) years - longer than what has been found in previous surveys of either low or high mass cores, and significantly longer than free fall or likely turbulent decay times. This implies some form of non-thermal support for cold cores during this early stage of star formation.

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) operated successfully during a 250-hour flight over Antarctica in December 2006 (BLAST06). As part of the calibration and pointing procedures, the red hypergiant star VY CMa was observed and used as the primary calibrator. Details of the overall BLAST06 calibration procedure are discussed. The 1-sigma absolute calibration is accurate to 10, 12, and 13% at the 250, 350, and 500 micron bands, respectively. The errors are highly correlated between bands resulting in much lower error for the derived shape of the 250-500 micron continuum. The overall pointing error is <5" rms for the 36, 42, and 60" beams. The performance of the optics and pointing systems is discussed.

Recently it has been proposed that there are two types of SN Ia progenitors -- short-lived and long-lived. On the basis of this idea, we develope a theory of a unified mechanism for the formation of the bimodal radial distribution of iron and oxygen in the Galactic disc. The underlying cause for the formation of the fine structure of the radial abundance pattern is the influence of spiral arms, specifically, the combined effect of the corotation resonance and turbulent diffusion. From our modelling we conclude that to explain the bimodal radial distributions simultaneously for oxygen and iron and to obtain approximately equal total iron output from different types of supernovae, the mean ejected iron mass per supernova event should be the same as quoted in literature if maximum mass of stars, that eject heavy elements, is $50 M_{\odot}$. For the upper mass limit of $70 M_{\odot}$ the production of iron by a supernova II explosion should be increased by about 1.5 times.

We have identified radio and/or mid-infrared counterparts to 198 out of 350 sources detected at >=5 sigma over ~ 9 square degrees centered on the Chandra Deep Field South (CDFS) by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) at 250, 350 and 500 um. We have matched 114 of these counterparts to optical sources with previously derived photometric redshifts and fitted SEDs to the BLAST fluxes and fluxes at 70 and 160 um acquired with the Spitzer Space Telescope. In this way, we have constrained dust temperatures, total far-infrared/sub-millimeter luminosities and star formation rates for each source. Our findings show that on average, the BLAST sources lie at significantly lower redshifts and have significantly lower rest-frame dust temperatures compared to submm sources detected in surveys conducted at 850 um. We demonstrate that an apparent increase in dust temperature with redshift in our sample arises as a result of selection effects. Finally, we provide the full multi-wavelength catalog of >= 5 sigma BLAST sources contained within the complete ~ 9 square degree survey area.

Submillimetre surveys during the past decade have discovered a population of luminous, high-redshift, dusty starburst galaxies. In the redshift range 1 <= z <= 4, these massive submillimetre galaxies go through a phase characterized by optically obscured star formation at rates several hundred times that in the local Universe. Half of the starlight from this highly energetic process is absorbed and thermally re-radiated by clouds of dust at temperatures near 30 K with spectral energy distributions peaking at 100 microns in the rest frame. At 1 <= z <= 4, the peak is redshifted to wavelengths between 200 and 500 microns. The cumulative effect of these galaxies is to yield extragalactic optical and far-infrared backgrounds with approximately equal energy densities. Since the initial detection of the far-infrared background (FIRB), higher-resolution experiments have sought to decompose this integrated radiation into the contributions from individual galaxies. Here we report the results of an extragalactic survey at 250, 350 and 500 microns. Combining our results at 500 microns with those at 24 microns, we determine that all of the FIRB comes from individual galaxies, with galaxies at z >= 1.2 accounting for 70 per cent of it. As expected, at the longest wavelengths the signal is dominated by ultraluminous galaxies at z > 1.

In this study we demonstrate for the first time that the unified Monte Carlo approach can be applied to model gas-grain chemistry in large reaction networks. Specifically, we build a time-dependent gas-grain chemical model of the interstellar medium, involving about 6000 gas-phase and 200 grain surface reactions. This model is used to test the validity of the standard and modified rate equation methods in models of dense and translucent molecular clouds and to specify under which conditions the use of the stochastic approach is desirable. We found that at temperatures 25--30 K gas-phase abundances of H$_2$O, NH$_3$, CO and many other gas-phase and surface species in the stochastic model differ from those in the deterministic models by more than an order of magnitude, at least, when tunneling is accounted for and/or diffusion energies are 3x lower than the binding energies. In this case, surface reactions, involving light species, proceed faster than accretion of the same species. In contrast, in the model without tunneling and with high binding energies, when the typical timescale of a surface recombination is greater than the timescale of accretion onto the grain, we obtain almost perfect agreement between results of Monte Carlo and deterministic calculations in the same temperature range. At lower temperatures ($\sim10$ K) gaseous and, in particular, surface abundances of most important molecules are not much affected by stochastic processes.

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 investigate general aspects of molecular line formation under conditions which are typical of prestellar cores. Focusing on simple linear molecules, we study formation of their rotational lines by radiative transfer simulations. We present a thermalization diagram to show the effects of collisions and radiation on the level excitation. We construct a detailed scheme (contribution chart) to illustrate the formation of emission line profiles. This chart can be used as an efficient tool to identify which parts of the cloud contribute to a specific line profile. We show how molecular line characteristics for uniform model clouds depend on hydrogen density, molecular column density, and kinetic temperature. The results are presented in a 2D plane to illustrate cooperative effects of the physical factors. We also use a core model with a non-uniform density distribution and chemical stratification to study the effects of cloud contraction and rotation on spectral line maps. We discuss the main issues that should be taken into account when dealing with interpretation and simulation of observed molecular lines.

We consider current state of star formation theory and requirements to observations in millimeter and submillimeter ranges which are necessary for resolution of the most actual problems of the physics of star formation. Two key features of star-forming regions which define observational requirements to their studies, are relatively low energy of processes that take place there and smallness of corresponding spatial scales. This is especially true for the objects in the latest stages of ``pre-stellar'' evolution, that is, hot cores, hyper- and ultracompact HII regions, and protoplanetary disks. Angular resolution, sensitivity, and spectral coverage in existing projects of ground-based and space telescopes of submillimeter and millimeter range are not completely adequate to necessary requirements. To obtain detailed information on star-forming regions as well as on individual protostars it is necessary to employ a space-based interferometer.

Spider is a balloon-borne experiment that will measure the polarization of the Cosmic Microwave Background over a large fraction of a sky at 1 degree resolution. Six monochromatic refracting millimeter-wave telescopes with large arrays of antenna-coupled transition-edge superconducting bolometers will provide system sensitivities of 4.2 and 3.1 micro K_cmb rt s at 100 and 150 GHz, respectively. A rotating half-wave plate will modulate the polarization sensitivity of each telescope, controlling systematics. Bolometer arrays operating at 225 GHz and 275 GHz will allow removal of polarized galactic foregrounds. In a 2-6 day first flight from Alice Springs, Australia in 2010, Spider will map 50% of the sky to a depth necessary to improve our knowledge of the reionization optical depth by a large factor.

We present a picture of star formation around the HII region Sh2-235 (S235) based upon data on the spatial distribution of young stellar clusters and the distribution and kinematics of molecular gas around S235. We observed 13CO(1-0) and CS(2-1) emission toward S235 with the Onsala Space Observatory 20-m telescope and analysed the star density distribution with archival data from the 2MASS survey. Dense molecular gas forms a shell-like structure at the south-eastern part of S235. The young clusters found with 2MASS data are embedded in this shell. The positional relationship of the clusters, the molecular shell and the HII region indicates that expansion of S235 is responsible for the formation of the clusters. The gas distribution in the S235 molecular complex is clumpy, which hampers interpretation exclusively on the basis of the morphology of the star forming region. We use data on kinematics of molecular gas to support the hypothesis of induced star formation, and distinguish three basic types of molecular gas components. The first type is primordial undisturbed gas of the giant molecular cloud, the second type is gas entrained in motion by expansion of the HII region (this is where the embedded clusters were formed), and the third type is a fast-moving gas, which might have been accelerated by winds from the newly formed clusters. The clumpy distribution of molecular gas and its kinematics around the HII region implies that the picture of triggered star formation around S235 can be a mixture of at least two possibilities: the "collect-and-collapse" scenario and the compression of pre-existing dense clumps by the shock wave.

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) operated successfully during a 100-hour flight from northern Sweden in June 2005 (BLAST05). As part of the calibration and pointing procedures, several compact sources were mapped, including solar system, Galactic, and extragalactic targets, specifically Pallas, CRL 2688, LDN 1014, IRAS 20126+4104, IRAS 21078+5211, IRAS 21307+5049, IRAS 22134+5834, IRAS 23011+6126, K3-50, W 75N, and Mrk 231. One additional source, Arp 220, was observed and used as our primary calibrator. Details of the overall BLAST05 calibration procedure are discussed here. The BLAST observations of each compact source are described, flux densities and spectral energy distributions are reported, and these are compared with previous measurements at other wavelengths. The 250, 350, and 500 um BLAST data can provide useful constraints to the amplitude and slope of the submillimeter continuum, which in turn may be useful for the improved calibration of other submillimeter instruments.

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) is a sub-orbital surveying experiment designed to study the evolutionary history and processes of star formation in local galaxies (including the Milky Way) and galaxies at cosmological distances. The BLAST continuum camera, which consists of 270 detectors distributed between 3 arrays, observes simultaneously in broad-band (30%) spectral-windows at 250, 350, and 500 microns. The optical design is based on a 2m diameter telescope, providing a diffraction-limited resolution of 30" at 250 microns. The gondola pointing system enables raster mapping of arbitrary geometry, with a repeatable positional accuracy of ~30"; post-flight pointing reconstruction to ~5" rms is achieved. The on-board telescope control software permits autonomous execution of a pre-selected set of maps, with the option of manual override. In this paper we describe the primary characteristics and measured in-flight performance of BLAST. BLAST performed a test-flight in 2003 and has since made two scientifically productive long-duration balloon flights: a 100-hour flight from ESRANGE (Kiruna), Sweden to Victoria Island, northern Canada in June 2005; and a 250-hour, circumpolar-flight from McMurdo Station, Antarctica, in December 2006.

We describe a map-making method which we have developed for the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) experiment, but which should have general application to data from other submillimeter arrays. Our method uses a Maximum Likelihood based approach, with several approximations, which allows images to be constructed using large amounts of data with fairly modest computer memory and processing requirements. This new approach, Signal And Noise Estimation Procedure Including Correlations (SANEPIC), builds upon several previous methods, but focuses specifically on the regime where there is a large number of detectors sampling the same map of the sky, and explicitly allowing for the the possibility of strong correlations between the detector timestreams. We provide real and simulated examples of how well this method performs compared with more simplistic map-makers based on filtering. We discuss two separate implementations of SANEPIC: a brute-force approach, in which the inverse pixel-pixel covariance matrix is computed; and an iterative approach, which is much more efficient for large maps. SANEPIC has been successfully used to produce maps using data from the 2005 BLAST flight.

We present the first results from a new 250, 350, and 500 micron Galactic Plane survey taken with the Balloon-borne Large-Aperture Submillimeter Telescope (BLAST) in 2005. This survey's primary goal is to identify and characterize high-mass proto-stellar objects (HMPOs). The region studied here covers 4 sq. deg near the open cluster NGC 6823 in the constellation Vulpecula (l=59). We find 60 compact sources (<60'' diameter) detected simultaneously in all three bands. Their spectral energy distributions (SEDs) are constrained through BLAST, IRAS, Spitzer MIPS, and MSX photometry, with inferred dust temperatures spanning ~12-40K assuming a dust emissivity index beta=1.5. The luminosity-to-mass ratio, a distance-independent quantity, spans ~0.2-130 L_\odot M_\odot^{-1}. Distances are estimated from coincident 13CO (1->0) velocities combined with a variety of other velocity and morphological data in the literature. In total, 49 sources are associated with a molecular cloud complex encompassing NGC 6823 (distance ~2.3kpc), 10 objects with the Perseus Arm (~8.5kpc) and one object is probably in the outer Galaxy (~14kpc). Near NGC 6823, the inferred luminosities and masses of BLAST sources span ~40-10^4 L_\odot, and ~15-700 M_\odot, respectively. The mass spectrum is compatible with molecular gas masses in other high-mass star forming regions. Several luminous sources appear to be Ultra Compact HII regions powered by early B stars. However, many of the objects are cool, massive gravitationally-bound clumps with no obvious internal radiation from a protostar, and hence excellent HMPO candidates.

Spider is a long-duration, balloon-borne polarimeter designed to measure large scale Cosmic Microwave Background (CMB) polarization with very high sensitivity and control of systematics. The instrument will map over half the sky with degree angular resolution in I, Q and U Stokes parameters, in four frequency bands from 96 to 275 GHz. Spider's ultimate goal is to detect the primordial gravity wave signal imprinted on the CMB B-mode polarization. One of the challenges in achieving this goal is the minimization of the contamination of B-modes by systematic effects. This paper explores a number of instrument systematics and observing strategies in order to optimize B-mode sensitivity. This is done by injecting realistic-amplitude, time-varying systematics in a set of simulated time-streams. Tests of the impact of detector noise characteristics, pointing jitter, payload pendulations, polarization angle offsets, beam systematics and receiver gain drifts are shown. Spider's default observing strategy is to spin continuously in azimuth, with polarization modulation achieved by either a rapidly spinning half-wave plate or a rapidly spinning gondola and a slowly stepped half-wave plate. Although the latter is more susceptible to systematics, results shown here indicate that either mode of operation can be used by Spider.

Molecular line observations of starless (prestellar) cores combined with a chemical evolution modeling and radiative transfer calculations are a powerful tool to study the earliest stages of star formation. However, conclusions drawn from such a modeling may noticeably depend on the assumed thermal structure of the cores. The assumption of isothermality, which may work well in chemo-dynamical studies, becomes a critical factor in molecular line formation simulations. We argue that even small temperature variations, which are likely to exist in starless cores, can have a non-negligible effect on the interpretation of molecular line data and derived core properties. In particular, ``chemically pristine'' isothermal cores (low depletion) can have centrally peaked C$^{18}$O and C$^{34}$S radial intensity profiles, while having ring-like intensity distributions in models with a colder center and/or warmer envelope assuming the same underlying chemical structure. Therefore, derived molecular abundances based on oversimplified thermal models may lead to a mis-interpretation of the line data.

The anisotropy of velocities in MHD turbulence is demonstrated explicitly by calculating the velocity gradients as a function of direction in representative simulations of decaying turbulence. It follows that the optical depths of spectral lines are anisotropic when there is MHD turbulence, and that this anisotropy influences the polarization characteristics of the emergent radiation. We calculate the linear polarization that results for the microwave lines of the CO molecule in star-forming gas and show that it is comparable to the polarization that is observed. This and our earlier result--that the anisotropy of MHD turbulence may be the cause for the absence of the Zeeman $\pi $-components in the spectra of OH mainline masers--are the first demonstrations of the occurrence of anisotropy in the optical depths caused by MHD turbulence. A non-local approximation is developed for the radiative transfer and the results are compared with those from a local (LVG) approximation.

This is the first paper in a series where we study the influence of turbulent diffusion and advective transport on the chemical evolution of protoplanetary disks, using a 2D flared disk model and a 2D mixing gas-grain chemical code with surface reactions. A first interesting result concerns the abundance of gas-phase CO in the outer regions of protoplanetary disks. In this Letter we argue that the gas-phase CO concentration in the disk regions, where the temperature is lower than $\sim25$ K, can be significantly enhanced due to the combined effect of vertical and radial mixing. This finding has a potential implication for the current observational data on the DM~Tau disk chemistry.

We have developed two redundant daytime star cameras to provide the fine pointing solution for the balloon-borne submillimeter telescope, BLAST. The cameras are capable of providing a reconstructed pointing solution with an absolute accuracy < 5 arcseconds. They are sensitive to stars down to magnitudes ~ 9 in daytime float conditions. Each camera combines a 1 megapixel CCD with a 200 mm f/2 lens to image a 2 degree x 2.5 degree field of the sky. The instruments are autonomous. An internal computer controls the temperature, adjusts the focus, and determines a real-time pointing solution at 1 Hz. The mechanical details and flight performance of these instruments are presented.

We present a detailed theoretical study of the isolated Bok globule CB17 (L1389) based on spectral maps of CS, HCO$^+$, C$^{18}$O, C$^{34}$S, and H$^{13}$CO$^+$ lines. A phenomenological model of prestellar core evolution, a time-dependent chemical model, and a radiative transfer simulation for molecular lines are combined to reconstruct the chemical and kinematical structure of this core. We developed a general criterion that allows to quantify the difference between observed and simulated spectral maps. By minimizing this difference, we find that very high and very low values of the effective sticking probability $S$ are not appropriate for the studied prestellar core. The most probable $S$ value for CB17 is 0.3--0.5. The spatial distribution of the intensities and self-absorption features of optically thick lines is indicative of UV irradiation of the core. By fitting simultaneously optically thin and optically thick transitions, we isolate the model that reproduces all the available spectral maps to a reasonable accuracy. The line asymmetry pattern in CB17 is reproduced by a combination of infall, rotation, and turbulent motions with velocities $\sim0.05$ km s$^{-1}$, $\sim0.1$ km s$^{-1}$, and $\sim0.1$ km s$^{-1}$, respectively. These parameters corresponds to energy ratios $E_{\rm rot}/E_{\rm grav}\approx0.03$, $E_{\rm therm}/E_{\rm grav}\approx0.8$, and $E_{\rm turb}/E_{\rm grav}\approx0.05$ (the rotation parameters are determined for $i=90^\circ$). The chemical age of the core is about 2 Myrs. In particular, this is indicated by the central depletion of CO, CS, and HCO$^+$. Based on the angular momentum value, we argue that the core is going to fragment, i.e., to form a binary (multiple) star. (abridged)

(abridged) The polarized thermal radiation at far infrared and submillimeter wavelengths from dust grains in interstellar clouds with irregular magnetic fields is simulated. We show that the polarization properties of quiescent clouds and of star-forming regions are consistent with Kolmogorov-like turbulent magnetic fields that are comparable in magnitude to the uniform component of the magnetic fields. If the beam size is much smaller than the correlation length $L_{\rm corr}$ of the fields, the calculated percentage polarization $p$ decreases to an asymptotic value when the number of correlation lengths $N_{\rm corr}$ through the cloud exceeds a few tens. For these values of $N_{\rm corr}$, the dispersion in the position angles $\sigma_\alpha$ is still appreciable--decreasing to only about $20^{\circ}$. However, when the finite size of a telescope beam is taken into account, the asymptotic value of $p$ is reached for fewer correlation lengths (smaller $N_{\rm corr}$) due to averaging over the beam; $\sigma_\alpha$ becames much smaller and consistent with the observational data. The smoothing of the polarization properties due to the combined effect of the thickness of the cloud and the finite size of the beam can be described by a single variable which we designate as the generalized number of correlation lengths. In addition, we study various factors that may contribute to the decrease in the linear polarization percentage with increasing intensity that is observed at submillimeter and far infrared wavelengths in many, though not in all, dark clouds (the "polarization hole" effect).

(abridged) We analyse the evolution of the fractional ionisation in a steady-state protoplanetary disc with a vertical temperature gradient and with gas-grain chemistry including surface reactions. The ionisation due to stellar X-rays, stellar and interstellar UV radiation, cosmic rays and radionuclide decay is taken into account. Using our reduction schemes as a tool for the analysis, we isolate small sets of chemical reactions that reproduce the evolution of the ionisation degree at representative disc locations with an accuracy of 50%-100%. Column densities of key molecules are calculated and compared to the results of other recent studies and observational data.

We consider astrophysical maser radiation that is created in the presence of mildly supersonic, magnetohydrodynamic (MHD) turbulence. The focus is on the OH masers for which the magnetic field is strong enough that the separations of the Zeeman components are greater than the spectral linebreadths. A longstanding puzzle has been the absence of the Zeeman $\pi $ components and the high circular polarization in the observed spectra of these masers. We first argue that the elongation of eddies along the field that has recently been recognized in MHD turbulence will enhance the optical depth parallel to the magnetic field in comparison with that perpendicular to the magnetic field. We then simulate maser emission with a numerical model of MHD turbulence to demonstrate quantitatively how the intensities of the linearly polarized $\pi $ components are suppressed and the intensities of the nearly circularly polarized $\sigma $ components are enhanced. This effect is also generic in the sense that most spectral lines in MHD turbulence with Mach number $M\sim 1$ should have larger optical depths parallel to the magnetic field than perpendicular. The effect is reduced considerably when $M<1$. The simulations also demonstrate that the velocity and magnetic field variations due to the turbulence can (but do not necessarily) cause one of the $\sigma $ components to be much more intense than the other, as is often observed for mainline OH masers.

With the chemical reaction rate database UMIST95 (Millar et al. 1997) we analyze how uncertainties in rate constants of gas-phase chemical reactions influence the modelling of molecular abundances in the interstellar medium. Random variations are introduced into the rate constants to estimate the scatter in theoretical abundances. Calculations are performed for dark and translucent molecular clouds where gas phase chemistry is adequate. Similar approach was used by Pineau des Forets & Roueff (2000) for the study of chemical bistability. All the species are divided into 6 sensitivity groups according to the value of the scatter in their model abundances computed with varied rate constants. It is shown that the distribution of species within these groups depends on the number of atoms in a molecule and on the adopted physical conditions. The simple method is suggested which allows to single out reactions that are most important for the evolution of a given species.

We present a new method to analyse and reduce chemical networks and apply this technique to the chemistry in molecular clouds. Using the technique, we investigated the possibility of reducing the number of chemical reactions and species in the UMIST 95 database simultaneously. In addition, we did the same reduction but with the ``objective technique'' in order to compare both methods. We found that it is possible to compute the abundance of carbon monoxide and fractional ionisation accurately with significantly reduced chemical networks in the case of pure gas-phase chemistry. For gas-grain chemistry involving surface reactions reduction is not worthwhile. Compared to the ``objective technique'' our reduction method is more effective but more time-consuming as well.

We develop a detailed chemical model for the starless cores of strongly magnetized molecular clouds, with the ambipolar diffusion-driven dynamic evolution of the clouds coupled to the chemistry through ion abundances. We concentrate on two representative model clouds in this initial study, one with magnetic fields and the other without. The model predictions on the peak values and spatial distributions of the column densities of CO, CCS, N$_2$H$^+$ and HCO$^+$ are compared with those observationally inferred for the well-studied starless core L1544, which is thought to be on the verge of star formation. We find that the magnetic model, in which the cloud is magnetically supported for several million years before collapsing dynamically, provides a reasonable overall fit to the available data on L1544; the fit is significantly worse for the non-magnetic model, in which the cloud collapses promptly. The observed large peak column density for N$_2$H$^+$ and clear central depression for CCS favor the magnetically-retarded collapse over the free-fall collapse. A relatively high abundance of CCS is found in the magnetic model, resulting most likely from an interplay of depletion and late-time hydrocarbon chemistry enhanced by CO depletion. These initial results lend some support to the standard picture of dense core formation in strongly magnetized clouds through ambipolar diffusion. They are at variance with those of Aikawa et al. (2001) who considered a set of models somewhat different from ours and preferred one in which the cloud collapses more or less freely for L1544.

Calculations are performed for the circular polarization of maser radiation from a turbulent, Keplerian disk that is intended to represent the sub-parsec disk at the nucleus of the galaxy NGC4258. The polarization in the calculations is a result of the Zeeman effect in the regime in which the Zeeman splitting is much less than the spectral linebreadth. Plausible configurations for turbulent magnetic and velocity fields in the disk are created by statistical methods. This turbulence, along with the Keplerian velocity gradients and the blending of the three hyperfine components to form the $6_{16} - 5_{23}$ masing transition of water, are key ingredients in determining the appearance of the polarized spectra that are calculated. These spectra are quite different from the polarized spectra that would be expected for a two-level transition where there is no hyperfine structure. The effect of the hyperfine structure on the polarization is most striking in the calculations for the maser emission that represents the central (or systemic) features of NGC4258. Information about magnetic fields is inferred from observations for polarized maser radiation and bears on the structure of accretion disks.

Calculations are performed for the linear polarization of starlight due to extinction by aligned dust grains when the starlight traverses a medium with irregular magnetic fields. This medium is intended to represent the optically thick components of interstellar clouds which are observed to make little, if any, contribution to the polarization of starlight. In agreement with the observations, we find that the polarization properties of the starlight---the average fractional polarization and the dispersion in position angles---can be essentially unchanged. For this, the rms of the irregular component must be greater than the average magnetic field, which in turn tends to imply that the turbulent velocities in these interstellar clouds are super-Alfvenic.

The strengths of magnetic fields in interstellar gas clouds are obtained through observations of the circular polarization of spectral line radiation. Irregularities in this magnetic field may be present due to turbulence, waves or perhaps other causes, and may play an essential role in the structure and evolution of the gas clouds. To infer information about these irregularities from the observational data, we develop statistical relationships between the rms values of the irregular component of the magnetic field and spatial variations in the circular polarization of the spectral line radiation. The irregularities are characterized in analogy with descriptions of turbulence---by a sum of Fourier waves having a power spectrum with a slope similar to that of Kolmogorov turbulence. For comparison, we also perform computations in which turbulent magnetic and velocity fields from representative MHD simulations by others are utilized. Although the effects of the variations about the mean value of the magnetic field along the path of a ray tend to cancel, a significant residual effect in the polarization of the emergent radiation remains for typical values of the relevant parameters. A map of observed spectra of the 21 cm line toward Orion A is analyzed and the results are compared with our calculations in order to infer the strength of the irregular component of the magnetic field. The rms of the irregular component is found to be comparable in magnitude to the mean magnetic field within the cloud. Hence, the turbulent and Alfven velocities should also be comparable.

Elliptical galaxies are often assumed to be the primary source of heavy elements in the intracluster medium (ICM), with the contribution of other morphological types being negligible. In this paper we argue that a role of spiral galaxies in the chemical evolution of the ICM is also important. This statement rests upon our recent calculations of the heavy element loss from a disk galaxy through the hot steady-state galactic wind and dust grains expulsion by stellar radiation pressure. We show that a typical disk galaxy is nearly as effective in enriching the ICM as an elliptical galaxy of the same mass. Having estimated the oxygen and iron loss from a single galaxy, we integrate them over the galactic mass spectrum. We show that the "effective" loss (per unit luminosity) from spiral galaxies is comparable to the loss from ellipticals. The dominant role of early-type galaxies in rich clusters is caused by that they outnumber spirals.

The linear polarization that is observed, together with likely changes in the orientation of the magnetic field along the line of sight and hence of the optical axes of the medium, can lead to the circular polarization that is observed in the radiation of the circumstellar SiO masers. A magnetic field greater than only about 30 mG is required, in contrast to 10-100 G that would be implied by the Zeeman interpretation. To assess quantitatively the likely changes in orientation of the magnetic field, calculations are performed with representative field configurations that are created by statistical sampling using a Kolmogorov-like power spectrum.