• ### Spatially resolved cold molecular outflows in ULIRGs(1805.03667)

May 9, 2018 astro-ph.CO, astro-ph.GA
We present new CO(2-1) observations of 3 low-z (~350 Mpc) ULIRG systems (6 nuclei) observed with ALMA at high-spatial resolution (~500 pc). We detect massive cold molecular gas outflows in 5 out of 6 nuclei (0.3-5)x10^8 Msun. These outflows are spatially resolved with deprojected radii of 0.25-1 kpc although high-velocity molecular gas is detected up to ~0.5-1.8 kpc (1-6 kpc deprojected). The mass outflow rates are 12-400 Msun/yr and the inclination corrected average velocity of the outflowing gas 350-550 km/s (v_max = 500-900 km/s). The origin of these outflows can be explained by the nuclear starbursts although the contribution of an obscured AGN can not be completely ruled out. The position angle (PA) of the outflowing gas along the kinematic minor axis of the nuclear molecular disk suggests that the outflow axis is perpendicular to the disk for three of these outflows. Only in one case, the outflow PA is clearly not along the kinematic minor axis. The outflow depletion times are 15-80 Myr which are slightly shorter than the star-formation (SF) depletion times (30-80 Myr). However, we estimate that only 15-30% of the outflowing gas will escape the gravitational potential of the nucleus. The majority of the outflowing gas will return to the disk after 5-10 Myr and become available to form new stars. Therefore, these outflows will not likely quench the nuclear starbursts. These outflows would be consistent with being driven by radiation pressure (momentum-driven) only if the coupling between radiation and dust increases with increasing SF rates. This can be achieved if the dust optical depth is higher in objects with higher SF. The relatively small sizes (<1 kpc) and dynamical times (<3 Myr) of the cold molecular outflows suggests that molecular gas cannot survive longer in the outflow environment or that it cannot form efficiently beyond these distances or times. (Abridged)
• ### Molecular gas masses of gamma-ray burst host galaxies(1804.06492)

April 17, 2018 astro-ph.GA, astro-ph.HE
The objectives of this paper are to analyse molecular gas properties of the first substantial sample of GRB hosts and test whether they are deficient in molecular gas. We obtained CO(2-1) observations of seven GRB hosts with the APEX and IRAM30m telescopes. We analyse these data together with all other hosts with previous CO observations. We obtained detections for three GRB hosts (980425, 080207 and 111005A) and upper limits for the remaining four (031203, 060505, 060814, 100316D). In our entire sample of twelve CO-observed GRB hosts, three are clearly deficient in molecular gas, even taking into account their metallicity (980425, 060814, and 080517). Four others are close to the best fit-line for other star-forming galaxies on the SFR-MH2 plot (051022, 060505, 080207, and 100316D). One host is clearly molecule-rich (111005A). Finally, for four GRB hosts the data is not deep enough to judge whether they are molecule-deficient (000418, 030329, 031203, 090423). The median value of the molecular gas depletion time, MH2/SFR, of GRB hosts is ~0.3 dex below that of other star-forming galaxies, but this result has low statistical significance. A Kolmogorov-Smirnov test performed on MH2/SFR shows only ~2 sigma difference between GRB hosts and other galaxies. This difference can partially be explained by metallicity effects, since the significance decreases to ~1 sigma for MH2/SFR vs. metallicity. Thus the sample of GRB hosts has consistent molecular properties to other galaxies, and can be treated as representative star-forming galaxies. Given the concentration of atomic gas recently found close to GRB and supernova sites, indicating recent gas inflow, our results imply that such inflow does not enhance the SFRs significantly, or that atomic gas converts efficiently into the molecular phase.
• ### The dust and [CII] morphologies of redshift ~4.5 sub-millimeter galaxies at ~200pc resolution: The absence of large clumps in the interstellar medium of high-redshift galaxies(1804.03663)

April 10, 2018 astro-ph.CO, astro-ph.GA
We present deep high resolution (0.03", 200pc) ALMA Band 7 observations covering the dust continuum and [CII] $\lambda157.7\mu$m emission in four $z\sim4.4-4.8$ sub-millimeter galaxies (SMGs) selected from the ALESS and AS2UDS surveys. The data show that the rest-frame 160$\mu$m (observed 345 GHz) dust emission is consistent with smooth morphologies on kpc scales for three of the sources. One source, UDS47.0, displays apparent substructure but this is also consistent with a smooth morphology, as indicated by simulations showing that smooth exponential disks can appear clumpy when observed at high angular resolution (0.03") and depth of these observations ($\sigma_{345\text{GHz}} \sim27-47\mu$Jy beam$^{-1}$). The four SMGs are bright [CII] emitters, and we extract [CII] spectra from the high resolution data, and recover $\sim20-100$% of the [CII] flux and $\sim40-80$% of the dust continuum emission, compared to the previous lower resolution observations. When tapered to 0.2" resolution our maps recover $\sim80-100$% of the continuum emission, indicating that $\sim60$% of the emission is resolved out on $\sim200$pc scales. We find that the [CII] emission in high-redshift galaxies is more spatially extended than the rest-frame 160$\mu$m dust continuum by a factor of $1.6\pm0.4$. By considering the $L_{\text{[CII]}}$/$L_{\text{FIR}}$ ratio as a function of the star-formation rate surface density ($\Sigma_{\text{SFR}}$) we revisit the [CII] deficit, and suggest that the decline in the $L_{\text{[CII]}}$/$L_{\text{FIR}}$ ratio as a function of $\Sigma_{\text{SFR}}$ is consistent with local processes. We also explore the physical drivers that may be responsible for these trends and can give rise to the properties found in the densest regions of SMGs.
• ### Molecular line emission in NGC 4945, imaged with ALMA(1802.09852)

Feb. 27, 2018 astro-ph.GA
NGC 4945 is one of the nearest (~3.8 Mpc; 1" ~ 19 pc) starburst galaxies. ALMA band 3 (3--4\,mm) observations of HCN, HCO+, CS, C3H2, SiO, HCO, and CH3C2H were carried out with ~2" resolution. The lines reveal a rotating nuclear disk of projected size 10" x 2" with position angle ~45 deg, inclination ~75 deg and an unresolved bright central core of size <2.5". The continuum source (mostly free-free radiation) is more compact than the nuclear disk by a linear factor of two but shows the same position angle and is centered 0.39" +_ 0.14" northeast of the nuclear accretion disk defined by H2O maser emission. Outside the nuclear disk, both HCN and CS delineate molecular arms on opposite sides of the dynamical center. These are connected by a (deprojected) 0.6 kpc sized molecular bridge, likely a dense gaseous bar seen almost ends-on, shifting gas from the front and back side into the nuclear disk. Modeling this nuclear disk located farther inside <100 pc) with tilted rings indicates a coplanar outflow reaching a characteristic deprojectd velocity of ~50 km/s. All our molecular lines, with the notable exception of CH3C2H, show significant absorption near the systemic velocity (~571 km/s), within a range of ~500-660 km/s. Apparently, only molecular transitions with low critical H2-density do not show absorption. The velocity field of the nuclear disk, derived from CH3C2H, provides evidence for rigid rotation in the inner few arcseconds and a dynamical mass of M = (2.1+_0.2) x 10^8 Mo inside a galactocentric radius of 2.45", with a significantly flattened rotation curve farther out. Velocity integrated line intensity maps with most pronounced absorption show molecular peak positions up to 1.5" southwest of the continuum peak, presumably due to absorption, which appears to be most severe slightly northeast of the nuclear maser disk.
• ### Herschel Spectroscopy of the Taffy Galaxies (UGC 12914/12915 = VV 254): Enhanced [C II] emission in the collisionally-formed bridge(1801.10230)

Jan. 30, 2018 astro-ph.GA
Using the PACS and SPIRE spectrometers on-board Herschel, we obtained observations of the Taffy galaxies (UGC 12914/12915) and bridge. The Taffy system is believed to be the result of a face-on collision between two gas-rich galaxies, in which the stellar disks passed through each other, but the gas was dispersed into a massive H I and molecular bridge between them. Emission is detected and mapped in both galaxies and the bridge in the [C II]157.7 $\mu$m and [O I]63.2 $\mu$m fine-structure lines. Additionally, SPIRE FTS spectroscopy detects the [C I] $^3$P$_2$$\rightarrow$$^3$P$_1$(809.3 GHz) and [C I] $^3$P$_1$$\rightarrow$$3$P$_0$(492.2 GHz) neutral carbon lines, and weakly detects high-J CO transitions in the bridge. These results indicate that the bridge is composed of a warm multi-phase medium consistent with shock and turbulent heating. Despite low star formation rates in the bridge, the [C II] emission appears to be enhanced, reaching [C II]/FIR ratios of 3.3% in parts of the bridge. Both the [C II] and [O I] lines show broad intrinsic multi-component profiles, similar to those seen in previous CO 1-0 and H I observations. The [C II] emission shares similar line profiles with both the double-peaked H I profiles and shares a high-velocity component with single-peaked CO profiles in the bridge, suggesting that the [C II] emission originates in both the neutral and molecular phases. We show that it is feasible that a combination of turbulently heated H$_2$ and high column-density H I, resulting from the galaxy collision, is responsible for the enhanced [C II] emission.
• ### The SCUBA-2 Cosmology Legacy Survey: The EGS deep field - II. Morphological transformation and multi-wavelength properties of faint submillimetre galaxies(1801.07718)

Jan. 23, 2018 astro-ph.GA
We present a multi-wavelength analysis of galaxies selected at 450 and 850um from the deepest SCUBA-2 observations in the Extended Groth Strip (EGS) field, which have an average depth of sigma_450=1.9 and sigma_850=0.46 mJy/beam over ~70 sq. arcmin. The final sample comprises 95 sources: 56 (59%) are detected at both wavelengths, 31 (33%) are detected only at 850um, and 8 (8%) are detected only at 450um. We identify counterparts for 75% of the whole sample. The redshift distributions of the 450 and 850um samples peak at different redshifts with median values of z=1.66 +/- 0.18 and z=2.30 +/- 0.20, respectively. However, the two populations have similar IR luminosities, SFRs, and stellar masses, with mean values of (1.5 +/- 0.2) x 10^12 L_sun, (150 +/- 20) M_sun/yr, and (9.0 +/- 0.6) x 10^10 M_sun, respectively. This places most of our sources (>85%) on the high-mass end of the main-sequence' of star-forming galaxies. Exploring the IR excess vs UV-slope (IRX-beta) relation we find that the most luminous galaxies are consistent with the Meurer law, while the less luminous galaxies lie below this relation. Using the results of a two-dimensional modelling of the HST H_160-band imaging, we derive a median Sersic index of n=1.4 +0.3 -0.1 and a median half-light radius of R_1/2 = 4.8 +/ 0.4 kpc. Based on a visual-like classification in the same band, we find that the dominant component for most of the galaxies at all redshifts is a disk-like structure, although there is a transition from irregular disks to disks with a spheroidal component at z~1.4, which morphologically supports the scenario of SMGs as progenitors of massive elliptical galaxies.
• ### The Herschel-ATLAS: magnifications and physical sizes of $500\,\mu$m-selected strongly lensed galaxies(1801.01831)

Jan. 5, 2018 astro-ph.GA
We perform lens modelling and source reconstruction of Submillimeter Array (SMA) data for a sample of 12 strongly lensed galaxies selected at 500$\mu$m in the Herschel Astrophysical Terahertz Large Area Survey H-ATLAS. A previous analysis of the same dataset used a single S\ersic profile to model the light distribution of each background galaxy. Here we model the source brightness distribution with an adaptive pixel scale scheme, extended to work in the Fourier visibility space of interferometry. We also present new SMA observations for seven other candidate lensed galaxies from the H-ATLAS sample. Our derived lens model parameters are in general consistent with previous findings. However, our estimated magnification factors, ranging from 3 to 10, are lower. The discrepancies are observed in particular where the reconstructed source hints at the presence of multiple knots of emission. We define an effective radius of the reconstructed sources based on the area in the source plane where emission is detected above 5$\sigma$. We also fit the reconstructed source surface brightness with an elliptical Gaussian model. We derive a median value $r_{eff}\,\sim 1.77\,$kpc and a median Gaussian full width at half maximum $\sim1.47\,$kpc. After correction for magnification, our sources have intrinsic star formation rates SFR$\,\sim900-3500\,M_{\odot}yr^{-1}$, resulting in a median star formation rate surface density $\Sigma_{SFR}\sim132\,M_{\odot}$ yr$^{-1}$ kpc$^{-2}$ (or $\sim 218\,M_{\odot}$ yr$^{-1}$ kpc$^{-2}$ for the Gaussian fit). This is consistent with what observed for other star forming galaxies at similar redshifts, and is significantly below the Eddington limit for a radiation pressure regulated starburst.
• ### Candidate high-z proto-clusters among the Planck compact sources, as revealed by Herschel-SPIRE(1712.07141)

Dec. 19, 2017 astro-ph.GA
By determining the nature of all the Planck compact sources within 808.4 deg^2 of large Herschel surveys, we have identified 27 candidate proto-clusters of dusty star forming galaxies (DSFGs) that are at least 3{\sigma} overdense in either 250, 350 or 500 $\mu$mm sources. We find roughly half of all the Planck compact sources are resolved by Herschel into multiple discrete objects, with the other half remaining unresolved by Herschel. We find a significant difference between versions of the Planck catalogues, with earlier releases hosting a larger fraction of candidate proto-clusters and Galactic Cirrus than later releases, which we ascribe to a difference in the filters used in the creation of the three catalogues. We find a surface density of DSFG candidate proto-clusters of (3.3 $\pm$ 0.7) x 10^(-2) sources deg^(-2), in good agreement with previous similar studies. We find that a Planck colour selection of S_{857}/S_{545} < 2 works well to select candidate proto-clusters, but can miss proto-clusters at z < 2. The Herschel colours of individual candidate proto-cluster members indicate our candidate proto-clusters all likely all lie at z > 1. Our candidate proto-clusters are a factor of 5 times brighter at 353 GHz than expected from simulations, even in the most conservative estimates. Further observations are needed to confirm whether these candidate proto-clusters are physical clusters, multiple proto-clusters along the line of sight, or chance alignments of unassociated sources.
• ### An extreme proto-cluster of luminous dusty starbursts in the early Universe(1709.02809)

Sept. 20, 2017 astro-ph.GA
We report the identification of an extreme proto-cluster of galaxies in the early Universe whose core (nicknamed Distant Red Core, DRC) is formed by at least ten dusty star-forming galaxies (DSFGs), confirmed to lie at $z_{\rm spec} = 4.002$ via detection of [CI](1-0), $^{12}$CO(6-5), $^{12}$CO(4-3), $^{12}$CO(2-1) and ${\rm H_2O} (2_{11} - 2_{02})$ emission lines, detected using ALMA and ATCA. The spectroscopically-confirmed components of the proto-cluster are distributed over a ${\rm 260\, kpc \times 310\, kpc}$ region and have a collective obscured star-formation rate (SFR) of $\sim 6500 \, M_\odot \, {\rm yr}^{-1}$, considerably higher than has been seen before in any proto-cluster of galaxies or over-densities of DSFGs at $z \gtrsim 4$. Most of the star formation is taking place in luminous DSFGs since no Ly$\alpha$ emitters are detected in the proto-cluster core, apart from a Ly$\alpha$ blob located next to one of the DRC dusty components and extending over $60\,{\rm kpc}$. The total obscured SFR of the proto-cluster could rise to ${\rm SFR} \sim 14,400 \, M_\odot \, {\rm yr}^{-1}$ if all the members of an over-density of bright DSFGs discovered around DRC in a wide-field LABOCA 870-$\mu$m image are part of the same structure. The total halo mass of DRC could be as high as $\sim 4.4 \times 10^{13}\,M_\odot$ and could be the progenitor of a Coma-like cluster at $z = 0$. The relatively short gas-depletion times of the DRC components suggest either the presence of a mechanism able to trigger extreme star formation simultaneously in galaxies spread over a few hundred kpc or the presence of gas flows from the cosmic web able to sustain star formation over several hundred million years.
• ### The SCUBA-2 Cosmology Legacy Survey: the nature of bright submm galaxies from 2 deg2 of 850-um imaging(1610.02409)

Sept. 16, 2017 astro-ph.GA
We present physical properties [redshifts (z), star-formation rates (SFRs) and stellar masses (Mstar)] of bright (S850>4mJy) submm galaxies in the ~2deg2 COSMOS and UDS fields selected with SCUBA-2/JCMT. We complete the galaxy identification process for all (~2000) S/N>3.5 850-um sources, but focus our scientific analysis on a high-quality sub-sample of 651 S/N>4 sources with complete multi-wavelength coverage including 1.1-mm imaging. We check the reliability of our identifications, and the robustness of the SCUBA-2 fluxes by revisiting the recent ALMA follow-up of 29 sources in our sample. Considering >4mJy ALMA sources, our identification method has a completeness of ~86 per cent with a reliability of ~92 per cent, and only ~15-20 per cent of sources are significantly affected by multiplicity (when a secondary component contributes >1/3 of the primary source flux). The impact of source blending on the 850-um source counts as determined with SCUBA-2 is modest; scaling the single-dish fluxes by ~0.9 reproduces the ALMA source counts. For our final SCUBA-2 sample we find median z=2.40+0.10-0.04, SFR=287+-6Moyr-1, and log(Mstar/Mo)=11.12+-0.02 (the latter for 349/651 sources with optical identifications). These properties clearly locate bright submm galaxies on the high-mass end of the 'main sequence' of star-forming galaxies out to z~6, suggesting that major mergers are not a dominant driver of the high-redshift submm-selected population. Their number densities are also consistent with the evolving galaxy stellar mass function. Hence, the submm galaxy population is as expected, albeit reproducing the evolution of the main sequence of star-forming galaxies remains a challenge for theoretical models/simulations.
• ### Witnessing the birth of the red sequence: the physical scale and morphology of dust emission in hyper-luminous starbursts in the early Universe(1709.04191)

Sept. 13, 2017 astro-ph.GA
We present high-spatial-resolution ($\sim 0.12''$ or $\approx 800 \, {\rm pc}$ at $z = 4.5$) ALMA $870\,\mu$m dust continuum observations of a sample of 44 ultrared dusty star-forming galaxies (DSFGs) selected from the H-ATLAS and HerMES far-infrared surveys because of their red colors from 250 to 500 $\mu$m: $S_{500} / S_{250} > 1.5$ and $S_{500} / S_{350} > 1.0$. With photometric redshifts in the range $z \sim 4$-6, our sample includes the most luminous starbursting systems in the early Universe known so far, with total obscured star-formation rates (SFRs) of up to $\sim 4,500 \, M_\odot \, {\rm yr}^{-1}$, as well as a population of lensed, less intrinsically luminous sources. The lower limit on the number of ultrared DSFGs at 870 $\mu$m (with flux densities measured from the ALMA maps and thus not affected by source confusion) derived in this work is in reasonable agreement with models of galaxy evolution, whereas there have been reports of conflicts at 500 $\mu$m (where flux densities are derived from SPIRE). Ultrared DSFGs have a variety of morphologies (from relatively extended disks with smooth radial profiles, to compact sources, both isolated and interacting) and an average size, $\theta_{\rm FWHM}$, of $1.46 \pm 0.41\, {\rm kpc}$, considerably smaller than the values reported in previous work for less-luminous DSFGs at lower redshifts. The size and the estimated gas-depletion times of our sources are compatible with their being the progenitors of the most massive, compact, red-and-dead galaxies at $z \sim 2$-3, and ultimately of local ultra-massive elliptical galaxies or massive galaxy clusters. We are witnessing the birth of the high-mass tail of the red sequence of galaxies.
• ### The Herschel Bright Sources (HerBS): Sample definition and SCUBA-2 observations(1709.01514)

Sept. 5, 2017 astro-ph.GA
We present the Herschel Bright Sources (HerBS) sample, a sample of bright, high-redshift Herschel sources detected in the 616.4 square degree H-ATLAS survey. The HerBS sample contains 209 galaxies, selected with a 500 {\mu}m flux density greater than 80 mJy and an estimated redshift greater than 2. The sample consists of a combination of HyLIRGs and lensed ULIRGs during the epoch of peak cosmic star formation. In this paper, we present SCUBA-2 observations at 850 ${\mu}$m of 189 galaxies of the HerBS sample, 152 of these sources were detected. We fit a spectral template to the Herschel-SPIRE and 850 ${\mu}$m SCUBA-2 flux densities of 22 sources with spectroscopically determined redshifts, using a two-component modified blackbody spectrum as a template. We find a cold- and hot-dust temperature of 21.29 K and 45.80 K, a cold-to-hot dust mass ratio of 26.62 and a $\beta$ of 1.83. The poor quality of the fit suggests that the sample of galaxies is too diverse to be explained by our simple model. Comparison of our sample to a galaxy evolution model indicates that the fraction of lenses is high. Out of the 152 SCUBA-2 detected galaxies, the model predicts 128.4 $\pm$ 2.1 of those galaxies to be lensed (84.5%). The SPIRE 500 ${\mu}$m flux suggests that out of all 209 HerBS sources, we expect 158.1 $\pm$ 1.7 lensed sources, giving a total lensing fraction of 76 per cent.
• ### The most distant, luminous, dusty star-forming galaxies: redshifts from NOEMA and ALMA spectral scans(1707.08967)

July 27, 2017 astro-ph.GA
We present 1.3- and/or 3-mm continuum images and 3-mm spectral scans, obtained using NOEMA and ALMA, of 21 distant, dusty, star-forming galaxies (DSFGs). Our sample is a subset of the galaxies selected by Ivison et al. (2016) on the basis of their extremely red far-infrared (far-IR) colours and low {\it Herschel} flux densities; most are thus expected to be unlensed, extraordinarily luminous starbursts at $z \gtrsim 4$, modulo the considerable cross-section to gravitational lensing implied by their redshift. We observed 17 of these galaxies with NOEMA and four with ALMA, scanning through the 3-mm atmospheric window. We have obtained secure redshifts for seven galaxies via detection of multiple CO lines, one of them a lensed system at $z=6.027$ (two others are also found to be lensed); a single emission line was detected in another four galaxies, one of which has been shown elsewhere to lie at $z=4.002$. Where we find no spectroscopic redshifts, the galaxies are generally less luminous by 0.3-0.4 dex, which goes some way to explaining our failure to detect line emission. We show that this sample contains amongst the most luminous known star-forming galaxies. Due to their extreme star-formation activity, these galaxies will consume their molecular gas in $\lesssim 100$ Myr, despite their high molecular gas masses, and are therefore plausible progenitors of the massive, red-and-dead' elliptical galaxies at $z \approx 3$.
• ### Low-redshift analogs of submm galaxies: a diverse population(1707.05329)

July 17, 2017 astro-ph.GA
We have combined the wide-area Herschel-ATLAS far-IR survey with spectroscopic redshifts from GAMA and SDSS to define a sample of 21 low--redshift ($z_{\rm spec} < 0.5$) analogs of submm galaxies (SMGs). These have been selected because their dust temperatures and total IR luminosities are similar to those for the classical high-redshift SMG population. As well as presenting the sample, in this paper we report $^{12}$CO(2-1) and $^{12}$CO(1-0) observations of 16 low-redshift analogs of SMGs taken with the IRAM-30m telescope. We have obtained that low-redshift analogs of SMGs represent a very diverse population, similar to what has been found for high-redshift SMGs. A large variety in the molecular gas excitation or $^{12}$CO(2-1)/$^{12}$CO(1-0) line ratio is seen, meaning that extrapolations from $J \geq 2$ CO lines can result in very uncertain molecular gas mass determinations. Our sources with $^{12}$CO(1-0) detections follow the dust--gas correlation found in previous work at different redshifts and luminosities. The molecular gas mass of low-redshift SMGs has an average value of $M_{\rm H_2} \sim 1.6 \times 10^{10}\,M_\odot$ and will be consumed in $\sim 100 \, {\rm Myr}$ . We also find a wide range of molecular gas fractions, with the highest values being compatible with those found in high-redshift SMGs with $^{12}$CO(1-0) detections, which are only the most luminous. Low-redshift SMGs offer a unique opportunity to study the properties of extreme star formation in a detail not possible at higher redshifts.
• ### HIFI Spectroscopy of ${\rm H_2O}$ submm Lines in Nuclei of Actively Star Forming Galaxies(1707.04914)

July 16, 2017 astro-ph.GA
We present a systematic survey of multiple velocity-resolved H$_2$O spectra using Herschel/HIFI towards nine nearby actively star forming galaxies. The ground-state and low-excitation lines (E$_{\rm up}\,\le 130\,{\rm K}$) show profiles with emission and absorption blended together, while absorption-free medium-excitation lines ($130\,{\rm K}\, \le\, E_{\rm up}\,\le\,350\,{\rm K}$) typically display line shapes similar to CO. We analyze the HIFI observation together with archival SPIRE/PACS H$_2$O data using a state-of-the-art 3D radiative transfer code which includes the interaction between continuum and line emission. The water excitation models are combined with information on the dust- and CO spectral line energy distribution to determine the physical structure of the interstellar medium (ISM). We identify two ISM components that are common to all galaxies: A warm ($T_{\rm dust}\,\sim\,40-70\,{\rm K}$), dense ($n({\rm H})\,\sim\,10^5-10^6\,{\rm cm^{-3}}$) phase which dominates the emission of medium-excitation H$_2$O lines. This gas phase also dominates the FIR emission and the CO intensities for $J_{\rm up} > 8$. In addition a cold ($T_{\rm dust}\,\sim\,20-30\,{\rm K}$), dense ($n({\rm H})\sim\,10^4- 10^5\,{\rm cm^{-3}}$) more extended phase is present. It outputs the emission in the low-excitation H$_2$O lines and typically also produces the prominent line absorption features. For the two ULIRGs in our sample (Arp 220 and Mrk 231) an even hotter and more compact (R$_s\,\le\,100$ pc) region is present which is possibly linked to AGN activity. We find that collisions dominate the water excitation in the cold gas and for lines with $E_{\rm up}\le300\,{\rm K}$ and $E_{\rm up}\le800\,{\rm K}$ in the warm and hot component, respectively. Higher energy levels are mainly excited by IR pumping.
• ### Recovering the properties of high redshift galaxies with different JWST broad-band filters(1706.05063)

June 15, 2017 astro-ph.GA
Imaging with the James Webb Space Telescope (JWST) will allow for observing the bulk of distant galaxies at the epoch of reionisation. The recovery of their properties, such as age, color excess E(B-V), specific star formation rate (sSFR) and stellar mass, will mostly rely on spectral energy distribution fitting, based on the data provided by JWST's two imager cameras, namely the Near Infrared Camera (NIRCam) and the Mid Infrared Imager (MIRI). In this work we analyze the effect of choosing different combinations of NIRCam and MIRI broad-band filters, from 0.6 {\mu}m to 7.7 {\mu}m, on the recovery of these galaxy properties. We performed our tests on a sample of 1542 simulated galaxies, with known input properties, at z=7-10. We found that, with only 8 NIRCam broad-bands, we can recover the galaxy age within 0.1 Gyr and the color excess within 0.06 mag for 70% of the galaxies. Besides, the stellar masses and sSFR are recovered within 0.2 and 0.3 dex, respectively, at z=7-9. Instead, at z=10, no NIRCam band traces purely the {\lambda}> 4000 {\AA} regime and the percentage of outliers in stellar mass (sSFR) increases by > 20% (> 90%), in comparison to z=9. The MIRI F560W and F770W bands are crucial to improve the stellar mass and the sSFR estimation at z=10. When nebular emission lines are present, deriving correct galaxy properties is challenging, at any redshift and with any band combination. In particular, the stellar mass is systematically overestimated in up to 0.3 dex on average with NIRCam data alone and including MIRI observations improves only marginally the estimation.
• ### VALES: I. The molecular gas content in star-forming dusty H-ATLAS galaxies up to z=0.35(1705.09826)

June 1, 2017 astro-ph.GA
We present an extragalactic survey using observations from the Atacama Large Millimeter/submillimeter Array (ALMA) to characterise galaxy populations up to $z=0.35$: the Valpara\'iso ALMA Line Emission Survey (VALES). We use ALMA Band-3 CO(1--0) observations to study the molecular gas content in a sample of 67 dusty normal star-forming galaxies selected from the $Herschel$ Astrophysical Terahertz Large Area Survey ($H$-ATLAS). We have spectrally detected 49 galaxies at $>5\sigma$ significance and 12 others are seen at low significance in stacked spectra. CO luminosities are in the range of $(0.03-1.31)\times10^{10}$ K km s$^{-1}$ pc$^2$, equivalent to $\log({\rm M_{gas}/M_{\odot}}) =8.9-10.9$ assuming an $\alpha_{\rm CO}$=4.6(K km s$^{-1}$ pc$^{2}$)$^{-1}$, which perfectly complements the parameter space previously explored with local and high-z normal galaxies. We compute the optical to CO size ratio for 21 galaxies resolved by ALMA at $\sim 3$."$5$ resolution (6.5 kpc), finding that the molecular gas is on average $\sim$ 0.6 times more compact than the stellar component. We obtain a global Schmidt-Kennicutt relation, given by $\log [\Sigma_{\rm SFR}/({\rm M_{\odot} yr^{-1}kpc^{-2}})]=(1.26 \pm 0.02) \times \log [\Sigma_{\rm M_{H2}}/({\rm M_{\odot}\,pc^{-2}})]-(3.6 \pm 0.2)$. We find a significant fraction of galaxies lying at intermediate efficiencies' between a long-standing mode of star-formation activity and a starburst, specially at $\rm L_{IR}=10^{11-12} L_{\odot}$. Combining our observations with data taken from the literature, we propose that star formation efficiencies can be parameterised by $\log [{\rm SFR/M_{H2}}]=0.19 \times {\rm (\log {L_{IR}}-11.45)}-8.26-0.41 \times \arctan[-4.84 (\log {\rm L_{IR}}-11.45) ]$. Within the redshift range we explore ($z<0.35$), we identify a rapid increase of the gas content as a function of redshift.
• ### VALES: III. The calibration between the dust continuum and interstellar gas content of star-forming galaxies(1702.07350)

Feb. 23, 2017 astro-ph.GA
We present the calibration between the dust continuum luminosity and interstellar gas content obtained from the Valpara\'{i}so ALMA Line Emission Survey (VALES) sample of 67 main-sequence star-forming galaxies at 0.02<$z$<0.35. We use CO(1-0) observations from the Atacama Large Millimetre/submillimetre Array (ALMA) to trace the molecular gas mass, $M_{\mathrm{H}_{2}}$, and estimate the rest-frame monochromatic luminosity at 850 $\mu$m, $L_{\nu_{850}}$, by extrapolating the dust continuum from MAGPHYS modelling of the far-ultraviolet to submillimetre spectral energy distribution sampled by the Galaxy And Mass Assembly (GAMA) survey. Adopting $\alpha_{\rm CO}$ = 6.5 (K km s$^{-1}$ pc$^{2}$)$^{-1}$, the average ratio of $L_{\nu_{850}}/M_{\mathrm{H}_{2}}$ = (6.4$\pm$1.4)$\times10^{19}$ erg s$^{-1}$ Hz$^{-1}$ $\mathrm{M}_{\odot}^{-1}$, in excellent agreement with literature values. We obtain a linear fit of $\log_{10}$ ($M_{\mathrm{H}_{2}}/\mathrm{M}_{\odot}$) = (0.92$\pm$0.02) $\log_{10}$ ($L_{\nu_{850}}$/erg s$^{-1}$ Hz$^{-1}$)-(17.31$\pm$0.59). We provide relations between $L_{\nu_{850}}$, $M_{\mathrm{H}_{2}}$ and $M_{\mathrm{ISM}}$ when combining the VALES and literature samples, and adopting a Galactic $\alpha_{\rm CO}$ value.
• ### Luminous, pc-scale CO 6-5 emission in the obscured nucleus of NGC1377(1702.05458)

Feb. 17, 2017 astro-ph.GA
High resolution submm observations are important in probing the morphology, column density and dynamics of obscured active galactic nuclei (AGNs). With high resolution (0.06 x 0.05) ALMA 690 GHz observations we have found bright (TB >80 K) and compact (FWHM 10x7 pc) CO 6-5 line emission in the nucleus of the extremely radio-quiet galaxy NGC1377. The integrated CO 6-5 intensity is aligned with the previously discovered jet/outflow of NGC1377 and is tracing the dense (n>1e4 cm-3), hot gas at the base of the outflow. The velocity structure is complex and shifts across the jet/outflow are discussed in terms of jet-rotation or separate, overlapping kinematical components. High velocity gas (deltaV +-145 km/s) is detected inside r<2-3 pc and we suggest that it is emerging from an inclined rotating disk or torus of position angle PA=140+-20 deg with a dynamical mass of approx 3e6 Msun. This mass is consistent with that of a supermassive black hole (SMBH), as inferred from the M-sigma relation. The gas mass of the proposed disk/torus constitutes <3% of the nuclear dynamical mass. In contrast to the intense CO 6-5 line emission, we do not detect dust continuum with an upper limit of S(690GHz)<2mJy. The corresponding, 5 pc, H2 column density is estimated to N(H2)<3e23 cm-2, which is inconsistent with a Compton Thick (CT) source. We discuss the possibility that CT obscuration may be occuring on small (subparsec) or larger scales. From SED fitting we suggest that half of the IR emission of NGC1377 is nuclear and the rest (mostly the far-infrared (FIR)) is more extended. The extreme radio quietness, and the lack of emission from other star formation tracers, raise questions on the origin of the FIR emission. We discuss the possibility that it is arising from the dissipation of shocks in the molecular jet/outflow or from irradiation by the nuclear source along the poles.
• ### The AT-LESS CO(1-0) survey of submillimetre galaxies in the Extended Chandra Deep Field South: First results on cold molecular gas in galaxies at z ~ 2(1701.05698)

Jan. 20, 2017 astro-ph.GA
We present the first results from our on-going Australia Telescope Compact Array survey of CO(1-0) in ALMA-identified submillimetre galaxies in the Extended Chandra Deep Field South. Strong detections of CO(1-0) emission from two submillimetre galaxies, ALESS 122.1 (z = 2.0232) and ALESS 67.1 (z = 2.1230), were obtained. We estimate gas masses of M_gas ~ 1.3 \times 10^{11} M_odot and M_gas ~ 1.0 \times 10^{11} M_\odot for ALESS 122.1 and ALESS 67.1, respectively, adopting alpha_CO = 1.0. Dynamical mass estimates from the kinematics of the CO(1-0) line yields M_dyn (sin i)^2 = 2.1 +- 1.1 \times 10^{11} M_odot and (3.2 +- 0.9) \times 10^{11} M_\odot for ALESS 122.1 and ALESS 67.1, respectively. This is consistent with the total baryonic mass estimates of these two systems. We examine star formation efficiency using the L_FIR versus L'_CO(1-0) relation for samples of local ULIRGs and LIRGs, and more distant star-forming galaxies, with CO(1-0) detections. We find some evidence of a shallower slope for ULIRGs and SMGs compared to less luminous systems, but a larger sample is required for definite conclusions. We determine gas-to-dust ratios of 170 +- 30 and 140 +- 30 for ALESS 122.1 and ALESS 67.1, respectively, showing ALESS 122.1 has an unusually large gas reservoir. By combining the 38.1 GHz continuum detection of ALESS 122.1 with 1.4 and 5.5 GHz data, we estimate that the free-free contribution to radio emission at 38.1 GHz is 34 +- 17 microJy, yielding a star formation rate (1400 +- 700 M_\odot yr^{-1}) consistent with that from the infrared luminosity.
• ### High dense gas fraction in intensely star-forming dusty galaxies at high redshift(1701.05901)

Jan. 20, 2017 astro-ph.GA
We present ALMA and VLA detections of the dense molecular gas tracers HCN, HCO$^+$ and HNC in two lensed, high-redshift starbursts selected from the {\it Herschel}-ATLAS survey: {\it H}-ATLAS\,J090740.0$-$004200 (SDP.9, $z \sim 1.6$) and {\it H}-ATLAS\,J091043.1$-$000321 (SDP.11, $z \sim 1.8$). ALMA observed the $J = 3-2$ transitions in both sources, while the VLA observed the $J = 1-0$ transitions in SDP.9. We have detected all observed HCN and HCO$^+$ lines in SDP.9 and SDP.11, and also HNC(3--2) in SDP.9. The amplification factors for both galaxies have been determined from sub-arcsec resolution CO and dust emission observations carried out with NOEMA and the SMA. The HNC(1--0)/HCN(1--0) line ratio in SDP.9 suggests the presence of photon-dominated regions, as it happens to most local (U)LIRGs. The CO, HCN and HCO$^+$ SLEDs of SDP.9 are compatible to those found for many local, infrared (IR) bright galaxies, indicating that the molecular gas in local and high-redshift dusty starbursts can have similar excitation conditions. We obtain that the correlation between total IR ($L_{\rm IR}$) and dense line ($L_{\rm dense}$) luminosity in SDP.9 and SDP.11 and local star-forming galaxies can be represented by a single relation. The scatter of the $L_{\rm IR} - L_{\rm dense}$ correlation, together with the lack of sensitive dense molecular gas tracer observations for a homogeneous sample of high-redshift galaxies, prevents us from distinguishing differential trends with redshift. Our results suggest that the intense star formation found in some high-redshift dusty, luminous starbursts is associated with more massive dense molecular gas reservoirs and higher dense molecular gas fractions.
• ### Evolution of cosmic star formation in the SCUBA-2 Cosmology Legacy Survey(1607.04283)

Jan. 5, 2017 astro-ph.GA, astro-ph.IM
We present a new exploration of the cosmic star-formation history and dust obscuration in massive galaxies at redshifts $0.5< z<6$. We utilize the deepest 450 and 850$\mu$m imaging from SCUBA-2 CLS, covering 230arcmin$^2$ in the AEGIS, COSMOS and UDS fields, together with 100-250$\mu$m imaging from Herschel. We demonstrate the capability of the T-PHOT deconfusion code to reach below the confusion limit, using multi-wavelength prior catalogues from CANDELS/3D-HST. By combining IR and UV data, we measure the relationship between total star-formation rate (SFR) and stellar mass up to $z\sim5$, indicating that UV-derived dust corrections underestimate the SFR in massive galaxies. We investigate the relationship between obscuration and the UV slope (the IRX-$\beta$ relation) in our sample, which is similar to that of low-redshift starburst galaxies, although it deviates at high stellar masses. Our data provide new measurements of the total SFR density (SFRD) in $M_\ast>10^{10}M_\odot$ galaxies at $0.5<z<6$. This is dominated by obscured star formation by a factor of $>10$. One third of this is accounted for by 450$\mu$m-detected sources, while one fifth is attributed to UV-luminous sources (brighter than $L^\ast_{UV}$), although even these are largely obscured. By extrapolating our results to include all stellar masses, we estimate a total SFRD that is in good agreement with previous results from IR and UV data at $z\lesssim3$, and from UV-only data at $z\sim5$. The cosmic star-formation history undergoes a transition at $z\sim3-4$, as predominantly unobscured growth in the early Universe is overtaken by obscured star formation, driven by the build-up of the most massive galaxies during the peak of cosmic assembly.
• ### The Dense Molecular Gas and Nuclear Activity in the ULIRG IRAS 13120-5453(1612.04401)

Dec. 21, 2016 astro-ph.GA
We present new ALMA Band 7 ($\sim340$ GHz) observations of the dense gas tracers HCN, HCO$^+$, and CS in the local, single-nucleus, ultraluminous infrared galaxy IRAS 13120-5453. We find centrally enhanced HCN (4-3) emission, relative to HCO$^+$ (4-3), but do not find evidence for radiative pumping of HCN. Considering the size of the starburst (0.5 kpc) and the estimated supernovae rate of $\sim1.2$ yr$^{-1}$, the high HCN/HCO$^+$ ratio can be explained by an enhanced HCN abundance as a result of mechanical heating by the supernovae, though the active galactic nucleus and winds may also contribute additional mechanical heating. The starburst size implies a high $\Sigma_{IR}$ of $4.7\times10^{12}$ $L_{\odot}$ kpc$^{-2}$, slightly below predictions of radiation-pressure limited starbursts. The HCN line profile has low-level wings, which we tentatively interpret as evidence for outflowing dense molecular gas. However, the dense molecular outflow seen in the HCN line wings is unlikely to escape the galaxy and is destined to return to the nucleus and fuel future star formation. We also present modeling of Herschel observations of the H$_2$O lines and find a nuclear dust temperature of $\sim40$ K. IRAS 13120-5453 has a lower dust temperature and $\Sigma_{IR}$ than is inferred for the systems termed "compact obscured nuclei" (such as Arp 220 and Mrk 231). If IRAS 13120-5453 has undergone a compact obscured nucleus phase, we are likely witnessing it at a time when the feedback has already inflated the nuclear ISM and diluted star formation in the starburst/AGN core.
• ### VALES: II. The physical conditions of interstellar gas in normal star-forming galaxies up to z=0.2 revealed by ALMA(1611.05867)

Nov. 17, 2016 astro-ph.GA
We use new Band-3 CO(1-0) observations taken with the Atacama Large Millimeter/submillimeter Array (ALMA) to study the physical conditions in the interstellar gas of a sample of 27 dusty main-sequence star-forming galaxies at 0.03<$z$<0.2 present in the Valpara\'iso ALMA Line Emission Survey (VALES). The sample is drawn from far-IR bright galaxies over $\sim$160 deg$^{2}$ in the Herschel Astrophysical Terahertz Large Area Survey (HATLAS), which is covered by Herschel [CII] 158 $\mu$m spectroscopy and far-infrared (FIR) photometry. The [CII] and CO lines are both detected at >5$\sigma$ in 26 sources. We find an average [CII] to CO(1-0) luminosity ratio of 3500$\pm$1200 for our sample that is consistent with previous studies. Using the [CII], CO and FIR measurements as diagnostics of the physical conditions of the interstellar medium, we compare these observations to the predictions of a photodissociation region (PDR) model to determine the gas density, surface temperature, pressure, and the strength of the incident far-ultraviolet (FUV) radiation field, $G_{0}$, normalised to the Habing Field. The majority of our sample exhibit hydrogen densities of 4 < $\log n/\mathrm{cm}^{3}$ < 5.5 and experience an incident FUV radiation field with strengths of 2 < $\log G_0$ < 3 when adopting standard adjustments. A comparison to galaxy samples at different redshifts indicates that the average strength of the FUV radiation field appears constant up to redshift $z\sim$6.4, yet the neutral gas density increases with redshift by a factor of $\sim$100, that persists regardless of various adjustments to our observable quantities. This evolution could provide an explanation for the observed evolution of the star formation rate density with cosmic time, yet could arise from a combination of observational biases when using different suites of emission lines as diagnostic tracers of PDR gas.
• ### The Herschel-ATLAS: a sample of 500{\mu}m-selected lensed galaxies over 600 square degrees(1611.03922)

Nov. 12, 2016 astro-ph.GA
We present a sample of 80 candidate strongly lensed galaxies with flux density above 100mJy at 500{\mu}m extracted from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS), over an area of 600 square degrees. Available imaging and spectroscopic data allow us to confirm the strong lensing in 20 cases and to reject it in one case. For other 8 objects the lensing scenario is strongly supported by the presence of two sources along the same line of sight with distinct photometric redshifts. The remaining objects await more follow-up observations to confirm their nature. The lenses and the background sources have median redshifts z_L = 0.6 and z_S = 2.5, respectively, and are observed out to z_L = 1.2 and z_S = 4.2. We measure the number counts of candidate lensed galaxies at 500{\mu}m and compare them with theoretical predictions, finding a good agreement for a maximum magnification of the background sources in the range 10-20. These values are consistent with the magnification factors derived from the lens modelling of individual systems. The catalogue presented here provides sub- mm bright targets for follow-up observations aimed at exploiting gravitational lensing to study with un-precedented details the morphological and dynamical properties of dusty star forming regions in z >~ 1.5 galaxies.