Christopher Semisch

Over half of the far-infrared background light comes from galaxies at z ≥ 1.2

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 μm in the rest frame. At 1 ≤ z ≤ 4, the peak is redshifted to wavelengths between 200 and 500 μm. 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 μm. Combining our results at 500 μm with those at 24 μm, we determine that all of the FIRB comes from individual galaxies, with galaxies at z ≥ 1.2 accounting for 70% of it. As expected, at the longest wavelengths the signal is dominated by ultraluminous galaxies at z > 1.

BLAST: RESOLVING THE COSMIC SUBMILLIMETER BACKGROUND

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) has made 1 deg2, 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 μm sources from the Far-Infrared Deep Extragalactic Legacy Survey, 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 μm, respectively. These numbers are more precise than previous estimates of the cosmic infrared background (CIB) and are consistent with 24 μm-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 μm. At all BLAST wavelengths, the relative intensity of high-z sources is higher for 24 μm-faint sources than that for 24 μm-bright sources. Galaxies identified as active galactic nuclei (AGNs) by their Infrared Array Camera colors are 1.6-2.6 times brighter than the average population at 250-500 μm, consistent with what is found for X-ray-selected AGNs. BzK-selected galaxies are found to be moderately brighter than typical 24 μm-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.

BLAST: the far‐infrared/radio correlation in distant 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 and the Giant Metre-wave Radio Telescope (GMRT) in the Extended Chandra Deep Field South (ECDFS). For a catalogue of BLAST 250-μm-selected galaxies, we remeasure the 70–870-μm flux densities at the positions of their most likely 24-μm 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_(250 μm)/S_(1400 MHz)]), and the bolometric equivalent, q_(IR). At z ≈ 0.6 , where our 250-μm filter probes rest-frame 160-μm 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-μm- and radio-selected galaxies. The difference between q_(IR) seen for 250-μm- and radio-selected galaxies suggests that star formation provides most of the IR luminosity in ≲100-μJy radio galaxies, but rather less for those in the mJy regime. For the 24-μm 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) ∝ (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.

Evidence for environmental changes in the submillimeter dust opacity

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

SANEPIC: a mapmaking method for time stream data from large arrays

We describe a mapmaking method that 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 on several previous methods but focuses specifically on the regime where there are a large number of detectors sampling the same map of the sky, and explicitly allowing for the possibility of strong correlations between the detector time streams. We provide real and simulated examples of how well this method performs compared with more simplistic mapmakers 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.

BLAST: A Far-Infrared Measurement of the History 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 Submillimeter 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 μm background is generated at z lsim 1 and the 500 μm background generated at z gsim 1. A significant growth is observed in the mean luminosity from ~109-1012 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.

Radio and mid-infrared identification of BLAST source counterparts in the Chandra Deep Field South

We have identified radio and/or mid-infrared counterparts to 198 out of 350 sources detected at ≥5σ over ~9 deg2 centered on the Chandra Deep Field South by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) at 250, 350, and 500 μm. We have matched 114 of these counterparts to optical sources with previously derived photometric redshifts and fitted spectral energy distributions to the BLAST fluxes and fluxes at 70 and 160 μm acquired with the Spitzer Space Telescope. In this way, we have constrained dust temperatures, total far-infrared/submillimeter 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 submillimeter sources detected in surveys conducted at 850 μm. 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 multiwavelength catalog of ≥5σ BLAST sources contained within the complete ~9 deg2 survey area.

BLAST: THE REDSHIFT SURVEY

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) has recently surveyed sime8.7 deg2 centered on Great Observatories Origins Deep Survey-South at 250, 350, and 500 μm. In Dye et al., we presented the catalog of sources detected at 5σ 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.

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) 2005: Calibration and Targeted Sources

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) operated successfully during a 100 hr flight from northern Sweden in 2005 June (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, W75N, 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 μm 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 Blast Survey of the Vela Molecular Cloud: Physical Properties of the Dense Cores in Vela-D

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