Felipe Navarrete

Evolution of the dust emission of massive galaxies up to z = 4 and constraints on their dominant mode of star formation

We aim to measure the average dust and molecular gas content of massive star-forming galaxies (>3 × 1010M⊙) up to z = 4 in the COSMOS field to determine if the intense star formation observed at high redshift is induced by major mergers or is caused by large gas reservoirs. Firstly, we measured the evolution of the average spectral energy distributions as a function of redshift using a stacking analysis of Spitzer, Herschel, LABOCA, and AzTEC data for two samples of galaxies: normal star-forming objects and strong starbursts, as defined by their distance to the main sequence. We found that the mean intensity of the radiation field ⟨ U ⟩ heating the dust (strongly correlated with dust temperature) increases with increasing redshift up to z = 4 in main-sequence galaxies. We can reproduce this evolution with simple models that account for the decrease in the gas metallicity with redshift. No evolution of ⟨ U ⟩ with redshift is found in strong starbursts. We then deduced the evolution of the molecular gas fraction (defined here as Mmol/ (Mmol + M⋆)) with redshift and found a similar, steeply increasing trend for both samples. At z ~ 4, this fraction reaches ~60%. The average position of the main-sequence galaxies is on the locus of the local, normal star-forming disks in the integrated Schmidt-Kennicutt diagram (star formation rate versus mass of molecular gas), suggesting that the bulk of the star formation up to z = 4 is dominated by secular processes.

Confirming the Existence of a Quiescent Galaxy Population out to z=3: A Stacking Analysis of Mid-, Far-Infrared and Radio Data

We performed a comprehensive stacking analysis on similar to 14,200 quiescent galaxy (QG) candidates at z = 0-3 across mid-, far-infrared (MIR and FIR), and radio wavelengths. Identified via their rest-frame NUV - r and r - J colors, the QG candidates (M-* = 10(9.8-12.2) M-circle dot) have drastically different IR and radio properties depending on their 24 mu m emission strength. The fraction of QG candidates with strong 24 mu m emission (equivalent to inferred star formation rates SFR24 \textgreater= 100 M-circle dot yr(-1), hereafter “IR-bright”) increases with redshift and peaks at 15%, and their stacked MIPS 24 mu m, Herschel (PACS and SPIRE) and VLA emissions are consistent with being star-forming galaxies (SFGs). In contrast, the majority of QG candidates are faint or undetected at 24 mu m individually (i.e., SFR24 \textless 100 M-circle dot yr(-1), hereafter “IR-faint”). Their low dust-obscured SFRs derived from Herschel stacking (SFRH less than or similar to 3, 15, 50 M-circle dot yr(-1) out to z similar to 1, 2, 3) are \textgreater 2.5-12.5x lower than compared to SFGs. This is consistent with the quiescence, as expected from their low unobscured SFRs, as inferred from modeling their ultraviolet-to-NIR photometry. The discrepancy between the L-IR derived from stacking Herschel and 24 mu m indicates that IR-faint QGs have dust SEDs that are different from those of SFGs. For the most massive (M-star \textgreater= 10(11) M-circle dot) IR-faint QGs at z \textless 1.5, the stacked 1.4 GHz emission is in excess of that expected from other SFR indicators, suggesting a widespread presence of low-luminosity active galactic nuclei. Our results reaffirm the existence of a significant population of QGs out to z. =. 3, thus corroborating the need to quench star formation in galaxies at early epochs.

New insights from deep VLA data on the potentially recoiling black hole CID-42 in the COSMOS field

We present deep 3-GHz Karl G. Jansky Very Large Array (VLA) observations of the potentially recoiling black hole CID-42 in the Cosmic Evolution Survey (COSMOS) field. This galaxy shows two optical nuclei in the Hubble Space Telescope/Advanced Camera for Surveys (HST/ACS) image and a large velocity offset of ≈1300 km s^(−1) between the broad and narrow Hβ emission line although the spectrum is not spacially resolved (Civano et al. 2010). The new 3 GHz VLA data have a bandwidth of 2 GHz and to correctly interpret the flux densities imaging was done with two different methods: multiscale multifrequency (MSMF) synthesis and spectral windows (SPWs) stacking. The final resolutions and sensitivities of these maps are 0.7 arcsec with rms = 4.6 μJy beam^(−1) and 0.9 arcsec with rms = 4.8 μJy beam^(−1), respectively. With a 7σ detection, we find that the entire observed 3-GHz radio emission can be associated with the south-eastern component of CID-42, coincident with the detected X-ray emission. We use our 3 GHz data combined with other radio data from the literature ranging from 320 MHz to 9 GHz, which include the VLA, Very Long Baseline Array (VLBA) and Giant Metrewave Radio Telescope (GMRT) data, to construct a radio synchrotron spectrum of CID-42. The radio spectrum suggests a type I unobscured radio-quiet flat-spectrum active galactic nucleus (AGN) in the south-eastern component which may be surrounded by a more extended region of old synchrotron electron population or shocks generated by the outflow from the supermassive black hole (SMBH). Our data are consistent with the recoiling black hole picture but cannot rule out the presence of an obscured and radio-quiet SMBH in the north-western component.

An ALMA survey of submillimetre galaxies in the COSMOS field: The extent of the radio-emitting region revealed by 3 GHz imaging with the Very Large Array

Context. The observed spatial scale of the radio continuum emission from star-forming galaxies can be used to investigate the spatial extent of active star formation, constrain the importance of cosmic-ray transport, and examine the effects of galaxy interactions. Aims. We determine the radio size distribution of a large sample of 152 submillimetre galaxies (SMGs) in the COSMOS field that were pre-selected at 1.1 mm, and later detected with the Atacama Large Millimetre/submillimetre Array (ALMA) in the observed-frame 1.3 mm dust continuum emission at a signal-to-noise ratio (S/N) of ≥5. Methods. We used the deep, subarcsecond-resolution (1σ = 2.3μJy beam^(-1); . 75) centimetre radio continuum observations taken by the Karl G. Jansky Very Large Array (VLA)-COSMOS 3 GHz Large Project. Results. One hundred and fifteen of the 152 target SMGs (76% ± 7%) were found to have a 3 GHz counterpart (≥ 4.2σ), which renders the radio detection rate notably high. The median value of the deconvolved major axis full width at half maximum (FWHM) size at 3 GHz is derived to be 0.59 ± 0.05 , or 4.6 ± 0.4 kpc in physical units, where the median redshift of the sources is z = 2.23 ± 0.13 (23% are spectroscopic and 77% are photometric values). The radio sizes are roughly log-normally distributed, and they show no evolutionary trend with redshift, or difference between different galaxy morphologies. We also derived the spectral indices between 1.4 and 3 GHz, and 3 GHz brightness temperatures for the sources, and the median values were found to be α_(1.4 GHz)^(3 GHz) = -0.67 (S_ν ∝ ν^α) and T_B = 12.6 ± 2 K. Three of the target SMGs, which are also detected with the Very Long Baseline Array (VLBA) at 1.4 GHz (AzTEC/C24b, 61, and 77a), show clearly higher brightness temperatures than the typical values, reaching T_B(3 GHz) > 10^(4.03) K for AzTEC/C61. Conclusions. The derived median radio spectral index agrees with a value expected for optically thin non-thermal synchrotron radiation, and the low median 3 GHz brightness temperature shows that the observed radio emission is predominantly powered by star formation and supernova activity. However, our results provide a strong indication of the presence of an active galactic nucleus in the VLBA and X-ray-detected SMG AzTEC/C61 (high TB and an inverted radio spectrum). The median radio-emitting size we have derived is ~ 1.5–3 times larger than the typical far-infrared dust-emitting sizes of SMGs, but similar to that of the SMGs’ molecular gas component traced through mid-J line emission of carbon monoxide. The physical conditions of SMGs probably render the diffusion of cosmic-ray electrons inefficient, and hence an unlikely process to lead to the observed extended radio sizes. Instead, our results point towards a scenario where SMGs are driven by galaxy interactions and mergers. Besides triggering vigorous starbursts, galaxy collisions can also pull out the magnetised fluids from the interacting disks, and give rise to a taffy-like synchrotron-emitting bridge. This provides an explanation for the spatially extended radio emission of SMGs, and can also cause a deviation from the well-known infrared-radio correlation owing to an excess radio emission. Nevertheless, further high-resolution observations are required to examine the other potential reasons for the very compact dust-emitting sizes of SMGs, such as the radial dust temperature and metallicity gradients.

An ALMA survey of submillimeter galaxies in the COSMOS field: Multiwavelength counterparts and redshift distribution

We carried out targeted ALMA observations of 129 fields in the COSMOS region at 1.25 mm, detecting 152 galaxies at S/N ≥ 5with an average continuum RMS of 150 μJy. These fields represent a S/N-limited sample of AzTEC/ASTE sources with 1.1 mm S/N≥ 4 over an area of 0.72 square degrees. Given ALMA’s fine resolution and the exceptional spectroscopic and multiwavelength photometric data available in COSMOS, this survey allows us unprecedented power in identifying submillimeter galaxy counterparts and determining their redshifts through spectroscopic or photometric means. In addition to 30 sources with prior spectroscopic redshifts, we identified redshifts for 113 galaxies through photometric methods and an additional nine sources with lower limits, which allowed a statistically robust determination of the redshift distribution. We have resolved 33 AzTEC sources into multi-component systems and our redshifts suggest that nine are likely to be physically associated. Our overall redshift distribution peaks at z ~ 2.0 with a high-redshift tail skewing the median redshift to z = 2.48 ± 0.05. We find that brighter millimeter sources are preferentially found at higher redshifts. Our faintestsources, with S_(1.25 mm) 1.8 mJy, have a median redshift of z = 3.08 ± 0.17. After accounting for spectral energy distribution shape and selection effects, these results are consistent with several previous submillimeter galaxy surveys, and moreover, support the conclusion that the submillimeter galaxy redshift distribution is sensitive to survey depth.

An ALMA survey of submillimetre galaxies in the COSMOS field: Physical properties derived from energy balance spectral energy distribution modelling

Context. Submillimetre galaxies (SMGs) represent an important source population in the origin and cosmic evolution of the most massive galaxies. Hence, it is imperative to place firm constraints on the fundamental physical properties of large samples of SMGs. Aims. We determine the physical properties of a sample of SMGs in the COSMOS field that were pre-selected at the observed-frame wavelength of λ_(obs) = 1.1 mm, and followed up at λ_(obs) = 1.3 mm with the Atacama Large Millimetre/submillimetre Array (ALMA). Methods. We used the MAGPHYS model package to fit the panchromatic (ultraviolet to radio) spectral energy distributions (SEDs) of 124 of the target SMGs, which lie at a median redshift of z = 2.30 (19.4% are spectroscopically confirmed). The SED analysis was complemented by estimating the gas masses of the SMGs by using the λ_(obs) = 1.3 mm dust emission as a tracer of the molecular gas component. Results. The sample median and 16th–84th percentile ranges of the stellar masses, obscured star formation rates, dust temperatures, and dust and gas masses were derived to be log(M⋆/M⊙) = 11.09^(+0.41)_(-0.53), SFR = 402^(+661)_(-233) M⊙ yr^(-1), T_(dust) = 39.7^(+9.7)_(-7.4) K, log(M_(dust)/M⊙) = 9.01^(+0.20)_(-0.31), and log(M_(gas)/M⊙ = 11.34^(+0.20)_(-0.23), respectively. The M_(dust)/M⋆ ratio was found to decrease as a function of redshift, while the M_(gas)/M_(dust) ratio shows the opposite, positive correlation with redshift. The derived median gas-to-dust ratio of 120^(+73)_(-30) agrees well with the canonical expectation. The gas fraction (M_(gas)/(M_(gas) + M⋆)) was found to range from 0.10 to 0.98 with a median of 0.62^(+0.27)_(-0.23). We found that 57.3% of our SMGs populate the main sequence (MS) of star-forming galaxies, while 41.9% of the sources lie above the MS by a factor of greater than three (one source lies below the MS). These super-MS objects, or starbursts, are preferentially found at z ≳ 3, which likely reflects the sensitivity limit of our source selection. We estimated that the median gas consumption timescale for our SMGs is ~535 Myr, and the super-MS sources appear to consume their gas reservoir faster than their MS counterparts. We found no obvious stellar mass–size correlations for our SMGs, where the sizes were measured in the observed-frame 3 GHz radio emission and rest-frame UV. However, the largest 3 GHz radio sizes are found among the MS sources. Those SMGs that appear irregular in the rest-frame UV are predominantly starbursts, while the MS SMGs are mostly disk-like. Conclusions. The physical parameter distributions of our SMGs and those of the equally bright, 870 μm selected SMGs in the ECDFS field (the so-called ALESS SMGs) are unlikely to be drawn from common parent distributions. This might reflect the difference in the pre-selection wavelength. Albeit being partly a selection bias, the abrupt jump in specific SFR and the offset from the MS of our SMGs at z ≳ 3 might also reflect a more efficient accretion from the cosmic gas streams, higher incidence of gas-rich major mergers, or higher star formation efficiency at z ≳ 3. We found a rather flat average trend between the SFR and dust mass, but a positive SFR−M_(gas) correlation. However, to address the questions of which star formation law(s) our SMGs follow, and how they compare with the Kennicutt-Schmidt law, the dust-emitting sizes of our sources need to be measured. Nonetheless, the larger radio-emitting sizes of the MS SMGs compared to starbursts is a likely indication of their more widespread, less intense star formation activity. The irregular rest-frame UV morphologies of the starburst SMGs are likely to echo their merger nature. The current stellar mass content of the studied SMGs is very high, so they must quench to form the so-called red-and-dead massive ellipticals. Our results suggest that the transition from high-z SMGs to local ellipticals via compact, quiescent galaxies (cQGs) at z ~ 2 might not be universal, and the latter population might also descend from the so-called blue nuggets. However, z ≳ 4 SMGs could be the progenitors of higher redshift, z ≳ 3 cQGs, while our results are also consistent with the possibility that ultra-massive early-type galaxies found at 1.2 ≲ z ≲ 2 experienced an SMG phase at z ≤ 3.