We present a new, publicly available library of dust spectral energy distributions (SEDs). These SEDs are characterized by only three parameters: the dust mass (Mdust), the dust temperature (Tdust), ...and the mid-to-total infrared color (IR8 ≡ LIR/L8). The latter measures the relative contribution of polycyclic aromatic hydrocarbon (PAH) molecules to the total infrared luminosity. We used this library to model star-forming galaxies at 0.5 < z < 4 in the deep CANDELS fields, using both individual detections and stacks of Herschel and ALMA imaging, and extending this sample to z = 0 using the Herschel Reference Survey. At first order, the dust SED of a galaxy was observed to be independent of stellar mass, but evolving with redshift. We found trends of increasing Tdust and IR8 with redshift and distance from the SFR–M∗ main sequence, and quantified for the first time their intrinsic scatter. Half of the observed variations of these parameters was captured by the above empirical relations, and after subtracting the measurement errors we found residual scatters of ΔTdust/Tdust = 12% and Δlog IR8 = 0.18 dex. We observed second order variations with stellar mass: massive galaxies (M∗ > 1011M⊙) at z ≤ 1 have slightly lower temperatures indicative of a reduced star formation efficiency, while low mass galaxies (M∗ < 1010M⊙) at z ≥ 1 showed reduced PAH emission, possibly linked to their lower metallicities. Building on these results, we constructed high-fidelity mock galaxy catalogs to predict the accuracy of infrared luminosities and dust masses determined using a single broadband measurement. Using a single James Webb Space Telescope (JWST) MIRI band, we found that LIR is typically uncertain by 0.15 dex, with a maximum of 0.25 dex when probing the rest-frame 8 μm, and this is not significantly impacted by typical redshift uncertainties. On the other hand, we found that ALMA bands 8 to 7 and 6 to 3 measured the dust mass at better than 0.2 and 0.15 dex, respectively, and independently of redshift, while bands 9 to 6 only measured LIR at better than 0.2 dex at z > 1, 3.2, 3.8, and 5.7, respectively. Starburst galaxies had their LIR significantly underestimated when measured by a single JWST or ALMA band, while their dust mass from a single ALMA band were moderately overestimated. This dust library and the results of this paper can be used immediately to improve the design of observing proposals, and interpret more accurately the large amount of archival data from Spitzer, Herschel and ALMA.
ABSTRACT
By compiling a comprehensive census of literature studies, we investigate the evolution of the main sequence (MS) of star-forming galaxies (SFGs) in the widest range of redshift (0 < z < 6) ...and stellar mass (108.5–1011.5 M⊙) ever probed. We convert all observations to a common calibration and find a remarkable consensus on the variation of the MS shape and normalization across cosmic time. The relation exhibits a curvature towards the high stellar masses at all redshifts. The best functional form is governed by two parameters: the evolution of the normalization and the turnover mass (M0(t)), which both evolve as a power law of the Universe age. The turn-over mass determines the MS shape. It marginally evolves with time, making the MS slightly steeper towards z ∼ 4–6. At stellar masses below M0(t), SFGs have a constant specific SFR (sSFR), while above M0(t) the sSFR is suppressed. We find that the MS is dominated by central galaxies. This allows to turn M0(t) into the corresponding host halo mass. This evolves as the halo mass threshold between cold and hot accretion regimes, as predicted by the theory of accretion, where the central galaxy is fed or starved of cold gas supply, respectively. We, thus, argue that the progressive MS bending as a function of the Universe age is caused by the lower availability of cold gas in haloes entering the hot accretion phase, in addition to black hole feedback. We also find qualitatively the same trend in the largest sample of star-forming galaxies provided by the IllustrisTNG simulation. Nevertheless, we still note large quantitative discrepancies with respect to observations, in particular at the high-mass end. These can not be easily ascribed to biases or systematics in the observed SFRs and the derived MS.
The star-formation history (SFH) of galaxies is a key assumption to derive their physical properties and can lead to strong biases. In this work, we derive the SFH of main sequence (MS) galaxies and ...show how the peak SFH of a galaxy depends on its seed mass at, for example, z = 5. This seed mass reflects the galaxy’s underlying dark matter (DM) halo environment. We show that, following the MS, galaxies undergo a drastic slow down of their stellar mass growth after reaching the peak of their SFH. According to abundance matching, these masses correspond to hot and massive DM halos which state could result in less efficient gas inflows on the galaxies and thus could be the origin of limited stellar mass growth. As a result, we show that galaxies, still on the MS, can enter the passive region of the UVJ diagram while still forming stars. The best fit to the MS SFH is provided by a right skew peak function for which we provide parameters depending on the seed mass of the galaxy. The ability of the classical analytical SFHs to retrieve the star-formation rate (SFR) of galaxies from spectral energy distribution (SED) fitting is studied. Due to mathematical limitations, the exponentially declining and delayed SFH struggle to model high SFR, which starts to be problematic at z > 2. The exponentially rising and log-normal SFHs exhibit the opposite behavior with the ability to reach very high SFR, and thus model starburst galaxies, but they are not able to model low values such as those expected at low redshift for massive galaxies. By simulating galaxies SED from the MS SFH, we show that these four analytical forms recover the SFR of MS galaxies with an error dependent on the model and the redshift. They are, however, sensitive enough to probe small variations of SFR within the MS, with an error ranging from 5 to 40% depending on the SFH assumption and redshift; but all the four fail to recover the SFR of rapidly quenched galaxies. However, these SFHs lead to an artificial gradient of age, parallel to the MS, which is not exhibited by the simulated sample. This gradient is also produced on real data as we show using a sample of real galaxies with redshifts between 1.5 and 2.5. Here, we propose an SFH composed of a delayed form to model the bulk of stellar population with the addition of a flexibility in the recent SFH. This SFH provides very good estimates of the SFR of MS, starbursts, and rapidly quenched galaxies at all redshift. Furthermore, when used on the real sample, the age gradient disappears which show its dependency on the SFH assumption made to perform the SED fitting.
We present ALMA Band 9 observations of the C II158 μm emission for a sample of 10 main-sequence galaxies at redshift z ˜ 2, with typical stellar masses (log M⋆/M⊙ ˜ 10.0-10.9) and star formation ...rates (˜35-115 M⊙ yr-1). Given the strong and well-understood evolution of the interstellar medium from the present to z = 2, we investigate the behaviour of the C II emission and empirically identify its primary driver. We detect C II from six galaxies (four secure and two tentative) and estimate ensemble averages including non-detections. The C II-to-infrared luminosity ratio (C II/LIR) of our sample is similar to that of local main-sequence galaxies (˜2 × 10-3), and ˜10 times higher than that of starbursts. The C II emission has an average spatial extent of 4-7 kpc, consistent with the optical size. Complementing our sample with literature data, we find that the C II luminosity correlates with galaxies' molecular gas mass, with a mean absolute deviation of 0.2 dex and without evident systematics: the C II-to-H2 conversion factor (α _C II ˜ 30 M⊙/L⊙) is largely independent of galaxies' depletion time, metallicity, and redshift. C II seems therefore a convenient tracer to estimate galaxies' molecular gas content regardless of their starburst or main-sequence nature, and extending to metal-poor galaxies at low and high redshifts. The dearth of C II emission reported for z > 6-7 galaxies might suggest either a high star formation efficiency or a small fraction of ultraviolet light from star formation reprocessed by dust.
Aims. We use high-resolution continuum images obtained with the Atacama Large Millimeter Array (ALMA) to probe the surface density of star formation in z ~ 2 galaxies and study the different physical ...properties between galaxies within and above the star-formation main sequence of galaxies. Methods. We use ALMA images at 870 μm with 0.2 arcsec resolution in order to resolve star formation in a sample of eight star-forming galaxies at z ~ 2 selected among the most massive Herschel galaxies in the GOODS-South field. This sample is supplemented with eleven galaxies from the public data of the 1.3 mm survey of the Hubble Ultra-Deep Field, HUDF. We derive dust and gas masses for the galaxies, compute their depletion times and gas fractions, and study the relative distributions of rest-frame ultraviolet (UV) and far-infrared (FIR) light. Results. ALMA reveals systematically dense concentrations of dusty star formation close to the center of the stellar component of the galaxies. We identify two different starburst regimes: (i) the classical population of starbursts located above the SFR-M⋆ main sequence, with enhanced gas fractions and short depletion times and (ii) a sub-population of galaxies located within the scatter of the main sequence that experience compact star formation with depletion timescales typical of starbursts of ~150 Myr. In both starburst populations, the FIR and UV are distributed in distinct regions and dust-corrected star formation rates (SFRs) estimated using UV-optical-near-infrared data alone underestimate the total SFR. Starbursts hidden in the main sequence show instead the lowest gas fractions of our sample and could represent the last stage of star formation prior to passivization. Being Herschel-selected, these main sequence galaxies are located in the high-mass end of the main sequence, hence we do not know whether these “starbursts hidden in the main sequence” also exist below 1011 M⊙. Active galactic nuclei (AGNs) are found to be ubiquitous in these compact starbursts, suggesting that the triggering mechanism also feeds the central black hole or that the active nucleus triggers star formation.
By using a set of different star formation rate (SFR) indicators, including Wide-field Infrared Survey Explorer (WISE) mid-infrared and H α emission, we study the slope of the main sequence (MS) of ...local star-forming galaxies at stellar masses larger than 10^{10} M_{⊙ }. The slope of the relation strongly depends on the SFR indicator used. In all cases, the local MS shows a bending at high stellar masses with respect to the slope obtained in the low-mass regime. While the distribution of galaxies in the upper envelope of the MS is consistent with a lognormal distribution, the lower envelope shows an excess of galaxies, which increases as a function of the stellar mass but varies as a function of the SFR indicator used. The scatter of the best lognormal distribution increases with stellar mass from ∼0.3 dex at 10^{10} M_{⊙ } to ∼0.45 at 10^{11} M_{⊙ }. The MS high-mass end is dominated by central galaxies of group-sized haloes with a red bulge and a disc redder than the lower mass counterparts. We argue that the MS bending in this region is due to two processes: (i) the formation of a bulge component as a consequence of the increased merger activity in groups, and (ii) the cold gas starvation induced by the hot halo environment, which cuts off the gas inflow on to the disc. Similarly, the increase of the MS scatter at high stellar masses would be explained by the larger spread of star formation histories of central group and cluster galaxies with respect to lower mass systems.
Submillimeter/millimeter observations of dusty star-forming galaxies with the Atacama Large Millimeter/submillimeter Array (ALMA) have shown that dust continuum emission generally occurs in compact ...regions smaller than the stellar distribution. However, it remains to be understood how systematic these findings are. Studies often lack homogeneity in the sample selection, target discontinuous areas with inhomogeneous sensitivities, and suffer from modest
u
v
coverage coming from single array configurations. GOODS-ALMA is a 1.1 mm galaxy survey over a continuous area of 72.42 arcmin
2
at a homogeneous sensitivity. In this version 2.0, we present a new low resolution dataset and its combination with the previous high resolution dataset from the survey, improving the
u
v
coverage and sensitivity reaching an average of
σ
= 68.4 μJy beam
−1
. A total of 88 galaxies are detected in a blind search (compared to 35 in the high resolution dataset alone), 50% at
S
/
N
peak
≥ 5 and 50% at 3.5 ≤
S
/
N
peak
≤ 5 aided by priors. Among them, 13 out of the 88 are optically dark or faint sources (
H
- or
K
-band dropouts). The sample dust continuum sizes at 1.1 mm are generally compact, with a median effective radius of
R
e
= 0
.
″
10 ± 0
.
″
05 (a physical size of
R
e
= 0.73 ± 0.29 kpc at the redshift of each source). Dust continuum sizes evolve with redshift and stellar mass resembling the trends of the stellar sizes measured at optical wavelengths, albeit a lower normalization compared to those of late-type galaxies. We conclude that for sources with flux densities
S
1.1 mm
> 1 mJy, compact dust continuum emission at 1.1 mm prevails, and sizes as extended as typical star-forming stellar disks are rare. The
S
1.1 mm
< 1 mJy sources appear slightly more extended at 1.1 mm, although they are still generally compact below the sizes of typical star-forming stellar disks.
Abstract
We investigate the fraction of close pairs and morphologically identified mergers on and above the star-forming main sequence (MS) at 0.2 ≤ z ≤2.0. The novelty of our work lies in the use of ...a non-parametric morphological classification performed on resolved stellar mass maps, reducing the contamination by non-interacting, high-redshift clumpy galaxies. We find that the merger fraction rapidly rises to ≥70 per cent above the MS, implying that – already at z ≳ 1 – starburst (SB) events (ΔMS ≥ 0.6) are almost always associated with a major merger (1:1 to 1:6 mass ratio). The majority of interacting galaxies in the SB region are morphologically disturbed, late-stage mergers. Pair fractions show little dependence on MS offset and pairs are more prevalent than late-stage mergers only in the lower half of the MS. In our sample, major mergers on the MS occur with a roughly equal frequency of ∼5–10 per cent at all masses ≳ 1010 M⊙. The MS major merger fraction roughly doubles between z = 0.2 and 2, with morphological mergers driving the overall increase at z ≳ 1. The differential redshift evolution of interacting pairs and morphologically classified mergers on the MS can be reconciled by evolving observability time-scales for both pairs and morphological disturbances. The observed variation of the late-stage merger fraction with ΔMS follows the perturbative 2-Star Formation Mode model, where any MS galaxy can experience a continuum of different star formation rate enhancements. This points to an SB–merger connection not only for extreme events, but also more moderate bursts which merely scatter galaxies upward within the MS, rather than fully elevating them above it.
ABSTRACT We present 0 4 resolution extinction-independent distributions of star formation and dust in 11 star-forming galaxies (SFGs) at z = 1.3-3.0. These galaxies are selected from sensitive ...blank-field surveys of the 2′ × 2′ Hubble Ultra-Deep Field at λ = 5 cm and 1.3 mm using the Karl G. Jansky Very Large Array and Atacama Large Millimeter/submillimeter Array. They have star formation rates (SFRs), stellar masses, and dust properties representative of massive main-sequence SFGs at z ∼ 2. Morphological classification performed on spatially resolved stellar mass maps indicates a mixture of disk and morphologically disturbed systems; half of the sample harbor X-ray active galactic nuclei (AGNs), thereby representing a diversity of z ∼ 2 SFGs undergoing vigorous mass assembly. We find that their intense star formation most frequently occurs at the location of stellar-mass concentration and extends over an area comparable to their stellar-mass distribution, with a median diameter of 4.2 1.8 kpc. This provides direct evidence of galaxy-wide star formation in distant blank-field-selected main-sequence SFGs. The typical galactic-average SFR surface density is 2.5 M yr−1 kpc−2, sufficiently high to drive outflows. In X-ray-selected AGN where radio emission is enhanced over the level associated with star formation, the radio excess pinpoints the AGNs, which are found to be cospatial with star formation. The median extinction-independent size of main-sequence SFGs is two times larger than those of bright submillimeter galaxies, whose SFRs are 3-8 times larger, providing a constraint on the characteristic SFR (∼300 M yr−1) above which a significant population of more compact SFGs appears to emerge.
We present high-fidelity, 30 mas (200 pc) resolution ALMA rest-frame 240 m observations of cold dust emission in three typical main-sequence star-forming galaxies (SFGs) at z ∼ 3 in the Hubble ...Ultra-Deep Field (HUDF). The cold dust is distributed within the smooth disklike central regions of star formation 1-3 kpc in diameter, despite their complex and disturbed rest-frame UV and optical morphologies. No dust substructures or clumps are seen down to 1-3 yr−1 (1 ) per 200 pc beam. No dust emission is observed at the locations of UV-emitting clumps, which lie 2-10 kpc from the bulk of star formation. Clumpy substructures can contribute no more than 1%-7% of the total star formation in these galaxies (3 upper limits). The lack of star-forming substructures in our HUDF galaxies is to be contrasted with the multiple substructures characteristic of submillimeter-selected galaxies (SMGs) at the same cosmic epoch, particularly the far-IR-bright SMGs with similarly high-fidelity ALMA observations of Hodge et al. Individual star-forming substructures in these SMGs contain ∼10%-30% of their total star formation. A substructure in these SMGs is often comparably bright in the far-infrared to (or in some cases brighter than) our typical SFGs, suggesting that these SMGs originate from a class of disruptive events involving multiple objects at the scale of our HUDF galaxies. The scale of the disruptive event found in our main-sequence SFGs, characterized by the lack of star-forming substructures at our resolution and sensitivity, could be less violent, e.g., gas-rich disk instability or minor mergers.