In this work, we present Alfnoor, a dedicated tool optimised for population studies of exoplanet atmospheres. Alfnoor combines the latest version of the retrieval algorithm TauREx 3, with the instrument noise simulator ArielRad and enables the simultaneous retrieval analysis of a large sample of exo-atmospheres. We applied this tool to the Ariel list of planetary candidates and focus on hydrogen dominated, cloudy atmospheres observed in transit with the Tier-2 mode (medium Ariel resolution). As a first experiment, we randomised the abundances - ranging from 10$^{-7}$ to 10$^{-2}$ - of the trace gases, which include H$_2$O, CH$_4$, CO, CO$_2$ and NH$_3$. This exercise allowed to estimate the detection limits for Ariel Tier-2 and Tier-3 modes when clouds are present. In a second experiment, we imposed an arbitrary trend between a chemical species and the effective temperature of the planet. A last experiment was run requiring molecular abundances being dictated by equilibrium chemistry at a certain temperature. Our results demonstrate the ability of Ariel Tier-2 and Tier-3 surveys to reveal trends between the chemistry and associated planetary parameters. Future work will focus on eclipse data, on atmospheres heavier than hydrogen and will be applied also to other observatories.

A new generation of exoplanet research beckons and with it the need for simulation tools that accurately predict signal and noise in transit spectroscopy observations. We developed ExoSim: an end-to-end simulator that models noise and systematics in a dynamical simulation. ExoSim improves on previous simulators in the complexity of its simulation, versatility of use and its ability to be generically applied to different instruments. It performs a dynamical simulation that can capture temporal effects, such as correlated noise and systematics on the light curve. It has also been extensively validated, including against real results from the Hubble WFC3 instrument. We find ExoSim is accurate to within 5% in most comparisons. ExoSim can interact with other models which simulate specific time-dependent processes. A dedicated star spot simulator allows ExoSim to produce simulated observations that include spot and facula contamination. ExoSim has been used extensively in the Phase A and B design studies of the ARIEL mission, and has many potential applications in the field of transit spectroscopy.

Although there exists a large sample of known exoplanets, little spectroscopic data exists that can be used to study their global atmospheric properties. This deficiency can be addressed by performing phase-resolved spectroscopy -- continuous spectroscopic observations of a planet's entire orbit about its host star -- of transiting exoplanets. Planets with characteristics suitable for atmospheric characterization have orbits of several days, thus phase curve observations are highly resource intensive, especially for shared use facilities. In this work, we show that an infrared spectrograph operating from a high altitude balloon platform can perform phase-resolved spectroscopy of hot Jupiter-type exoplanets with performance comparable to a space-based telescope. Using the EXoplanet Climate Infrared TElescope (EXCITE) experiment as an example, we quantify the impact of the most important systematic effects that we expect to encounter from a balloon platform. We show an instrument like EXCITE will have the stability and sensitivity to significantly advance our understanding of exoplanet atmospheres. Such an instrument will both complement and serve as a critical bridge between current and future space-based near infrared spectroscopic instruments.

Ariel has been selected as ESA's M4 mission for launch in 2028 and is designed for the characterisation of a large and diverse population of exoplanetary atmospheres to provide insights into planetary formation and evolution within our Galaxy. Here we present a study of Ariel's capability to observe currently-known exoplanets and predicted TESS discoveries. We use the Ariel Radiometric model (ArielRad) to simulate the instrument performance and find that ~2000 of these planets have atmospheric signals which could be characterised by Ariel. This list of potential planets contains a diverse range of planetary and stellar parameters. From these we select an example Mission Reference Sample (MRS), comprised of 1000 diverse planets to be completed within the primary mission life, which is consistent with previous studies. We also explore the mission capability to perform an in-depth survey into the atmospheres of smaller planets, which may be enriched or secondary. Earth-sized planets and Super-Earths with atmospheres heavier than H/He will be more challenging to observe spectroscopically. However, by studying the time required to observe ~110 Earth-sized/Super-Earths, we find that Ariel could have substantial capability for providing in-depth observations of smaller planets. Trade-offs between the number and type of planets observed will form a key part of the selection process and this list of planets will continually evolve with new exoplanet discoveries replacing predicted detections. The Ariel target list will be constantly updated and the MRS re-selected to ensure maximum diversity in the population of planets studied during the primary mission life.

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

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

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

The Twinkle space telescope has been designed for the characterisation of exoplanets and Solar System objects. Operating in a low Earth, Sun-synchronous orbit, Twinkle is equipped with a 45 cm telescope and visible (0.4 - 1um) and infrared (1.3 - 4.5um) spectrometers which can be operated simultaneously. Twinkle is a general observatory which will provide on-demand observations of a wide variety of targets within wavelength ranges that are currently not accessible using other space telescopes or accessible only to oversubscribed observatories in the short-term future. Here we explore the ability of Twinkle's spectrometers to characterise the currently-known exoplanets. We study the spectral resolution achievable by combining multiple observations for various planetary and stellar types. We also simulate spectral retrievals for some well-known planets (HD 209458 b, GJ 3470 b and 55 Cnc e). From the exoplanets known today, we find that with a single transit or eclipse, Twinkle could probe 89 planets at low spectral resolution (R < 20) as well as 12 planets at higher resolution (R > 20) in channel 1 (1.3 - 4.5um). With 10 observations, the atmospheres of 144 planets could be characterised with R < 20 and 81 at higher resolutions. Upcoming surveys will reveal thousands of new exoplanets, many of which will be located within Twinkle's field of regard. TESS in particular is predicted to discover many targets around bright stars which will be suitable for follow-up observations. We include these anticipated planets and find that the number of planets Twinkle could observe in the near infrared in a single transit or eclipse increases to 558 for R > 20 and 41 at lower resolutions. By stacking 10 transits or eclipses, there are 1185 potential targets for study at R < 20 as well as 388 planets at higher resolutions.

Stellar variability from pulsations and granulation presents a source of correlated noise that can impact the accuracy and precision of multi-band photometric transit observations of exoplanets. This can potentially cause biased measurements in the transmission or emission spectrum or underestimation of the final error bars on the spectrum. ARIEL is a future space telescope and instrument designed to perform a transit spectroscopic survey of a large sample of exoplanets. In this paper we perform simulations to assess the impact of stellar variability arising from pulsations and granulation on ARIEL observations of GJ 1214b and HD 209458b. We take into account the correlated nature of stellar noise, quantify it, and compare it to photon noise. In the range 1.95-7.8 \textmu m, stellar pulsation and granulation noise has insignificant impact compared to photon noise for both targets. In the visual range the contribution increases significantly but remains small in absolute terms and will have minimal impact on the transmission spectra of the targets studied. The impact of pulsation and granulation will be greatest for planets with low scale height atmospheres and long transit times around bright stars.

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

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

The Mexico-UK Sub-millimetre Camera for AsTronomy (MUSCAT) is a large-format, millimetre-wave camera consisting of 1,500 background-limited lumped-element kinetic inductance detectors (LEKIDs) scheduled for deployment on the Large Millimeter Telescope (Volc\'an Sierra Negra, Mexico) in 2018. MUSCAT is designed for observing at 1.1 mm and will utilise the full 40' field of view of the LMTs upgraded 50-m primary mirror. In its primary role, MUSCAT is designed for high-resolution follow-up surveys of both galactic and extra-galactic sub-mm sources identified by Herschel. MUSCAT is also designed to be a technology demonstrator that will provide the first on-sky demonstrations of novel design concepts such as horn-coupled LEKID arrays and closed continuous cycle miniature dilution refrigeration. Here we describe some of the key design elements of the MUSCAT instrument such as the novel use of continuous sorption refrigerators and a miniature dilutor for continuous 100-mK cooling of the focal plane, broadband optical coupling to Aluminium LEKID arrays using waveguide chokes and anti-reflection coating materials as well as with the general mechanical and optical design of MUSCAT. We explain how MUSCAT is designed to be simple to upgrade and the possibilities for changing the focal plane unit that allows MUSCAT to act as a demonstrator for other novel technologies such as multi-chroic polarisation sensitive pixels and on-chip spectrometry in the future. Finally, we will report on the current status of MUSCAT's commissioning.

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

MUSCAT is a large format mm-wave camera scheduled for installation on the Large Millimeter Telescope Alfonso Serrano (LMT) in 2018. The MUSCAT focal plane is based on an array of horn coupled lumped-element kinetic inductance detectors optimised for coupling to the 1.1mm atmospheric window. The detectors are fed with fully baffled reflective optics to minimize stray-light contamination. This combination will enable background-limited performance at 1.1 mm across the full 4 arcminute field-of-view of the LMT. The easily accessible focal plane will be cooled to 100 mK with a new closed cycle miniature dilution refrigerator that permits fully continuous operation. The MUSCAT instrument will demonstrate the science capabilities of the LMT through two relatively short science programmes to provide high resolution follow-up surveys of Galactic and extra-galactic sources previously observed with the Herschel space observatory, after the initial observing campaigns. In this paper, we will provide an overview of the overall instrument design as well as an update on progress and scheduled installation on the LMT.

The Crab nebula is a supernova remnant exhibiting a highly polarized synchrotron radiation at radio and millimeter wavelengths. It is the brightest source in the microwave sky with an extension of 7 by 5 arcminutes and commonly used as a standard candle for any experiment which aims at measuring the polarization of the sky. Though its spectral energy distribution has been well characterized in total intensity, polarization data are still lacking at millimetre wavelengths. We report in this paper high resolution (18 arcsec FWHM) observations of the Crab nebula in total intensity and linear polarization at 150 GHz with the NIKA camera. NIKA, operated at the IRAM 30 m telescope from 2012 to 2015, is a camera made of Lumped Element Kinetic Inductance Detectors (LEKIDs) observing the sky at 150 and 260 GHz. From these observations we are able to reconstruct the spatial distribution of the polarization degree and angle of the Crab nebula, which is found to be compatible with previous observations at lower and higher frequencies. Averaging across the source and using other existing data sets we find that the Crab nebula polarization angle is consistent with being constant over a wide range of frequencies with a value of -87.7$^\circ$ +- 0.3 in Galactic coordinates. We also present the first estimation of the Crab nebula spectral energy distribution polarized flux in a wide frequency range: 30-353 GHz. Assuming a single power law emission model we find that the polarization spectral index $\beta_{pol}$ = - 0.347 +- 0.026 is compatible with the intensity spectral index $\beta$ = - 0.323 +- 0.001.

Since the discovery of the first extrasolar planet more than twenty years ago, we have discovered nearly four thousand planets orbiting stars other than the Sun. Current observational instruments (on board the Hubble Space Telescope, Spitzer, and on ground-based facilities) have allowed the scientific community to obtain important information on the physical and chemical properties of these planets. However, for a more in-depth characterisation of these worlds, more powerful telescopes are needed. Thanks to the high sensitivity of their instruments, the next generation of space observatories (e.g. JWST, ARIEL) will provide observations of unprecedented quality, allowing us to extract far more information than what was previously possible. Such high quality observations will provide constraints on theoretical models of exoplanet atmospheres and lead to a greater understanding of their physics and chemistry. Important modelling efforts have been carried out during the past few years, showing that numerous parameters and processes (such as the elemental abundances, temperature, mixing, etc.) are likely to affect the atmospheric composition of exoplanets and subsequently the observable spectra. In this manuscript, we review the different parameters that can influence the molecular composition of exoplanet atmospheres. We show that the high-precision of ARIEL observations will improve our view and characterisation of exoplanet atmospheres. We also consider future developments that are necessary to improve atmospheric models, driven by the need to interpret the available observations.

We discuss a systematic effect associated with measuring polarization with a continuously rotating half-wave plate. The effect was identified with the data from the E and B Experiment (EBEX), which was a balloon-borne instrument designed to measure the polarization of the CMB as well as that from Galactic dust. The data show polarization fraction larger than 10\% while less than 3\% were expected from instrumental polarization. We give evidence that the excess polarization is due to detector non-linearity in the presence of a continuously rotating HWP. The non-linearity couples intensity signals into polarization. We develop a map-based method to remove the excess polarization. Applying this method for the 150 (250) GHz bands data we find that 81\% (92\%) of the excess polarization was removed. Characterization and mitigation of this effect is important for future experiments aiming to measure the CMB B-modes with a continuously rotating HWP.

The thermal emission of dust grains is a powerful tool for probing cold, dense regions of molecular gas in the ISM, and so constraining dust properties is key to obtaining accurate measurements of dust mass and temperature. By placing constraints on the dust emissivity spectral index, beta, towards two star-forming infrared dark clouds, SDC18.888 and SDC24.489, we evaluate the role of mass concentration in the associated star-formation activity. We exploit the simultaneous 1.2mm and 2.0mm imaging capability of NIKA on the IRAM 30m telescope to construct maps of beta for both clouds, and by incorporating Herschel observations, we create H2 column density maps with 13" resolution. While we find no significant systematic radial variations around the most massive clumps in either cloud on >0.1 pc scales, their mean beta values are significantly different, with beta = 2.07 +/- 0.09 (rand) +/- 0.25 (syst) for SDC18.888 and beta = 1.71 +/- 0.09 (rand) +/- 0.25 (syst) for SDC24.489. These differences could be a consequence of the very different environments in which both clouds lie, and we suggest that the proximity of SDC18.888 to the W39 HII region may raise beta on scales of 1 pc. We also find that the mass in SDC24.489 is more centrally concentrated and circularly symmetric than in SDC18.888, and is consistent with a scenario in which spherical globally-collapsing clouds concentrate a higher fraction of their mass into a single core than elongated clouds that will more easily fragment, distributing their mass into many cores. We demonstrate that beta variations towards interstellar clouds can be robustly constrained with high-SNR NIKA observations, providing more accurate estimates of their masses. The methods presented here will be applied to the Galactic Star Formation with NIKA2 (GASTON) large programme, extending our analysis to a statistically significant sample of star-forming clouds.

The New IRAM KID Arrays 2 (NIKA2) consortium has just finished installing and commissioning a millimetre camera on the IRAM 30 m telescope. It is a dual-band camera operating with three frequency multiplexed kilo-pixels arrays of Lumped Element Kinetic Inductance Detectors (LEKID) cooled at 150 mK, designed to observe the intensity and polarisation of the sky at 260 and 150 GHz (1.15 and 2 mm). NIKA2 is today an IRAM resident instrument for millimetre astronomy, such as Intra Cluster Medium from intermediate to distant clusters and so for the follow-up of Planck satellite detected clusters, high redshift sources and quasars, early stages of star formation and nearby galaxies emission. We present an overview of the instrument performance as it has been evaluated at the end of the commissioning phase.

Using conventional superconductor theory we discuss and validate a model that describes the energy-resolving performance of an aluminium LEKID to single-photon absorption events. While aluminium is not the optimum material for single-photon counting applications, this material is well understood and is used to understand the underlying device physics of these detectors. We also discuss data analysis techniques used to extract single-photon detections from noisy data.

Sub-structures in the hot gas of galaxy clusters are related to their formation history and to the astrophysical processes at play in the intracluster medium (ICM). The thermal Sunyaev-Zel'dovich (tSZ) effect is directly sensitive to the line-of-sight integrated ICM pressure, and is thus particularly adapted to study ICM sub-structures. We apply structure-enhancement filtering algorithms to high resolution tSZ observations of distant clusters, in order to search for pressure discontinuities, compressions, as well as secondary peaks in the ICM. The same filters are applied to synthetic tSZ images extracted from RHAPSODY-G hydrodynamic simulations, in order to better interpret the extracted features. We also study the noise propagation trough the filters and quantify the impact of systematic effects, point source residuals being identified as the dominant potential contaminant. In 3 of our 6 NIKA clusters we identify features at high S/N that show clear evidence for merger events. In MACSJ0717 (z=0.55), three strong pressure gradients are observed on the E, SE and W sectors, and two main peaks in the pressure distribution are identified. We observe a lack of tSZ compact structure in the cool-core cluster PSZ1G045.85 (z=0.61), and a tSZ gradient ridge dominates in the SE. In the highest z cluster, CLJ1227 (z=0.89), we detect a ~45" (360 kpc) long ridge pressure gradient associated with a secondary pressure peak in the W region. Our results show that current tSZ facilities have now reached the angular resolution and sensitivity to allow an exploration of the details of pressure sub-structures in clusters, even at high z. This opens the possibility to quantify the impact of the dynamical state on the relation between the tSZ signal and the mass of clusters, which is important when using tSZ clusters to test cosmological models. This work also marks the first NIKA cluster sample data release.

Millimeter-wave continuum astronomy is today an indispensable tool for both general Astrophysics studies and Cosmology. General purpose, large field-of-view instruments are needed to map the sky at intermediate angular scales not accessible by the high-resolution interferometers and by the coarse angular resolution space-borne or ground-based surveys. These instruments have to be installed at the focal plane of the largest single-dish telescopes. In this context, we have constructed and deployed a multi-thousands pixels dual-band (150 and 260 GHz, respectively 2mm and 1.15mm wavelengths) camera to image an instantaneous field-of-view of 6.5arc-min and configurable to map the linear polarization at 260GHz. We are providing a detailed description of this instrument, named NIKA2 (New IRAM KID Arrays 2), in particular focusing on the cryogenics, the optics, the focal plane arrays based on Kinetic Inductance Detectors (KID) and the readout electronics. We are presenting the performance measured on the sky during the commissioning runs that took place between October 2015 and April 2017 at the 30-meter IRAM (Institut of Millimetric Radio Astronomy) telescope at Pico Veleta. NIKA2 has been successfully deployed and commissioned, performing in-line with the ambitious expectations. In particular, NIKA2 exhibits FWHM angular resolutions of around 11 and 17.5 arc-seconds at respectively 260 and 150GHz. The NEFD (Noise Equivalent Flux Densities) demonstrated on the maps are, at these two respective frequencies, 33 and 8 mJy*sqrt(s). A first successful science verification run has been achieved in April 2017. The instrument is currently offered to the astronomical community during the coming winter and will remain available for at least the next ten years.

CASTAway is a mission concept to explore our Solar System's main asteroid belt. Asteroids and comets provide a window into the formation and evolution of our Solar System and the composition of these objects can be inferred from space-based remote sensing using spectroscopic techniques. Variations in composition across the asteroid populations provide a tracer for the dynamical evolution of the Solar System. The mission combines a long-range (point source) telescopic survey of over 10,000 objects, targeted close encounters with 10 to 20 asteroids and serendipitous searches to constrain the distribution of smaller (e.g. 10 m) size objects into a single concept. With a carefully targeted trajectory that loops through the asteroid belt, CASTAway would provide a comprehensive survey of the main belt at multiple scales. The scientific payload comprises a 50 cm diameter telescope that includes an integrated low-resolution (R = 30 to 100) spectrometer and visible context imager, a thermal (e.g. 6 to 16 microns) imager for use during the flybys, and modified star tracker cameras to detect small (approx. 10 m) asteroids. The CASTAway spacecraft and payload have high levels of technology readiness and are designed to fit within the programmatic and cost caps for a European Space Agency medium class mission, whilst delivering a significant increase in knowledge of our Solar System.

The development of precision cosmology with clusters of galaxies requires high-angular resolution Sunyaev-Zel'dovich (SZ) observations. As for now, arcmin resolution SZ observations (e.g. SPT, ACT and Planck) only allowed detailed studies of the intra cluster medium for low redshift clusters (z<0.2). With both a wide field of view (6.5 arcmin) and a high angular resolution (17.7 and 11.2 arcsec at 150 and 260 GHz), the NIKA2 camera installed at the IRAM 30-m telescope (Pico Veleta, Spain), will bring valuable information in the field of SZ imaging of clusters of galaxies. The NIKA2 SZ observation program will allow us to observe a large sample of clusters (50) at redshifts between 0.4 and 0.9. As a pilot study for NIKA2, several clusters of galaxies have been observed with the pathfinder, NIKA, at the IRAM 30-m telescope to cover the various configurations and observation conditions expected for NIKA2.

We propose a method to map the temperature distribution of the hot gas in galaxy clusters that uses resolved images of the thermal Sunyaev-Zel'dovich (tSZ) effect in combination with X-ray data. Application to images from the New IRAM KIDs Array (NIKA) and XMM-Newton allows us to measure and determine the spatial distribution of the gas temperature in the merging cluster MACS J0717.5+3745, at $z=0.55$. Despite the complexity of the target object, we find a good morphological agreement between the temperature maps derived from X-ray spectroscopy only -- using XMM-Newton ($T_{\rm XMM}$) and Chandra ($T_{\rm CXO}$) -- and the new gas-mass-weighted tSZ+X-ray imaging method ($T_{\rm SZ}$). We correlate the temperatures from tSZ+X-ray imaging and those from X-ray spectroscopy alone and find that $T_{\rm SZ}$ is higher than $T_{\rm XMM}$ and lower than $T_{\rm CXO}$ by $\sim 10\%$ in both cases. Our results are limited by uncertainties in the geometry of the cluster gas, contamination from kinetic SZ ($\sim 10\%$), and the absolute calibration of the tSZ map ($7\%$). Investigation using a larger sample of clusters would help minimise these effects.

Context: In the past decade, sensitive, resolved Sunyaev-Zel'dovich (SZ) studies of galaxy clusters have become common. Whereas many previous SZ studies have parameterized the pressure profiles of galaxy clusters, non-parametric reconstructions will provide insights into the thermodynamic state of the intracluster medium (ICM). Aims: We seek to recover the non-parametric pressure profiles of the high redshift ($z=0.89$) galaxy cluster CLJ 1226.9+3332 as inferred from SZ data from the MUSTANG, NIKA, Bolocam, and Planck instruments, which all probe different angular scales. Methods: Our non-parametric algorithm makes use of logarithmic interpolation, which under the assumption of ellipsoidal symmetry is analytically integrable. For MUSTANG, NIKA, and Bolocam we derive a non-parametric pressure profile independently and find good agreement among the instruments. In particular, we find that the non-parametric profiles are consistent with a fitted gNFW profile. Given the ability of Planck to constrain the total signal, we include a prior on the integrated Compton Y parameter as determined by Planck. Results: For a given instrument, constraints on the pressure profile diminish rapidly beyond the field of view. The overlap in spatial scales probed by these four datasets is therefore critical in checking for consistency between instruments. By using multiple instruments, our analysis of CLJ 1226.9+3332 covers a large radial range, from the central regions to the cluster outskirts: $0.05 R_{500} < r < 1.1 R_{500}$. This is a wider range of spatial scales than is typical recovered by SZ instruments. Similar analyses will be possible with the new generation of SZ instruments such as NIKA2 and MUSTANG2.

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

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

The characterization of dust properties in the interstellar medium (ISM) is key for star formation. Mass estimates are crucial to determine gravitational collapse conditions for the birth of new stellar objects in molecular clouds. However, most of these estimates rely on dust models that need further observational constraints from clouds to prestellar and protostellar cores. We present results of a study of dust emissivity changes based on mm-continuum data obtained with the NIKA camera at the IRAM-30m telescope. Observing dust emission at 1.15 mm and 2 mm allows us to constrain the dust emissivity index ($\beta$) in the Rayleigh-Jeans tail of the dust spectral energy distribution (SED) far from its peak emission, where the contribution of other parameters (i.e. dust temperature) is important. Focusing on the Taurus molecular cloud, a low-mass star-forming regions in the Gould Belt, we analyze the emission properties of several distinct objects in the B213 filament: three prestellar cores, two Class-0/I protostellar cores and one Class-II object. By means of the ratio of the two NIKA channel-maps, we show that in the Rayleigh-Jeans approximation the dust emissivity index varies among the objects. For one prestellar and two protostellar cores, we produce a robust study using Herschel data to constrain the dust temperature of the sources. By using the Abel transform inversion technique we get accurate radial $\beta$ profiles. We find systematic spatial variations of $\beta$ in the protostellar cores that is not observed in the prestellar core. While in the former case $\beta$ decreases toward the center, in the latter it remains constant. Moreover, $\beta$ appears anticorrelated with the dust temperature. We discuss the implication of these results in terms of dust grain evolution between pre- and protostellar cores.

The Atmospheric Remote-sensing Infrared Exoplanet Large-survey mission (ARIEL) is one of the three present candidates for the ESA M4 (the fourth medium mission) launch opportunity. The proposed Payload will perform a large unbiased spectroscopic survey from space concerning the nature of exoplanets atmospheres and their interiors to determine the key factors affecting the formation and evolution of planetary systems. ARIEL will observe a large number (>500) of warm and hot transiting gas giants, Neptunes and super-Earths around a wide range of host star types, targeting planets hotter than 600 K to take advantage of their well-mixed atmospheres. It will exploit primary and secondary transits spectroscopy in the 1.2-8 um spectral range and broad-band photometry in the optical and Near IR (NIR). The main instrument of the ARIEL Payload is the IR Spectrometer (AIRS) providing low-resolution spectroscopy in two IR channels: Channel 0 (CH0) for the 1.95-3.90 um band and Channel 1 (CH1) for the 3.90-7.80 um range. It is located at the intermediate focal plane of the telescope and common optical system and it hosts two IR sensors and two cold front-end electronics (CFEE) for detectors readout, a well defined process calibrated for the selected target brightness and driven by the Payload's Instrument Control Unit (ICU).

The E and B Experiment (EBEX) was a long-duration balloon-borne instrument designed to measure the polarization of the cosmic microwave background (CMB) radiation. EBEX was the first balloon-borne instrument to implement a kilo-pixel array of transition edge sensor (TES) bolometric detectors and the first CMB experiment to use the digital version of the frequency domain multiplexing system for readout of the TES array. The scan strategy relied on 40 s peak-to-peak constant velocity azimuthal scans. We discuss the unique demands on the design and operation of the payload that resulted from these new technologies and the scan strategy. We describe the solutions implemented including the development of a power system designed to provide a total of at least 2.3 kW, a cooling system to dissipate 590 W consumed by the detectors' readout system, software to manage and handle the data of the kilo-pixel array, and specialized attitude reconstruction software. We present flight performance data showing faultless management of the TES array, adequate powering and cooling of the readout electronics, and constraint of attitude reconstruction errors such that the spurious B-modes they induced were less than 10% of CMB B-mode power spectrum with $r=0.05$.

Magnetic fields, which play a major role in a large number of astrophysical processes from galactic to cosmological scales, can be traced via observations of dust polarization as demonstrated by the Planck satellite results. In particular, low-resolution observations of dust polarization have demonstrated that Galactic filamentary structures, where star formation takes place, are associated to well organized magnetic fields. A better understanding of this process requires detailed observations of galactic dust polarization on scales of 0.01 to 0.1 pc. Such high-resolution polarization observations can be carried out at the IRAM 30 m telescope using the recently installed NIKA2 camera, which features two frequency bands at 260 and 150 GHz (respectively 1.15 and 2.05 mm), the 260 GHz band being polarization sensitive. NIKA2 so far in commissioning phase, has its focal plane filled with ~3300 detectors to cover a Field of View (FoV) of 6.5 arcminutes diameter. The NIKA camera, which consisted of two arrays of 132 and 224 Lumped Element Kinetic Inductance Detectors (LEKIDs) and a FWHM (Full-Width-Half-Maximum) of 12 and 18.2 arcsecond at 1.15 and 2.05 mm respectively, has been operated at the IRAM 30 m telescope from 2012 to 2015 as a test-bench for NIKA2. NIKA was equipped of a room temperature polarization system (a half wave plate (HWP) and a grid polarizer facing the NIKA cryostat window). The fast and continuous rotation of the HWP permits the quasi simultaneous reconstruction of the three Stokes parameters, I, Q and U at 150 and 260 GHz. This paper presents the first polarization measurements with KIDs and reports the polarization performance of the NIKA camera and the pertinence of the choice of the polarization setup in the perspective of NIKA2. (abridged)

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

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

Measurement of the gas velocity distribution in galaxy clusters provides insight into the physics of mergers, through which large scale structures form in the Universe. Velocity estimates within the intracluster medium (ICM) can be obtained via the Sunyaev-Zel'dovich (SZ) effect, but its observation is challenging both in term of sensitivity requirement and control of systematic effects, including the removal of contaminants. In this paper we report resolved observations, at 150 and 260 GHz, of the SZ effect toward the triple merger MACS J0717.5+3745 (z=0.55), using data obtained with the NIKA camera at the IRAM 30m telescope. Assuming that the SZ signal is the sum of a thermal (tSZ) and a kinetic (kSZ) component and by combining the two NIKA bands, we extract for the first time a resolved map of the kSZ signal in a cluster. The kSZ signal is dominated by a dipolar structure that peaks at -5.1 and +3.4 sigma, corresponding to two subclusters moving respectively away and toward us and coincident with the cold dense X-ray core and a hot region undergoing a major merging event. We model the gas electron density and line-of-sight velocity of MACS J0717.5+3745 as four subclusters. Combining NIKA data with X-ray observations from XMM-Newton and Chandra, we fit this model to constrain the gas line-of-sight velocity of each component, and we also derive, for the first time, a velocity map from kSZ data (i.e. that is model-dependent). Our results are consistent with previous constraints on the merger velocities, and thanks to the high angular resolution of our data, we are able to resolve the structure of the gas velocity. Finally, we investigate possible contamination and systematic effects with a special care given to radio and submillimeter galaxies. Among the sources that we detect with NIKA, we find one which is likely to be a high redshift lensed submillimeter galaxy.

The determination of the thermodynamic properties of clusters of galaxies at intermediate and high redshift can bring new insights into the formation of large-scale structures. It is essential for a robust calibration of the mass-observable scaling relations and their scatter, which are key ingredients for precise cosmology using cluster statistics. Here we illustrate an application of high resolution $(< 20$ arcsec) thermal Sunyaev-Zel'dovich (tSZ) observations by probing the intracluster medium (ICM) of the \planck-discovered galaxy cluster \psz\ at redshift $z = 0.61$, using tSZ data obtained with the NIKA camera, which is a dual-band (150 and 260~GHz) instrument operated at the IRAM 30-meter telescope. We deproject jointly NIKA and \planck\ data to extract the electronic pressure distribution from the cluster core ($R \sim 0.02\, R_{500}$) to its outskirts ($R \sim 3\, R_{500}$) non-parametrically for the first time at intermediate redshift. The constraints on the resulting pressure profile allow us to reduce the relative uncertainty on the integrated Compton parameter by a factor of two compared to the \planck\ value. Combining the tSZ data and the deprojected electronic density profile from \xmm\ allows us to undertake a hydrostatic mass analysis, for which we study the impact of a spherical model assumption on the total mass estimate. We also investigate the radial temperature and entropy distributions. These data indicate that \psz\ is a massive ($M_{500} \sim 5.5 \times 10^{14}$ M$_{\odot}$) cool-core cluster. This work is part of a pilot study aiming at optimizing the treatment of the NIKA2 tSZ large program dedicated to the follow-up of SZ-discovered clusters at intermediate and high redshifts. (abridged)

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

We present the first major data release of the largest single key-project in area carried out in open time with the Herschel Space Observatory. The Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) is a survey of 600 deg^2 in five photometric bands - 100, 160, 250, 350 and 500 um - with the PACS and SPIRE cameras. In this paper and a companion paper (Bourne et al. 2016) we present the survey of three fields on the celestial equator, covering a total area of 161.6 deg^2 and previously observed in the Galaxy and Mass Assembly (GAMA) spectroscopic survey. This paper describes the Herschel images and catalogues of the sources detected on the SPIRE 250 um images. The 1-sigma noise for source detection, including both confusion and instrumental noise, is 7.4, 9.4 and 10.2 mJy at 250, 350 and 500 um. Our catalogue includes 120230 sources in total, with 113995, 46209 and 11011 sources detected at >4-sigma at 250, 350 and 500 um. The catalogue contains detections at >3-sigma at 100 and 160 um for 4650 and 5685 sources, and the typical noise at these wavelengths is 44 and 49 mJy. We include estimates of the completeness of the survey and of the effects of flux bias and also describe a novel method for determining the true source counts. The H-ATLAS source counts are very similar to the source counts from the deeper HerMES survey at 250 and 350 um, with a small difference at 500 um. Appendix A provides a quick start in using the released datasets, including instructions and cautions on how to use them.

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

The New IRAM KID Array 2 (NIKA 2) is a dual-band camera operating with three frequency-multiplexed kilopixels arrays of Lumped Element Kinetic Inductance Detectors (LEKID) cooled at 150 mK. NIKA 2 is designed to observe the intensity and polarisation of the sky at 1.15 and 2.0 mm wavelength from the IRAM 30 m telescope. The NIKA 2 instrument represents a huge step in performance as compared to the NIKA pathfinder instrument, which has already shown state-of-the-art detector and photometric performance. After the commissioning planned to be accomplished at the end of 2016, NIKA 2 will be an IRAM resident instrument for the next ten years or more. NIKA 2 should allow the astrophysical community to tackle a large number of open questions reaching from the role of the Galactic magnetic field in star formation to the discrepancy between cluster-based and CMB-based cosmology possibly induced by the unknown cluster physics. We present an overview of the commissioning phase together with some first results.

NIKA2 (New IRAM KID Arrays) is a dual band (150 and 260 GHz) imaging camera based on Kinetic Inductance Detectors (KIDs) and designed to work at the IRAM 30 m telescope (Pico Veleta, Spain). Built on the experience of the NIKA prototype, NIKA2 has been installed at the 30 m focal plane in October 2015 and the commissioning phase is now ongoing. Through the thermal Sunyaev-Zeldovich (tSZ) effect, NIKA2 will image the ionized gas residing in clusters of galaxies with a resolution of 12 and 18 arcsec FWHM (at 150 and 260 GHz, respectively). We report on the recent tSZ measurements with the NIKA camera and discuss the future objectives for the NIKA2 SZ large Program, 300h of observation dedicated to SZ science. With this program we intend to perform a high angular resolution follow-up of a cosmologically-representative sample of clusters belonging to SZ catalogues, with redshift greater than 0.5. The main output of the program will be the study of the redshift evolution of the cluster pressure profile as well as that of the scaling laws relating the cluster global properties.

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

The prototype of the NIKA2 camera, NIKA, is an instrument operating at the IRAM 30-m telescope, which can observe simultaneously at 150 and 260GHz. One of the main goals of NIKA2 is to measure the pressure distribution in galaxy clusters at high resolution using the thermal SZ (tSZ) effect. Such observations have already proved to be an excellent probe of cluster pressure distributions even at high redshifts. However, an important fraction of clusters host submm and/or radio point sources, which can significantly affect the reconstructed signal. Here we report on <20" resolution observations at 150 and 260GHz of the cluster MACSJ1424, which hosts both radio and submm point sources. We examine the morphology of the tSZ signal and compare it to other datasets. The NIKA data are combined with Herschel satellite data to study the SED of the submm point source contaminants. We then perform a joint reconstruction of the intracluster medium (ICM) electronic pressure and density by combining NIKA, Planck, XMM-Newton, and Chandra data, focusing on the impact of the radio and submm sources on the reconstructed pressure profile. We find that large-scale pressure distribution is unaffected by the point sources because of the resolved nature of the NIKA observations. The reconstructed pressure in the inner region is slightly higher when the contribution of point sources are removed. We show that it is not possible to set strong constraints on the central pressure distribution without accurately removing these contaminants. The comparison with X-ray only data shows good agreement for the pressure, temperature, and entropy profiles, which all indicate that MACSJ1424 is a dynamically relaxed cool core system. The present observations illustrate the possibility of measuring these quantities with a relatively small integration time, even at high redshift and without X-ray spectroscopy.

NIKA 2 (New Instrument of Kids Array) is a next generation continuum and polarized instrument successfully installed in October 2015 at the IRAM 30 m telescope on Pico-Veleta (Granada, Spain). NIKA 2 is a high resolution dual-band camera, operating with frequency multiplexed LEKIDs (Lumped Element Kinetic Inductance Detectors) cooled at 100 mK. Dual color images are obtained thanks to the simultaneous readout of a 1020 pixels array at 2 mm and 1140 x 2 pixels arrays at 1.15 mm with a final resolution of 18 and 12 arcsec respectively, and 6.5 arcmin of Field of View (FoV). The two arrays at 1.15 mm allow us to measure the linear polarization of the incoming light. This will place NIKA 2 as an instrument of choice to study the role of magnetic fields in the star formation process. The NIKA experiment, a prototype for NIKA 2 with a reduced number of detectors (about 400 LEKIDs) and FoV (1.8 arcmin), has been successfully operated at the IRAM 30 telescope in several open observational campaigns. The performance of the NIKA 2 polarization setup has been successfully validated with the NIKA prototype.

The E and B Experiment (EBEX) is a balloon-borne polarimeter designed to measure the polarization of the cosmic microwave background radiation and to characterize the polarization of galactic dust. EBEX was launched December 29, 2012 and circumnavigated Antarctica observing $\sim$6,000 square degrees of sky during 11 days at three frequency bands centered around 150, 250 and 410 GHz. EBEX was the first experiment to operate a kilo-pixel array of transition-edge sensor bolometers and a continuously rotating achromatic half-wave plate aboard a balloon platform. It also pioneered the use of detector readout based on digital frequency domain multiplexing. We describe the temperature calibration of the experiment. The gain response of the experiment is calibrated using a two-step iterative process. We use signals measured on passes across the Galactic plane to convert from readout-system counts to power. The effective smoothing scale of the EBEX optics and the star camera-to-detector offset angles are determined through \c{hi}2 minimization using the compact HII region RCW 38. This two-step process is initially performed with parameters measured before the EBEX 2013 flight and then repeated until the calibration factor and parameters converge.

NIKA2 (New IRAM KID Array 2) is a camera dedicated to millimeter wave astronomy based upon kilopixel arrays of Kinetic Inductance Detectors (KID). The pathfinder instrument, NIKA, has already shown state-of-the-art detector performance. NIKA2 builds upon this experience but goes one step further, increasing the total pixel count by a factor $\sim$10 while maintaining the same per pixel performance. For the next decade, this camera will be the resident photometric instrument of the Institut de Radio Astronomie Millimetrique (IRAM) 30m telescope in Sierra Nevada (Spain). In this paper we give an overview of the main components of NIKA2, and describe the achieved detector performance. The camera has been permanently installed at the IRAM 30m telescope in October 2015. It will be made accessible to the scientific community at the end of 2016, after a one-year commissioning period. When this happens, NIKA2 will become a fundamental tool for astronomers worldwide.

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

We have developed a passive 350 GHz (850 {\mu}m) video-camera to demonstrate lumped element kinetic inductance detectors (LEKIDs) -- designed originally for far-infrared astronomy -- as an option for general purpose terrestrial terahertz imaging applications. The camera currently operates at a quasi-video frame rate of 2 Hz with a noise equivalent temperature difference per frame of $\sim$0.1 K, which is close to the background limit. The 152 element superconducting LEKID array is fabricated from a simple 40 nm aluminum film on a silicon dielectric substrate and is read out through a single microwave feedline with a cryogenic low noise amplifier and room temperature frequency domain multiplexing electronics.

NIKA is a dual-band camera operating with 315 frequency multiplexed LEKIDs cooled at 100 mK. NIKA is designed to observe the sky in intensity and polarisation at 150 and 260 GHz from the IRAM 30-m telescope. It is a test-bench for the final NIKA2 camera. The incoming linear polarisation is modulated at four times the mechanical rotation frequency by a warm rotating multi-layer Half Wave Plate. Then, the signal is analysed by a wire grid and finally absorbed by the LEKIDs. The small time constant (< 1ms ) of the LEKID detectors combined with the modulation of the HWP enables the quasi-simultaneous measurement of the three Stokes parameters I, Q, U, representing linear polarisation. In this paper we present results of recent observational campaigns demonstrating the good performance of NIKA in detecting polarisation at mm wavelength.

The thermal Sunyaev-Zel'dovich (tSZ) effect is expected to provide a low scatter mass proxy for galaxy clusters since it is directly proportional to the cluster thermal energy. The tSZ observations have proven to be a powerful tool for detecting and studying them, but high angular resolution observations are now needed to push their investigation to a higher redshift. In this paper, we report high angular (< 20 arcsec) resolution tSZ observations of the high-redshift cluster CL J1226.9+3332 (z=0.89). It was imaged at 150 and 260 GHz using the NIKA camera at the IRAM 30-meter telescope. The 150 GHz map shows that CL J1226.9+3332 is morphologically relaxed on large scales with evidence of a disturbed core, while the 260 GHz channel is used mostly to identify point source contamination. NIKA data are combined with those of Planck and X-ray from Chandra to infer the cluster's radial pressure, density, temperature, and entropy distributions. The total mass profile of the cluster is derived, and we find $M_{500} = 5.96^{+1.02}_{-0.79} $ x $10^{14} M_{\odot}$ within the radius $R_{500} = 930^{+50}_{-43}$ kpc, at a 68% confidence level. ($R_{500}$ is the radius within which the average density is 500 times the critical density at the cluster's redshift.) NIKA is the prototype camera of NIKA2, a KIDs (kinetic inductance detectors) based instrument to be installed at the end of 2015. This work is, therefore, part of a pilot study aiming at optimizing tSZ NIKA2 large programs.

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

The Exoplanet Characterisation Observatory (EChO) has been studied as a space mission concept by the European Space Agency in the context of the M3 selection process. Through direct measurement of the atmospheric chemical composition of hundreds of exoplanets, EChO would address fundamental questions such as: What are exoplanets made of? How do planets form and evolve? What is the origin of exoplanet diversity? More specifically, EChO is a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large, diverse and well-defined planetary sample within its four to six year mission lifetime. In this paper we use the end-to-end instrument simulator EChOSim to model the currently discovered targets, to gauge which targets are observable and assess the EChO performances obtainable for each observing tier and time. We show that EChO would be capable of observing over 170 relativity diverse planets if it were launched today, and the wealth of optimal targets for EChO expected to be discovered in the next 10 years by space and ground-based facilities is simply overwhelming. In addition, we build on previous molecular detectability studies to show what molecules and abundances will be detectable by EChO for a selection of real targets with various molecular compositions and abundances. EChO's unique contribution to exoplanetary science will be in identifying the main constituents of hundreds of exoplanets in various mass/temperature regimes, meaning that we will be looking no longer at individual cases but at populations. Such a universal view is critical if we truly want to understand the processes of planet formation and evolution in various environments. In this paper we present a selection of key results. The full results are available online (http://www.ucl.ac.uk/exoplanets/echotargetlist/).

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

NIKA (New IRAM KID Arrays) is a dual-band imaging instrument installed at the IRAM (Institut de RadioAstronomie Millimetrique) 30-meter telescope at Pico Veleta (Spain). Two distinct Kinetic Inductance Detectors (KID) focal planes allow the camera to simultaneously image a field-of-view of about 2 arc-min in the bands 125 to 175 GHz (150 GHz) and 200 to 280 GHz (240 GHz). The sensitivity and stability achieved during the last commissioning Run in June 2013 allows opening the instrument to general observers. We report here the latest results, in particular in terms of sensitivity, now comparable to the state-of-the-art Transition Edge Sensors (TES) bolometers, relative and absolute photometry. We describe briefly the next generation NIKA-2 instrument, selected by IRAM to occupy, from 2015, the continuum imager/polarimeter slot at the 30-m telescope.

Clusters of galaxies provide valuable information on the evolution of the Universe and large scale structures. Recent cluster observations via the thermal Sunyaev-Zel'dovich (tSZ) effect have proven to be a powerful tool to detect and study them. In this context, high resolution tSZ observations (~ tens of arcsec) are of particular interest to probe intermediate and high redshift clusters. Observations of the tSZ effect will be carried out with the millimeter dual-band NIKA2 camera, based on Kinetic Inductance Detectors (KIDs) to be installed at the IRAM 30-meter telescope in 2015. To demonstrate the potential of such an instrument, we present tSZ observations with the NIKA camera prototype, consisting of two arrays of 132 and 224 detectors that observe at 140 and 240 GHz with a 18.5 and 12.5 arcsec angular resolution, respectively. The cluster RX J1347.5-1145 was observed simultaneously at 140 and 240 GHz. We used a spectral decorrelation technique to remove the atmospheric noise and obtain a map of the cluster at 140 GHz. The efficiency of this procedure has been characterized through realistic simulations of the observations. The observed 140 GHz map presents a decrement at the cluster position consistent with the tSZ nature of the signal. We used this map to study the pressure distribution of the cluster by fitting a gNFW model to the data. Subtracting this model from the map, we confirm that RX J1347.5-1145 is an ongoing merger, which confirms and complements previous tSZ and X-ray observations. For the first time, we demonstrate the tSZ capability of KID based instruments. The NIKA2 camera with ~ 5000 detectors and a 6.5 arcmin field of view will be well-suited for in-depth studies of the intra cluster medium in intermediate to high redshifts, which enables the characterization of recently detected clusters by the Planck satellite.

High resolution observations of the thermal Sunyaev-Zel'dovich (tSZ) effect are necessary to allow the use of clusters of galaxies as a probe for large scale structures at high redshifts. With its high resolution and dual-band capability at millimeter wavelengths, the NIKA camera can play a significant role in this context. NIKA is based on newly developed Kinetic Inductance Detectors (KIDs) and operates at the IRAM 30m telescope, Pico Veleta, Spain. In this paper, we give the status of the NIKA camera, focussing on the KID technology. We then present observations of three galaxy clusters: RX J1347.5-1145 as a demonstrator of the NIKA capabilities and the recent observations of CL J1226.9+3332 (z = 0.89) and MACS J0717.5+3745 (z = 0.55). We also discuss prospects for the final NIKA2 camera, which will have a 6.5 arcminute field of view with about 5000 detectors in two bands at 150 and 260 GHz.

The New IRAM KID Array (NIKA) instrument is a dual-band imaging camera operating with Kinetic Inductance Detectors (KID) cooled at 100 mK. NIKA is designed to observe the sky at wavelengths of 1.25 and 2.14 mm from the IRAM 30 m telescope at Pico Veleta with an estimated resolution of 13\,arcsec and 18 arcsec, respectively. This work presents the performance of the NIKA camera prior to its opening to the astrophysical community as an IRAM common-user facility in early 2014. NIKA is a test bench for the final NIKA2 instrument to be installed at the end of 2015. The last NIKA observation campaigns on November 2012 and June 2013 have been used to evaluate this performance and to improve the control of systematic effects. We discuss here the dynamical tuning of the readout electronics to optimize the KID working point with respect to background changes and the new technique of atmospheric absorption correction. These modifications significantly improve the overall linearity, sensitivity, and absolute calibration performance of NIKA. This is proved on observations of point-like sources for which we obtain a best sensitivity (averaged over all valid detectors) of 40 and 14 mJy.s$^{1/2}$ for optimal weather conditions for the 1.25 and 2.14 mm arrays, respectively. NIKA observations of well known extended sources (DR21 complex and the Horsehead nebula) are presented. This performance makes the NIKA camera a competitive astrophysical instrument.

The New IRAM KID Array (NIKA) is a dual-band camera operating with frequency multiplexed arrays of Lumped Element Kinetic Inductance Detectors (LEKIDs) cooled to 100 mK. NIKA is designed to observe the intensity and polarisation of the sky at 1.25 and 2.14 mm from the IRAM 30 m telescope. We present the improvements on the control of systematic effects and astrophysical results made during the last observation campaigns between 2012 and 2014.

EBEX is a balloon-borne telescope designed to measure the polarization of the cosmic microwave background radiation. During its eleven day science flight in the Austral Summer of 2012, it operated 955 spider-web transition edge sensor (TES) bolometers separated into bands at 150, 250 and 410 GHz. This is the first time that an array of TES bolometers has been used on a balloon platform to conduct science observations. Polarization sensitivity was provided by a wire grid and continuously rotating half-wave plate. The balloon implementation of the bolometer array and readout electronics presented unique development requirements. Here we present an outline of the readout system, the remote tuning of the bolometers and Superconducting QUantum Interference Device (SQUID) amplifiers, and preliminary current noise of the bolometer array and readout system.

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

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

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

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

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

EChOSim is the end-to-end time-domain simulator of the Exoplanet Characterisation Observatory (EChO) space mission. EChOSim has been developed to assess the capability EChO has to detect and characterize the atmospheres of transiting exoplanets, and through this revolutionize the knowledge we have of the Milky Way and of our place in the Galaxy. Here we discuss the details of the EChOSim implementation and describe the models used to represent the instrument and to simulate the detection. Software simulators have assumed a central role in the design of new instrumentation and in assessing the level of systematics affecting the measurements of existing experiments. Thanks to its high modularity, EChOSim can simulate basic aspects of several existing and proposed spectrometers for exoplanet transits, including instruments on the Hubble Space Telescope and Spitzer, or ground-based and balloon borne experiments. A discussion of different uses of EChOSim is given, including examples of simulations performed to assess the EChO mission.

We demonstrate the effectiveness of the Exoplanet Characterisation Observatory mission concept for constraining the atmospheric properties of hot and warm gas giants and super Earths. Synthetic primary and secondary transit spectra for a range of planets are passed through EChOSim (Waldmann & Pascale 2014) to obtain the expected level of noise for different observational scenarios; these are then used as inputs for the NEMESIS atmospheric retrieval code and the retrieved atmospheric properties (temperature structure, composition and cloud properties) compared with the known input values, following the method of Barstow et al. (2013a). To correctly retrieve the temperature structure and composition of the atmosphere to within 2 {\sigma}, we find that we require: a single transit or eclipse of a hot Jupiter orbiting a sun-like (G2) star at 35 pc to constrain the terminator and dayside atmospheres; 20 transits or eclipses of a warm Jupiter orbiting a similar star; 10 transits/eclipses of a hot Neptune orbiting an M dwarf at 6 pc; and 30 transits or eclipses of a GJ1214b-like planet.

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

Atmospheric spectroscopy of extrasolar planets is an intricate business. Atmospheric signatures typically require a photometric precision of $1 \times 10^{-4}$ in flux over several hours. Such precision demands high instrument stability as well as an understanding of stellar variability and an optimal data reduction and removal of systematic noise. In the context of the EChO mission concept, we here discuss the data reduction and analysis pipeline developed for the EChO end-to-end simulator EChOSim. We present and discuss the step by step procedures required in order to obtain the final exoplanetary spectrum from the EChOSim`raw data' using a simulated observation of the secondary eclipse of the hot-Neptune 55 Cnc e.

The New IRAM KIDs Array (NIKA) is a pathfinder instrument devoted to millimetric astronomy. In 2009 it was the first multiplexed KID camera on the sky; currently it is installed at the focal plane of the IRAM 30-meters telescope at Pico Veleta (Spain). We present preliminary data from the last observational run and the ongoing developments devoted to the next NIKA-2 kilopixels camera, to be commissioned in 2015. We also report on the latest laboratory measurements, and recent improvements in detector cosmetics and read-out electronics. Furthermore, we describe a new acquisition strategy allowing us to improve the photometric accuracy, and the related automatic tuning procedure.

We present the first detection of the thermal Sunyaev-Zel'dovich (tSZ) effect from a cluster of galaxies performed with a KIDs (Kinetic Inductance Detectors) based instrument. The tSZ effect is a distortion of the black body CMB (Cosmic Microwave Background) spectrum produced by the inverse Compton interaction of CMB photons with the hot electrons of the ionized intra-cluster medium. The massive, intermediate redshift cluster RX J1347.5-1145 has been observed using NIKA (New IRAM KIDs arrays), a dual-band (140 and 240 GHz) mm-wave imaging camera, which exploits two arrays of hundreds of KIDs: the resonant frequencies of the superconducting resonators are shifted by mm-wave photons absorption. This tSZ cluster observation demonstrates the potential of the next generation NIKA2 instrument, being developed for the 30m telescope of IRAM, at Pico Veleta (Spain). NIKA2 will have 1000 detectors at 140GHz and 2x2000 detectors at 240GHz, providing in that band also a measurement of the linear polarization. NIKA2 will be commissioned in 2015.

We report on deep near-infrared observations obtained with the Wide Field Camera 3 (WFC3) onboard the Hubble Space Telescope (HST) of the first five confirmed gravitational lensing events discovered by the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). We succeed in disentangling the background galaxy from the lens to gain separate photometry of the two components. The HST data allow us to significantly improve on previous constraints of the mass in stars of the lensed galaxy and to perform accurate lens modelling of these systems, as described in the accompanying paper by Dye et al. We fit the spectral energy distributions of the background sources from near-IR to millimetre wavelengths and use the magnification factors estimated by Dye et al. to derive the intrinsic properties of the lensed galaxies. We find these galaxies to have star-formation rates of approximately 400 to 2000 M_sol/yr, with approximately (6-25)x10^10 M_sol of their baryonic mass already turned into stars. At these rates of star formation, all remaining molecular gas will be exhausted in less than 100 Myr, reaching a final mass in stars of a few 10^11 M_sol. These galaxies are thus proto-ellipticals caught during their major episode of star formation, and observed at the peak epoch z=1.5-3 of the cosmic star formation history of the Universe.

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

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

Photometric stability is a key requirement for time-resolved spectroscopic observations of transiting extrasolar planets. In the context of the Exoplanet Characterisation Observatory (EChO) mission design, we here present and investigate means of translating spacecraft pointing instabilities as well as temperature fluctuation of its optical chain into an overall error budget of the exoplanetary spectrum to be retrieved. Given the instrument specifications as of date, we investigate the magnitudes of these photometric instabilities in the context of simulated observations of the exoplanet HD189733b secondary eclipse.

We present a measurement of the angular power spectrum of the cosmic far-infrared background (CFIRB) anisotropies in one of the extragalactic fields of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) at 250, 350 and 500 \mu m bands. Consistent with recent measurements of the CFIRB power spectrum in Herschel-SPIRE maps, we confirm the existence of a clear one-halo term of galaxy clustering on arcminute angular scales with large-scale two-halo term of clustering at 30 arcminutes to angular scales of a few degrees. The power spectrum at the largest angular scales, especially at 250 \mu m, is contaminated by the Galactic cirrus. The angular power spectrum is modeled using a conditional luminosity function approach to describe the spatial distribution of unresolved galaxies that make up the bulk of the CFIRB. Integrating over the dusty galaxy population responsible for the background anisotropies, we find that the cosmic abundance of dust, relative to the critical density, to be between \Omega_dust=10^{-6} and 8 x 10^{-6} in the redshift range z ~ 0-3. This dust abundance is consistent with estimates of the dust content in the Universe using quasar reddening and magnification measurements in the SDSS.

We have measured the clustering properties of low-redshift (z < 0.3) sub-mm galaxies detected at 250 micron in the Herschel-ATLAS Science Demonstration Phase (SDP) field. We selected a sample for which we have high-quality spectroscopic redshifts, obtained from reliably matching the 250-micron sources to a complete (for r < 19.4) sample of galaxies from the GAMA database. Both the angular and spatial clustering strength are measured for all z < 0.3 sources as well as for five redshift slices with thickness delta z=0.05 in the range 0.05 < z < 0.3. Our measured spatial clustering length r_0 is comparable to that of optically-selected, moderately star-forming (blue) galaxies: we find values around 5 Mpc. One of the redshift bins contains an interesting structure, at z = 0.164.

We present new observations from Z-Spec, a broadband 185-305 GHz spectrometer, of five sub-millimeter bright lensed sources selected from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) science demonstration phase (SDP) catalog. We construct a redshift finding algorithm using combinations of the signal-to-noise of all the lines falling in the Z-Spec bandpass to determine redshifts with high confidence, even in cases where the signal-to-noise in individual lines is low. We measure the dust continuum in all sources and secure CO redshifts for four out of five (z~1.5-3). In one source, SDP.17, we tentatively identify two independent redshifts and a water line, confirmed at z=2.308. Our sources have properties characteristic of dusty starburst galaxies, with magnification-corrected star formation rates of 10^(2-3) M_sun/yr. Lower limits for the dust masses (~a few 10^8 M_sun) and spatial extents (~1 kpc equivalent radius) are derived from the continuum spectral energy distributions, corresponding to dust temperatures between 54 and 69K. In the LTE approximation, we derive relatively low CO excitation temperatures (< 100 K) and optical depths (tau<1). Performing a non-LTE excitation analysis using RADEX, we find that the CO lines measured by Z-Spec (from J=4->3 to 10->9, depending on the galaxy) localize the best solutions to either a high-temperature / low-density region, or a low-temperature / high-density region near the LTE solution, with the optical depth varying accordingly. Observations of additional CO lines, CO(1-0) in particular, are needed to constrain the non-LTE models.

We present high-spatial resolution imaging obtained with the Submillimeter Array (SMA) at 880um and the Keck Adaptive Optics (AO) system at Ks-band of a gravitationally lensed sub-millimeter galaxy (SMG) at z=4.243 discovered in the Herschel-Astrophysical Terahertz Large Area Survey. The SMA data (angular resolution ~0.6") resolve the dust emission into multiple lensed images, while the Keck AO Ks-band data (angular resolution ~0.1") resolve the lens into a pair of galaxies separated by 0.3". We present an optical spectrum of the foreground lens obtained with the Gemini-South telescope that provides a lens redshift of z_lens = 0.595 +/- 0.005. We develop and apply a new lens modeling technique in the visibility plane that shows that the SMG is magnified by a factor of mu = 4.1 +/- 0.2 and has an intrinsic infrared (IR) luminosity of L_IR = (2.1 +/- 0.2) x 10^13 Lsun. We measure a half-light radius of the background source of r_s = 4.4 +/- 0.5 kpc which implies an IR luminosity surface density of Sigma_IR = (3.4 +/- 0.9) x 10^11 Lsun kpc^-2, a value that is typical of z > 2 SMGs but significantly lower than IR luminous galaxies at z~0. The two lens galaxies are compact (r_lens ~ 0.9 kpc) early-types with Einstein radii of theta_E1 = 0.57 +/- 0.01 and theta_E2 = 0.40 +/- 0.01 that imply masses of M_lens1 = (7.4 +/- 0.5) x 10^10 Msun and M_lens2 = (3.7 +/- 0.3) x 10^10 Msun. The two lensing galaxies are likely about to undergo a dissipationless merger, and the mass and size of the resultant system should be similar to other early-type galaxies at z~0.6. This work highlights the importance of high spatial resolution imaging in developing models of strongly lensed galaxies discovered by Herschel.

We use spitzer-IRAC data to identify near-infrared counterparts to submillimeter galaxies detected with Herschel-SPIRE at 250um in the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). Using a likelihood ratio analysis we identify 146 reliable IRAC counterparts to 123 SPIRE sources out of the 159. We find that, compared to the field population, the SPIRE counterparts occupy a distinct region of 3.6 and 4.5um color-magnitude space, and we use this property to identify a further 23 counterparts to 13 SPIRE sources. The IRAC identification rate of 86% is significantly higher than those that have been demonstrated with wide-field ground-based optical and near-IR imaging of Herschel fields. We estimate a false identification rate of 3.6%, corresponding to 4 to 5 sources. Among the 73 counterparts that are undetected in SDSS, 57 have both 3.6 and 4.5um coverage. Of these 43 have [3.6] - [4.5]> 0 indicating that they are likely to be at z > 1.4. Thus, ~ 40% of identified SPIRE galaxies are likely to be high redshift (z > 1.4) sources. We discuss the statistical properties of the IRAC-identified SPIRE galaxy sample including far-IR luminosities, dust temperatures, star-formation rates, and stellar masses. The majority of our detected galaxies have 10^10 to 10^11 L_sun total IR luminosities and are not intense starbursting galaxies as those found at z ~ 2, but they have a factor of 2 to 3 above average specific star-formation rates compared to near-IR selected galaxy samples.

We identify near-infrared Ks band counterparts to Herschel-ATLAS sub-mm sources, using a preliminary object catalogue from the VISTA VIKING survey. The sub-mm sources are selected from the H-ATLAS Phase 1 catalogue of the GAMA 9h field, which includes all objects detected at 250, 350 or 500 um with the SPIRE instrument. We apply and discuss a likelihood ratio (LR) method for VIKING candidates within a search radius of 10" of the 22,000 SPIRE sources with a 5 sigma detection at 250 um. We find that 11,294(51%) of the SPIRE sources have a best VIKING counterpart with a reliability $R\ge 0.8$, and the false identification rate of these is estimated to be 4.2%. We expect to miss ~5% of true VIKING counterparts. There is evidence from Z-J and J-Ks colours that the reliable counterparts to SPIRE galaxies are marginally redder than the field population. We obtain photometric redshifts for ~68% of all (non-stellar) VIKING candidates with a median redshift of 0.405. Comparing to the results of the optical identifications supplied with the Phase I catalogue, we find that the use of medium-deep near-infrared data improves the identification rate of reliable counterparts from 36% to 51%.

A laboratory prototype spectral/spatial interferometer has been constructed to demonstrate the feasibility of the double Fourier technique at Far Infrared (FIR) wavelengths (0.15 - 1 THz). It is planned to use this demonstrator to investigate and validate important design features and data processing methods for future astronomical FIR interferometer instruments. In building this prototype we have had to address several key technologies to provide an end-end system demonstration of this double Fourier interferometer. We report on the first results taken when viewing single slit and double slit sources at the focus of a large collimator used to simulate real sources at infinity. The performance of the prototype instrument for these specific field geometries is analyzed to compare with the observed interferometric fringes and to demonstrate image reconstruction capabilities.

While the selection of strongly lensed galaxies with 500{\mu}m flux density S_500>100 mJy has proven to be rather straightforward (Negrello et al. 2010), for many applications it is important to analyze samples larger than the ones obtained when confining ourselves to such a bright limit. Moreover, only by probing to fainter flux densities is possible to exploit strong lensing to investigate the bulk of the high-z star-forming galaxy population. We describe HALOS (the Herschel-ATLAS Lensed Objects Selection), a method for efficiently selecting fainter candidate strongly lensed galaxies, reaching a surface density of ~1.5-2 deg^-2, i.e. a factor of about 4 to 6 higher than that at the 100 mJy flux limit. HALOS will allow the selection of up to ~1000 candidate strongly lensed galaxies (with amplifications \mu>2) over the full H-ATLAS survey area. Applying HALOS to the H-ATLAS Science Demonstration Phase field (~14.4 deg^2) we find 31 candidate strongly lensed galaxies, whose candidate lenses are identified in the VIKING near-infrared catalog. Using the available information on candidate sources and candidate lenses we tentatively estimate a ~72% purity of the sample. The redshift distribution of the candidate lensed sources is close to that reported for most previous surveys for lensed galaxies, while that of candidate lenses extends to substantially higher redshifts than found in the other surveys. The counts of candidate strongly lensed galaxies are also in good agreement with model predictions (Lapi et al. 2011). Even though a key ingredient of the method is the deep near-infrared VIKING photometry, we show that H-ATLAS data alone allow the selection of a similarly deep sample of candidate strongly lensed galaxies with an efficiency close to 50%; a slightly lower surface density (~1.45 deg^-2) can be reached with a ~70% efficiency.

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.

A dedicated mission to investigate exoplanetary atmospheres represents a major milestone in our quest to understand our place in the universe by placing our Solar System in context and by addressing the suitability of planets for the presence of life. EChO -the Exoplanet Characterisation Observatory- is a mission concept specifically geared for this purpose. EChO will provide simultaneous, multi-wavelength spectroscopic observations on a stable platform that will allow very long exposures. EChO will build on observations by Hubble, Spitzer and groundbased telescopes, which discovered the first molecules and atoms in exoplanetary atmospheres. EChO will simultaneously observe a broad enough spectral region -from the visible to the mid-IR- to constrain from one single spectrum the temperature structure of the atmosphere and the abundances of the major molecular species. The spectral range and resolution are tailored to separate bands belonging to up to 30 molecules to retrieve the composition and temperature structure of planetary atmospheres. The target list for EChO includes planets ranging from Jupiter-sized with equilibrium temperatures Teq up to 2000 K, to those of a few Earth masses, with Teq ~300 K. We have baselined a dispersive spectrograph design covering continuously the 0.4-16 micron spectral range in 6 channels (1 in the VIS, 5 in the IR), which allows the spectral resolution to be adapted from several tens to several hundreds, depending on the target brightness. The instrument will be mounted behind a 1.5 m class telescope, passively cooled to 50 K, with the instrument structure and optics passively cooled to ~45 K. EChO will be placed in a grand halo orbit around L2. We have also undertaken a first-order cost and development plan analysis and find that EChO is easily compatible with the ESA M-class mission framework.

Exploiting the Herschel-ATLAS Science Demonstration Phase (SDP) survey data, we have determined the luminosity functions (LFs) at rest-frame wavelengths of 100 and 250 micron and at several redshifts z>1, for bright sub-mm galaxies with star formation rates (SFR) >100 M_sun/yr. We find that the evolution of the comoving LF is strong up to z~2.5, and slows down at higher redshifts. From the LFs and the information on halo masses inferred from clustering analysis, we derived an average relation between SFR and halo mass (and its scatter). We also infer that the timescale of the main episode of dust-enshrouded star formation in massive halos (M_H>3*10^12 M_sun) amounts to ~7*10^8 yr. Given the SFRs, which are in the range 10^2-10^3 M_sun/yr, this timescale implies final stellar masses of order of 10^11-10^12 M_sun. The corresponding stellar mass function matches the observed mass function of passively evolving galaxies at z>1. The comparison of the statistics for sub-mm and UV selected galaxies suggests that the dust-free, UV bright phase, is >10^2 times shorter than the sub-mm bright phase, implying that the dust must form soon after the onset of star formation. Using a single reference Spectral Energy Distribution (SED; the one of the z~2.3 galaxy SMM J2135-0102), our simple physical model is able to reproduce not only the LFs at different redshifts > 1 but also the counts at wavelengths ranging from 250 micron to ~1 mm. Owing to the steepness of the counts and their relatively broad frequency range, this result suggests that the dispersion of sub-mm SEDs of z>1 galaxies around the reference one is rather small.

The Herschel survey, H-ATLAS, with its large areal coverage, has recently discovered a number of bright, strongly lensed high-z submillimeter galaxies. The strong magnification makes it possible to study molecular species other than CO, which are otherwise difficult to observe in high-z galaxies. Among the lensed galaxies already identified by H-ATLAS, the source J090302.9-014127B (SDP.17b) at z = 2.305 is remarkable due to its excitation conditions and a tentative detection of the H2O 202-111 emission line (Lupu et al. 2010). We report observations of this line in SDP.17b using the IRAM interferometer equipped with its new 277- 371GHz receivers. The H2O line is detected at a redshift of z = 2.3049+/-0.0006, with a flux of 7.8+/-0.5 Jy km s-1 and a FWHM of 250+/-60 km s-1. The new flux is 2.4 times weaker than the previous tentative detection, although both remain marginally consistent within 1.6-sigma. The intrinsic line luminosity and ratio of H2O(202-111)/CO8-7 seem comparable with those of the nearby starburst/enshrouded-AGN Mrk 231, suggesting that SDP.17b could also host a luminous AGN. The detection of a strong H2O 202-111 line in SDP.17b implies an efficient excitation mechanism of the water levels that must occur in very dense and warm interstellar gas.

We present Herschel-SPIRE Fourier Transform Spectrometer (FTS) and radio follow-up observations of two Herschel-ATLAS (H-ATLAS) detected strongly lensed distant galaxies. In one of the targeted galaxies H-ATLAS J090311.6+003906 (SDP.81) we detect [OIII] 88\mum and [CII] 158\mum lines at a signal-to-noise ratio of ~5. We do not have any positive line identification in the other fainter target H-ATLAS J091305.0-005343 (SDP.130). Currently SDP.81 is the faintest sub-mm galaxy with positive line detections with the FTS, with continuum flux just below 200 mJy in the 200-600 \mum wavelength range. The derived redshift of SDP.81 from the two detections is z=3.043 +/-0.012, in agreement with ground-based CO measurements. This is the first detection by Herschel of the [OIII] 88\mum line in a galaxy at redshift higher than 0.05. Comparing the observed lines and line ratios with a grid of photo-dissociation region (PDR) models with different physical conditions, we derive the PDR cloud density n ~ 2000 cm^{-3} and the far-UV ionizing radiation field G_0 ~ 200 (in units of the Habing field -- the local Galactic interstellar radiation field of 1.6x10^{-6} W/m^2). Using the CO derived molecular mass and the PDR properties we estimate the effective radius of the emitting region to be 500-700 pc. These characteristics are typical for star-forming, high redshift galaxies. The radio observations indicate that SDP.81 deviates significantly from the local FIR/radio correlation, which hints that some fraction of the radio emission is coming from an AGN. The constraints on the source size from millimiter-wave observations put a very conservative upper limit of the possible AGN contribution to less than 33%. These indications, together with the high [OIII]/FIR ratio and the upper limit of [OI] 63\mum/[CII] 158\mum suggest that some fraction of the ionizing radiation is likely to originate from an AGN.

We present the first direct and unbiased measurement of the evolution of the dust mass function of galaxies over the past 5 billion years of cosmic history using data from the Science Demonstration Phase of the Herschel-ATLAS. The sample consists of galaxies selected at 250{\mu}m which have reliable counterparts from SDSS at z < 0.5, and contains 1867 sources. Dust masses are calculated using both a single temperature grey-body model for the spectral energy distribution and also using a model with multiple temperature components. The dust temperature for either model shows no trend with redshift. Splitting the sample into bins of redshift reveals a strong evolution in the dust properties of the most massive galaxies. At z = 0.4 - 0.5, massive galaxies had dust masses about five times larger than in the local Universe. At the same time, the dust-to-stellar mass ratio was about 3-4 times larger, and the optical depth derived from fitting the UV-sub-mm data with an energy balance model was also higher. This increase in the dust content of massive galaxies at high redshift is difficult to explain using standard dust evolution models and requires a rapid gas consumption timescale together with either a more top-heavy IMF, efficient mantle growth, less dust destruction or combinations of all three. This evolution in dust mass is likely to be associated with a change in overall ISM mass, and points to an enhanced supply of fuel for star formation at earlier cosmic epochs.

We report ground-based follow-up observations of the exceptional source, ID141, one the brightest sources detected so far in the H-ATLAS cosmological survey. ID141 was observed using the IRAM 30-meter telescope and Plateau de Bure interferometer (PdBI), the Submillimeter Array (SMA) and the Atacama Pathfinder Experiment (APEX) submillimeter telescope to measure the dust continuum and emission lines of the main isotope of carbon monoxide and carbon ([C I] and [C II]). The detection of strong CO emission lines with the PdBI confirms that ID141 is at high redshift (z=4.243 +/- 0.001). The strength of the continuum and emission lines suggests that ID141 is gravitationally lensed. The width (Delta V (FWHM) ~ 800 km/s}) and asymmetric profiles of the CO and carbon lines indicate orbital motion in a disc or a merger. The properties derived for ID141 are compatible with a ultraluminous (L_FIR ~ 8.5 +/- 0.3 x 10^13/mu_L Lsun, where mu_L is the amplification factor, dense (n ~ 10^4 cm^-3) and warm (T_kin ~ 40K) starburst galaxy, with an estimated star-formation rate of (0.7 to 1.7) x 10^4/mu_L Msun/yr. The carbon emission lines indicate a dense (n ~ 10^4 cm^-3) Photo-Dominated Region, illuminated by a far-UV radiation field a few thousand times more intense than that in our Galaxy. In conclusion, the physical properties of the high-z galaxy, ID141, are remarkably similar to those of local ultraluminous infrared galaxies.

We present a technique to identify optical counterparts of 250 um-selected sources from the Herschel-ATLAS survey. Of the 6621 250 um > 32 mJy sources in our science demonstration catalogue we find that ~60 percent have counterparts brighter than r=22.4 mag in the Sloan Digital Sky Survey. Applying a likelihood ratio technique we are able to identify 2423 of the counterparts with a reliability R > 0.8. This is approximately 37 percent of the full 250 micron catalogue. We have estimated photometric redshifts for each of these 2423 reliable counterparts, while 1099 also have spectroscopic redshifts collated from several different sources, including the GAMA survey. We estimate the completeness of identifying counterparts as a function of redshift, and present evidence that 250 um-selected Herschel-ATLAS galaxies have a bimodal redshift distribution. Those with reliable optical identifications have a redshift distribution peaking at z ~ 0.25 +/- 0.05, while sub-mm colours suggest that a significant fraction with no counterpart above the r-band limit have z > 1. We also suggest a method for selecting populations of strongly-lensed high redshift galaxies. Our identifications are matched to UV--NIR photometry from the GAMA survey, and these data are available as part of the Herschel-ATLAS public data release.

We investigate the ultraviolet and optical properties and environment of low redshift galaxies detected in the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) science demonstration data. We use the Sloan Digital Sky Survey seventh release and the Galaxy And Mass Assembly database to select galaxies with r_Petro < 19.0 mag in the redshift range 0.02 < z < 0.2 and look for their submillimeter counterparts in H-ATLAS. Our results show that at low redshift, H-ATLAS detects mainly blue/star-forming galaxies with a minor contribution from red systems which are highly obscured by dust. In addition we find that the colour of a galaxy rather than the local density of its environment determines whether it is detectable by H-ATLAS. The average dust temperature of galaxies that are simultaneously detected by both PACS and SPIRE is 25K \pm 4K, independent of environment. This analysis provides a glimpse of the potential of the H-ATLAS data to investigate the submillimeter properties of galaxies in the local universe.

We report the detection of a significant excess in the surface density of far-infrared sources from the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) within ~1 Mpc of the centres of 66 optically-selected clusters of galaxies in the SDSS with <z>~0.25. From the analysis of the multiwavelength properties of their counterparts we conclude that the far-infrared emission is associated with dust-obscured star formation and/or active galactic nuclei within galaxies in the clusters themselves. The excess reaches a maximum at a radius of ~0.8 Mpc, where we find 1.0\pm0.3 S_250um>34 mJy sources on average per cluster above what would be expected for random field locations. If the far-infrared emission is dominated by star formation (as opposed to AGN) then this corresponds to an average star formation rate of ~7 M_sun/yr per cluster in sources with L_IR>5d10 L_sun. Although lensed sources make a negligible contribution to the excess signal, a fraction of the sources around the clusters could be gravitationally lensed, and we have identified a sample of potential cases of cluster-lensed Herschel sources that could be targeted in follow-up studies.

The Herschel-ATLAS is a survey of 550 square degrees with the Herschel Space Observatory in five far--infrared and submillimetre bands. The first data for the survey, observations of a field 4x4 sq. degrees in size, were taken during the Science Demonstration Phase, and reach a 5 sigma noise level of 33 mJy/beam at 250 microns. This paper describes the source extraction methods used to create the corresponding Science Demonstration Phase catalogue, which contains 6876 sources, selected at 250 microns, within ~14 sq. degrees. SPIRE sources are extracted using a new method specifically developed for Herschel data; PACS counterparts of these sources are identified using circular apertures placed at the SPIRE positions. Aperture flux densities are measured for sources identified as extended after matching to optical wavelengths. The reliability of this catalogue is also discussed, using full simulated maps at the three SPIRE bands. These show that a significant number of sources at 350 and 500 microns have undergone flux density enhancements of up to a factor of ~2, due mainly to source confusion. Correction factors are determined for these effects. The SDP dataset and corresponding catalogue will be available from http://www.h-atlas.org/.

We use the science demonstration field data of the Herschel-ATLAS to study how star formation, traced by the far-infrared Herschel data, is related to both the accretion luminosity and redshift of quasars selected from the Sloan Digital Sky Survey and the 2SLAQ survey. By developing a maximum likelihood estimator to investigate the presence of correlations between the far-infrared and optical luminosities we find evidence that the star-formation in quasar hosts is correlated with both redshift and quasar accretion luminosity. Assuming a relationship of the form L_IR \propto L_QSO^{\theta} (1 + z)^{\zeta}, we find {\theta} = 0.22 +/- 0.08 and {\zeta} = 1.6 +/- 0.4, although there is substantial additional uncertainty in {\zeta} of order +/- 1, due to uncertainties in the host galaxy dust temperature. We find evidence for a large intrinsic dispersion in the redshift dependence, but no evidence for intrinsic dispersion in the correlation between L_QSO and L_IR, suggesting that the latter may be due to a direct physical connection between star formation and black hole accretion. This is consistent with the idea that both the quasar activity and star formation are dependent on the same reservoir of cold gas, so that they are both affected by the influx of cold gas during mergers or heating of gas via feedback processes.

We use multiwavelength data from the Galaxy And Mass Assembly (GAMA) and Herschel ATLAS (H-ATLAS) surveys to compare the relationship between various dust obscuration measures in galaxies. We explore the connections between the ultraviolet (UV) spectral slope, $\beta$, the Balmer decrement, and the far infrared (IR) to $150\,$nm far ultraviolet (FUV) luminosity ratio. We explore trends with galaxy mass, star formation rate (SFR) and redshift in order to identify possible systematics in these various measures. We reiterate the finding of other authors that there is a large scatter between the Balmer decrement and the $\beta$ parameter, and that $\beta$ may be poorly constrained when derived from only two broad passbands in the UV. We also emphasise that FUV derived SFRs, corrected for dust obscuration using $\beta$, will be overestimated unless a modified relation between $\beta$ and the attenuation factor is used. Even in the optimum case, the resulting SFRs have a significant scatter, well over an order of magnitude. While there is a stronger correlation between the IR to FUV luminosity ratio and $\beta$ parameter than with the Balmer decrement, neither of these correlations are particularly tight, and dust corrections based on $\beta$ for high redshift galaxy SFRs must be treated with caution. We conclude with a description of the extent to which the different obscuration measures are consistent with each other as well as the effects of including other galactic properties on these correlations.

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.

We present physical properties of two submillimeter selected gravitationally lensed sources, identified in the Herschel Astrophysical Terahertz Large Area Survey. These submillimeter galaxies (SMGs) have flux densities > 100 mJy at 500 um, but are not visible in existing optical imaging. We fit light profiles to each component of the lensing systems in Spitzer IRAC 3.6 and 4.5 um data and successfully disentangle the foreground lens from the background source in each case, providing important constraints on the spectral energy distributions (SEDs) of the background SMG at rest-frame optical-near-infrared wavelengths. The SED fits show that these two SMGs have high dust obscuration with Av ~4 to 5 and star formation rates of ~100 M_sun/yr. They have low gas fractions and low dynamical masses compared to 850 um selected galaxies.

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.

We measure the projected cross-correlation between low redshift (z < 0.5) far-IR selected galaxies in the SDP field of the Herschel-ATLAS (H-ATLAS) survey and optically selected galaxies from the Galaxy and Mass Assembly (GAMA) redshift survey. In order to obtain robust correlation functions, we restrict the analysis to a subset of 969 out of 6900 H-ATLAS galaxies, which have reliable optical counterparts with r<19.4 mag and well-determined spectroscopic redshifts. The overlap region between the two surveys is 12.6 sq. deg; the matched sample has a median redshift of z ~ 0.2. The cross-correlation of GAMA and H-ATLAS galaxies within this region can be fitted by a power law, with correlation length r_0 ~ 4.63 +/- 0.51 Mpc. Comparing with the corresponding auto-correlation function of GAMA galaxies within the SDP field yields a relative bias (averaged over 2-8 Mpc) of H-ATLAS and GAMA galaxies of b_H/b_G ~ 0.6. Combined with clustering measurements from previous optical studies, this indicates that most of the low redshift H-ATLAS sources are hosted by halos with masses comparable to that of the Milky Way. The correlation function appears to depend on the 250 um luminosity, L_250, with bright (median luminosity \nu L_250 ~ 1.6 x 10^10 L_sun) objects being somewhat more strongly clustered than faint (\nu L_250 ~ 4.0 x 10^9 L_sun) objects. This implies that galaxies with higher dust-obscured star formation rates are hosted by more massive halos.

Gravitational lensing is a powerful astrophysical and cosmological probe and is particularly valuable at submillimeter wavelengths for the study of the statistical and individual properties of dusty starforming galaxies. However the identification of gravitational lenses is often time-intensive, involving the sifting of large volumes of imaging or spectroscopic data to find few candidates. We used early data from the Herschel Astrophysical Terahertz Large Area Survey to demonstrate that wide-area submillimeter surveys can simply and easily detect strong gravitational lensing events, with close to 100% efficiency.

We present a joint analysis of the overlapping BLAST 250, 350, 500um, and LABOCA 870um observations (from the LESS survey) of the ECDF-S. Out to z~3, the BLAST filters sample near the peak wavelength of far-infrared (FIR) emission from galaxies (rest-frame wavelengths ~60--200um), primarily produced by dust heated through absorption in star-forming clouds. However, identifying counterparts to individual BLAST peaks is very challenging, given the large beams (FWHM 36--60"). In contrast, the 870um observations have a significantly smaller 19" FWHM beam, and are sensitive to higher redshifts (z~1--5, and potentially beyond) due to the more favourable negative K-correction. We use the LESS data, as well as deep Spitzer and VLA imaging, to identify 118 individual sources that produce significant emission in the BLAST bands. We characterize the temperatures and FIR luminosities for a subset of 69 sources which have well-measured submm SEDs and redshift measurements out to z~3. For flux-limited sub-samples in each BLAST band, and a dust emissivity index \beta=2.0, we find a median temperature T=30K (all bands) as well as median redshifts: z=1.1 (interquartile range 0.2--1.9) for S250 > 40mJy; z=1.3 (interquartile range 0.6--2.1) for S350 > 30mJy; and z=1.6 (interquartile range 1.3--2.3) for S500 > 20mJy. Taking into account selection effects for our survey (a bias toward detecting lower-temperature galaxies), we find no evidence for evolution in the local FIR-temperature correlation out to z~2.5. Comparing with star-forming galaxy SED templates, about 8% of our sample appears to exhibit significant excesses in the radio and/or mid-IR, consistent with those sources harbouring an AGN. We describe the following techniques in two appendices: our `matched filter' for identifying sources in the presence of point-source confusion; and our approach for identifying counterparts using likelihood ratios.

We have reduced the data taken with the Spectral and Photometric Imaging Receiver (SPIRE) photometer on board the Herschel Space Observatory in the Science Demonstration Phase (SDP) of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). We describe the data reduction, which poses specific challenges, both because of the sheer size of the data, and because only two scans are made for each region. We implement effective solutions to process the bolometric timelines into maps, and show that correlations among detectors are negligible, and that the photometer is stable on time scales up to 250 s. This is longer than the time the telescope takes to cross the observed sky region, and it allows us to use naive binning methods for an optimal reconstruction of the sky emission. The maps have equal contribution of confusion and white instrumental noise, and the latter is estimated to 5.3, 6.4, and 6.7 mJy/beam (1-{\sigma}), at 250, 350, and 500 \mu{m}, respectively. This pipeline is used to reduce other H-ATLAS observations, as they became available, and we discuss how it can be used with the optimal map maker implemented in the Herschel Interactive Processing Environment (HIPE), to improve computational efficiency and stability. The SDP dataset is available from http://www.h-atlas.org/.

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 reduction of data taken with the PACS instrument on board the Herschel Space Observatory in the Science Demonstration Phase of the Herschel-ATLAS (H-ATLAS) survey, specifically data obtained for a 4x4-deg^2 region using Herschel's fast-scan (60 arcsec/s) parallel mode. We describe in detail a pipeline for data reduction using customised procedures within HIPE from data retrieval to the production of science-quality images. We found that the standard procedure for removing Cosmic-Ray glitches also removed parts of bright sources and so implemented an effective two-stage process to minimise these problems. The pronounced 1/f noise is removed from the timelines using 3.4- and 2.5-arcmin boxcar high-pass filters at 100 and 160-um. Empirical measurements of the point-spread function (PSF) are used to determine the encircled energy fraction as a function of aperture size. For the 100- and 160-um bands, the effective PSFs are ~9 and ~13 arcsec (FWHM), and the 90-per-cent encircled energy radii are 13 and 18 arcsec. Astrometric accuracy is good to ~<2 arcsec. The noise in the final maps is correlated between neighbouring pixels and rather higher than advertised prior to launch. For a pair of cross-scans, the 5-sigma point-source sensitivities are 125-165 mJy for 9-13-arcsec radius apertures at 100-um and 150-240 mJy for 13-18-arcsec radius apertures at 160-um.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

EBEX is a NASA-funded balloon-borne experiment designed to measure the polarization of the cosmic microwave background (CMB). Observations will be made using 1432 transition edge sensor (TES) bolometric detectors read out with frequency multiplexed SQuIDs. EBEX will observe in three frequency bands centered at 150, 250, and 410 GHz, with 768, 384, and 280 detectors in each band, respectively. This broad frequency coverage is designed to provide valuable information about polarized foreground signals from dust. The polarized sky signals will be modulated with an achromatic half wave plate (AHWP) rotating on a superconducting magnetic bearing (SMB) and analyzed with a fixed wire grid polarizer. EBEX will observe a patch covering ~1% of the sky with 8' resolution, allowing for observation of the angular power spectrum from \ell = 20 to 1000. This will allow EBEX to search for both the primordial B-mode signal predicted by inflation and the anticipated lensing B-mode signal. Calculations to predict EBEX constraints on r using expected noise levels show that, for a likelihood centered around zero and with negligible foregrounds, 99% of the area falls below r = 0.035. This value increases by a factor of 1.6 after a process of foreground subtraction. This estimate does not include systematic uncertainties. An engineering flight was launched in June, 2009, from Ft. Sumner, NM, and the long duration science flight in Antarctica is planned for 2011. These proceedings describe the EBEX instrument and the North American engineering flight.

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 present the hardware and software systems implementing autonomous operation, distributed real-time monitoring, and control for the EBEX instrument. EBEX is a NASA-funded balloon-borne microwave polarimeter designed for a 14 day Antarctic flight that circumnavigates the pole. To meet its science goals the EBEX instrument autonomously executes several tasks in parallel: it collects attitude data and maintains pointing control in order to adhere to an observing schedule; tunes and operates up to 1920 TES bolometers and 120 SQUID amplifiers controlled by as many as 30 embedded computers; coordinates and dispatches jobs across an onboard computer network to manage this detector readout system; logs over 3~GiB/hour of science and housekeeping data to an onboard disk storage array; responds to a variety of commands and exogenous events; and downlinks multiple heterogeneous data streams representing a selected subset of the total logged data. Most of the systems implementing these functions have been tested during a recent engineering flight of the payload, and have proven to meet the target requirements. The EBEX ground segment couples uplink and downlink hardware to a client-server software stack, enabling real-time monitoring and command responsibility to be distributed across the public internet or other standard computer networks. Using the emerging dirfile standard as a uniform intermediate data format, a variety of front end programs provide access to different components and views of the downlinked data products. This distributed architecture was demonstrated operating across multiple widely dispersed sites prior to and during the EBEX engineering flight.

We have measured star formation histories (SFHs) and stellar masses of galaxies detected by the Balloon-borne Large Aperture Sub-millimetre Telescope (BLAST) over approximately 9 square degrees centred on the Chandra Deep Field South. We have applied the recently developed SFH reconstruction method of Dye et al. to optical, near-infrared and mid-infrared photometry of 92 BLAST galaxies. We find significant differences between the SFHs of low mass (<10^11 M_sol) and high mass (>10^11 M_sol) galaxies. On average, low mass galaxies exhibit a dominant late burst of star formation which creates a large fraction of their stellar mass. Conversely, high mass systems tend to have a significant amount of stellar mass that formed much earlier. We also find that the high mass SFHs evolve more strongly than the low mass SFHs. These findings are consistent with the phenomenon of downsizing observed in optically selected samples of galaxies.

We describe the cryostat and supporting electronics for the EBEX experiment. EBEX is a balloon-borne polarimeter designed to measure the B-mode polarization of the cosmic microwave background radiation. The instrument includes a 1.5 meter Gregorian-type telescope and 1432 bolometric transition edge sensor detectors operating at 0.3 K. Electronics for monitoring temperatures and controlling cryostat refrigerators is read out over CANbus. A timing system ensures the data from all subsystems is accurately synchronized. EBEX completed an engineering test flight in June 2009 during which the cryogenics and supporting electronics performed according to predictions. The temperatures of the cryostat were stable, and an analysis of a subset of the data finds no scan synchronous signal in the cryostat temperatures. Preparations are underway for an Antarctic flight.

Aims: We aim to demonstrate that the Herschel ATLAS (H-ATLAS) is suitable for a blind and unbiased survey for debris disks by identifying candidate debris disks associated with main sequence stars in the initial science demonstration field of the survey. We show that H-ATLAS reveals a population of far-infrared/sub-mm sources that are associated with stars or star-like objects on the SDSS main-sequence locus. We validate our approach by comparing the properties of the most likely candidate disks to those of the known population. Methods: We use a photometric selection technique to identify main sequence stars in the SDSS DR7 catalogue and a Bayesian Likelihood Ratio method to identify H-ATLAS catalogue sources associated with these main sequence stars. Following this photometric selection we apply distance cuts to identify the most likely candidate debris disks and rule out the presence of contaminating galaxies using UKIDSS LAS K-band images. Results: We identify 78 H-ATLAS sources associated with SDSS point sources on the main-sequence locus, of which two are the most likely debris disk candidates: H-ATLAS J090315.8 and H-ATLAS J090240.2. We show that they are plausible candidates by comparing their properties to the known population of debris disks. Our initial results indicate that bright debris disks are rare, with only 2 candidates identified in a search sample of 851 stars. We also show that H-ATLAS can derive useful upper limits for debris disks associated with Hipparcos stars in the field and outline the future prospects for our debris disk search programme.

We have determined the luminosity function of 250um-selected galaxies detected in the ~14 sq.deg science demonstration region of the Herschel-ATLAS project out to a redshift of z=0.5. Our findings very clearly show that the luminosity function evolves steadily out to this redshift. By selecting a sub-group of sources within a fixed luminosity interval where incompleteness effects are minimal, we have measured a smooth increase in the comoving 250um luminosity density out to z=0.2 where it is 3.6+1.4-0.9 times higher than the local value.

To investigate the poorly constrained sub-mm counts and spectral properties of blazars we searched for these in the Herschel-ATLAS (H-ATLAS) science demostration phase (SDP) survey catalog. We cross-matched 500$\mu$m sources brighter than 50 mJy with the FIRST radio catalogue. We found two blazars, both previously known. Our study is among the first blind blazar searches at sub-mm wavelengths, i.e., in the spectral regime where little is still known about the blazar SEDs, but where the synchrotron peak of the most luminous blazars is expected to occur. Our early results are consistent with educated extrapolations of lower frequency counts and question indications of substantial spectral curvature downwards and of spectral upturns at mm wavelengths. One of the two blazars is identified with a Fermi/LAT $\gamma$-ray source and a WMAP source. The physical parameters of the two blazars are briefly discussed.These observations demonstrate that the H-ATLAS survey will provide key information about the physics of blazars and their contribution to sub-mm counts.

We present colour-colour diagrams of detected sources in the Herschel-ATLAS Science Demonstration Field from 100 to 500 microns using both PACS and SPIRE. We fit isothermal modified black bodies to the spectral energy distribution (SED) to extract the dust temperature of sources with counterparts in Galaxy And Mass Assembly (GAMA) or SDSS surveys with either a spectroscopic or a photometric redshift. For a subsample of 330 sources detected in at least three FIR bands with a significance greater than 3 $\sigma$, we find an average dust temperature of $(28 \pm 8)$K. For sources with no known redshift, we populate the colour-colour diagram with a large number of SEDs generated with a broad range of dust temperatures and emissivity parameters, and compare to colours of observed sources to establish the redshift distribution of this sample. For another subsample of 1686 sources with fluxes above 35 mJy at 350 microns and detected at 250 and 500 microns with a significance greater than 3$\sigma$, we find an average redshift of $2.2 \pm 0.6$.

We present a derivation of the star formation rate per comoving volume of quasar host galaxies, derived from stacking analyses of far-infrared to mm-wave photometry of quasars with redshifts 0<z<6 and absolute I-band magnitudes -22>I_AB>-32. We use the science demonstration observations of the first ~16 deg^2 from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) in which there are 240 quasars from the Sloan Digital Sky Survey (SDSS) and a further 171 from the 2dF-SDSS LRG and QSO (2SLAQ) survey. We supplement this data with a compilation of data from IRAS, ISO, Spitzer, SCUBA and MAMBO. H-ATLAS alone statistically detects the quasars in its survey area at >5sigma at 250, 350 and 500um. From the compilation as a whole we find striking evidence of downsizing in quasar host galaxy formation: low-luminosity quasars with absolute magnitudes in the range -22>I_AB>-24 have a comoving star formation rate (derived from 100um rest-frame luminosities) peaking between redshifts of 1 and 2, while high-luminosity quasars with I_AB<-26 have a maximum contribution to the star formation density at z~3. The volume-averaged star formation rate of -22>I_AB>-24 quasars evolves as (1+z)^{2.3 +/- 0.7} at z<2, but the evolution at higher luminosities is much faster reaching (1+z)^{10 +/- 1} at -26>I_AB>-28. We tentatively interpret this as a combination of a declining major merger rate with time and gas consumption reducing fuel for both black hole accretion and star formation.

Aims. The Herschel-ATLAS survey (H-ATLAS) will be the largest area survey to be undertaken by the Herschel Space Observatory. It will cover 550 sq. deg. of extragalactic sky at wavelengths of 100, 160, 250, 350 and 500 microns when completed, reaching flux limits (5 sigma) from 32 to 145mJy. We here present galaxy number counts obtained for SPIRE observations of the first ~14 sq. deg. observed at 250, 350 and 500 microns. Methods. Number counts are a fundamental tool in constraining models of galaxy evolution. We use source catalogs extracted from the H-ATLAS maps as the basis for such an analysis. Correction factors for completeness and flux boosting are derived by applying our extraction method to model catalogs and then applied to the raw observational counts. Results. We find a steep rise in the number counts at flux levels of 100-200mJy in all three SPIRE bands, consistent with results from BLAST. The counts are compared to a range of galaxy evolution models. None of the current models is an ideal fit to the data but all ascribe the steep rise to a population of luminous, rapidly evolving dusty galaxies at moderate to high redshift.

We present measurements of the angular correlation function of galaxies selected from the first field of the H-ATLAS survey. Careful removal of the background from galactic cirrus is essential, and currently dominates the uncertainty in our measurements. For our 250 micron-selected sample we detect no significant clustering, consistent with the expectation that the 250 micron-selected sources are mostly normal galaxies at z<~ 1. For our 350 micron and 500 micron-selected samples we detect relatively strong clustering with correlation amplitudes A of 0.2 and 1.2 at 1', but with relatively large uncertainties. For samples which preferentially select high redshift galaxies at z~2-3 we detect significant strong clustering, leading to an estimate of r_0 ~ 7-11 h^{-1} Mpc. The slope of our clustering measurements is very steep, delta~2. The measurements are consistent with the idea that sub-mm sources consist of a low redshift population of normal galaxies and a high redshift population of highly clustered star-bursting galaxies.

We use Herschel PACS and SPIRE observations of the edge-on spiral galaxy UGC 4754, taken as part of the H-ATLAS SDP observations, to investigate the dust energy balance in this galaxy. We build detailed SKIRT radiative models based on SDSS and UKIDSS maps and use these models to predict the far-infrared emission. We find that our radiative transfer model underestimates the observed FIR emission by a factor of two to three. Similar discrepancies have been found for other edge-on spiral galaxies based on IRAS, ISO, and SCUBA data. Thanks to the good sampling of the SED at FIR wavelengths, we can rule out an underestimation of the FIR emissivity as the cause for this discrepancy. Instead we support highly obscured star formation that contributes little to the optical extinction as a more probable explanation.

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.

The physical properties of galactic cirrus emission are not well characterized. BOOMERanG is a balloon-borne experiment designed to study the cosmic microwave background at high angular resolution in the millimeter range. The BOOMERanG 245 and 345GHz channels are sensitive to interstellar signals, in a spectral range intermediate between FIR and microwave frequencies. We look for physical characteristics of cirrus structures in a region at high galactic latitudes (b~-40{\deg}) where BOOMERanG performed its deepest integration, combining the BOOMERanG data with other available datasets at different wavelengths. We have detected eight emission patches in the 345 GHz map, consistent with cirrus dust in the Infrared Astronomical Satellite maps. The analysis technique we have developed allows to identify the location and the shape of cirrus clouds, and to extract the flux from observations with different instruments at different wavelengths and angular resolutions. We study the integrated flux emitted from these cirrus clouds using data from Infrared Astronomical Satellite (IRAS), DIRBE, BOOMERanG and Wilkinson Microwave Anisotropy Probe (WMAP) in the frequency range 23-3000 GHz (13 mm 100 micron wavelength). We fit the measured spectral energy distributions with a combination of a grey body and a power-law spectra considering two models for the thermal emission. The temperature of the thermal dust component varies in the 7 - 20 K range and its emissivity spectral index is in the 1 - 5 range. We identified a physical relation between temperature and spectral index as had been proposed in previous works. This technique can be proficiently used for the forthcoming Planck and Herschel missions data.

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) has recently surveyed ~8.7 deg^2 centered on GOODS-South at 250, 350, and 500 microns. In Dye et al. (2009) we presented the catalogue of sources detected at 5-sigma in at least one band in this field and the probable counterparts to these sources in other wavebands. In this paper, we present the results of a redshift survey in which we succeeded in measuring redshifts for 82 of these counterparts. The spectra show that the BLAST counterparts are mostly star-forming galaxies but not extreme ones when compared to those found in the Sloan Digital Sky Survey. Roughly one quarter of the BLAST counterparts contain an active nucleus. We have used the spectroscopic redshifts to carry out a test of the ability of photometric redshift methods to estimate the redshifts of dusty galaxies, showing that the standard methods work well even when a galaxy contains a large amount of dust. We have also investigated the cases where there are two possible counterparts to the BLAST source, finding that in at least half of these there is evidence that the two galaxies are physically associated, either because they are interacting or because they are in the same large-scale structure. Finally, we have made the first direct measurements of the luminosity function in the three BLAST bands. We find strong evolution out to z=1, in the sense that there is a large increase in the space-density of the most luminous galaxies. We have also investigated the evolution of the dust-mass function, finding similar strong evolution in the space-density of the galaxies with the largest dust masses, showing that the luminosity evolution seen in many wavebands is associated with an increase in the reservoir of interstellar matter in galaxies.

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.

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.

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.

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.

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.

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

The Sunyaev-Zel'dovich (SZ) effect is the inverse Compton-scattering of cosmic microwave background (CMB) photons by hot electrons in the intervening gas throughout the universe. The effect has a distinct spectral signature that allows its separation from other signals in multifrequency CMB datasets. Using CMB anisotropies measured at three frequencies by the BOOMERanG 2003 flight we constrain SZ fluctuations in the 10 arcmin to 1 deg angular range. Propagating errors and potential systematic effects through simulations, we obtain an overall upper limit of 15.3 uK (2 sigma) for rms SZ fluctuations in a broad bin between multipoles of of 250 and 1200 at the Rayleigh-Jeans (RJ) end of the spectrum. When combined with other CMB anisotropy and SZ measurements, we find that the local universe normalization of the density perturbations is sigma-8(SZ) < 0.96 at the 95% confidence level, consistent with sigma-8 determined from primordial perturbations.

We use Minkowski Functionals (MF) to constrain a primordial non-Gaussian contribution to the CMB intensity field as observed in the 150 GHz and 145 GHz BOOMERanG maps from the 1998 and 2003 flights, respectively, performing for the first time a joint analysis of the two datasets. A perturbative expansion of the MF formulae in the limit of a weakly non-Gaussian field yields analytical formulae, derived by Hikage et al. (2006), which can be used to constrain the coupling parameter f_NL without the need for non-Gaussian simulations. We find -1020<f_NL<390 at 95% CL, significantly improving the previous constraints by De Troia et al. (2007) on the BOOMERanG 2003 dataset. These are the best f_NL limits to date for suborbital probes.

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.

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.

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.

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.

We analyze the BOOMERanG 2003 (B03) 145 GHz temperature map to constrain the amplitude of a non Gaussian, primordial contribution to CMB fluctuations. We perform a pixel space analysis restricted to a portion of the map chosen in view of high sensitivity, very low foreground contamination and tight control of systematic effects. We set up an estimator based on the three Minkowski functionals which relies on high quality simulated data, including non Gaussian CMB maps. We find good agreement with the Gaussian hypothesis and derive the first limits based on BOOMERanG data for the non linear coupling parameter f_NL as -300<f_NL<650 at 68% CL and -800<f_NL<1050 at 95% CL.

The trispectrum of the cosmic microwave background can be used to assess the level of non-Gaussianity on cosmological scales. It probes the fourth order moment, as a function of angular scale, of the probability distribution function of fluctuations and has been shown to be sensitive to primordial non-gaussianity, secondary anisotropies (such as the Ostriker-Vishniac effect) and systematic effects (such as astrophysical foregrounds). In this paper we develop a formalism for estimating the trispectrum from high resolution sky maps which incorporates the impact of finite sky coverage. This leads to a series of operations applied to the data set to minimize the effects of contamination due to the Gaussian component and correlations between estimates at different scales. To illustrate the effect of the estimation process, we apply our procedure to the BOOMERanG data set and show that it is consistent with Gaussianity. This work presents the first estimation of the CMB trispectrum on sub-degree scales.

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 the cosmological parameters from the CMB intensity and polarization power spectra of the 2003 Antarctic flight of the BOOMERANG telescope. The BOOMERANG data alone constrains the parameters of the $\Lambda$CDM model remarkably well and is consistent with constraints from a multi-experiment combined CMB data set. We add LSS data from the 2dF and SDSS redshift surveys to the combined CMB data set and test several extensions to the standard model including: running of the spectral index, curvature, tensor modes, the effect of massive neutrinos, and an effective equation of state for dark energy. We also include an analysis of constraints to a model which allows a CDM isocurvature admixture.

We present the BOOMERanG-03 experiment and maps of the Stokes parameters I, Q, U of the microwave sky obtained during a 14 day balloon flight in 2003. Three regions of the southern sky were surveyed: a deep survey (~ 90 square degrees) and a shallow survey (~ 750 square degrees) at high Galactic latitudes (both centered at RA ~ 5.5 h, dec ~ -45 deg) and a survey of ~ 300 square degrees across the Galactic plane at RA ~ 9.1 h, dec ~ -47 deg. All three surveys were carried out in three wide frequency bands centered at 145, 245 and 345 GHz, with an angular resolution of ~ 10'. The 145 GHz maps of Stokes I are dominated by Cosmic Microwave Background (CMB) temperature anisotropy, which is mapped with high signal to noise ratio. The measured anisotropy pattern is consistent with the pattern measured in the same region by BOOMERanG-98 and by WMAP. The 145 GHz maps of Stokes Q and U provide a robust statistical detection of polarization of the CMB when subjected to a power spectrum analysis. This amplitude of the polarization is consistent with that of the CMB in the $\Lambda$CDM cosmological scenario. At 145 GHz, in the CMB surveys, the intensity and polarization of the astrophysical foregrounds are found to be negligible with respect to the cosmological signal. At 245 and 345 GHz we detect ISD emission correlated to the 3000 GHz IRAS/DIRBE maps, and give upper limits for any other non-CMB component. We also present intensity maps of the surveyed section of the Galactic plane. These are compared to monitors of different interstellar components, showing that a variety of emission mechanisms is present in that region.

We report measurements of the CMB polarization power spectra from the January 2003 Antarctic flight of BOOMERANG. The primary results come from six days of observation of a patch covering 0.22% of the sky centered near R.A. = 82.5 deg., Dec= -45 deg. The observations were made using four pairs of polarization sensitive bolometers operating in bands centered at 145 GHz. Using two independent analysis pipelines, we measure a non-zero <EE> signal in the range 100< l <1000 with a significance 4.8-sigma, a 2-sigma upper limit of 8.6 uK^2 for any <BB> contribution, and a 2-sigma upper limit of 7.0 uK^2 for the <EB> spectrum. Estimates of foreground intensity fluctuations and the non-detection of <BB> and <EB> signals rule out any significant contribution from galactic foregrounds. The results are consistent with a Lambda-CDM cosmology seeded by adiabatic perturbations. We note that this is the first detection of CMB polarization with bolometric detectors.

We present a measurement of the temperature-polarization angular cross power spectrum, <TE>, of the Cosmic Microwave Background. The result is based on $\sim 200$ hours of data from 8 polarization sensitive bolometers operating at 145 GHz during the 2003 flight of BOOMERANG. We detect a significant <TE> correlation in the $\ell$-range between 50 and 950 with a statistical significance > 3.5 sigma. Contamination by polarized foreground emission and systematic effects are negligible in comparison with statistical uncertainty. The spectrum is consistent with previous detections and with the "concordance model" that assumes adiabatic initial conditions. This is the first measurement of <TE> using bolometric detectors.

We report on observations of the Cosmic Microwave Background (CMB) obtained during the January 2003 flight of Boomerang . These results are derived from 195 hours of observation with four 145 GHz Polarization Sensitive Bolometer (PSB) pairs, identical in design to the four 143 GHz Planck HFI polarized pixels. The data include 75 hours of observations distributed over 1.84% of the sky with an additional 120 hours concentrated on the central portion of the field, itself representing 0.22% of the full sky. From these data we derive an estimate of the angular power spectrum of temperature fluctuations of the CMB in 24 bands over the multipole range (50 < l < 1500). A series of features, consistent with those expected from acoustic oscillations in the primordial photon-baryon fluid, are clearly evident in the power spectrum, as is the exponential damping of power on scales smaller than the photon mean free path at the epoch of last scattering (l > 900). As a consistency check, the collaboration has performed two fully independent analyses of the time ordered data, which are found to be in excellent agreement.

The recent study of BOOMERanG 150 GHz Cosmic Microwave Background (CMB) radiation maps have detected ellipticity of the temperature anisotropy spots independent on the temperature threshold. The effect has been found for spots up to several degrees in size, where the biases of the ellipticity estimator and of the noise are small. To check the effect, now we have studied, with the same algorithm and in the same sky region, the WMAP maps. We find ellipticity of the same average value also in WMAP maps, despite of the different sensitivity of the two experiments to low multipoles. Large spot elongations had been detected also for the COBE-DMR maps. If this effect is due to geodesic mixing and hence due to non precisely zero curvature of the hyperbolic Universe, it can be linked to the origin of WMAP low multipoles anomaly.

We have measured the ellipticity of several degree scale anisotropies in the BOOMERanG maps of the Cosmic Microwave Background (CMB) at 150 GHz. The average ellipticity is around 2.6-2.7. The biases of the estimator of the ellipticity and for the noise are small in this case. Large spot elongation had been detected also for COBE-DMR maps. If this effect is due to geodesic mixing, it would indicate a non precisely zero curvature of the Universe which is among the discussed reasons of the WMAP low multipole anomaly. Both effects are related to the diameter of the Universe: the geodesics mixing through hyperbolic geometry, low multipoles through boundary conditions.This common reason can also be related with the origin of the the cosmological constant: the modes of vacuum fluctuations conditioned by the boundary conditions lead to a value of the cosmological constant being in remarkable agreement with the supernovae observations.

In the 1998-99 flight, BOOMERanG has produced maps of $\sim 4 %$ of the sky at high Galactic latitudes, at frequencies of 90, 150, 240 and 410 GHz, with resolution $\simgt 10'$. The faint structure of the Cosmic Microwave Background at horizon and sub-horizon scales is evident in these maps. These maps compare well to the maps recently obtained at lower frequencies by the WMAP experiment. Here we compare the amplitude and morphology of the structures observed in the two sets of maps. We also outline the polarization sensitive version of BOOMERanG, which was flown early this year to measure the linear polarization of the microwave sky at 150, 240 and 350 GHz.

The properties of the Cosmic Microwave Background (CMB) maps carry valuable cosmological information. Here we report the results of the analysis hot and cold CMB anisotropy spots in the BOOMERanG 150 GHz map in terms of number, area, ellipticity, vs. temperature threshold. We carried out this analysis for the map obtained by summing independent measurement channels (signal plus noise map) and for a comparison map (noise only map) obtained by differencing the same channels. The anisotropy areas (spots) have been identified for both maps for various temperature thresholds and a catalog of the spots has been produced. The orientation (obliquity) of the spots is random for both maps. We computed the mean elongation of spots obtained from the maps at a given temperature threshold using a simple estimator. We found that for the sum map there is a region of temperature thresholds where the average elongation is not dependent on the threshold. Its value is ~ 2.3 for cold areas and ~ 2.2 for hot areas. This is a non-trivial result. The bias of the estimator is less than 0.4 for areas of size less than 30', and smaller for larger areas. The presence of noise also biases the ellipticity by less than 0.3. These biases have not been subtracted in the results quoted above. The threshold independent and random obliquity behaviour in the sum map is stable against pointing reconstruction accuracy and noise level of the data, thus confirming that these are actual properties of the dataset. The data used here give a hint of high ellipticity for the largest spots. Analogous elongation properties of CMB anisotropies had been detected for COBE-DMR 4-year data. If this is due to geodesics mixing, it would point to a non zero curvature of the Universe.

BOOMERANG is a balloon-borne telescope designed for long duration (LDB) flights around Antarctica. The second LDB Flight of BOOMERANG took place in January 2003. The primary goal of this flight was to measure the polarization of the CMB. The receiver uses polarization sensitive bolometers at 145 GHz. Polarizing grids provide polarization sensitivity at 245 and 345 GHz. We describe the BOOMERANG telescope noting changes made for 2003 LDB flight, and discuss some of the issues involved in the measurement of polarization with bolometers. Lastly, we report on the 2003 flight and provide an estimate of the expected results.

We describe the Millimeter wave Anisotropy eXperiment IMaging Array (MAXIMA), a balloon-borne experiment designed to measure the temperature anisotropy of the Cosmic Microwave Background (CMB) on angular scales of 10' to 5 degrees . MAXIMA mapped the CMB using 16 bolometric detectors in spectral bands centered at 150 GHz, 240 GHz, and 410 GHz, with 10' resolution at all frequencies. The combined receiver sensitivity to CMB anisotropy was ~40 microK/rt(sec). Systematic parasitic contributions were minimized by using four uncorrelated spatial modulations, thorough crosslinking, multiple independent CMB observations, heavily baffled optics, and strong spectral discrimination. Pointing reconstruction was accurate to 1', and absolute calibration was better than 4%. Two MAXIMA flights with more than 8.5 hours of CMB observations have mapped a total of 300 deg^2 of the sky in regions of negligible known foreground emission. MAXIMA results have been released in previous publications. MAXIMA maps, power spectra and correlation matrices are publicly available at http://cosmology.berkeley.edu/maxima