High-redshift Dusty Galaxies Research Articles

Simba-EoR: early galaxy formation in the simba simulation including a new sub-grid interstellar medium model

ABSTRACT We update the dust model present within the simba galaxy simulations with a self-consistent framework for the co-evolution of dust and molecular hydrogen populations in the interstellar medium, and use this to explore $z \ge 6$ galaxy evolution. In addition to tracking the evolution of dust and molecular hydrogen abundances, our model fully integrates these species into the simba simulation, explicitly modelling their impact on physical processes such as star formation and cooling through the inclusion of a novel two-phase sub-grid model for interstellar gas. Running two cosmological simulations down to $z \sim 6$ we find that our simba-EoR model displays a generally tighter concordance with observational data than fiducial simba. Additionally we observe that our simba-EoR models increase star formation activity at early epochs, producing larger dust-to-gas ratios consequently. Finally, we discover a significant population of hot dust at $\sim 100$ K, aligning with contemporaneous observations of high-redshift dusty galaxies, alongside the large $\sim 20$ K population typically identified.

Metal factories in the early Universe

ABSTRACT We have estimated the mass of metals in the molecular gas in 13 dusty star-forming galaxies at $z \sim 4$ in which the gas, based on previous observations, lies in a cold rotating disc. We estimated the metal masses using either the submillimetre line or continuum emission from three tracers of the overall metal content – carbon atoms, carbon monoxide molecules, and dust grains – using the first simultaneous calibration of all three tracers. We obtain very similar mass estimates from the different tracers, which are similar to the entire metal content of a present-day massive early-type galaxy. We used the dynamical masses of these galaxies to estimate an upper limit on the mass of the molecular gas in each galaxy, allowing us to estimate a lower limit on the metal abundance of the gas, finding values for many of the galaxies well above the solar value. We show that the high metal masses and metal abundances are what is expected shortly after the formation of a galaxy for a top-heavy IMF. We suggest a scenario for galaxy evolution in which massive galaxies reach a high metal abundance during their formation phase, which is then gradually reduced by dry mergers with lower mass galaxies. We show that the metals in the outflows from high-redshift dusty star-forming galaxies can quantitatively explain the long-standing puzzle that such a large fraction of the metals in galaxy clusters ($\simeq$0.75) is in the intracluster gas rather than in the galaxies themselves.

Faint mm NIKA2 dusty star-forming galaxies: Finding the high-redshift population

Aims. High-redshift dusty star-forming galaxies (DSFGs) are proposed to be the progenitors of massive quiescent galaxies arising at cosmic noon, providing a crucial insight into the formation, assembly, and early quenching of massive galaxies in the early Universe. However, their high redshift combined with high dust obscuration adds significant difficulties to their redshift measurement, which is mandatory for detailed studies of their physical properties. Blind mm spectral scans are the most unbiased way in prinicple for obtaining accurate spectroscopic redshifts for these sources, but identifying faint molecular and atomic lines within limited telescope time for faint DSFGs is also difficult with these scans. Methods. We developed a new framework to constrain the source redshift. The method jointly accounts for the detection and/or nondetection of spectral lines and the prior information from the photometric redshift and total infrared luminosity from spectral energy distribution analysis. The method uses the estimated total infrared luminosity to predict the line fluxes at given redshifts and generates model spectra. The redshift-dependent spectral models were then compared with the observed spectra to determine the redshift. Results. We applied this joint redshift analysis method to four high-z dusty star-forming galaxy candidates selected from the NIKA2 observations of the HLSJ091828.6+514223 (HLS) field that were further observed by NOEMA with blind spectral scans. These sources only have Herschel SPIRE photometry as ancillary data. They were selected because SPIRE counterparts are faint or entirely lacking and thus favor to select the highest-redshift candidates. The method finds a spectroscopic redshift of 4 in the five NOEMA-counterpart detected sources, with z > 3. Based on these measurements, we derived the CO and [CI] lines and mm continuum fluxes from the NOEMA data and studied the properties of their interstellar medium and star formation. We find cold dust temperatures in some of the HLS sources compared to the general population of submm galaxies, which might be related to the bias introduced by the SPIRE-dropout selection. All sources except for one have a short gas-depletion time of a few hundred million years, which is typical of high-z submm galaxies. The only exception shows a longer gas-depletion time of up to a few billion years. This is comparable to the gas-depletion times of main-sequence galaxies at the same redshift. Furthermore, we identify a possible overdensity of dusty star-forming galaxies at z = 5.2 that is traced by two sources in our sample, as well as a lensed galaxy HLSJ091828.6+514223. Conclusions. We demonstrate that our method when applied to mm-selected DSFGs is able to determine the redshift accurately. This accuracy with only multiple emission lines with a low signal-to-noise ratio shows promising potential for the blind redshift search in large samples of high-z DSFGs, even in the absence of optical to near infrared photometric redshifts.

Metal and dust evolution in ALMA REBELS galaxies: insights for future JWST observations

ABSTRACT Atacama Large Millimeter/submillimeter Array (ALMA) observations revealed the presence of significant amounts of dust in the first Gyr of Cosmic time. However, the metal and dust build-up picture remains very uncertain due to the lack of constraints on metallicity. JWST has started to reveal the metal content of high-redshift targets, which may lead to firmer constraints on high-redshift dusty galaxies evolution. In this work, we use detailed chemical and dust evolution models to explore the evolution of galaxies within the ALMA Reionization Era Bright Emission Line Survey (REBELS) survey, testing different metallicity scenarios that could be inferred from JWST observations. In the models, we track the build-up of stellar mass using non-parametric star formation histories for REBELS galaxies. Different scenarios for metal and dust evolution are simulated by allowing different prescriptions for gas flows and dust processes. The model outputs are compared with measured dust scaling relations, by employing metallicity-dependent calibrations for the gas mass based on the [C ii] 158 μm line. Independently of the galaxies metal content, we found no need for extreme dust prescriptions to explain the dust masses revealed by ALMA. However, different levels of metal enrichment will lead to different dominant dust production mechanisms, with stardust production dominant over other interstellar medium dust processes only in the metal-poor case. This points out how metallicity measurements from JWST will significantly improve our understanding of the dust build-up in high-redshift galaxies. We also show that models struggle to reproduce observables such as dust-to-gas and dust-to-stellar ratios simultaneously, possibly indicating an overestimation of the gas mass through current calibrations, especially at high metallicities.

SUNRISE: The rich molecular inventory of high-redshift dusty galaxies revealed by broadband spectral line surveys

Understanding the nature of high-redshift dusty galaxies requires a comprehensive view of their interstellar medium (ISM) and molecular complexity. However, the molecular ISM at high redshifts is commonly studied using only a few species beyond 12C16O, limiting our understanding. In this paper, we present the results of deep 3 mm spectral line surveys using the NOrthern Extended Millimeter Array (NOEMA) targeting two strongly lensed dusty galaxies observed when the Universe was less than 1.8 Gyr old: APM 08279+5255, a quasar at redshift z = 3.911, and NCv1.143 (H-ATLAS J125632.7+233625), a z = 3.565 starburst galaxy. The spectral line surveys cover rest-frame frequencies from about 330 to 550 GHz for both galaxies. We report the detection of 38 and 25 emission lines in APM 08279+5255 and NCv1.143, respectively. These lines originate from 17 species, namely CO, 13CO, C18O, CN, CCH, HCN, HCO+, HNC, CS, C34S, H2O, H3O+, NO, N2H+, CH, c-C3H2, and the vibrationally excited HCN and neutral carbon. The spectra reveal the chemical richness and the complexity of the physical properties of the ISM. By comparing the spectra of the two sources and combining the analysis of the molecular gas excitation, we find that the physical properties and the chemical imprints of the ISM are different: the molecular gas is more excited in APM 08279+5255, which exhibits higher molecular gas temperatures and densities compared to NCv1.143; the molecular abundances in APM 08279+5255 are akin to the values of local active galactic nuclei (AGN), showing boosted relative abundances of the dense gas tracers that might be related to high-temperature chemistry and/or the X-ray-dominated regions, while NCv1.143 more closely resembles local starburst galaxies. The most significant differences between the two sources are found in H2O: the 448 GHz ortho-H2O(423 − 330) line is significantly brighter in APM 08279+5255, which is likely linked to the intense far-infrared radiation from the dust powered by AGN. Our astrochemical model suggests that, at such high column densities, far-ultraviolet radiation is less important in regulating the ISM, while cosmic rays (and/or X-rays and shocks) are the key players in shaping the molecular abundances and the initial conditions of star formation. Both our observed CO isotopologs line ratios and the derived extreme ISM conditions (high gas temperatures, densities, and cosmic-ray ionization rates) suggest the presence of a top-heavy stellar initial mass function. From the ∼330–550 GHz continuum, we also find evidence of nonthermal millimeter flux excess in APM 08279+5255 that might be related to the central supermassive black hole. Such deep spectral line surveys open a new window into the physics and chemistry of the ISM and the radiation field of galaxies in the early Universe.

PRUSSIC

We present deep ALMA Band 3 observations of the HCN, HCO+, and HNC(4–3) emission in SDP.81, a well-studied z = 3.042, strongly lensed galaxy. These lines trace the high-density gas, which remains almost entirely unexplored in z ≥ 1 galaxies. Additionally, these dense-gas tracers are potentially powerful diagnostics of the mechanical heating of the interstellar medium. While the HCN(4–3) and HNC(4–3) lines are not detected, the HCO+(4–3) emission is clearly detected and resolved. This is the third detection of this line in a high-redshift star-forming galaxy. We find an unusually high HCO+/HCN intensity ratio of ≥2.2. Based on the modelling of the photodissociation region, the most likely explanation for the elevated HCO+/HCN ratio is that SDP.81 has low mechanical heating, making up less than 10% of the total energy budget, along with a sub-solar metallicity of Z ≈ 0.5 Z⊙. While such conditions might not be representative of the general population of high-redshift dusty galaxies, a lower-than-solar metallicity might significantly impact gas masses inferred from CO observations. In addition, we report the detection of CO(0–1) absorption from the foreground lensing galaxy and CO(1–0) emission from a massive companion to the lensing galaxy, approximately 50 kpc to the south-east.

Observations of neutral carbon in 29 high-z lensed dusty star-forming galaxies and the comparison of gas mass tracers

Context. The nature and evolution of high-redshift dusty star-forming galaxies (high-z DSFGs) remain open questions. Their massive gas reservoirs play an important role in driving the intense star formation rates hosted in these galaxies. Aims. Our aim is to estimate the molecular gas content of high-z DSFGs by using various gas mass tracers, such as the [CI], CO, [CII] emission lines and the dust content. These tracers need to be well calibrated as they are all limited by uncertainties on factors such as αCO, XCI, α[CII], and δGDR, thereby affecting the accurate determination of the gas mass. The main goal of our work is to check the consistency between the gas mass tracers and to cross-calibrate the uncertain factors. Methods. We observed the two [CI] line transitions for 29 South Pole telescope Submillimeter Galaxies (SPT-SMGs) with the Atacama Large Millimeter/submillimeter Array–Atacama Compact Array (ALMA-ACA). Additionally, we also present new Atacama Pathfinder EXperiment (APEX) observations of [CII] line for nine of these galaxies. We combined our observations with the rich ancillary data of low- and mid-J CO lines, ancillary [CII] line, and dust mass estimations for these galaxies. Results. We find a nearly linear relation between the infrared luminosity and [CI] luminosity if we fit the starbursts and main-sequence galaxies separately. We measure a median [CI]-derived excitation temperature of 34.5 ± 2.1 K. We probed the properties of the interstellar medium (ISM), such as density and radiation field intensity, using [CI] to mid- or high-J CO lines and the [CI]-to-infrared luminosity ratio, and find similar values to the SMG populations in the literature. Finally, the gas masses estimated from [CI], CO, dust, and [CII] do not exhibit any significant trend with the infrared luminosity or the dust temperature. We provide the various cross-calibrations between these tracers. Conclusions. Our study confirms that [CI] is a suitable tracer of the molecular gas content, and shows an overall agreement between all the classical gas tracers used at high redshift. However, their absolute calibration, and thus the gas depletion timescale measurements, remain uncertain.

The nature of 500 micron risers – II. Multiplicities and environments of sub-mm faint dusty star-forming galaxies

ABSTRACT We present Submillimeter Common User Bolometer Array 2 (SCUBA-2) and Submillimeter Array (SMA) follow-up observations of four candidate high-redshift dusty star-forming galaxies, selected as sources with rising spectral energy distributions in the 250, 350, and 500 μm Herschel SPIRE bands. Previous SMA observations showed no counterparts to these sources, but in our deeper sub-mm observations we detect counterparts to all four 500 μm risers, with three resolving into multiple systems. For these three multiple systems, the SMA 345 GHz (≈870 μm) observations recover $123 \pm 73{{\ \rm per\ cent}}$, $60 \pm 15{{\ \rm per\ cent}}$, and $19 \pm 4{{\ \rm per\ cent}}$ of the integrated 850 μm flux density from SCUBA-2, indicating that there may be additional sources below our SMA detection limit making up a dense, protocluster core. The fourth 500 μm riser was observed at a lower frequency and so we cannot make a similar comparison. We estimate photometric redshifts based on FIR/sub-mm colours, finding that three fourths likely lie at z ≥ 2. This fits with the interpretation that the 500 μm riser selection criterion selects both intrinsically red, individual galaxies at z > 4, and multiple systems at more moderate redshifts, artificially reddened by the effects of blending. We use the SCUBA-2 850 μm maps to investigate the environments of these 500 μm risers. By constructing cumulative number counts and estimating photometric redshifts for surrounding SCUBA-2 detections, we find that one of our 500 μm risers could plausibly reside in a z ≥ 2 protocluster. We infer that bright 500 μm risers with faint 850 μm flux densities are typically multiple systems at z ≥ 2 that may reside in overdensities of bright sub-mm galaxies.

High resolution spectral imaging of CO(7–6), [CI](2–1), and continuum of three high-z lensed dusty star-forming galaxies using ALMA

High-redshift dusty star-forming galaxies with very high star formation rates (500−3000 M⊙ yr−1) are key to understanding the formation of the most extreme galaxies in the early Universe. Characterising the gas reservoir of these systems can reveal the driving factor behind the high star formation. Using molecular gas tracers such as, high-J CO lines, neutral carbon lines, and the dust continuum, we can estimate the gas density and radiation field intensity in their interstellar media. In this paper, we present high resolution (∼0.4″) observations of CO(7−6), [CI](2−1), and dust continuum of three lensed galaxies from the South pole telescope – sub-millimetre galaxies (SPT-SMG) sample at z ∼ 3 with the Atacama Large Millimetre/submillimetre Array. Our sources have high intrinsic star formation rates (> 850 M⊙ yr−1) and rather short depletion timescales (< 100 Myr). Based on the L[CI](2−1)/LCO(7 − 6) and L[CI](2−1)/LIR ratios, our galaxy sample has similar radiation field intensities and gas densities compared to other submillimetre galaxies. We performed visibility-based lens modelling on these objects to reconstruct the kinematics in the source plane. We find that the cold gas masses of the sources are compatible with simple dynamical mass estimates using ULIRG-like values of the CO-H2 conversion factor αCO, but not Milky Way-like values. We find diverse source kinematics in our sample: SPT0103−45 and SPT2147−50 are likely rotating disks, while SPT2357−51 is possibly a major merger. The analysis presented in the paper could be extended to a larger sample to determine better statistics of morphologies and interstellar medium properties of high-z dusty star-forming galaxies.

Measurements of the Dust Properties in z ≃ 1–3 Submillimeter Galaxies with ALMA

Abstract We present Atacama Large Millimeter/submillimeter Array (ALMA) 2 mm continuum observations of a complete and unbiased sample of 99 870 μm selected submillimeter galaxies (SMGs) in the Extended Chandra Deep Field South (ALESS). Our observations of each SMG reach average sensitivities of 53 μJy beam−1. We measure the flux densities for 70 sources, for which we obtain a typical 870 μm-to-2 mm flux ratio of 14 ± 5. We do not find a redshift dependence of this flux ratio, which would be expected if the dust emission properties of our SMGs were the same at all redshifts. By combining our ALMA measurements with existing Herschel/SPIRE observations, we construct a (biased) subset of 27 galaxies for which the cool dust emission is sufficiently well sampled to obtain precise constraints on their dust properties using simple isothermal models. Thanks to our new 2 mm observations, the dust emissivity index is well constrained and robust against different dust opacity assumptions. The median dust emissivity index of our SMGs is β ≃ 1.9 ± 0.4, consistent with the emissivity index of dust in the Milky Way and other local and high-redshift galaxies, as well as classical dust-grain model predictions. We also find a negative correlation between the dust temperature and β, similar to low-redshift observational and theoretical studies. Our results indicate that β ≃ 2 in high-redshift dusty star-forming galaxies, implying little evolution in dust-grain properties between our SMGs and local dusty galaxy samples, and suggesting that these high-mass and high-metallicity galaxies have dust reservoirs driven by grain growth in their interstellar medium.

An ALMA Survey of the SCUBA-2 Cosmology Legacy Survey UKIDSS/UDS Field: The Far-infrared/Radio Correlation for High-redshift Dusty Star-forming Galaxies

Abstract We study the radio properties of 706 submillimeter galaxies (SMGs) selected at 870 μm with the Atacama Large Millimeter Array from the SCUBA-2 Cosmology Legacy Survey map of the Ultra Deep Survey field. We detect 273 SMGs at >4σ in deep Karl G. Jansky Very Large Array 1.4 GHz observations, of which a subset of 45 SMGs are additionally detected in 610 MHz Giant Metre-Wave Radio Telescope imaging. We quantify the far-infrared/radio correlation (FIRRC) through parameter q IR, defined as the logarithmic ratio of the far-infrared and radio luminosity, and include the radio-undetected SMGs through a stacking analysis. We determine a median q IR = 2.20 ± 0.03 for the full sample, independent of redshift, which places these z ∼ 2.5 dusty star-forming galaxies 0.44 ± 0.04 dex below the local correlation for both normal star-forming galaxies and local ultra-luminous infrared galaxies (ULIRGs). Both the lack of redshift evolution and the offset from the local correlation are likely the result of the different physical conditions in high-redshift starburst galaxies, compared to local star-forming sources. We explain the offset through a combination of strong magnetic fields (B ≳ 0.2 mG), high interstellar medium (ISM) densities and additional radio emission generated by secondary cosmic rays. While local ULIRGs are likely to have similar magnetic field strengths, we find that their compactness, in combination with a higher ISM density compared to SMGs, naturally explains why local and high-redshift dusty star-forming galaxies follow a different FIRRC. Overall, our findings paint SMGs as a homogeneous population of galaxies, as illustrated by their tight and nonevolving far-infrared/radio correlation.

The assembly of dusty galaxies at z ≥ 4: statistical properties

ABSTRACT The recent discovery of high-redshift dusty galaxies implies a rapid dust enrichment of their interstellar medium (ISM). To interpret these observations, we run a cosmological simulation in a 30 h−1 cMpc/size volume down to z ≈ 4. We use the hydrodynamical code dustygadget, which accounts for the production of dust by stellar populations and its evolution in the ISM. We find that the cosmic dust density parameter (Ωd) is mainly driven by stellar dust at z ≳ 10, so that mass- and metallicity-dependent yields are required to assess the dust content in the first galaxies. At z ≲ 9, the growth of grains in the ISM of evolved systems [log(M⋆/M⊙) > 8.5] significantly increases their dust mass, in agreement with observations in the redshift range 4 ≲ z < 8. Our simulation shows that the variety of high-redshift galaxies observed with the Atacama Large Millimeter Array can naturally be accounted for by modelling the grain growth time-scale as a function of the physical conditions in the gas cold phase. In addition, the trends of dust-to-metal and dust-to-gas (${\cal D}$) ratios are compatible with the available data. A qualitative investigation of the inhomogeneous dust distribution in a representative massive halo at z ≈ 4 shows that dust is found from the central galaxy up to the closest satellites along polluted filaments with $\rm log({\cal D}) \le -2.4$, but sharply declines at distances d ≳ 30 kpc along many lines of sight, where $\rm log({\cal D}) \lesssim -4.0$.

Dense-gas tracers and carbon isotopes in five 2.5 < z < 4 lensed dusty star-forming galaxies from the SPT SMG sample

The origin of the high star formation rates (SFR) observed in high-redshift dusty star-forming galaxies is still unknown. Large fractions of dense molecular gas might provide part of the explanation, but there are few observational constraints on the amount of dense gas in high-redshift systems dominated by star formation. In this paper, we present the results of our Atacama large millimeter array (ALMA) program targeting dense-gas tracers (HCN(5-4), HCO+(5-4), and HNC(5-4)) in five strongly lensed galaxies from the South Pole Telescope (SPT) submillimeter galaxy sample. We detected two of these lines (S/N > 5) in SPT-125-47 at z = 2.51 and tentatively detected all three (S/N ∼ 3) in SPT0551-50 at z = 3.16. Since a significant fraction of our target lines is not detected, we developed a statistical method to derive unbiased mean properties of our sample taking into account both detections and non-detections. On average, the HCN(5-4) and HCO+(5-4) luminosities of our sources are a factor of ∼1.7 fainter than expected, based on the local L′HCN(5-4) − LIR relation, but this offset corresponds to only ∼2σ if we consider sample variance. We find that both the HCO+/HCN and HNC/HCN flux ratios are compatible with unity. The first ratio is expected for photo-dominated regions (PDRs) while the second is consistent with PDRs or X-ray dominated regions (XDRs) and/or mid-infrared (IR) pumping of HNC. Our sources are at the high end of the local relation between the star formation efficiency, determined using the LIR/[CI] and LIR/CO ratios, and the dense-gas fraction, estimated using the HCN/[CI] and HCN/CO ratios. Finally, in SPT0125-47, which has the highest signal-to-noise ratio, we found that the velocity profiles of the lines tracing dense (HCN, HCO+) and lower-density (CO, [CI]) molecular gas are similar. In addition to these lines, we obtained one robust and one tentative detection of 13CO(4-3) and found an average I 12CO(4-3)/I13CO(4-3) flux ratio of 26.1−3.5+4.5, indicating a young but not pristine interstellar medium. We argue that the combination of large and slightly enriched gas reservoirs and high dense-gas fractions could explain the prodigious star formation in these systems.

An Analysis of ALMA Deep Fields and the Perceived Dearth of High-z Galaxies

Abstract Deep, pencil-beam surveys from ALMA at 1.1–1.3 mm have uncovered an apparent absence of high-redshift dusty galaxies, with existing redshift distributions peaking around z ∼ 1.5–2.5. This has led to a perceived dearth of dusty systems at z ≳ 4 and the conclusion, according to some models, that the early universe was relatively dust-poor. In this paper, we extend the backward-evolution galaxy model described by Casey et al. to the ALMA regime (in depth and area) and determine that the measured number counts and redshift distributions from ALMA deep field surveys are fully consistent with constraints of the infrared luminosity function (IRLF) at z < 2.5 determined by single-dish submillimeter and millimeter surveys conducted on much larger angular scales (∼1–10 deg2). We find that measured 1.1–1.3 mm number counts are most constraining for the measurement of the faint-end slope of the IRLF at z ≲ 2.5 instead of the prevalence of dusty galaxies at z ≳ 4. Recent studies have suggested that UV-selected galaxies at z > 4 may be particularly dust-poor, but we find that their millimeter-wave emission cannot rule out consistency with the Calzetti dust attenuation law, even by assuming relatively typical, cold-dust (T dust ≈ 30 K) spectral energy distributions. Our models suggest that the design of ALMA deep fields requires substantial revision to constrain the prevalence of z > 4 early universe obscured starbursts. The most promising avenue for detection and characterization of such early dusty galaxies will come from future ALMA 2 mm blank-field surveys covering a few hundred arcmin2 and the combination of existing and future dual-purpose 3 mm data sets.

Towards a census of high-redshift dusty galaxies with Herschel

Context. Over the last decade a large number of dusty star-forming galaxies has been discovered up to redshift z = 2 − 3 and recent studies have attempted to push the highly confused Herschel SPIRE surveys beyond that distance. To search for z ≥ 4 galaxies they often consider the sources with fluxes rising from 250 μm to 500 μm (so-called “500 μm-risers”). Herschel surveys offer a unique opportunity to efficiently select a large number of these rare objects, and thus gain insight into the prodigious star-forming activity that takes place in the very distant Universe. Aims. We aim to implement a novel method to obtain a statistical sample of 500 μm-risers and fully evaluate our selection inspecting different models of galaxy evolution. Methods. We consider one of the largest and deepest Herschel surveys, the Herschel Virgo Cluster Survey. We develop a novel selection algorithm which links the source extraction and spectral energy distribution fitting. To fully quantify selection biases we make end-to-end simulations including clustering and lensing. Results. We select 133 500 μm-risers over 55 deg2, imposing the criteria: S500 > S350 > S250, S250 > 13.2 mJy and S500 > 30 mJy. Differential number counts are in fairly good agreement with models, displaying a better match than other existing samples. The estimated fraction of strongly lensed sources is 24+6-5% based on models. Conclusions. We present the faintest sample of 500 μm-risers down to S250 = 13.2 mJy. We show that noise and strong lensing have an important impact on measured counts and redshift distribution of selected sources. We estimate the flux-corrected star formation rate density at 4 < z < 5 with the 500 μm-risers and find it to be close to the total value measured in far-infrared. This indicates that colour selection is not a limiting effect to search for the most massive, dusty z > 4 sources.

SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES): Faint-end Counts at 450 μm

Abstract The SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES) is a three-year JCMT Large Program aiming to reach the 450 μm confusion limit in the COSMOS-CANDELS region to study a representative sample of the high-redshift far-infrared galaxy population that gives rise to the bulk of the far-infrared background. We present the first-year data from STUDIES. We reached a 450 μm noise level of 0.91 mJy for point sources at the map center, covered an area of 151 arcmin2, and detected 98 and 141 sources at 4.0σ and 3.5σ, respectively. Our derived counts are best constrained in the 3.5–25 mJy regime using directly detected sources. Below the detection limits, our fluctuation analysis further constrains the slope of the counts down to 1 mJy. The resulting counts at 1–25 mJy are consistent with a power law having a slope of −2.59 (±0.10 for 3.5–25 mJy, and − 0.7 + 0.4 for 1–3.5 mJy). There is no evidence of a faint-end termination or turnover of the counts in this flux density range. Our counts are also consistent with previous SCUBA-2 blank-field and lensing-cluster surveys. The integrated surface brightness from our counts down to 1 mJy is 90.0 ± 17.2 Jy deg−2, which can account for up to 83 − 16 + 15 % of the COBE 450 μm background. We show that Herschel counts at 350 and 500 μm are significantly higher than our 450 μm counts, likely caused by its large beam and source clustering. High angular resolution instruments like SCUBA-2 at 450 μm are therefore highly beneficial for measuring the luminosity and spatial density of high-redshift dusty galaxies.

Insights on star-formation histories and physical properties of 1.2 ≤z ≲ 4 Herschel-detected galaxies

We test the impact of using variable star forming histories (SFHs) and the use of the IR luminosity (LIR) as a constrain on the physical parameters of high redshift dusty star-forming galaxies. We explore in particular the stellar properties of galaxies in relation with their location on the SFR-M* diagram. We perform SED fitting of the UV-NIR and FIR emissions of a large sample of GOODS-Herschel galaxies, for which rich multi-wavelength observations are available. We test different SFHs and imposing energy conservation in the SED fitting process, to face issues like the age-extinction degeneracy and produce SEDs consistent with observations. Our models work well for the majority of the sample, with the notable exception of the high LIR end, for which we have indications that our simple energy conservation approach cannot hold true. We find trends in the SFHs fitting our sources depending on stellar mass M* and z. Trends also emerge in the characteristic timescales of the SED models depending on the location on the SFR-M* diagram. We show that whilst using the same available observational data, we can produce galaxies less star-forming than usually inferred, if we allow declining SFHs, while properly reproducing their observables. These sources can be post-starbursts undergoing quenching, and their SFRs are potentially overestimated if inferred from their LIR. Fitting without the IR constrain leads to a strong preference for declining SFHs, while its inclusion increases the preference of rising SFHs, more so at high z, in tentative agreement with the cosmic star formation history. Keeping in mind that the sample is biased towards high LIR, the evolution shaped by our model appears as both bursty (initially) and steady-lasting (later on). The global SFH of the sample follows the cosmic SFH with a small scatter, and is compatible with the "downsizing" scenario of galaxy evolution.

HerMES: a search for high-redshift dusty galaxies in the HerMES Large Mode Survey – catalogue, number counts and early results

Selecting sources with rising flux densities towards longer wavelengths from Herschel/SPIRE maps is an efficient way to produce a catalogue rich in high-redshift (z > 4) dusty star-forming galaxies. The effectiveness of this approach has already been confirmed by spectroscopic follow-up observations, but the previously available catalogues made this way are limited by small survey areas. Here we apply a map-based search method to 274 deg$^2$ of the HerMES Large Mode Survey (HeLMS) and create a catalogue of 477 objects with SPIRE flux densities $S_{500} > S_{350} >S_{250}$ and a $5 \sigma$ cut-off $S_{500} > $ 52 mJy. From this catalogue we determine that the total number of these "red" sources is at least an order of magnitude higher than predicted by galaxy evolution models. These results are in agreement with previous findings in smaller HerMES fields; however, due to our significantly larger sample size we are also able to investigate the shape of the red source counts for the first time. We have obtained spectroscopic redshift measurements for two of our sources using the Atacama Large Millimeter/submillimeter Array (ALMA). The redshifts z = 5.1 and z = 3.8 confirm that with our selection method we can indeed find high-redshift dusty star-forming galaxies.

Herschel protocluster survey: a search for dusty star-forming galaxies in protoclusters atz= 2–3

We present a Herschel/SPIRE survey of three protoclusters at z=2-3 (2QZCluster, HS1700, SSA22). Based on the SPIRE colours (S350/S250 and S500/S350) of 250 $\mu$m sources, we selected high redshift dusty star-forming galaxies potentially associated with the protoclusters. In the 2QZCluster field, we found a 4-sigma overdensity of six SPIRE sources around 4.5' (~2.2 Mpc) from a density peak of H$\alpha$ emitters at z=2.2. In the HS1700 field, we found a 5-sigma overdensity of eight SPIRE sources around 2.1' (~1.0 Mpc) from a density peak of LBGs at z=2.3. We did not find any significant overdensities in SSA22 field, but we found three 500 $\mu$m sources are concentrated 3' (~1.4 Mpc) east to the LAEs overdensity. If all the SPIRE sources in these three overdensities are associated with protoclusters, the inferred star-formation rate densities are 10$^3$-10$^4$ times higher than the average value at the same redshifts. This suggests that dusty star-formation activity could be very strongly enhanced in z~2-3 protoclusters. Further observations are needed to confirm the redshifts of the SPIRE sources and to investigate what processes enhance the dusty star-formation activity in z~2-3 protoclusters.

THE REDSHIFT DISTRIBUTION OF DUSTY STAR-FORMING GALAXIES FROM THE SPT SURVEY

ABSTRACT We use the Atacama Large Millimeter/submillimeter Array (ALMA) in Cycle 1 to determine spectroscopic redshifts of high-redshift dusty star-forming galaxies (DSFGs) selected by their 1.4 mm continuum emission in the South Pole Telescope (SPT) survey. We present ALMA 3 mm spectral scans between 84 and 114 GHz for 15 galaxies and targeted ALMA 1 mm observations for an additional eight sources. Our observations yield 30 new line detections from CO, [C i], [N ii], H2O and NH3. We further present Atacama Pathfinder Experiment [C ii] and CO mid-J observations for seven sources for which only a single line was detected in spectral-scan data from ALMA Cycle 0 or Cycle 1. We combine the new observations with previously published and new millimeter/submillimeter line and photometric data of the SPT-selected DSFGs to study their redshift distribution. The combined data yield 39 spectroscopic redshifts from molecular lines, a success rate of >85%. Our sample represents the largest data set of its kind today and has the highest spectroscopic completeness among all redshift surveys of high-z DSFGs. The median of the redshift distribution is z = 3.9 ± 0.4, and the highest-redshift source in our sample is at z = 5.8. We discuss how the selection of our sources affects the redshift distribution, focusing on source brightness, selection wavelength, and strong gravitational lensing. We correct for the effect of gravitational lensing and find the redshift distribution for 1.4 mm selected sources with a median redshift of z = 3.1 ± 0.3. Comparing to redshift distributions selected at shorter wavelengths from the literature, we show that selection wavelength affects the shape of the redshift distribution.