Intensity mapping provides a unique means to probe the epoch of reionization (EoR), when the neutral intergalactic medium was ionized by energetic photons emitted from the first galaxies. The C ii ...158 m fine-structure line is typically one of the brightest emission lines of star-forming galaxies and thus a promising tracer of the global EoR star formation activity. However, C ii intensity maps at 6 z 8 are contaminated by interloping CO rotational line emission (3 ≤ Jupp ≤ 6) from lower-redshift galaxies. Here we present a strategy to remove the foreground contamination in upcoming C ii intensity mapping experiments, guided by a model of CO emission from foreground galaxies. The model is based on empirical measurements of the mean and scatter of the total infrared luminosities of galaxies at z < 3 and with stellar masses M * > 10 8 M selected in the K-band from the COSMOS/UltraVISTA survey, which can be converted to CO line strengths. For a mock field of the Tomographic Ionized-carbon Mapping Experiment, we find that masking out the "voxels" (spectral-spatial elements) containing foreground galaxies identified using an optimized CO flux threshold results in a z-dependent criterion m K AB 22 (or M * 10 9 M ) at z < 1 and makes a C ii/COtot power ratio of 10 at k = 0.1 h/Mpc achievable, at the cost of a moderate 8% loss of total survey volume.
We quantify the fraction of the cosmic infrared background (CIB) that originates from galaxies identified in the UV/optical/near-infrared by stacking 81,250 (~35.7 arcmin super(-2)) K-selected ...sources (K sub(AB) < 24.0) split according to their rest-frame U - V versus V - J colors into 72,216 star-forming and 9034 quiescent galaxies, on maps from Spitzer/MIPS (24 mu m), Herschel/PACS (100, 160 mu m), Herschel/SPIRE (250, 350, 500 mu m), and AzTEC (1100 mu m). Of that total, about 95% originates from star-forming galaxies, while the remaining 5% is from apparently quiescent galaxies. The CIB at lambda lap 200 mu m appears to be sourced predominantly from galaxies at z lap 1, while at lambda gap 200 mu the bulk originates from 1 lap z lap 2. Stacking analyses were performed using simstack, a novel algorithm designed to account for possible biases in the stacked flux density due to clustering. It is made available to the public at www.astro.caltech.edu/~viero/viero_homepage/toolbox.html.
The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) is a suborbital 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 three arrays, observes simultaneously in broadband (30%) spectral windows at 250, 350, and 500 mum. The optical design is based on a 2 m diameter telescope, providing a diffraction-limited resolution of 30 super(image ) at 250 mum. The gondola pointing system enables raster mapping of arbitrary geometry, with a repeatable positional accuracy of image30 super(image ); postflight pointing reconstruction to image5 super(image ) rms is achieved. The onboard telescope control software permits autonomous execution of a preselected 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 hr flight from ESRANGE (Kiruna), Sweden to Victoria Island, northern Canada in 2005 June; and a 250 hr, circumpolar flight from McMurdo Station, Antarctica, in 2006 December.
We develop a framework for joint constraints on the CO luminosity function based on power spectra (PS) and voxel intensity distributions (VID) and apply this to simulations of CO Mapping Array ...Pathfinder (COMAP), a CO intensity mapping experiment. This Bayesian framework is based on a Markov chain Monte Carlo (MCMC) sampler coupled to a Gaussian likelihood with a joint PS + VID covariance matrix computed from a large number of fiducial simulations and re-calibrated with a small number of simulations per MCMC step. The simulations are based on dark matter halos from fast peak patch simulations combined with the LCO(Mhalo) model of Li et al. We find that the relative power to constrain the CO luminosity function depends on the luminosity range of interest. In particular, the VID is more sensitive at large luminosities, while the PS and the VID are both competitive at small and intermediate luminosities. The joint analysis is superior to using either observable separately. When averaging over CO luminosities ranging between , and over 10 cosmological realizations of COMAP Phase 2, the uncertainties (in dex) are larger by 58% and 30% for the PS and VID, respectively, when compared to the joint analysis (PS + VID). This method is generally applicable to any other random field, with a complicated likelihood, as long a fast simulation procedure is available.
We present multi-wavelength detections of nine candidate gravitationally-lensed dusty star-forming galaxies (DSFGs) selected at 218GHz (1.4mm) from the ACT equatorial survey. Among the brightest ACT ...sources, these represent the subset of the total ACT sample lying in
SPIRE fields, and all nine of the 218GHz detections were found to have bright
counterparts. By fitting their spectral energy distributions (SEDs) with a modified blackbody model with power-law temperature distribution, we find the sample has a median redshift of
(68 per cent confidence interval), as expected for 218GHz selection, and an apparent total infrared luminosity of
, which suggests that they are either strongly lensed sources or unresolved collections of unlensed DSFGs. The effective apparent diameter of the sample is
, further evidence of strong lensing or multiplicity, since the typical diameter of dusty star-forming galaxies is 1.0-2.5 kpc. We emphasize that the effective apparent diameter derives from SED modelling without the assumption of optically thin dust (as opposed to image morphology). We find that the sources have substantial optical depth.
to dust around the peak in the modified blackbody spectrum (
⩽ 500
m), a result that is robust to model choice.
We present measurements of the auto- and cross-frequency power spectra of the cosmic infrared background (CIB) at 250, 350, and 500 mu m (1200, 860, and 600 GHz) from observations totaling ~70 deg ...sub(2) made with the SPIRE instrument aboard the Herschel Space Observatory. We measure a fractional anisotropy delta I/I = 14% + or - 4% , detecting signatures arising from the clustering of dusty star-forming galaxies in both the linear (2-halo) and nonlinear (1-halo) regimes; and that the transition from the 2- to 1-halo terms, below which power originates predominantly from multiple galaxies within dark matter halos, occurs at k sub(theta) ~ 0.10-0.12 aremin super(-1) (l ~ 2160-2380), from 250 to 500 mu m. We measure the cross-correlation power spectra between bands, finding that bands which are farthest apart are the least correlated, as well as hints of a reduction in the correlation between bands when resolved sources are more aggressively masked.
We present measurements of the power spectra of cosmic infrared background (CIB) and cosmic microwave background (CMB) fluctuations in six frequency bands. Maps at the lower three frequency bands, ...95, 150, and 220 GHz (3330, 2000, and 1360 m) are from the South Pole Telescope, while the upper three frequency bands, 600, 857, and 1200 GHz (500, 350, 250 m) are observed with Herschel/SPIRE. From these data, we produce 21 angular power spectra (6 auto- and 15 cross-frequency) spanning the multipole range 600 ≤ ≤ 11,000. Our measurements are the first to cross-correlate measurements near the peak of the CIB spectrum with maps at 95 GHz, complementing and extending the measurements from Planck Collaboration et al. at 143-857 GHz. The observed fluctuations originate largely from clustered, infrared-emitting, dusty star-forming galaxies, the CMB, and to a lesser extent radio galaxies, active galactic nuclei, and the Sunyaev-Zel'dovich effect.
The Herschel Extragalactic Legacy Project (HELP) brings together a vast range of data from many astronomical observatories. Its main focus is on the Herschel data, which maps dust-obscured star ...formation over 1300 deg2. With this unprecedented combination of data sets, it is possible to investigate how the star formation versus stellar mass relation (main sequence) of star-forming galaxies depends on environment. In this pilot study, we explore this question within 0.1 < z < 3.2 using data in the COSMOS field. We estimate the local environment from a smoothed galaxy density field using the full photometric redshift probability distribution. We estimate star formation rates by stacking the SPIRE data from the Herschel Multi-tiered Extragalactic Survey. Our analysis rules out the hypothesis that the main sequence for star-forming systems is independent of environment at 1.5 < z < 2, while a simple model in which the mean specific star formation rate declines with increasing environmental density gives a better description. However, we cannot exclude a simple hypothesis in which the main sequence for star-forming systems is independent of environment at z < 1.5 and z > 2. We also estimate the evolution of the star formation rate density in the COSMOS field, and our results are consistent with previous measurements at z < 1.5 and z > 2 but we find a
$1.4^{+0.3}_{-0.2}$
times higher peak value of the star formation rate density at z ∼ 1.9.
We present measurements of the mean mid-infrared to submillimetre flux densities of massive (M
★≳ 1011 M⊙) galaxies at redshifts 1.7 < z < 2.9, obtained by stacking positions of known objects taken ...from the GOODS NICMOS Survey (GNS) catalogue on maps at 24
m (Spitzer/MIPS); 70, 100 and 160
m (Herschel/PACS); 250, 350 and 500
m (BLAST); and 870
m (LABOCA). A modified blackbody spectrum fit to the stacked flux densities indicates a median interquartile star formation rate (SFR) of SFR = 6348, 81 M⊙ 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, disc-like and spheroid-like, according to their Sérsic indices, n. We find evidence that most of the star formation is occurring in n≤ 2 (disc-like) galaxies, with median interquartile SFR = 122100, 150 M⊙ yr−1, while there are indications that the n > 2 (spheroid-like) population may be forming stars at a median interquartile SFR = 149, 20 M⊙ 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.
We use new large area far infrared maps ranging from 65 to 500 mum obtained with the AKARI and the Balloon-borne Large Aperture Submillimeter Telescope 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 ~35 K and with an estimated mass of 0.06 M sun. 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 mum 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.