Search

We present new JWST/NIRSpec IFU observations of the J1000+0234 system at $z=4.54$, the dense core of a galaxy protocluster hosting a massive, dusty star forming galaxy (DSFG) with a low luminosity radio counterpart. The new data reveals two extended, high equivalent width (EW$_0 > 1000$ {\AA}) nebulae at each side of the DSFG disk along its minor axis (namely O3-N and O3-S). On one hand, O3-N's spectrum shows a prominent FWHM $\sim1300$ km s$^{-1}$ broad and blueshifted component, suggesting an outflow origin. On the other hand, O3-S stretches over parsec and has a velocity gradient that spans $800$ km s$^{-1}$ but no evidence of a broad component. Both sources, however, seem to be powered at least partially by an active galactic nucleus (AGN), so we classify them as extended emission-line regions (EELRs). The strongest evidence comes from the detection of the high-ionization [Ne V] $\lambda3427$ line toward O3-N, which paired with the non-detection of hard X-rays implies an obscuring column density above the Compton-thick regime. In O3-S, the [Ne V] line is not detected, but we measure a He II well above the expectation for star formation. We interpret this as O3-S being externally irradiated by the AGN, akin to the famous Hanny's Voorwerp object in the local Universe. In addition, more classical line ratio diagnostics (e.g. [O III]/H$\beta$ vs [N II]/H$\alpha$) put the DSFG itself in the AGN region of the diagrams, and hence the most probable host of the AGN. These results showcase the ability of JWST of unveiling highly obscured AGN at high redshifts., Comment: 5 pages, 4 figures plus 5 appendices (incl. 3 extra figures and one table). Submitted to A&A on July 17th 2024

We present new $\lambda_{\rm rest}=77$ $\mu$m dust continuum observations from the ALMA of HZ10 (CRISTAL-22), a dusty main-sequence galaxy at $z$=5.66 as part of the [CII] Resolved Ism in STar-forming Alma Large program, CRISTAL. The high angular resolution of the ALMA Band 7 and new Band 9 data($\sim{0}''.4$) reveals the complex structure of HZ10, which comprises two main components (HZ10-C and HZ10-W) and a bridge-like dusty emission between them (the Bridge). We model the dust spectral energy distribution (SED) to constrain the physical conditions of the interstellar medium (ISM) and its variations among the different components identified in HZ10. We find that HZ10-W (the more UV-obscured component) has an SED dust temperature of $T_{\rm SED}$$\sim$51.2$\pm13.1$ K; this is $\sim$5 K higher (although still consistent) than that of the central component and previous global estimations for HZ10. Our new ALMA data allow us to reduce by a factor of $\sim$2.3 the uncertainties of global $T_{\rm SED}$ measurements compared to previous studies. Interestingly, HZ10-W shows a lower [CII]/FIR ratio compared to the other two components (although still within the uncertainties), suggesting a harder radiation field destroying polycyclic aromatic hydrocarbon associated with [CII] emission (e.g., active galactic nuclei or young stellar populations). While HZ10-C appears to follow the tight IRX-$\beta_{\rm UV}$ relation seen in local UV-selected starburst galaxies and high-$z$ star-forming galaxies, we find that both HZ10-W and the Bridge depart from this relation and are well described by dust-screen models with holes in front of a hard UV radiation field. This suggests that the UV emission (likely from young stellar populations) is strongly attenuated in the more dusty components of the HZ10 system.

The ALMA [CII] Resolved Ism in STar-forming gALaxies (CRISTAL) survey is a Cycle 8 ALMA Large Programme that studies the cold gas component of high-redshift galaxies. Its sub-arcsecond resolution observations are key to disentangling physical mechanisms that shape galaxies during cosmic dawn. In this paper, we explore the morphology and kinematics of the cold gas, star-forming, and stellar components in the star-forming main-sequence galaxy CRISTAL-05/HZ3, at z = 5.54. Our analysis includes 0.3" spatial resolution (~2 kpc) ALMA observations of the [CII] line. While CRISTAL-05 was previously classified as a single source, our observations reveal that the system is a close interacting pair surrounded by an extended component of carbon-enriched gas. This is imprinted in the disturbed elongated [CII] morphology and the separation of the two components in the position-velocity diagram (~100 km/s). The central region is composed of two components, named C05-NW and C05-SE, with the former being the dominant one. A significant fraction of the [CII] arises beyond the close pair up to 10 kpc, while the regions forming new massive stars and the stellar component seem compact (r_[CII] ~ 4 r_UV), as traced by rest-frame UV and optical imaging obtained with the Hubble Space Telescope and the James Webb Space Telescope. Our kinematic model, using the DYSMALpy software, yields a minor contribution of dark matter of C05-NW within a radius of ~2x Reff. Finally, we explore the resolved [CII]/FIR ratios as a proxy for shock-heating produced by this merger. We argue that the extended [CII] emission is mainly caused by the merger, which could not be discerned with lower-resolution observations. Our work emphasizes the need for high-resolution observations to fully characterize the dynamic stages of infant galaxies and the physical mechanisms that drive the metal enrichment of the circumgalactic medium., Comment: Submitted to A&A - comments are welcome! - 19 pages, 13 figures

We present a new method for modelling the kinematics of galaxies from interferometric observations by performing the optimization of the kinematic model parameters directly in visibility-space instead of the conventional approach of fitting velocity fields produced with the CLEAN algorithm in real-space. We demonstrate our method on ALMA observations of $^{12}$CO (2$-$1), (3$-$2) or (4$-$3) emission lines from an initial sample of 30 massive 850$\mu$m-selected dusty star-forming galaxies with far-infrared luminosities $\gtrsim$$\,10^{12}\,$L$_{\odot}$ in the redshift range $z \sim\,$1.2$-$4.7. Using the results from our modelling analysis for the 12 sources with the highest signal-to-noise emission lines and disk-like kinematics, we conclude the following: (i) Our sample prefers a CO-to-$H_2$ conversion factor, of $\alpha_{\rm CO} = 0.92 \pm 0.36$; (ii) These far-infrared luminous galaxies follow a similar Tully$-$Fisher relation between the circularized velocity, $V_{\rm circ}$, and baryonic mass, $M_{\rm b}$, as more typical star-forming samples at high redshift, but extend this relation to much higher masses $-$ showing that these are some of the most massive disk-like galaxies in the Universe; (iii) Finally, we demonstrate support for an evolutionary link between massive high-redshift dusty star-forming galaxies and the formation of local early-type galaxies using the both the distributions of the baryonic and kinematic masses of these two populations on the $M_{\rm b}\,-\,\sigma$ plane and their relative space densities.

We investigate the degree of dust obscured star formation in 49 massive (${\rm log}_{10}(M_{\star}/{\rm M}_{\odot})>9$) Lyman-break galaxies (LBGs) at $z = 6.5$-$8$ observed as part of the ALMA Reionization Era Bright Emission Line Survey (REBELS) large program. By creating deep stacks of the photometric data and the REBELS ALMA measurements we determine the average rest-frame UV, optical and far-infrared (FIR) properties which reveal a significant fraction ($f_{\rm obs} = 0.4$-$0.7$) of obscured star formation, consistent with previous studies. From measurements of the rest-frame UV slope, we find that the brightest LBGs at these redshifts show bluer ($\beta \simeq -2.2$) colours than expected from an extrapolation of the colour-magnitude relation found at fainter magnitudes. Assuming a modified blackbody spectral-energy distribution (SED) in the FIR (with dust temperature of $T_{\rm d} = 46\,{\rm K}$ and $\beta_{\rm d} = 2.0$), we find that the REBELS sources are in agreement with the local ''Calzetti-like'' starburst Infrared-excess (IRX)-$\beta$ relation. By reanalysing the data available for 108 galaxies at $z \simeq 4$-$6$ from the ALPINE ALMA large program using a consistent methodology and assumed FIR SED, we show that from $z \simeq 4$-$8$, massive galaxies selected in the rest-frame UV have no appreciable evolution in their derived IRX-$\beta$ relation. When comparing the IRX-$M_{\star}$ relation derived from the combined ALPINE and REBELS sample to relations established at $z < 4$, we find a deficit in the IRX, indicating that at $z > 4$ the proportion of obscured star formation is lower by a factor of $\gtrsim 3$ at a given a $M_{\star}$. Our IRX-$\beta$ results are in good agreement with the high-redshift predictions of simulations and semi-analytic models for $z \simeq 7$ galaxies with similar stellar masses and SFRs., Comment: 18 pages, 10 figures, 2 tables (plus 1 figure and 2 tables in the appendix). Submitted to MNRAS, comments welcome

A key to understanding the formation of the first galaxies is to quantify the content of the molecular gas as the fuel for star formation activity through the epoch of reionization. In this paper, we use the 158$\mu$m [CII] fine-structure emission line as a tracer of the molecular gas in the interstellar medium (ISM) in a sample of $z=6.5-7.5$ galaxies recently unveiled by the Reionization Era Bright Line Emission Survey, REBELS, with the Atacama Large Millimeter/submillimeter Array. We find substantial amounts of molecular gas ($\sim10^{10.5}\ M_\odot$) comparable to those found in lower redshift galaxies for similar stellar masses ($\sim10^{10}\ M_\odot$). The REBELS galaxies appear to follow the standard scaling relations of molecular gas to stellar mass ratio ($\mu_{\rm mol}$) and gas depletion timescale ($t_{\rm dep}$) with distance to the star-forming main-sequence expected from extrapolations of $z\sim1-4$ observations. We find median values at $z\sim7$ of $\mu_{\rm mol}=2.6_{-1.4}^{4.1}$ and $t_{\rm dep}=0.5_{-0.14}^{+0.26}$ Gyr, indicating that the baryonic content of these galaxies is gas-phase dominated and little evolution from $z\sim7$ to 4. Our measurements of the cosmic density of molecular gas, log$(\rho_{\rm mol}/(M_\odot {\rm Mpc}^{-3}))=6.34^{+0.34}_{-0.31}$, indicate a steady increase by an order of magnitude from $z\sim7$ to 4., Comment: Accepted for publication in Astronomy & Astrophysics. 11 pages, 7 figures. Typo fixed in an author's name

We present Atacama Large Millimeter/Submillimeter Array (ALMA) [CII] and $\sim158$ $\rm\mu m$ continuum observations of REBELS-25, a massive, morphologically complex ultra-luminous infrared galaxy (ULIRG; $L_{\rm IR}=1.5^{+0.8}_{-0.5}\times10^{12}$ L$_\odot$) at $z=7.31$, spectroscopically confirmed by the Reionization Era Bright Emission Line Survey (REBELS) ALMA Large Programme. REBELS-25 has a significant stellar mass of $M_{*}=8^{+4}_{-2}\times10^{9}$ M$_\odot$. From dust-continuum and ultraviolet observations, we determine a total obscured + unobscured star formation rate of SFR $=199^{+101}_{-63}$ M$_\odot$ yr$^{-1}$. This is about four times the SFR estimated from an extrapolated main-sequence. We also infer a [CII]-based molecular gas mass of $M_{\rm H_2}=5.1^{+5.1}_{-2.6}\times10^{10}$ $M_\odot$, implying a molecular gas depletion time of $ t_{\rm depl, H_2}=0.3^{+0.3}_{-0.2}$ Gyr. We observe a [CII] velocity gradient consistent with disc rotation, but given the current resolution we cannot rule out a more complex velocity structure such as a merger. The spectrum exhibits excess [CII] emission at large positive velocities ($\sim500$ km s$^{-1}$), which we interpret as either a merging companion or an outflow. In the outflow scenario, we derive a lower limit of the mass outflow rate of 200 M$_\odot$ yr$^{-1}$, which is consistent with expectations for a star formation-driven outflow. Given its large stellar mass, SFR and molecular gas reservoir $\sim700$ Myr after the Big Bang, we explore the future evolution of REBELS-25. Considering a simple, conservative model assuming an exponentially declining star formation history, constant star formation efficiency, and no additional gas inflow, we find that REBELS-25 has the potential to evolve into a galaxy consistent with the properties of high-mass quiescent galaxies recently observed at $z\sim4$., Comment: Accepted for publication in MNRAS. 21 pages, 8 figures

We present the first observational infrared luminosity function (IRLF) measurement in the Epoch of Reionization (EoR) based on a UV-selected galaxy sample with ALMA spectroscopic observations. Our analysis is based on the ALMA large program Reionization Era Bright Emission Line Survey (REBELS), which targets 42 galaxies at $\mathrm{z=6.4-7.7}$ with [CII] 158$\micron$ line scans. 16 sources exhibit a dust detection, 15 of which are also spectroscopically confirmed through the [CII] line. The IR luminosities of the sample range from $\log L_{IR}/L_\odot=11.4$ to 12.2. Using the UVLF as a proxy to derive the effective volume for each of our target sources, we derive IRLF estimates, both for detections and for the full sample including IR luminosity upper limits. The resulting IRLFs are well reproduced by a Schechter function with the characteristic luminosity of $\log L_{*}/L_\odot=11.6^{+0.2}_{-0.1}$. Our observational results are in broad agreement with the average of predicted IRLFs from simulations at $z\sim7$. Conversely, our IRLFs lie significantly below lower redshift estimates, suggesting a rapid evolution from $z\sim4$ to $z\sim7$, into the reionization epoch. The inferred obscured contribution to the cosmic star-formation rate density at $z\sim7$ amounts to $\mathrm{log(SFRD/M_{\odot}/yr/Mpc^{3}) = -2.66^{+0.17}_{-0.14} }$ which is at least $\sim$10\% of UV-based estimates. We conclude that the presence of dust is already abundant in the EoR and discuss the possibility of unveiling larger samples of dusty galaxies with future ALMA and JWST observations., Comment: 9 pages, 5 figures

We present the average [CII] $158\,\rm{\mu m}$ emission line sizes of UV-bright star-forming galaxies at $z\sim7$. Our results are derived from a stacking analysis of [CII] $158\,\rm{\mu m}$ emission lines and dust continua observed by ALMA, taking advantage of the large program Reionization Era Bright Emission Line Survey (REBELS). We find that the average [CII] emission at $z\sim7$ has an effective radius $r_e$ of $2.2\pm0.2\,\rm{kpc}$. It is $\gtrsim2\times$ larger than the dust continuum and the rest-frame UV emission, in agreement with recently reported measurements for $z\lesssim6$ galaxies. Additionally, we compared the average [CII] size with $4

We include a fully coupled treatment of metal and dust enrichment into the Delphi semi-analytic model of galaxy formation to explain the dust content of 13 Lyman Break Galaxies (LBGs) detected by the Atacama Large millimetre Array (ALMA) REBELS Large Program at $z\simeq 7$. We find that the galaxy dust mass, $M_d$, is regulated by the combination of SNII dust production, astration, shock destruction, and ejection in outflows; grain growth (with a standard timescale $\tau_0= 30$ Myr) plays a negligible role. The model predicts a dust-to-stellar mass ratio of $\sim 0.07-0.1\%$ and a UV-to-total star formation rate relation such that $log (\psi_{\rm UV}) = -0.05 ~[log (\psi)]^{2} + 0.86 ~log(\psi) -0.05$ (implying that 55-80\% of the star formation is obscured) for REBELS galaxies with stellar mass $M_* = 10^{9-10} M_\odot$. This relation reconciles the intrinsic UV luminosity of LBGs with their observed luminosity function at $z=7$. However, 2 out of the 13 systems show dust-to-stellar mass ratios ($\sim 0.94-1.1\%$) that are up to $18\times$ larger than expected from the fiducial relation. Due to the physical coupling between dust and metal enrichment, even decreasing $\tau_0$ to very low values (0.3 Myr) only increases the dust-to-stellar mass ratio by a factor $ \sim 2$. Given that grain growth is not a viable explanation for such high observed ratios of the dust-to-stellar mass, we propose alternative solutions., Comment: Accepted to MNRAS

We analyse FIR dust continuum measurements for 14 galaxies ($z\approx 7$) in the ALMA REBELS LP to derive their physical properties. Our model uses three input data: (a) the UV spectral slope, $\beta$, (b) the observed UV continuum flux at $1500$A, $F_{\rm UV}$, (c) the observed continuum flux at $\approx 158\mu$m, $F_{158}$, and considers Milky Way (MW) and SMC extinction curves, along with different dust geometries. We find that REBELS galaxies have (28-90.5)% of their star formation obscured; the total (UV+IR) star formation rates are in the range $31.5 < {\rm SFR}/ (M_\odot {\rm yr}^{-1}) < 129.5$. The sample-averaged dust mass and temperature are $(1.3\pm 1.1)\times 10^7 M_\odot$ and $52 \pm 11$ K, respectively. In some galaxies dust is abundant (REBELS-14, $M'_d \approx 3.4 \times 10^7 M_\odot$), or hot (REBELS-18, $T'_d \approx 67$ K). The dust distribution is compact ($<0.3$ kpc for 70% of the galaxies). The dust yield per supernova is $0.1 \le y_d/M_\odot \le 3.3$, with 70% of the galaxies requiring $y_d < 0.25 M_\odot$. Three galaxies (REBELS-12, 14, 39) require $y_d > 1 M_\odot$. With the SFR predicted by the model and a MW extinction curve, REBELS galaxies detected in [CII] nicely follow the local $L_{\rm CII}-$SFR relation, and are approximately located on the Kennicutt-Schmidt relation. The sample-averaged gas depletion time is of $0.11\, y_P^{-2}$ Gyr, where $y_P$ is the ratio of the gas-to-stellar distribution radius. For some systems a solution simultaneously matching the observed ($\beta, F_{\rm UV}, F_{158}$) values cannot be found. This occurs when the index $I_m = (F_{158}/F_{\rm UV})/(\beta-\beta_{\rm int})$, where $\beta_{\rm int}$ is the intrinsic UV slope, exceeds $I_m^*\approx 1120$ for a MW curve. For these objects we argue that the FIR and UV emitting regions are not co-spatial, questioning the use of the IRX-$\beta$ relation., Comment: 22 pages, 7 figures, accepted by MNRAS. Comments welcome

We study the members of the protocluster around AzTEC3 submillimeter galaxy at z=5.3. We analyzed the data from the MUSE instrument in an area of 1.4x1.4 arcmin^2 around AzTEC3 and derived information on the Lya line in emission. We compared the Lya profile of various regions of the environment with the zELDA radiative transfer model, revealing the neutral gas distribution and kinematics. We identified 10 Lya emitting sources, including 2 regions with extended emission: one embedding AzTEC3 and LBG3, a star-forming galaxy located 12 kpc north of the SMG and another toward LBG-1, a star-forming galaxy located 90 kpc to the southeast. The sources appear distributed in an elongated configuration of about 70'' in extent. The number of sources confirms the overdensity around AzTEC3. For the AzTEC3+LBG3 system, the Lya emission appears redshifted and more spatially extended than the [CII] line emission. Similarly, the Lya line spectrum is broader in velocity than [CII] for LBG1. In the former spectrum, the Lya emission is elongated to the north of LBG3 and to the south of AzTEC3, where a faint Lya emitting galaxy is also located. The elongated structures could resemble tidal features due to the interaction of the two galaxies with AzTEC3. Also, we find a bridge of gas, revealed by the Lya emission between AzTEC3 and LBG3. The Lya emission toward LBG1 embeds its three components. The HI kinematics support the idea of a merger of the three components. Given the availability of CO and [CII] observations from previous campaigns, and our Lya information, we find evidence of starburst-driven phenomena and interactions around AzTEC-3. The stellar mass of the galaxies of the overdensity and the Lya luminosity of the HI nebula associated with AzTEC-3 imply a dark matter halo of 10^12 Msun at z=5.3 that could evolve into a cluster of 2x10^14 Msun at z=0., Comment: Submitted to A&A, Accepted on Feb 4th, 2022

ALMA observations have revealed the presence of dust in the first generations of galaxies in the Universe. However, the dust temperature $T_d$ remains mostly unconstrained due to the few available FIR continuum data at redshift $z>5$. This introduces large uncertainties in several properties of high-$z$ galaxies, namely their dust masses, infrared luminosities, and obscured fraction of star formation. Using a new method based on simultaneous [CII] 158$\mu$m line and underlying dust continuum measurements, we derive $T_ d$ in the continuum and [CII] detected $z\approx 7$ galaxies in the ALMA Large Project REBELS sample. We find $39\ \mathrm{K} < T_d < 58\ \mathrm{K}$, and dust masses in the narrow range $M_d = (0.9-3.6)\times 10^7 M_{\odot}$. These results allow us to extend for the first time the reported $T_d(z)$ relation into the Epoch of Reionization. We produce a new physical model that explains the increasing $T_ d(z)$ trend with the decrease of gas depletion time, $t_{dep}=M_g/\mathrm{SFR}$, induced by the higher cosmological accretion rate at early times; this hypothesis yields $T_d \propto (1+z)^{0.4}$. The model also explains the observed $T_d$ scatter at a fixed redshift. We find that dust is warmer in obscured sources, as a larger obscuration results in more efficient dust heating. For UV-transparent (obscured) galaxies, $T_d$ only depends on the gas column density (metallicity), $T_d \propto N_H^{1/6}$ ($T_d \propto Z^{-1/6}$). REBELS galaxies are on average relatively transparent, with effective gas column densities around $N_H \simeq (0.03-1)\times 10^{21} \mathrm{cm}^{-2}$. We predict that other high-$z$ galaxies (e.g. MACS0416-Y1, A2744-YD4), with estimated $T_d \gg 60$ K, are significantly obscured, low-metallicity systems. In fact $T_d$ is higher in metal-poor systems due to their smaller dust content, which for fixed $L_{ IR}$ results in warmer temperatures., Comment: Accepted for publication in MNRAS

Over the past decades, rest-frame ultraviolet (UV) observations have provided large samples of UV luminous galaxies at redshift (z) greater than 6, during the so-called epoch of reionization. While a few of these UV identified galaxies revealed significant dust reservoirs, very heavily dust-obscured sources at these early times have remained elusive. They are limited to a rare population of extreme starburst galaxies, and companions of rare quasars. These studies conclude that the contribution of dust-obscured galaxies to the cosmic star formation rate density at $z>6$ is sub-dominant. Recent ALMA and Spitzer observations have identified a more abundant, less extreme population of obscured galaxies at $z=3-6$. However, this population has not been confirmed in the reionization epoch so far. Here, we report the discovery of two dust-obscured star forming galaxies at $z=6.6813\pm0.0005$ and $z=7.3521\pm0.0005$. These objects are not detected in existing rest-frame UV data, and were only discovered through their far-infrared [CII] lines and dust continuum emission as companions to typical UV-luminous galaxies at the same redshift. The two galaxies exhibit lower infrared luminosities and star-formation rates than extreme starbursts, in line with typical star-forming galaxies at $z\sim7$. This population of heavily dust-obscured galaxies appears to contribute 10-25 per cent to the $z>6$ cosmic star formation rate density., Comment: 16 pages, 7 figures, 4 tables, published in Nature

The Reionization Era Bright Emission Line Survey (REBELS) is a cycle-7 ALMA Large Program (LP) that is identifying and performing a first characterization of many of the most luminous star-forming galaxies known in the z>6.5 universe. REBELS is providing this probe by systematically scanning 40 of the brightest UV-selected galaxies identified over a 7-deg**2 area for bright 158-micron [CII] and 88-micron [OIII] lines and dust-continuum emission. Selection of the 40 REBELS targets was done by combining our own and other photometric selections, each of which is subject to extensive vetting using three completely independent sets of photometry and template-fitting codes. Building on the observational strategy deployed in two pilot programs, we are increasing the number of massive interstellar medium (ISM) reservoirs known at z>6.5 by ~4-5x to >30. In this manuscript, we motivate the observational strategy deployed in the REBELS program and present initial results. Based on the 60.6 hours of ALMA observations taken in the first year of the program (November 2019 to January 2020), 18 highly significant >~7sigma [CII] lines have already been discovered, the bulk of which (13/18) also show >~3.3 sigma dust-continuum emission. These newly discovered lines more than triple the number of bright ISM-cooling lines known in the z>6.5 universe, such that the number of ALMA-derived redshifts at z>6.5 already rival Lya redshift discoveries. An analysis of the completeness of our search results vs. star formation rate (SFR) suggests an ~79% efficiency in scanning for [CII] when the SFR(UV+IR) is in excess of 28 M_sol/yr. These new LP results further demonstrate ALMA's efficiency as a "redshift machine", particularly in the epoch of reionization., Comment: 30 pages, 4 tables, 18 figures, accepted for publication in ApJ

We present Atacama Large Millimetre Array (ALMA) 2mm continuum observations of a complete and unbiased sample of 99 870micron-selected sub-millimeter galaxies (SMGs) in the Extended Chandra Deep Field South (ALESS). Our observations of each SMG reach average sensitivities of 53 microJy/beam. We measure the flux densities for 70 sources, for which we obtain a typical 870micron-to-2mm flux ratio of 14 +/- 5. We do not find a redshift dependence of this flux ratio, which would be expected if the dust emission properties of our SMGs were the same at all redshifts. By combining our ALMA measurements with existing Herschel/SPIRE observations, we construct a (biased) subset of 27 galaxies for which the cool dust emission is sufficiently well sampled to obtain precise constraints on their dust properties using simple isothermal models. Thanks to our new 2mm observations, the dust emissivity index is well-constrained and robust against different dust opacity assumptions. The median dust emissivity index of our SMGs is $\beta\simeq1.9\pm0.4$, consistent with the emissivity index of dust in the Milky Way and other local and high-redshift galaxies, as well as classical dust grain model predictions. We also find a negative correlation between the dust temperature and $\beta$, similar to low-redshift observational and theoretical studies. Our results indicate that $\beta\simeq2$ in high-redshift dusty star-forming galaxies, implying little evolution in dust grain properties between our SMGs and local dusty galaxy samples, and suggesting these high-mass and high-metallicity galaxies have dust reservoirs driven by grain growth in their ISM., Comment: 31 pages, 20 figures, accepted for publication in the Astrophysical Journal

We present the first results from the 2mm Mapping Obscuration to Reionization (MORA) survey, the largest ALMA contiguous blank-field survey to-date with a total area of 184 sq. arcmin and the only at 2mm to search for dusty star-forming galaxies (DSFGs). We use the 13 sources detected above 5sigma to estimate the first ALMA galaxy number counts at this wavelength. These number counts are then combined with the state-of-the-art galaxy number counts at 1.2mm and 3mm and with a backward evolution model to place constraints on the evolution of the IR luminosity function and dust-obscured star formation in the last 13 billion years. Our results suggest a steep redshift evolution on the space density of DSFGs and confirm the flattening of the IR luminosity function at faint luminosities, with a slope of $\alpha_{LF} = -0.42^{+0.02}_{-0.04}$. We conclude that the dust-obscured component, which peaks at z=2-2.5, has dominated the cosmic history of star formation for the past ~12 billion years, back to z~4. At z=5, the dust-obscured star formation is estimated to be ~35% of the total star formation rate density and decreases to 25%-20% at z=6-7, implying a minor contribution of dust-enshrouded star formation in the first billion years of the Universe. With the dust-obscured star formation history constrained up to the end of the epoch of reionization, our results provide a benchmark to test galaxy formation models, to study the galaxy mass assembly history, and to understand the dust and metal enrichment of the Universe at early times., Comment: Accepted for publication in The Astrophysical Journal

Extreme emission-line galaxies (EELGs) at redshift z=1-2 provide a unique view of metal-poor, starburst sources that are the likely drivers of the cosmic reionization at z$\geq6$. However, the molecular gas reservoirs of EELGs - the fuel for their intense star-formation - remain beyond the reach of current facilities. We present ALMA [CII] and PdBI CO(2-1) observations of a z=1.8, strongly lensed EELG SL2S 0217, a bright Lyman-$\alpha$ emitter with a metallicity 0.05 $Z_\odot$. We obtain a tentative (3-4$\sigma$) detection of the [CII] line and set an upper limit on the [CII]/SFR ratio of $\leq1\times10^6$ $L_\odot$/($M_\odot$ yr$^{-1}$), based on the synthesized images and visibility-plane analysis. The CO(2-1) emission is not detected. Photoionization modelling indicates that up to 80% of the [CII] emission originates from neutral or molecular gas, although we can not rule out that the gas is fully ionized. The very faint [CII] emission is in line with both nearby metal-poor dwarfs and high-redshift Lyman-$\alpha$ emitters, and predictions from hydrodynamical simulations. However, the [CII] line is 30$\times$ fainter than predicted by the De Looze et al. [CII]-SFR relation for local dwarfs, illustrating the danger of extrapolating locally-calibrated relations to high-redshift, metal-poor galaxies., Comment: Corrected to match the published version. We thank Y. Harikane and I. Smail for their comments on the previous arXiv version

We study the radio properties of 706 sub-millimeter galaxies (SMGs) selected at 870$\mu$m with the Atacama Large Millimeter Array from the SCUBA-2 Cosmology Legacy Survey map of the Ultra Deep Survey field. We detect 273 SMGs at $>4\sigma$ in deep Karl G. Jansky Very Large Array 1.4 GHz observations, of which a subset of 45 SMGs are additionally detected in 610 MHz Giant Metre-Wave Radio Telescope imaging. We quantify the far-infrared/radio correlation through parameter $q_\text{IR}$, defined as the logarithmic ratio of the far-infrared and radio luminosity, and include the radio-undetected SMGs through a stacking analysis. We determine a median $q_\text{IR} = 2.20\pm0.03$ for the full sample, independent of redshift, which places these $z\sim2.5$ dusty star-forming galaxies $0.44\pm0.04$ dex below the local correlation for both normal star-forming galaxies and local ultra-luminous infrared galaxies (ULIRGs). Both the lack of redshift-evolution and the offset from the local correlation are likely the result of the different physical conditions in high-redshift starburst galaxies, compared to local star-forming sources. We explain the offset through a combination of strong magnetic fields ($B\gtrsim0.2$mG), high interstellar medium (ISM) densities and additional radio emission generated by secondary cosmic rays. While local ULIRGs are likely to have similar magnetic field strengths, we find that their compactness, in combination with a higher ISM density compared to SMGs, naturally explains why local and high-redshift dusty star-forming galaxies follow a different far-infrared/radio correlation. Overall, our findings paint SMGs as a homogeneous population of galaxies, as illustrated by their tight and non-evolving far-infrared/radio correlation., Comment: Accepted for publication in ApJ; 18 pages + appendices; 7 figures in main text

We present a survey of the molecular gas in 61 submillimetre galaxies (SMGs) selected from 870$\mu$m continuum surveys of the COSMOS, UDS and ECDFS fields, using the Atacama Large Millimeter Array (ALMA) and the Northern Extended Millimeter Array (NOEMA). 46 $^{12}$CO ($J=$2-5) emission lines are detected in 45 of the targets at $z=$1.2-4.8, with redshifts indicating that those which are submillimetre bright and undetected/faint in the optical/near-infrared typically lie at higher redshifts, with a gradient of $\Delta z/\Delta S_{870}=$0.11$\pm$0.04mJy$^{-1}$. We also supplement our data with literature sources to construct a statistical CO spectral line energy distribution and find the $^{12}$CO line luminosities in SMGs peak at $J_{\rm up}\sim$6, consistent with the Cosmic Eyelash, among similar studies. Our SMGs lie mostly on or just above the main sequence, displaying a decrease in their gas depletion timescales $t_{\rm dep} = M_{\rm gas}/{\rm SFR}$ with redshift in the range $z\sim$1-5 and a median of 200$\pm$50Myr at $z\sim$2.8. This coincides with an increase in molecular gas fraction $\mu_{\rm gas} = M_{\rm gas}/M_\ast$ across the same redshift range. Finally we demonstrate that the $M_{\rm baryon}$-$\sigma$ distribution of our SMGs is consistent with that followed by early-type galaxies in the Coma cluster, providing strong support to the suggestion that SMGs are progenitors of massive local spheroidal galaxies. On the basis of this we suggest that the SMG populations above and below an 870-$\mu$m flux limit of $S_{870}\sim$5mJy may correspond to the division between slow- and fast-rotators seen in local early-type galaxies., Comment: Submitted to MNRAS

We update the spectral modeling code MAGPHYS to include a 2175\AA\ absorption feature in its UV-to-near-IR dust attenuation prescription. This allows us to determine the strength of this feature and the shape of the dust attenuation curve in ~5000 star-forming galaxies at 0.1

We present a new investigation of the dependence of H2O maser detection rates and properties on the mid-IR AGN luminosity, L_AGN, and the obscuring column density, N_H, based on mid-IR and hard X-ray photometry. Based on spectral energy distribution fitting that allows for decomposition of the black hole accretion and star-formation components in the mid-infrared, we show that the megamaser (disk maser) detection rate increases sharply for galaxies with 12 micron AGN luminosity L^{AGN}_{12 micron} greater than 10^42 erg/s, from ~<3%(~<2%) to ~12%(~5%). By using the ratio of the observed X-ray to mid-IR AGN luminosity as an indicator of N_H, we also find that the megamaser (disk maser) detection rates are boosted to 15%(7%) and 20%(9%) for galaxies with N_H >= 10^23 cm^{-2} and N_H >= 10^{24} cm^{-2}, respectively. Combining these column density cuts with a constraint for high L^{AGN}_{12 micron} (>=10^42 erg/s) predicts further increases in the megamaser (disk maser) detection rates to 19%(8%) and 27%(14%), revealing unprecedented potential boosts of the megamaser and disk maser detection rates by a factor of 7-15 relative to the current rates, depending on the chosen sample selection criteria. A noteworthy aspect of these new predictions is that the completeness rates are only compromised mildly, with the rates remaining at the level of ~95%(~50%) for sources with N_H >= 10^{23} cm^{-2} (N_H >= 10^24 cm^-2). Applying these selection methods to the current X-ray AGN surveys predicts the detection of >~15 new megamaser disks., Comment: 22 pages, 11 figures, 4 tables, Accepted by ApJ

The rich information on (sub)millimeter dust continuum emission from distant galaxies in the public Atacama Large Millimeter/submillimeter Array (ALMA) archive is contained in thousands of inhomogeneous observations from individual PI-led programs. To increase the usability of these data for studies deepening our understanding of galaxy evolution, we have developed automated mining pipelines for the ALMA archive in the COSMOS field (A3COSMOS) that efficiently exploit the available information for large numbers of galaxies across cosmic time and keep the data products in sync with the increasing public ALMA archive: (a) a dedicated ALMA continuum imaging pipeline, (b) two complementary photometry pipelines for both blind source extraction and prior source fitting, (c) a counterpart association pipeline utilizing the multiwavelength data available (including quality assessment based on machine-learning techniques), (d) an assessment of potential (sub)millimeter line contribution to the measured ALMA continuum, and (e) extensive simulations to provide statistical corrections to biases and uncertainties in the ALMA continuum measurements. Application of these tools yields photometry catalogs with ~1000 (sub)millimeter detections (spurious fraction ~8%-12%) from over 1500 individual ALMA continuum images. Combined with ancillary photometric and redshift catalogs and the above quality assessments, we provide robust information on redshift, stellar mass, and star formation rate for ~700 galaxies at redshifts 0.5-6 in the COSMOS field (with undetermined selection function). The ALMA photometric measurements and galaxy properties are released publicly within our blind extraction, prior fitting, and galaxy property catalogs, plus the images. These products will be updated on a regular basis in the future., Comment: 56 pages, 36 figures, 5 tables, published on ApJS. Catalogs are available at journal website in FITS format. Current and future catalogs and continuum images will be released at COSMOS IRSA server https://irsa.ipac.caltech.edu/data/COSMOS/overview.html and A3COSMOS website https://sites.google.com/view/a3cosmos/data

We analyse the physical properties of a large, homogeneously selected sample of ALMA-located sub-mm galaxies (SMGs). This survey, AS2UDS, identified 707 SMGs across the ~1 sq.deg. field, including ~17 per cent, which are undetected at $K$>~25.7 mag. We interpret their ultraviolet-to-radio data using MAGPHYS and determine a median redshift of z=2.61+-0.08 (1$\sigma$ range of z=1.8-3.4) with just ~6 per cent at z>4. Our survey provides a sample of massive dusty galaxies at z>~1, with median dust and stellar masses of $M_d$=(6.8+-0.3)x10$^{8}$M$_\odot$ (thus, gas masses of ~10$^{11}$M$_\odot$) and $M_\ast=$(1.26+-0.05)x10$^{11}$M$_\odot$. We find no evolution in dust temperature at a constant far-infrared luminosity across z~1.5-4. The gas mass function of our sample increases to z~2-3 and then declines at z>3. The space density and masses of SMGs suggest that almost all galaxies with $M_\ast$>~3x10$^{11}$M$_\odot$ have passed through an SMG-like phase. The redshift distribution is well fit by a model combining evolution of the gas fraction in halos with the growth of halo mass past a threshold of $M_h$~6x10$^{12}$M$_\odot$, thus SMGs may represent the highly efficient collapse of gas-rich massive halos. We show that SMGs are broadly consistent with simple homologous systems in the far-infrared, consistent with a centrally illuminated starburst. Our study provides strong support for an evolutionary link between the active, gas-rich SMG population at z>1 and the formation of massive, bulge-dominated galaxies across the history of the Universe., Comment: 33 pages, 16 figures, MNRAS accepted version

We present an enhanced version of the multiwavelength spectral modeling code MAGPHYS that allows the estimation of galaxy photometric redshift and physical properties (e.g., stellar mass, star formation rate, dust attenuation) simultaneously, together with robust characterization of their uncertainties. The self-consistent modeling over ultraviolet to radio wavelengths in MAGPHYS+photo-z is unique compared to standard photometric redshift codes. The broader wavelength consideration is particularly useful for breaking certain degeneracies in color vs. redshift for dusty galaxies with limited observer-frame ultraviolet and optical data (or upper limits). We demonstrate the success of the code in estimating redshifts and physical properties for over 4,000 infrared-detected galaxies at 0.4

We present the catalogue and properties of sources in AS2UDS, an 870-$\mu$m continuum survey with the Atacama Large Millimetre/sub-millimetre Array (ALMA) of 716 single-dish sub-millimetre sources detected in the UKIDSS/UDS field by the SCUBA-2 Cosmology Legacy Survey. In our sensitive ALMA follow-up observations we detect 708 sub-millimetre galaxies (SMGs) at $>$\,4.3$\sigma$ significance across the $\sim$\,1-degree diameter field. We combine our precise ALMA positions with the extensive multi-wavelength coverage in the UDS field to fit the spectral energy distributions of our SMGs to derive a median redshift of $z_{\rm phot}=$\,2.61$\pm$0.09. This large sample reveals a statistically significant trend of increasing sub-millimetre flux with redshift suggestive of galaxy downsizing. 101 ALMA maps do not show a $>$\,4.3$\sigma$ SMG, but we demonstrate from stacking {\it Herschel} SPIRE observations at these positions, that the vast majority of these blank maps correspond to real single-dish sub-millimetre sources. We further show that these blank maps contain an excess of galaxies at $z_{\rm phot}=$\,1.5--4 compared to random fields, similar to the redshift range of the ALMA-detected SMGs. In addition, we combine X-ray and mid-infrared active galaxy nuclei activity (AGN) indicators to yield a likely range for the AGN fraction of 8--28\,\% in our sample. Finally, we compare the redshifts of this population of high-redshift, strongly star-forming galaxies with the inferred formation redshifts of massive, passive galaxies being found out to $z\sim$\,2, finding reasonable agreement -- in support of an evolutionary connection between these two classes of massive galaxy., Comment: Submitted to MNRAS; Comments welcome. Full catalogue will be made publicly available on acceptance of paper

We present 0.15-arcsec (1 kpc) resolution ALMA observations of the [CII] 157.74 um line and rest-frame 160-um continuum emission in two z~3 dusty, star-forming galaxies - ALESS 49.1 and ALESS 57.1, combined with resolved CO(3-2) observations. In both sources, the [CII] surface brightness distribution is dominated by a compact core $\leq$1 kpc in radius, a factor of 2-3 smaller than the extent of the CO(3-2) emission. In ALESS 49.1, we find an additional extended (8-kpc radius), low surface-brightness [CII] component. Based on an analysis of mock ALMA observations, the [CII] and 160-um continuum surface brightness distributions are inconsistent with a single-Gaussian surface brightness distribution with the same size as the CO(3-2) emission. The [CII] rotation curves flatten at $\simeq$2 kpc radius, suggesting the kinematics of the central regions are dominated by a baryonic disc. Both galaxies exhibit a strong [CII]/FIR deficit on 1-kpc scales, with FIR-surface-brightness to [CII]/FIR slope steeper than in local star-forming galaxies. A comparison of the [CII]/CO(3-2) observations with PDR models suggests a strong FUV radiation field ($G_0\sim10^4$) and high gas density ($n\mathrm{(H)}\sim10^4-10^5$ cm$^{-3}$) in the central regions of ALESS 49.1 and 57.1. The most direct interpretation of the pronounced [CII]/FIR deficit is a thermal saturation of the C+ fine-structure levels at temperatures $\geq$500 K, driven by the strong FUV field., Comment: Submitted to ApJ

We carry out a blind search of 3mm continuum sources using the ALMA Science Archive to derive the first galaxy number counts at this wavelength. The analyzed data are drawn from observations towards three extragalactic legacy fields: COSMOS, CDF-S, and the UDS comprising more than 130 individual ALMA Band 3 pointings and an effective survey area of ~200 sq. arcmin with a continuum sensitivity that allows for the direct detection of unlensed Dusty Star-Forming Galaxies (DSFGs) dust emission beyond the epoch of reionization. We present a catalog of 16 sources detected at >5sigma with flux densities S_3mm = 60 - 600 uJy from which number counts are derived. These number counts are then used to place constraints on the volume density of DSFGs with an empirical backward evolution model. Our measured 3mm number counts indicate that the contribution of DSFGs to the cosmic star formation rate density at z >~ 4 is non-negligible. This is contrary to the generally adopted assumption of a sharply decreasing contribution of obscured galaxies at z > 4 as inferred by optical and near-infrared surveys. This work demonstrates the power of ALMA 3mm observations which can reach outstanding continuum sensitivities during typical spectral line science programs. Further constraints on 3mm-selected galaxies will be essential to refine models of galaxy formation and evolution as well as models of early Universe dust production mechanisms., Comment: Accepted for in The Astrophysical Journal

We present sub-kpc-scale mapping of the 870 $\mu$m ALMA continuum emission in six luminous ($L_{\rm IR}~\sim~5~\times10^{12}$ L$_{\odot}$) submillimeter galaxies (SMGs) from the ALESS survey of the Extended Chandra Deep Field South. Our high-fidelity 0.07$''$-resolution imaging ($\sim$500 pc) reveals robust evidence for structures with deconvolved sizes of $\lesssim$0.5-1 kpc embedded within (dominant) exponential dust disks. The large-scale morphologies of the structures within some of the galaxies show clear curvature and/or clump-like structures bracketing elongated nuclear emission, suggestive of bars, star-forming rings, and spiral arms. In this interpretation, the ratio of the `ring' and `bar' radii (1.9$\pm$0.3) agrees with that measured for such features in local galaxies. These potential spiral/ring/bar structures would be consistent with the idea of tidal disturbances, with their detailed properties implying flat inner rotation curves and Toomre-unstable disks (Q<1). The inferred one-dimensional velocity dispersions ($\sigma_{\rm r}\lesssim$ 70-160 km s$^{-1}$) are marginally consistent with the limits implied if the sizes of the largest structures are comparable to the Jeans length. We create maps of the star formation rate density ($\Sigma_{\rm SFR}$) on $\sim$500 pc scales and show that the SMGs are able to sustain a given (galaxy-averaged) $\Sigma_{\rm SFR}$ over much larger physical scales than local (ultra-)luminous infrared galaxies. However, on 500 pc scales, they do not exceed the Eddington limit set by radiation pressure on dust. If confirmed by kinematics, the potential presence of non-axisymmetric structures would provide a means for net angular momentum loss and efficient star formation, helping to explain the very high star formation rates measured in SMGs., Comment: 18 pages, 6 figures, Accepted to ApJ

We present a multiwavelength study of an atypical submillimeter galaxy in the GOODS-North field, with the aim to understand its physical properties of stellar and dust emission, as well as the central AGN activity. Although it is shown that the source is likely an extremely dusty galaxy at high redshift, its exact position of submillimeter emission is unknown. With the new NOEMA interferometric imaging, we confirm that the source is a unique dusty galaxy. It has no obvious counterpart in the optical and even NIR images observed with HST at lambda~<1.4um. Photometric-redshift analyses from both stellar and dust SED suggest it to likely be at z~>4, though a lower redshift at z~>3.1 cannot be fully ruled out (at 90% confidence interval). Explaining its unusual optical-to-NIR properties requires an old stellar population (~0.67 Gyr), coexisting with a very dusty ongoing starburst component. The latter is contributing to the FIR emission, with its rest-frame UV and optical light being largely obscured along our line of sight. If the observed fluxes at the rest-frame optical/NIR wavelengths were mainly contributed by old stars, a total stellar mass of ~3.5x10^11Msun would be obtained. An X-ray spectral analysis suggests that this galaxy harbors a heavily obscured AGN with N_H=3.3x10^23 cm^-2 and an intrinsic 2-10 keV luminosity of L_X~2.6x10^44 erg/s, which places this object among distant type 2 quasars. The radio emission of the source is extremely bright, which is an order of magnitude higher than the star-formation-powered emission, making it one of the most distant radio-luminous dusty galaxies. The combined characteristics of the galaxy suggest that the source appears to have been caught in a rare but critical transition stage in the evolution of submillimeter galaxies, where we are witnessing the birth of a young AGN and possibly the earliest stage of its jet formation and feedback., Comment: 13 pages in printer format, 10 figures, 1 table, accepted for publication in the A&A

The thermodynamic aspects of the polymorphic phase transition from {\alpha}-indomethacin to {\gamma}-indomethacin are the fundamental key to find the most bioavailable phase of indomethacin. In the present work, varying the temperature and solvent permittivity changes the polymorphic transitions. Hence, the thermodynamic properties such as enthalpy, Gibbs free energy, and entropy of both indomethacin polymorphs are determined in terms of the solvent permittivity as functions of indomethacin's temperature in a vacuum, which are crucially related to the stability, spontaneity, and reversibility of the polymorphic transformation., Comment: 17 pages, 5 figures

We use ALMA observations of four sub-millimetre galaxies (SMGs) at $z\sim2-3$ to investigate the spatially resolved properties of the inter-stellar medium (ISM) at scales of 1--5 kpc (0.1--0.6$''$). The velocity fields of our sources, traced by the $^{12}$CO($J$=3-2) emission, are consistent with disk rotation to first order, implying average dynamical masses of $\sim$3$\times10^{11}$M$_{\odot}$ within two half-light radii. Through a Bayesian approach we investigate the uncertainties inherent to dynamically constraining total gas masses. We explore the covariance between the stellar mass-to-light ratio and CO-to-H$_{2}$ conversion factor, $\alpha_{\rm CO}$, finding values of $\alpha_{\rm CO}=1.1^{+0.8}_{-0.7}$ for dark matter fractions of 15 \%. We show that the resolved spatial distribution of the gas and dust continuum can be uncorrelated to the stellar emission, challenging energy balance assumptions in global SED fitting. Through a stacking analysis of the resolved radial profiles of the CO(3-2), stellar and dust continuum emission in SMG samples, we find that the cool molecular gas emission in these sources (radii $\sim$5--14 kpc) is clearly more extended than the rest-frame $\sim$250 $\mu$m dust continuum by a factor $>2$. We propose that assuming a constant dust-to-gas ratio, this apparent difference in sizes can be explained by temperature and optical-depth gradients alone. Our results suggest that caution must be exercised when extrapolating morphological properties of dust continuum observations to conclusions about the molecular gas phase of the ISM., Comment: Submitted to ApJ, 16 pages, 6 figures

Aims. We aim to investigate if the environment (characterised by the host dark matter halo mass) plays any role in shaping the galaxy star formation main sequence (MS). Methods. The Galaxy and Mass Assembly project (GAMA) combines a spectroscopic survey with photometric information in 21 bands from the far-ultraviolet (FUV) to the far-infrared (FIR). Stellar masses and dust-corrected star-formation rates (SFR) are derived from spectral energy distribution (SED) modelling using MAGPHYS. We use the GAMA galaxy group catalogue to examine the variation of the fraction of star-forming galaxies (SFG) and properties of the MS with respect to the environment. Results. We examine the environmental dependence for stellar mass selected samples without preselecting star-forming galaxies and study any dependence on the host halo mass separately for centrals and satellites out to z ~ 0.3. We find the SFR distribution at fixed stellar mass can be described by the combination of two Gaussians (referred to as the star-forming Gaussian and the quiescent Gaussian). Using the observed bimodality to define SFG, we investigate how the fraction of SFG F(SFG) and properties of the MS change with environment. For centrals, the position of the MS is similar to the field but with a larger scatter. No significant dependence on halo mass is observed. For satellites, the position of the MS is almost always lower (by ~0.2 dex) compared to the field and the width is almost always larger. F(SFG) is similar between centrals (in different halo mass bins) and field galaxies. However, for satellites F(SFG) decreases with increasing halo mass and this dependence is stronger towards lower redshift., Comment: 6 pages, 6 figures, accepted for publication in A&A