Search

Exploiting IRAM 30 m CO spectroscopy, we find that SDSS post-merger galaxies display gas fractions and depletion times enhanced by 25-50%, a mildly higher CO excitation, and standard molecular-to-atomic gas ratios, compared to non-interacting galaxies with similar redshift, stellar mass ($M_{\star}$) and star-formation rate (SFR). To place these results in context, we compile further samples of interacting or starbursting galaxies, from pre-coalescence kinematic pairs to post-starbursts, carefully homogenising gas mass, $M_{\star}$ and SFR measurements in the process. We explore systematics by duplicating our analysis for different SFR and $M_{\star}$ estimators, finding good qualitative agreement in general. Molecular gas fractions and depletion times are enhanced in interacting pairs, albeit less than for post-mergers. Among all samples studied, gas fraction and depletion time enhancements appear largest in young (a few 100 Myr) post-starbursts. While there is only partial overlap between post-mergers and post-starbursts, this suggests that molecular gas reservoirs are boosted throughout most stages of galaxy interactions, plausibly due to torque-driven inflows of halo gas and gas compression. The gas fraction and depletion time offsets of mergers and post-starbursts anti-correlate with their distance from the galaxy main sequence $\Delta({\rm MS})$, evidencing the role of SFE in driving the high SFRs of the strongest starbursts. Post-starbursts display the steepest dependency of gas fraction and SFE-offsets on $\Delta({\rm MS})$, with an evolving normalisation that reflects gas reservoir depletion over time. Our multi-sample analysis paints a coherent picture of the starburst-merger throughout the low-z merger sequence. It reconciles contradictory literature findings by highlighting that gas fraction enhancements and SFE variations both play their part in merger-driven star formation., Comment: 46 pages, 21 figures. Submitted to the Open Journal of Astrophysics

Recent deep and wide radio surveys extend the studies for radio-excess active galactic nuclei (radio-AGNs) to lower luminosities and higher redshifts, providing new insights into the abundance and physical origin of radio-AGNs. Here we focus on the cosmic evolution, physical properties and AGN-host galaxy connections of radio-AGNs selected from a sample of ~ 500,000 galaxies at 0 < z < 4 in GOODS-N, GOODS-S, and COSMOS fields. Combining deep radio data with multi-band, de-blended far-infrared (FIR) and sub-millimeter data, we identify 1162 radio-AGNs through radio excess relative to the FIR-radio relation. We study the cosmic evolution of 1.4 GHz radio luminosity functions (RLFs) for star-forming galaxies (SFGs) and radio-AGNs, which are well described by a pure luminosity evolution of $L_*\propto (1+z)^{-0.31z+3.41}$ and a pure density evolution of $\Phi_*\propto (1+z)^{-0.80z+2.88}$, respectively. We derive the turnover luminosity above which the number density of radio-AGNs surpasses that of SFGs. This crossover luminosity increases as increasing redshift, from $10^{22.9}$ W Hz$^{-1}$ at z ~ 0 to $10^{25.2}$ W Hz$^{-1}$ at z ~ 4. At full redshift range (0 < z < 4), we further derive the probability ($p_{radio}$) of SFGs and quiescent galaxies (QGs) hosting a radio-AGN as a function of stellar mass ($M_*$), radio luminosity ($L_R$), and redshift (z), which yields $p_{radio}\propto (1+z)^{3.54}M_*^{1.02}L_R^{-0.90}$ for SFGs, and $p_{radio}\propto (1+z)^{2.38}M_*^{1.39}L_R^{-0.60}$ for QGs, respectively. It indicates that radio-AGNs in QGs prefer to reside in more massive galaxies with larger $L_R$ than those in SFGs, and radio-AGN fraction increases towards higher redshift in both SFGs and QGs with a more rapid increase in SFGs. Further, we find that the radio-AGN fraction depends on accretion states of BHs and redshift in SFGs, while in QGs it also depends on BH (or galaxy) mass., Comment: Accepted by Astronomy & Astrophysics (A&A), 25 pages, 17 figures, 5 tables

Sources that leak Lyman-continuum (LyC) photons and lead to the reionisation of the universe are intensely studied using multiple observing facilities. Recently, the Low-redshift LyC Survey (LzLCS) has found the first large sample of LyC emitting galaxies at low redshift ($z\sim 0.3$) with the Hubble Space Telescope/Cosmic Origins Spectrograph. The LzLCS sample contains a robust estimate of the LyC escape fraction ($f_\mathrm{esc}^\mathrm{LyC}$) for 66 galaxies spanning a wide range of $f_\mathrm{esc}^\mathrm{LyC}$. Here we, for the first time, aim to study the radio continuum (RC) properties of LzLCS sources and their dependence on $f_\mathrm{esc}^\mathrm{LyC}$. We present Karl G. Jansky Very Large Array RC observations at C (4-8 GHz), S (2-4 GHz) and L (1-2 GHz) bands for a sub-sample of the LzLCS sources. The radio spectral index ($\alpha^{\mathrm{3GHz}}_\mathrm{6GHz}$) spans a wide range from being flat ( $\geq -0.1$) to very steep ($\leq -1.0$). The strongest leakers in our sample show flat $\alpha^{\mathrm{3GHz}}_\mathrm{6GHz}$, weak leakers have $\alpha^{\mathrm{3GHz}}_\mathrm{6GHz}$ close to normal star-forming galaxies, and non-leakers are characterized by steep $\alpha^{\mathrm{3GHz}}_\mathrm{6GHz}$. We argue that a combination of young ages, free-free absorption, and a flat cosmic-ray energy spectrum can altogether lead to a flat $\alpha^{\mathrm{3GHz}}_\mathrm{6GHz}$ for strong leakers. Non-leakers are characterized by steep spectra which can arise due to break/cutoff at high frequencies. Such a cutoff in the spectrum can arise in a single injection model of CRs characteristic of galaxies which have recently stopped star formation. Such a relation between $\alpha^{\mathrm{3GHz}}_\mathrm{6GHz}$ and $f_\mathrm{esc}^\mathrm{LyC}$ hints at the interesting role of supernovae, CRs, and magnetic fields in facilitating the escape ( and/or the lack) of LyC photons. (Abridged), Comment: 20 pages, 17 figures, 4 tables, 8 appendix pages, Astronomy & Astrophysics in press

SKA-MID surveys will be the first in the radio domain to achieve clearly sub-arcsecond resolution at high sensitivity over large areas, opening new science applications for galaxy evolution. To investigate the potential of these surveys, we create simulated SKA-MID images of a $\sim$0.04 deg$^{2}$ region of GOODS-North, constructed using multi-band HST imaging of 1723 real galaxies containing significant substructure at $0

We present NOEMA CO (2-1) line and ALMA 870 $\mu$m continuum observations of a main-sequence galaxy at $z=1.38$. The galaxy was initially selected as a "gas-giant", based on the gas mass derived from sub-mm continuum (log$(M_{\rm gas}/M_{\odot})=11.20\pm0.20$), however the gas mass derived from CO (2-1) luminosity brings down the gas mass to a value consistent with typical star-forming galaxies at that redshift (log$(M_{\rm gas}/M_{\odot})=10.84\pm0.03$). Despite that the dust-to-stellar mass ratio remains elevated above the scaling relations by a factor of 5. We explore the potential physical picture and consider an underestimated stellar mass and optically thick dust as possible causes. Based on the updated gas-to-stellar mass ratio we rule out the former, and while the latter can contribute to the dust mass overestimate it is still not sufficient to explain the observed physical picture. Instead, possible explanations include enhanced HI reservoirs, CO-dark H$_2$ gas, an unusually high metallicity, or the presence of an optically dark, dusty contaminant. Using the ALMA data at 870 $\mu$m coupled with $HST$/ACS imaging, we find extended morphology in dust continuum and clumpy star-formation in rest-frame UV in this galaxy, and a tentative $\sim 10$ kpc dusty arm is found bridging the galaxy center and a clump in F814W image. The galaxy shows levels of dust obscuration similar to the so-called $HST$-dark galaxies at higher redshifts, and would fall into the optically faint/dark $JWST$ color-color selection at $z>2$. It is therefore possible that our object could serve as low-$z$ analog of the $HST$-dark populations. This galaxy serves as a caveat to the gas masses based on the continuum alone, with a larger sample required to unveil the full picture., Comment: 9 pages, 5 figures, 1 table. Submitted to A&A Letters

Our current understanding of the cosmic star formation history at z>3 is primarily based on UV-selected galaxies (i.e., LBGs). Recent studies of H-dropouts have revealed that we may be missing a large proportion of star formation that is taking place in massive galaxies at z>3. In this work, we extend the H-dropout criterion to lower masses to select optically dark/faint galaxies (OFGs), in order to complete the census between LBGs and H-dropouts. Our criterion (H> 26.5 mag & [4.5] < 25 mag) combined with a de-blending technique is designed to select not only extremely dust-obscured massive galaxies but also normal star-forming galaxies. In total, we identified 27 OFGs at z_phot > 3 (z_med=4.1) in the GOODS-ALMA field, covering a wide distribution of stellar masses with log($M_{\star}$/$M_{\odot}$) = 9.4-11.1. We find that up to 75% of the OFGs with log($M_{\star}$/$M_{\odot}$) = 9.5-10.5 were neglected by previous LBGs and H-dropout selection techniques. After performing stacking analyses, the OFGs exhibit shorter gas depletion timescales, slightly lower gas fractions, and lower dust temperatures than typical star-forming galaxies. Their SFR_tot (SFR_ IR+SFR_UV) is much larger than SFR_UVcorr (corrected for dust extinction), with SFR_tot/SFR_UVcorr = $8\pm1$, suggesting the presence of hidden dust regions in the OFGs that absorb all UV photons. The average dust size measured by a circular Gaussian model fit is R_e(1.13 mm)=1.01$\pm$0.05 kpc. We find that the cosmic SFRD at z>3 contributed by massive OFGs is at least two orders of magnitude higher than the one contributed by equivalently massive LBGs. Finally, we calculate the combined contribution of OFGs and LBGs to the cosmic SFRD at z=4-5 to be 4 $\times$ 10$^{-2}$ $M_{\odot}$ yr$^{-1}$Mpc$^{-3}$, which is about 0.15 dex (43%) higher than the SFRD derived from UV-selected samples alone at the same redshift., Comment: 26 pages, 14 figures, 8 tables. Accepted for publication in A&A

Future deep HI surveys will be essential for understanding the nature of galaxies and the content of the Universe. However, the large volume of these data will require distributed and automated processing techniques. We introduce LiSA, a set of python modules for the denoising, detection and characterization of HI sources in 3D spectral data. LiSA was developed and tested on the Square Kilometer Array Science Data Challenge 2 dataset, and contains modules and pipelines for easy domain decomposition and parallel execution. LiSA contains algorithms for 2D-1D wavelet denoising using the starlet transform and flexible source finding using null-hypothesis testing. These algorithms are lightweight and portable, needing only a few user-defined parameters reflecting the resolution of the data. LiSA also includes two convolutional neural networks developed to analyse data cubes which separate HI sources from artifacts and predict the HI source properties. All of these components are designed to be as modular as possible, allowing users to mix and match different components to create their ideal pipeline. We demonstrate the performance of the different components of LiSA on the SDC2 dataset, which is able to find 95% of HI sources with SNR > 3 and accurately predict their properties., Comment: 18 pages, 16 figures

Radio free-free emission is considered to be one of the most reliable tracers of star formation in galaxies. However, as it constitutes the faintest part of the radio spectrum -- being roughly an order of magnitude less luminous than radio synchrotron emission at the GHz frequencies typically targeted in radio surveys -- the usage of free-free emission as a star formation rate tracer has mostly remained limited to the local Universe. Here we perform a multi-frequency radio stacking analysis using deep Karl G. Jansky Very Large Array observations at 1.4, 3, 5, 10 and 34 GHz in the COSMOS and GOODS-North fields to probe free-free emission in typical galaxies at the peak of cosmic star formation. We find that $z \sim 0.5 - 3$ star-forming galaxies exhibit radio emission at rest-frame frequencies of $\sim 65 - 90$ GHz that is $\sim 1.5 - 2\times$ fainter than would be expected from a simple combination of free-free and synchrotron emission, as in the prototypical starburst galaxy M82. We interpret this as a deficit in high-frequency synchrotron emission, while the level of free-free emission is as expected from M82. We additionally provide the first constraints on the cosmic star formation history using free-free emission at $0.5 \lesssim z \lesssim 3$, which are in good agreement with more established tracers at high redshift. In the future, deep multi-frequency radio surveys will be crucial in order to accurately determine the shape of the radio spectrum of faint star-forming galaxies, and to further establish radio free-free emission as a tracer of high-redshift star formation., Comment: Accepted for publication in ApJ; 22 pages + appendices, 11 figures in main text

We study the interstellar medium (ISM) properties as a function of the molecular gas size for 77 infrared-selected galaxies at $z \sim 1.3$. Molecular gas sizes are measured on ALMA images that combine CO(2-1), CO(5-4) and underlying continuum observations, and include CO(4-3), CO(7-6)+[CI]($^3 P_2-^3P_1$), [CI]($^3 P_1-^3P_0$) observations for a subset of the sample. The $\gtrsim 46 \%$ of our galaxies have a compact molecular gas reservoir, and lie below the optical disks mass-size relation. Compact galaxies on and above the main sequence have higher CO excitation and star formation efficiency than galaxies with extended molecular gas reservoirs, as traced by CO(5-4)/CO(2-1) and CO(2-1)/$L_{\rm IR, SF}$ ratios. Average CO+[CI] spectral line energy distributions indicate higher excitation in compacts relative to extended sources. Using CO(2-1) and dust masses as molecular gas mass tracers, and conversion factors tailored to their ISM conditions, we measure lower gas fractions in compact main-sequence galaxies compared to extended sources. We suggest that the sub-millimetre compactness, defined as the ratio between the molecular gas and the stellar size, is an unavoidable information to be used with the main sequence offset to describe the ISM properties of galaxies, at least above $M_{\star} \geqslant 10^{10.6}$ M$_{\odot}$, where our observations fully probe the main sequence scatter. Our results are consistent with mergers driving the gas in the nuclear regions, enhancing the CO excitation and star formation efficiency. Compact main-sequence galaxies are consistent with being an early post-starburst population following a merger-driven starburst episode, stressing the important role of mergers in the evolution of massive galaxies., Comment: Accepted on MNRAS on October 11th 2021. Added Figure 10, 11, 12 following referee's comments. Abstract slightly modified to fit to arXiv's requirements

The infrared-radio correlation (IRRC) underpins many commonly used radio luminosity-star formation rate (SFR) calibrations. In preparation for the new generation of radio surveys we revisit the IRRC of low-$z$ galaxies by (a) drawing on the best currently available IR and 1.4 GHz radio photometry, plus ancillary data over the widest possible area, and (b) carefully assessing potential systematics. We compile a catalogue of $\sim$9,500 z $<$ 0.2 galaxies and derive their 1.4 GHz radio ($L_{\mathrm{1.4}}$), total IR, and monochromatic IR luminosities in up to seven bands, allowing us to parameterize the wavelength-dependence of monochromatic IRRCs from 22-500 $\mu$m. For the first time for low-$z$ samples, we quantify how poorly matched IR and radio survey depths bias measured median IR/radio ratios, $\overline{q}_{\mathrm{TIR}}$, and discuss the level of biasing expected for low-z IRRC studies in ASKAP/MeerKAT fields. For our subset of $\sim$2,000 high-confidence star-forming galaxies we find a median $\overline{q}_{\mathrm{TIR}}$ of 2.54 (scatter: 0.17 dex). We show that $\overline{q}_{\mathrm{TIR}}$ correlates with $L_{\mathrm{1.4}}$, implying a non-linear IRRC with slope 1.11$\pm$0.01. Our new $L_{\mathrm{1.4}}$-SFR calibration, which incorporates this non-linearity, reproduces SFRs from panchromatic SED fits substantially better than previous IRRC-based recipes. Finally, we match the evolutionary slope of recently measured $\overline{q}_{\mathrm{TIR}}$-redshift trends without having to invoke redshift evolution of the IRRC. In this framework, the redshift evolution of $\overline{q}_{\mathrm{TIR}}$ reported at GHz frequencies in the literature is the consequence of a partial, redshift-dependent sampling of a non-linear IRRC obeyed by low-$z$ {\it and} distant galaxies., Comment: Accepted for publication in MNRAS. 29 pages, 20 figures, 9 tables including appendices

Using the deepest 1.2 mm continuum map to date in the Hubble Ultra Deep Field obtained as part of the ALMA Spectroscopic Survey (ASPECS) large program, we measure the cosmic density of dust and implied gas (H$_{2}+$H I) mass in galaxies as a function of look-back time. We do so by stacking the contribution from all $H$-band selected galaxies above a given stellar mass in distinct redshift bins, $\rho_{\rm dust}(M_\ast>M,z)$ and $\rho_{\rm gas}(M_\ast>M,z)$. At all redshifts, $\rho_{\rm dust}(M_\ast>M,z)$ and $\rho_{\rm gas}(M_\ast>M,z)$ grow rapidly as $M$ decreases down to $10^{10}\,M_\odot$, but this growth slows down towards lower stellar masses. This flattening implies that at our stellar mass-completeness limits ($10^8\,M_\odot$ and $10^{8.9}\,M_\odot$ at $z\sim0.4$ and $z\sim3$), both quantities converge towards the total cosmic dust and gas mass densities in galaxies. The cosmic dust and gas mass densities increase at early cosmic time, peak around $z\sim2$, and decrease by a factor $\sim4$ and 7, compared to the density of dust and molecular gas in the local universe, respectively. The contribution of quiescent galaxies -- i.e., with little on-going star-formation -- to the cosmic dust and gas mass densities is minor ($\lesssim10\%$). The redshift evolution of the cosmic gas mass density resembles that of the star-formation rate density, as previously found by CO-based measurements. This confirms that galaxies have relatively constant star-formation efficiencies (within a factor $\sim2$) across cosmic time. Our results also imply that by $z\sim0$, a large fraction ($\sim90\%$) of dust formed in galaxies across cosmic time has been destroyed or ejected to the intergalactic medium., Comment: 32 pages, 13 figures, Accepted for publication in ApJ

We use the results from the ALMA large program ASPECS, the spectroscopic survey in the Hubble Ultra Deep Field (HUDF), to constrain CO luminosity functions of galaxies and the resulting redshift evolution of $\rho$(H$_2$). The broad frequency range covered enables us to identify CO emission lines of different rotational transitions in the HUDF at $z>1$. We find strong evidence that the CO luminosity function evolves with redshift, with the knee of the CO luminosity function decreasing in luminosity by an order of magnitude from $\sim$2 to the local universe. Based on Schechter fits, we estimate that our observations recover the majority (up to $\sim$90%, depending on the assumptions on the faint end) of the total cosmic CO luminosity at $z$=1.0-3.1. After correcting for CO excitation, and adopting a Galactic CO-to-H$_2$ conversion factor, we constrain the evolution of the cosmic molecular gas density $\rho$(H$_2$): this cosmic gas density peaks at $z\sim1.5$ and drops by factor of $6.5_{-1.4}^{+1.8}$ to the value measured locally. The observed evolution in $\rho$(H$_2$) therefore closely matches the evolution of the cosmic star formation rate density $\rho_{\rm SFR}$. We verify the robustness of our result with respect to assumptions on source inclusion and/or CO excitation. As the cosmic star formation history can be expressed as the product of the star formation efficiency and the cosmic density of molecular gas, the similar evolution of $\rho$(H$_2$) and $\rho_{\rm SFR}$ leaves only little room for a significant evolution of the average star formation efficiency in galaxies since $z\sim 3$ (85% of cosmic history)., Comment: 22 pages, 11 figures. Paper re-submitted to ApJ after addressing the first round of comments by the referee

Galaxy interactions are often accompanied by an enhanced star formation rate (SFR). Since molecular gas is essential for star formation, it is vital to establish whether, and by how much, galaxy interactions affect the molecular gas properties. We investigate the effect of interactions on global molecular gas properties by studying a sample of 58 galaxies in pairs and 154 control galaxies. Molecular gas properties are determined from observations with the JCMT, PMO, CSO telescopes, and supplemented with data from the xCOLD GASS and JINGLE surveys at $^{12}$CO(1-0) and $^{12}$CO(2-1). The SFR, gas mass ($M_\mathrm{H_{2}}$), and gas fraction ($f_{gas}$) are all enhanced in galaxies in pairs by $\sim$ 2.5 times compared to the controls matched in redshift, mass, and effective radius, while the enhancement of star formation efficiency (SFE $\equiv$ SFR/$M_{H_{2}}$) is less than a factor of 2. We also find that the enhancements in SFR, $M_{H_{2}}$ and $f_{gas}$ increase with decreasing pair separation and are larger in systems with smaller stellar mass ratio. Conversely, the SFE is only enhanced in close pairs (separation $<$ 20 kpc) and equal-mass systems; therefore most galaxies in pairs lie in the same parameter space on the SFR-$M_{H_{2}}$ plane as controls. This is the first time that the dependence of molecular gas properties on merger configurations is probed statistically with a relatively large sample and with a carefully-selected control sample for individual galaxies. We conclude that galaxy interactions do modify the molecular gas properties, although the strength of the effect is merger configuration dependent., Comment: 26 pages, 11 figures, accepted for publication in ApJ

Deep continuum surveys at radio wavelengths are able to cover large areas, yield high angular resolution, and do not suffer from dust extinction, thus providing a robust way to measure the star formation history of the universe. However, with the current sensitivities of existing radio telescopes, it remains challenging to detect galaxies that dominate the cosmic star formation history even with extremely long integrations. With the ngVLA, a new portion of parameter space will be opened up for radio continuum surveys: deep (~200 nJy/bm), large-area (~1 deg^2), sub-arcsecond surveys at high frequencies (~8 GHz), where the observed radio emission from high-z galaxies should be dominated by free-free emission, providing a robust measurement for the star formation history of the universe. By being able to image the star formation activity, unbiased by dust, for a large, homogeneous sample of galaxies with a wide range of luminosities into the epoch of reionization, such surveys with the ngVLA will be highly complementary to those conducted by JWST, which will only be sensitive to un-obscured star formation., Comment: 6 pages, 2 figures. To be published in the ASP Monograph Series, "Science with a Next-Generation VLA", ed. E. J. Murphy (ASP, San Francisco, CA)

Radio continuum observations have proven to be a workhorse in our understanding of the star formation process (i.e., stellar birth and death) from galaxies both in the nearby universe and out to the highest redshifts. In this article we focus on how the ngVLA will transform our understanding of star formation by enabling one to map and decompose the radio continuum emission from large, heterogeneous samples of nearby galaxies on $\gtrsim 10$\,pc scales to conduct a proper accounting of the energetic processes powering it. At the discussed sensitivity and angular resolution, the ngVLA will simultaneously be able to create maps of current star formation activity at $\sim$100\,pc scales, as well as detect and characterize (e.g., size, spectral shape, density, etc.) discrete H{\sc ii} regions and supernova remnants on 10\,pc scales in galaxies out to the distance of the Virgo cluster. Their properties can then be used to see how they relate to the local and global ISM and star formation conditions. Such investigations are essential for understanding the astrophysics of high-$z$ measurements of galaxies, allowing for proper modeling of galaxy formation and evolution., Comment: 9 pages, 3 figures; To be published in the ASP Monograph Series, "Science with a Next-Generation VLA", ed. E. J. Murphy (ASP, San Francisco, CA). arXiv admin note: substantial text overlap with arXiv:1510.06431

In this work we update the L-Galaxies semi-analytic model (SAM) to better follow the physical processes responsible for the growth of bulges via disc instabilities (leading to pseudo-bulges) and mergers (leading to classical bulges). We address the former by considering the contribution of both stellar and gaseous discs in the stability of the galaxy, and we update the latter by including dissipation of energy in gas-rich mergers. Furthermore, we introduce angular momentum losses during cooling and find that an accurate match to the observed correlation between stellar disc scale length and mass at z ~ 0.0 requires that the gas loses 20% of its initial specific angular momentum to the corresponding dark matter halo during the formation of the cold gas disc. We reproduce the observed trends between the stellar mass and specific angular momentum for both disc- and bulge-dominated galaxies, with the former rotating faster than the latter of the same mass. We conclude that a two-component instability recipe provides a morphologically diverse galaxy sample which matches the observed fractional breakdown of galaxies into different morphological types. This recipe also enables us to obtain an excellent fit to the morphology-mass relation and stellar mass function of different galactic types. Finally, we find that energy dissipation during mergers reduces the merger remnant sizes and allows us to match the observed mass-size relation for bulge-dominated systems.

Using infrared data from the Herschel Space Observatory and Karl G. Jansky Very Large Array (VLA) 3 GHz observations in the COSMOS field, we investigate the redshift evolution of the infrared-radio correlation (IRRC) for star-forming galaxies (SFGs) we classify as either spheroid- or disc-dominated based on their morphology. The sample predominantly consists of disc galaxies with stellar mass ${\gtrsim}10^{10}\,M_{\odot}$, and residing on the star-forming main sequence (MS). After the removal of AGN using standard approaches, we observe a significant difference between the redshift-evolution of the median IR/radio ratio $\overline{q}_{\mathrm{TIR}}$ of (i) a sample of ellipticals, plus discs with a substantial bulge component (`spheroid-dominated' SFGs) and, (ii) virtually pure discs and irregular systems (`disc-dominated' SFGs). The spheroid-dominated population follows a declining $\overline{q}_{\mathrm{TIR}}$ vs. $z$ trend similar to that measured in recent evolutionary studies of the IRRC. However, for disc-dominated galaxies, where radio and IR emission should be linked to star formation in the most straightforward way, we measure very little change in $\overline{q}_{\mathrm{TIR}}$. This suggests that low-redshift calibrations of radio emission as an SFR-tracer may remain valid out to at least $z\,{\simeq}\,1\,{-}\,1.5$ for pure star-forming systems. We find that the different redshift-evolution of $q_{\rm TIR}$ for the spheroid- and disc-dominated sample is mainly due to an increasing radio excess for spheroid-dominated galaxies at $z\,{\gtrsim}\,$0.8, hinting at some residual AGN activity in these systems. This finding demonstrates that in the absence of AGN the IRRC is independent of redshift, and that radio observations can therefore be used to estimate SFRs at all redshifts for genuinely star-forming galaxies., Comment: Resubmitted to MNRAS after addressing reviewer feedback

We present results on the dust attenuation of galaxies at redshift $\sim3-6$ by studying the relationship between the UV spectral slope ($\beta_{\rm UV}$) and the infrared excess (IRX; $L_{\rm IR}$/$L_{\rm UV}$) using ALMA far-infrared continuum observations. Our study is based on a sample of 67 massive, star-forming galaxies with a median mass of $M_{\ast}\sim 10^{10.7}\,M_{\rm \odot}$ spanning a redshift range $z=2.6-3.7$ (median $z=3.2$) that were observed with ALMA at $\lambda_{rest}=300\,{\rm \mu m}$. Both the individual ALMA detections (41 sources) and stacks including all galaxies show the IRX-$\beta_{\rm UV}$ relationship at $z\sim3$ is mostly consistent with that of local starburst galaxies on average. However, we find evidence for a large dispersion around the mean relationship by up to $\pm0.5$ dex. Nevertheless, the locally calibrated dust correction factors based on the IRX-$\beta_{\rm UV}$ relation are on average applicable to main-sequence $z\sim3$ galaxies. This does not appear to be the case at even higher redshifts, however. Using public ALMA observations of $z\sim4-6$ galaxies we find evidence for a significant evolution in the IRX-$\beta_{\rm UV}$ and the IRX-$M_{\ast}$ relations beyond $z\sim3$ toward lower IRX values. We discuss several caveats that could affect these results, including the assumed dust temperature. ALMA observations of larger $z>3$ galaxy samples will be required to confirm this intriguing redshift evolution., Comment: 8 pages, 4 figures, submitted to MNRAS; v2 - Small changes to text

HE0450-2958, an interacting quasar-starburst galaxy pair at $z$ = 0.285, is one of the best known examples of strong star formation activity in the presence of a quasar-driven jet. We present new multi-band JVLA-imaging covering 1 to 6 GHz and reaching an angular resolution of up to $0{}_{{}^.}^{"}6$ (a 6-fold improvement over existing radio data). We confirm the previous detection of a spatially extended radio component around the quasar indicating that there is on-going star formation activity in the quasar host galaxy. For the first time, we directly detect a jet-like bipolar outflow from the quasar aligned with its companion star-forming galaxy (SFG) and several blobs of ionized gas in its vicinity identified in previous studies. Within the companion SFG we find evidence for a flattening of the synchrotron spectral index towards the point of intersection with the jet axis, further suggesting that the outflow may actually be impacting its interstellar medium (ISM). We discuss two possible mechanisms that could have triggered the starburst in the companion SFG: a wet-dry merger with the quasar and jet-induced star formation. While triggering through interaction-driven gas dynamics cannot be excluded with current data, our new observations make HE0450-2958 a strong candidate for jet-induced star formation, and one of the rare links between local systems (like Minkowski's Object or Centaurus A) and the high-z regime where radio-optical alignments suggest that this phenomenon could be more common., Comment: 11 pages including 5 figures, published in MNRAS

We investigate the multiplicity of extragalactic sources detected by the Herschel Space Observatory in the COSMOS field. Using 3.6- and 24-$\mu$m catalogues, in conjunction with 250-$\mu$m data from Herschel, we seek to determine if a significant fraction of Herschel sources are composed of multiple components emitting at 250 $\mu$m. We use the XID+ code, using Bayesian inference methods to produce probability distributions of the possible contributions to the observed 250-$\mu$m flux for each potential component. The fraction of Herschel flux assigned to the brightest component is highest for sources with total 250-$\mu$m fluxes < 45 mJy; however, the flux in the brightest component is still highest in the brightest Herschel sources. The faintest 250-$\mu$m sources (30-45 mJy) have the majority of their flux assigned to a single bright component; the second brightest component is typically significantly weaker, and contains the remainder of the 250-$\mu$m source flux. At the highest 250-$\mu$m fluxes (45-110 mJy), the brightest and second brightest components are assigned roughly equal fluxes, and together are insufficient to reach 100 per cent of the 250-$\mu$m source flux. This indicates that additional components are required, beyond the brightest two components, to reproduce the observed flux. 95 per cent of the sources in our sample have a second component that contains more than 10 per cent of the total source flux. Particularly for the brightest Herschel sources, assigning the total flux to a single source may overestimate the flux contributed by around 150 per cent., Comment: 12 pages, 7 figures. MNRAS accepted

Gas and dust in star-forming galaxies at the peak epoch of galaxy assembly are presently the topic of intense study, but little is known about the interstellar medium (ISM) of distant, passively evolving galaxies. We report on a deep 3 mm-band search with IRAM/PdBI for molecular (H$_2$) gas in a massive ($M_{\star}{\sim}6{\times}10^{11}M_{\odot}$) elliptical galaxy at z=1.4277, the first observation of this kind ever attempted. We place a 3$\sigma$ upper limit of 0.32 Jy km/s on the flux of the CO($J$=2$\rightarrow$1) line or $L'_{\rm CO}$$<$8.8$\times$10$^{9}$ K km/s pc$^2$, assuming a disk-like CO-morphology and a circular velocity scaling with the stellar velocity dispersion as in local early-type galaxies (ETGs). This translates to an H$_2$ mass of $<$3.9$\times$10$^{10}$($\alpha_{\rm CO}$/4.4)$M_{\odot}$ or a gas fraction of $\lesssim$6% assuming a Salpeter initial mass function (IMF) and an ISM dominated by H$_2$, as observed in many local, high-mass ellipticals. This low value approaches that of local ETGs, suggesting that the low star formation activity in massive, high-z passive galaxies reflects a true dearth of gas and a lesser role for inhibitive mechanisms like morphological quenching., Comment: 6 pages, 3 figures & 1 table. Accepted for publication in ApJL

In this chapter, we highlight a number of science investigations that are enabled by the inclusion of Band~5 ($4.6-13.8$ GHz) for SKA1-MID science operations, while focusing on the astrophysics of star formation over cosmic time. For studying the detailed astrophysics of star formation at high-redshift, surveys at frequencies $\gtrsim$10 GHz have the distinct advantage over traditional $\sim$1.4 GHz surveys as they are able to yield higher angular resolution imaging while probing higher rest frame frequencies of galaxies with increasing redshift, where emission of star-forming galaxies becomes dominated by thermal (free-free) radiation. In doing so, surveys carried out at $\gtrsim$10 GHz provide a robust, dust-unbiased measurement of the massive star formation rate by being highly sensitive to the number of ionizing photons that are produced. To access this powerful star formation rate diagnostic requires that Band~5 be available for SKA1-MID. We additionally present a detailed science case for frequency coverage extending up to 30 GHz during full SKA2 operations, as this allows for highly diverse science while additionally providing contiguous frequency coverage between the SKA and ALMA, which will likely be the two most powerful interferometers for the coming decades. To enable this synergy, it is crucial that the dish design of the SKA be flexible enough to include the possibility of being fit with receivers operating up to 30 GHz., Comment: 15 pages, 3 figures, to appear as part of 'Continuum Science' in Proceedings "Advancing Astrophysics with the Square Kilometre Array" PoS(AASKA14)085; "Kilometer" replaced by "Kilometre" in updated version

We present an analysis of the deepest Herschel images in four major extragalactic fields GOODS-North, GOODS-South, UDS and COSMOS obtained within the GOODS-Herschel and CANDELS-Herschel key programs. The picture provided by 10497 individual far-infrared detections is supplemented by the stacking analysis of a mass-complete sample of 62361 star-forming galaxies from the CANDELS-HST H band-selected catalogs and from two deep ground-based Ks band-selected catalogs in the GOODS-North and the COSMOS-wide fields, in order to obtain one of the most accurate and unbiased understanding to date of the stellar mass growth over the cosmic history. We show, for the first time, that stacking also provides a powerful tool to determine the dispersion of a physical correlation and describe our method called "scatter stacking" that may be easily generalized to other experiments. We demonstrate that galaxies of all masses from z=4 to 0 follow a universal scaling law, the so-called main sequence of star-forming galaxies. We find a universal close-to-linear slope of the logSFR-logM* relation with evidence for a flattening of the main sequence at high masses (log(M*/Msun) > 10.5) that becomes less prominent with increasing redshift and almost vanishes by z~2. This flattening may be due to the parallel stellar growth of quiescent bulges in star-forming galaxies. Within the main sequence, we measure a non varying SFR dispersion of 0.3 dex. The specific SFR (sSFR=SFR/M*) of star-forming galaxies is found to continuously increase from z=0 to 4. Finally we discuss the implications of our findings on the cosmic SFR history and show that more than 2/3 of present-day stars must have formed in a regime dominated by the main sequence mode. As a consequence we conclude that, although omnipresent in the distant Universe, galaxy mergers had little impact in shaping the global star formation history over the last 12.5 Gyr., Comment: 30 pages, 26 figures, accepted for publication in A&A

Limited studies have been performed on the radio-loud fraction in X-ray selected type 1 AGN samples. The consistency between various radio-loudness definitions also needs to be checked. We measure the radio-loudness of the 407 type 1 AGNs in the XMM-COSMOS quasar sample using nine criteria from the literature (six defined in the rest-frame and three defined in the observed frame): $R_L=\log(L_{5GHz}/L_B)$, $q_{24}=\log(L_{24\mu m}/L_{1.4GHz})$, $R_{uv}=\log(L_{5GHz}/L_{2500\AA})$, $R_{i}=\log(L_{1.4GHz}/L_i)$, $R_X=\log(\nu L_{\nu}(5GHz)/L_X)$, $P_{5GHz}=\log(P_{5GHz}(W/Hz/Sr))$, $R_{L,obs}=\log(f_{1.4GHz}/f_B)$ (observed frame), $R_{i,obs}=\log(f_{1.4GHz}/f_i)$ (observed frame), and $q_{24, obs}=\log(f_{24\mu m}/f_{1.4GHz})$ (observed frame). Using any single criterion defined in the rest-frame, we find a low radio-loud fraction of $\lesssim 5\%$ in the XMM-COSMOS type 1 AGN sample, except for $R_{uv}$. Requiring that any two criteria agree reduces the radio-loud fraction to $\lesssim 2\%$ for about 3/4 of the cases. The low radio-loud fraction cannot be simply explained by the contribution of the host galaxy luminosity and reddening. The $P_{5GHz}=\log(P_{5GHz}(W/Hz/Sr))$ gives the smallest radio-loud fraction. Two of the three radio-loud fractions from the criteria defined in the observed frame without k-correction ($R_{L,obs}$ and $R_{i,obs}$) are much larger than the radio-loud fractions from other criteria., Comment: 12 pages, 7 figures, MNRAS submitted

We use deep panchromatic datasets in the GOODS-N field, from GALEX to the deepest Herschel far-infrared and VLA radio continuum imaging, to explore, using mass-complete samples, the evolution of the star formation activity and dust attenuation of star-forming galaxies to z~4. Our main results can be summarized as follows: i) the slope of the SFR-M correlation is consistent with being constant, and equal to ~0.8 at least up to z~1.5, while its normalization keeps increasing with redshift; ii) for the first time here we are able to explore the FIR-radio correlation for a mass-selected sample of star-forming galaxies: the correlation does not evolve up to z~4; iii) we confirm that galaxy stellar mass is a robust proxy for UV dust attenuation in star-forming galaxies, with more massive galaxies being more dust attenuated, strikingly we find that this attenuation relation evolves very weakly with redshift, the amount of dust attenuation increasing by less than 0.3 magnitudes over the redshift range [0.5-4] for a fixed stellar mass, as opposed to a tenfold increase of star formation rate; iv) the correlation between dust attenuation and the UV spectral slope evolves in redshift, with the median UV spectral slope of star-forming galaxies becoming bluer with redshift. By z~3, typical UV slopes are inconsistent, given the measured dust attenuation, with the predictions of commonly used empirical laws. Finally, building on existing results, we show that gas reddening is marginally larger (by a factor of around 1.3) than stellar reddening at all redshifts probed, and also that the amount of dust attenuation at a fixed ISM metallicity increases with redshift. We speculate that our results support evolving ISM conditions of typical star-forming galaxies such that at z~1.5 Main Sequence galaxies have ISM conditions getting closer to those of local starbursts., Comment: 19 pages, 12 figure, accepted for publication in ApJ

We present deep IRAM Plateau de Bure Interferometer (PdBI) observations, searching for CO-emission toward two massive, non-lensed Lyman break galaxies (LBGs) at z=3.216 and 4.058. With one low significance CO detection (3.5 sigma) and one sensitive upper limit, we find that the CO lines are >~ 3-4 times weaker than expected based on the relation between IR and CO luminosities followed by similarly, massive galaxies at z=0-2.5. This is consistent with a scenario in which these galaxies have low metallicity, causing an increased CO-to-H_2 conversion factor, i.e., weaker CO-emission for a given molecular (H_2) mass. The required metallicities at z>3 are lower than predicted by the fundamental metallicity relation (FMR) at these redshifts, consistent with independent evidence. Unless our galaxies are atypical in this respect, detecting molecular gas in normal galaxies at z>3 may thus remain challenging even with ALMA., Comment: 9 pages, 3 figures, accepted for publication in Astrophysical Journal Letters

Star-forming galaxies (SFGs) display a continuous distribution of specific star formation rates (sSFR) which can be approximated by the superposition of two log-normal distributions. The 1st of these encompasses the main sequence (MS) of SFGs, the 2nd one a rarer population of starbursts (SB). We show that the sSFR-distribution of SBs can be regarded as the result of a physical process (plausibly merging) taking the mathematical form of a log-normal boosting kernel and enhancing star formation activity. We explore the utility of splitting the star-forming population into MS and SB galaxies - an approach we term "2-Star Formation Mode" (2-SFM) framework - for understanding their molecular gas properties. Variations of star formation efficiency (SFE) and gas fractions among SFGs take a simple, redshift-independent form, once these quantities are normalized to the value of an average MS galaxy. The change in SFE for galaxies undergoing a starburst event scales supra-linearly with the SFR-increase, as expected for mergers. This implies a continuous distribution of galaxies in the Schmidt-Kennicutt plane that separates more clearly into loci for SBs and normal galaxies than observed in SFR vs. M* space. SBs with the largest deviations (>10-fold) from the MS, like many local ULIRGs, are not average SBs, but even rarer events with progenitors having larger gas fractions than typical MS galaxies. We statistically infer that the gas fractions of typical SBs decrease by a factor of 2 to 3 with respect to their direct MS progenitors, as expected to occur in short-term SFR-boosts during which internal gas reservoirs are drained more quickly than gas is accreted from the cosmic web. We predict variations of the CO-to-H2 conversion factor in the SFR-M* plane and provide evidence that the higher sSFR of distant galaxies is a direct consequence of larger gas fraction in these systems., Comment: 37 pages, 18 figures. Published in ApJ, 793, 19

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 for 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...) is available online at http://irfu.cea.fr/Sap/Phocea/Page/index.php?id=537., Comment: 7 pages, 4 figures, 1 table, ApJL in press

We present a parameter-less approach to predict the shape of the infrared (IR) luminosity function (LF) at redshifts z < 2. It requires no tuning and relies on only three observables: (1) the redshift evolution of the stellar mass function for star-forming galaxies, (2) the evolution of the specific star formation rate (sSFR) of main-sequence galaxies, and (3) the double-Gaussian decomposition of the sSFR-distribution at fixed stellar mass into a contribution (assumed redshift- and mass-invariant) from main-sequence and starburst activity. This self-consistent and simple framework provides a powerful tool for predicting cosmological observables: observed IR LFs are successfully matched at all z < 2, suggesting a constant or only weakly redshift-dependent contribution (8-14%) of starbursts to the star formation rate density. We separate the contributions of main-sequence and starburst activity to the global IR LF at all redshifts. The luminosity threshold above which the starburst component dominates the IR LF rises from log(LIR/Lsun) = 11.4 to 12.8 over 0 < z < 2, reflecting our assumed (1+z)^2.8-evolution of sSFR in main-sequence galaxies., Comment: 7 pages, 4 figures & 1 table. Accepted for publication in ApJL. Minor typos corrected in v2 following receipt of proofs

We explore the evolution of the specific star formation rate (SSFR) for 3.6um-selected galaxies of different M_* in the COSMOS field. The average SFR for sub-sets of these galaxies is estimated with stacked 1.4GHz radio continuum emission. We separately consider the total sample and a subset of galaxies (SF) that shows evidence for substantive recent star formation in the rest-frame optical SED. At 0.22, at least above 4x10^10M_Sun where our conclusions are most robust. We find a tight correlation with power-law dependence, SSFR (M_*)^beta, between SSFR and M_* at all z. It tends to flatten below ~10^10M_Sun if quiescent galaxies are included; if they are excluded a shallow index beta_SFG -0.4 fits the correlation. On average, higher M_* objects always have lower SSFRs, also among SF galaxies. At z>1.5 there is tentative evidence for an upper SSFR-limit that an average galaxy cannot exceed. It is suggested by a flattening of the SSFR-M_* relation (also for SF sources), but affects massive (>10^10M_Sun) galaxies only at the highest z. Below z=1.5 there thus is no direct evidence that galaxies of higher M_* experience a more rapid waning of their SSFR than lower M_* SF systems. In this sense, the data rule out any strong 'downsizing'. We combine our results with recent measurements of the galaxy (stellar) mass function in order to determine the characteristic mass of a SF galaxy (M_*=10^(10.6\pm0.4)M_Sun). In this sense, too, there is no 'downsizing'. Our analysis constitutes the most extensive SFR density determination with a single technique to z=3. Recent Herschel results are consistent with our results, but rely on far smaller samples., Comment: 37 pages, 14 figures, 7 tables; accepted for publication in the Astrophysical Journal; High resolution versions of all figures available at www.mpia-hd.mpg.de/homes/karim/research.html

Previous observational studies of the infrared (IR)-radio relation out to high redshift employed any detectable star forming systems at a given redshift within the restricted area of cosmological survey fields. Consequently, the evolution inferred relies on a comparison between the average IR/radio properties of (i) very IR-luminous high-z sources and (ii) more heterogeneous low(er)-z samples that often lack the strongest IR emitters. In this report we consider populations of objects with comparable luminosities over the last 10 Gyr by taking advantage of deep IR (esp. Spitzer 24 micron) and VLA 1.4 GHz observations of the COSMOS field. Consistent with recent model predictions, both Ultra Luminous Infrared Galaxies (ULIRGs) and galaxies on the bright end of the evolving IR luminosity function do not display any change in their average IR/radio ratios out to z~2 when corrected for bias. Uncorrected data suggested ~0.3 dex of positive evolution., Comment: 9 pages, 4 figures. Accepted for publication in ApJL.

We compare the surface brightness-inclination relation for a sample of COSMOS pure disk galaxies at z~0.7 with an artificially redshifted sample of SDSS disks well matched to the COSMOS sample in terms of rest-frame photometry and morphology, as well as their selection and analysis. The offset between the average surface brightness of face-on and edge-on disks in the redshifted SDSS sample matches that predicted by measurements of the optical depth of galactic disks in the nearby universe. In contrast, large disks at z~0.7 have a virtually flat surface brightness-inclination relation, suggesting that they are more opaque than their local counterparts. This could be explained by either an increased amount of optically thick material in disks at higher redshift, or a different spatial distribution of the dust., Comment: 7 pages, 3 figures. Accepted for publication in ApJL.

VLA 1.4 GHz (rms noise ~0.012 mJy) and MIPS 24 and 70 micron (rms noise ~0.02 and ~1.7 mJy, respectively) observations covering the 2 square degree COSMOS field are combined with an extensive multi-wavelength data set to study the evolution of the IR-radio relation at intermediate and high redshift. With ~4500 sources -- of which ~30% have spectroscopic redshifts -- the current sample is significantly larger than previous ones used for the same purpose. Both monochromatic IR/radio flux ratios (q24 & q70), as well as the ratio of the total IR and the 1.4 GHz luminosity (qTIR) are used as indicators for the IR/radio properties of star forming galaxies and AGN. Using a sample jointly selected at IR and radio wavelengths in order to reduce selection biases, we provide firm support for previous findings that the IR-radio relation remains unchanged out to at least z~1.4. Moreover, based on data from ~150 objects we also find that the local relation likely still holds at 2.5