We present a new technique to measure multi-wavelength "super-deblended" photometry from highly confused images, which we apply to Herschel and ground-based far-infrared (FIR) and (sub-)millimeter ...(mm) data in the northern field of the Great Observatories Origins Deep Survey. There are two key novelties. First, starting with a large database of deep Spitzer 24 m and VLA 20 cm detections that are used to define prior positions for fitting the FIR/submm data, we perform an active selection of useful priors independently at each frequency band, moving from less to more confused bands. Exploiting knowledge of redshift and all available photometry, we identify hopelessly faint priors that we remove from the fitting pool. This approach significantly reduces blending degeneracies and allows reliable photometry to be obtained for galaxies in FIR+mm bands. Second, we obtain well-behaved, nearly Gaussian flux density uncertainties, individually tailored to all fitted priors for each band. This is done by exploiting extensive simulations that allow us to calibrate the conversion of formal fitting uncertainties to realistic uncertainties, depending on directly measurable quantities. We achieve deeper detection limits with high fidelity measurements and uncertainties at FIR+mm bands. As an illustration of the utility of these measurements, we identify 70 galaxies with and reliable FIR+mm detections. We present new constraints on the cosmic star formation rate density at , finding a significant contribution from dusty galaxies that are missed by optical-to-near-infrared color selection. Photometric measurements for 3306 priors, including more than 1000 FIR+mm detections, are released publicly with our catalog.
Follow-up observations at high-angular resolution of bright submillimeter galaxies selected from deep extragalactic surveys have shown that the single-dish sources are comprised of a blend of several ...galaxies. Consequently, number counts derived from low- and high-angular-resolution observations are in tension. This demonstrates the importance of resolution effects at these wavelengths and the need for realistic simulations to explore them. We built a new 2 deg2 simulation of the extragalactic sky from the far-infrared to the submillimeter. It is based on an updated version of the 2SFM (two star-formation modes) galaxy evolution model. Using global galaxy properties generated by this model, we used an abundance-matching technique to populate a dark-matter lightcone and thus simulate the clustering. We produced maps from this simulation and extracted the sources, and we show that the limited angular resolution of single-dish instruments has a strong impact on (sub)millimeter continuum observations. Taking into account these resolution effects, we are reproducing a large set of observables, as number counts and their evolution with redshift and cosmic infrared background power spectra. Our simulation consistently describes the number counts from single-dish telescopes and interferometers. In particular, at 350 and 500 μm, we find that the number counts measured by Herschel between 5 and 50 mJy are biased towards high values by a factor ~2, and that the redshift distributions are biased towards low redshifts. We also show that the clustering has an important impact on the Herschel pixel histogram used to derive number counts from P(D) analysis. We find that the brightest galaxy in the beam of a 500 μm Herschel source contributes on average to only ~60% of the Herschel flux density, but that this number will rise to ~95% for future millimeter surveys on 30 m-class telescopes (e.g., NIKA2 at IRAM). Finally, we show that the large number density of red Herschel sources found in observations but not in models might be an observational artifact caused by the combination of noise, resolution effects, and the steepness of color- and flux density distributions. Our simulation, called Simulated Infrared Dusty Extragalactic Sky (SIDES), is publicly available.
We interpret the wide variety of redshift distributions for galaxies found by far-infrared and (sub) millimeter deep surveys depending on their depth and wavelength using our phenomenological model ...of galaxy evolution. This model reproduces the observed redshift distributions from 100 μm to 1.4 mm and, especially, the increase in the median redshift with survey wavelength without any new parameter tuning. This median redshift also varies significantly with the depth of the surveys, and deeper surveys do not necessarily probe higher redshifts. Paradoxically, at a fixed wavelength and flux limit, the lensed sources are not always at higher redshift. We found that the higher redshift of 1.4 mm-selected South Pole Telescope (SPT) sources compared to other SMG surveys is caused not only by the lensing selection, but also by the longer wavelength. This SPT sample is expected to be dominated by a population of lensed main-sequence galaxies and a minor contribution (~10%) of unlensed extreme starbursts.
Abstract
We present the characteristics of 2 mm selected sources from the largest Atacama Large Millimeter/submillimeter Array (ALMA) blank-field contiguous survey conducted to date, the Mapping ...Obscuration to Reionization with ALMA (MORA) survey covering 184 arcmin
2
at 2 mm. Twelve of 13 detections above 5
σ
are attributed to emission from galaxies, 11 of which are dominated by cold dust emission. These sources have a median redshift of
〈
z
2
mm
〉
=
3.6
−
0.3
+
0.4
primarily based on optical/near-infrared photometric redshifts with some spectroscopic redshifts, with 77% ± 11% of sources at
z
> 3 and 38% ± 12% of sources at
z
> 4. This implies that 2 mm selection is an efficient method for identifying the highest-redshift dusty star-forming galaxies (DSFGs). Lower-redshift DSFGs (
z
< 3) are far more numerous than those at
z
> 3 yet are likely to drop out at 2 mm. MORA shows that DSFGs with star formation rates in excess of 300
M
⊙
yr
−1
and a relative rarity of ∼10
−5
Mpc
−3
contribute ∼30% to the integrated star formation rate density at 3 <
z
< 6. The volume density of 2 mm selected DSFGs is consistent with predictions from some cosmological simulations and is similar to the volume density of their hypothesized descendants: massive, quiescent galaxies at
z
> 2. Analysis of MORA sources’ spectral energy distributions hint at steeper empirically measured dust emissivity indices than reported in typical literature studies, with
〈
β
〉
=
2.2
−
0.4
+
0.5
. The MORA survey represents an important step in taking census of obscured star formation in the universe’s first few billion years, but larger area 2 mm surveys are needed to more fully characterize this rare population and push to the detection of the universe’s first dusty galaxies.
We reproduce the mid-infrared to radio galaxy counts with a new empirical model based on our current understanding of the evolution of main-sequence (MS) and starburst (SB) galaxies. We rely on a ...simple spectral energy distribution (SED) library based on Herschel observations: a single SED for the MS and another one for SB, getting warmer with redshift. Our model is able to reproduce recent measurements of galaxy counts performed with Herschel, including counts per redshift slice. This agreement demonstrates the power of our 2-Star-Formation Modes (2SFM) decomposition in describing the statistical properties of infrared sources and their evolution with cosmic time. We discuss the relative contribution of MS and SB galaxies to the number counts at various wavelengths and flux densities. We also show that MS galaxies are responsible for a bump in the 1.4 GHz radio counts around 50 mu Jy. Material of the model (predictions, SED library, mock catalogs, etc.) is available online.
We present a physical characterization of MM J100026.36+021527.9 (a.k.a. "Mambo-9"), a dusty star-forming galaxy (DSFG) at z = 5.850 0.001. This is the highest-redshift unlensed DSFG (and fourth most ...distant overall) found to date and is the first source identified in a new 2 mm blank-field map in the COSMOS field. Though identified in prior samples of DSFGs at 850 m to 1.2 mm with unknown redshift, the detection at 2 mm prompted further follow-up as it indicated a much higher probability that the source was likely to sit at z > 4. Deep observations from the Atacama Large Millimeter and submillimeter Array (ALMA) presented here confirm the redshift through the secure detection of 12CO(J = 6→5) and p-H2O (21,1 → 20,2). Mambo-9 is composed of a pair of galaxies separated by 6 kpc with corresponding star formation rates of 590 M yr−1 and 220 M yr−1, total molecular hydrogen gas mass of (1.7 0.4) × 1011M , dust mass of (1.3 0.3) × 109M , and stellar mass of ( ) × 109M . The total halo mass, (3.3 0.8) × 1012M , is predicted to exceed 1015M by z = 0. The system is undergoing a merger-driven starburst that will increase the stellar mass of the system tenfold in τdepl = 40−80 Myr, converting its large molecular gas reservoir (gas fraction of ) into stars. Mambo-9 evaded firm spectroscopic identification for a decade, following a pattern that has emerged for some of the highest-redshift DSFGs found. And yet, the systematic identification of unlensed DSFGs like Mambo-9 is key to measuring the global contribution of obscured star formation to the star formation rate density at z 4, the formation of the first massive galaxies, and the formation of interstellar dust at early times ( 1 Gyr).
Abstract
The 2 mm Mapping Obscuration to Reionization with ALMA (MORA) Survey was designed to detect high-redshift (
z
≳ 4), massive, dusty star-forming galaxies (DSFGs). Here we present two likely ...high-redshift sources, identified in the survey, whose physical characteristics are consistent with a class of optical/near-infrared (OIR)-invisible DSFGs found elsewhere in the literature. We first perform a rigorous analysis of all available photometric data to fit spectral energy distributions and estimate redshifts before deriving physical properties based on our findings. Our results suggest the two galaxies, called MORA-5 and MORA-9, represent two extremes of the “OIR-dark” class of DSFGs. MORA-5 (
z
phot
=
4.3
−
1.3
+
1.5
) is a significantly more active starburst with a star formation rate (SFR) of
830
−
190
+
340
M
⊙
yr
−1
compared to MORA-9 (
z
phot
=
4.3
−
1.0
+
1.3
), whose SFR is a modest
200
−
60
+
250
M
⊙
yr
−1
. Based on the stellar masses (
M
⋆
≈ 10
10−11
M
⊙
), space density (
n
∼ (5 ± 2) × 10
−6
Mpc
−3
, which incorporates two other spectroscopically confirmed OIR-dark DSFGs in the MORA sample at
z
= 4.6 and
z
= 5.9), and gas depletion timescales (<1 Gyr) of these sources, we find evidence supporting the theory that OIR-dark DSFGs are the progenitors of recently discovered 3 <
z
< 4 massive quiescent galaxies.
Recent studies have revealed a strong correlation between the star formation rate (SFR) and stellar mass of the majority of star-forming galaxies, the so-called star-forming main sequence. An ...empirical modeling approach (the 2-SFM framework) that distinguishes between the main sequence and rarer starburst galaxies is capable of reproducing most statistical properties of infrared galaxies, such as number counts, luminosity functions, and redshift distributions. In this paper, we extend this approach by establishing a connection between stellar mass and halo mass with the technique of abundance matching. Based on a few simple assumptions and a physically motivated formalism, our model successfully predicts the (cross-)power spectra of the cosmic infrared background (CIB), the cross-correlation between CIB and cosmic microwave background (CMB) lensing, and the correlation functions of bright, resolved infrared galaxies measured by Herschel, Planck, ACT, and SPT. We use this model to infer the redshift distribution of CIB-anisotropies and of the CIB × CMB lensing signal, as well as the level of correlation between CIB-anisotropies at different wavelengths. We study the link between dark matter halos and star-forming galaxies in the framework of our model. We predict that more than 90% of cosmic star formation activity occurs in halos with masses between 1011.5 and 1013.5 M⊙. If taking subsequent mass growth of halos into account, this implies that the majority of stars were initially (at z > 3) formed in the progenitors of clusters (Mh(z = 0) > 1013.5 M⊙), then in groups (1012.5 < Mh(z = 0) < 1013.5 M⊙) at 0.5 < z < 3, and finally in Milky-Way-like halos (1011.5 < Mh(z = 0) < 1012.5 M⊙) at z < 0.5. At all redshifts, the dominant contribution to the SFR density stems from halos of mass ~1012 M⊙, in which the instantaneous star formation efficiency – defined here as the ratio between SFR and baryonic accretion rate – is maximal (~70%). The strong redshift-evolution of SFR in the galaxies that dominate the CIB is thus plausibly driven by increased accretion from the cosmic web onto halos of this characteristic mass scale. Material (effective spectral energy distributions, differential emissivities of halos, relations between Mh and SFR) associated to this model is available at http://irfu.cea.fr/Sap/Phocea/Page/index.php?id=537.
ABSTRACT
We present an extensive ALMA spectroscopic follow-up programme of the $z\, {=}\, 4.3$ structure SPT2349–56, one of the most actively star-forming protocluster cores known, to identify ...additional members using their C ii 158 μm and CO(4–3) lines. In addition to robustly detecting the 14 previously published galaxies in this structure, we identify a further 15 associated galaxies at $z\, {=}\, 4.3$, resolving 55$\, {\pm }\,$5 per cent of the 870 μm flux density at 0.5 arcsec resolution compared to 21 arcsec single-dish data. These galaxies are distributed into a central core containing 23 galaxies extending out to 300 kpc in diameter, and a northern extension, offset from the core by 400 kpc, containing three galaxies. We discovered three additional galaxies in a red Herschel-SPIRE source 1.5 Mpc from the main structure, suggesting the existence of many other sources at the same redshift as SPT2349–56 that are not yet detected in the limited coverage of our data. An analysis of the velocity distribution of the central galaxies indicates that this region may be virialized with a mass of (9$\pm 5)\, {\times }\, 10^{12}$ M⊙, while the two offset galaxy groups are about 30 and 60 per cent less massive and show significant velocity offsets from the central group. We calculate the C ii and far-infrared number counts, and find evidence for a break in the C ii luminosity function. We estimate the average SFR density within the region of SPT2349–56 containing single-dish emission (a proper diameter of 720 kpc), assuming spherical symmetry, to be roughly 4$\, {\times }\, 10^4$ M⊙ yr−1 Mpc−3; this may be an order of magnitude greater than the most extreme examples seen in simulations.
Massive galaxy-scale outflows of gas are one of the most commonly invoked mechanisms to regulate the growth and evolution of galaxies throughout the universe. While the gas in outflows spans a large ...range of temperatures and densities, the cold molecular phase is of particular interest because molecular outflows may be capable of suppressing star formation in galaxies by removing the star-forming gas. We have conducted the first survey of molecular outflows at z > 4, targeting 11 strongly lensed dusty, star-forming galaxies (DSFGs) with high-resolution Atacama Large Millimeter/submillimeter Array observations of OH 119 m absorption as an outflow tracer. In this first paper, we give an overview of the survey, focusing on the detection rate and structure of molecular outflows. We find unambiguous evidence for outflows in 8/11 galaxies (73%), more than tripling the number known at z > 4. This implies that molecular winds in z > 4 DSFGs must have both a near-unity occurrence rate and large opening angles to be detectable in absorption. Lensing reconstructions reveal that 500 pc scale clumpy structures in the outflows are common. The individual clumps are not directly resolved, but from optical depth arguments we expect that future observations will require 50-200 pc spatial resolution to resolve them. We do not detect high-velocity C ii wings in any of the sources with clear OH outflows, indicating that C ii is not a reliable tracer of molecular outflows. Our results represent a first step toward characterizing molecular outflows at z > 4 at the population level, demonstrating that large-scale outflows are ubiquitous among early massive, dusty galaxies.