Search

We present the results of a {\it James Webb Space Telescope} NIRCam investigation into the young massive star cluster (YMC) population in the luminous infrared galaxy VV 114. We identify 374 compact YMC candidates with a $S/N \geq 3$, 5, and 5 at F150W, F200W, and F356W respectively. A direct comparison with our {\it HST} cluster catalog reveals that $\sim 20\%$ of these sources are undetected at optical wavelengths. Based on {\it yggdrasil} stellar population models, we identify 17 YMC candidates in our {\it JWST} imaging alone with F150W-F200W and F200W-F356W colors suggesting they are all very young, dusty ($A_{V} = 5 - 15$), and massive ($10^{5.8} < M_{\odot} < 10^{6.1}$). The discovery of these `hidden' sources, many of which are found in the `overlap' region between the two nuclei, quadruples the number of $t < 3$ Myr clusters, and nearly doubles the number of $t < 6$ Myr clusters detected in VV 114. Now extending the cluster age distribution ($dN/d\tau \propto \tau^{\gamma}$) to the youngest ages, we find a slope of $\gamma = -1.30 \pm 0.39$ for $10^{6} < \tau (\mathrm{yr}) < 10^{7}$, which is consistent with the previously determined value from $10^{7} < \tau (\mathrm{yr}) < 10^{8.5}$, and confirms that VV 114 has a steep age distribution slope for all massive star clusters across the entire range of cluster ages observed. Finally, the consistency between our {\it JWST}- and {\it HST}-derived age distribution slopes indicates that the balance between cluster formation and destruction has not been significantly altered in VV 114 over the last 0.5 Gyr., Comment: 13 pages, 4 figures, 1 table, published in the Astrophysical Journal

Observations and simulations of interacting galaxies and mergers in the local universe have shown that interactions can significantly enhance the star formation rates (SFR) and fueling of Active Galactic Nuclei (AGN). However, at higher redshift, some simulations suggest that the level of star formation enhancement induced by interactions is lower due to the higher gas fractions and already increased SFRs in these galaxies. To test this, we measure the SFR enhancement in a total of 2351 (1327) massive ($M_*>10^{10}M_\odot$) major ($1, Comment: 23 pages, 13 figures, Accepted for publication in ApJ

We present the results of a stacking analysis performed on Spitzer/Infrared Spectrograph high-resolution mid-infrared spectra of luminous infrared galaxies (LIRGs) in the Great Observatories All-Sky LIRG Survey (GOALS). By binning on mid-infrared active galactic nucleus (AGN) fraction and stacking spectra, we detect bright emission lines [Ne II] and [Ne III], which trace star formation, and fainter emission lines [Ne V] and [O IV], which trace AGN activity, throughout the sample. We find the [Ne II] luminosity is fairly constant across all AGN fraction bins, while the [O IV] and [Ne V] luminosities increase by over an order of magnitude. Our measured average line ratios, [Ne V]/[Ne II] and [O IV]/[Ne II], at low AGN fraction are similar to H II galaxies while the line ratios at high AGN fraction are similar to LINERs and Seyferts. We decompose the [O IV] luminosity into star-formation and AGN components by fitting the [O IV] luminosity as a function of the [Ne II] luminosity and the mid-infrared AGN fraction. The [O IV] luminosity in LIRGs is dominated by star formation for mid-infrared AGN fractions $\lesssim0.3$. With the corrected [O IV] luminosity, we calculate black hole accretion rates ranging from $10^{-5}$ M$_{\odot}$/yr at low AGN fractions to 0.2 M$_{\odot}$/yr at the highest AGN fractions. We find that using the [O IV] luminosity, without correcting for star formation, can lead to an overestimate of the BHAR by up to a factor of 30 in starburst dominated LIRGs. Finally, we show the BHAR/SFR ratio increases by more than three orders of magnitude as a function of mid-infrared AGN fraction in LIRGs., 17 pages, 15 figures. Accepted to ApJ

While low-luminosity galaxies dominate number counts at all redshifts, their contribution to cosmic reionization is poorly understood due to a lack of knowledge of their physical properties. We isolate a sample of 35 z ≈ 4–5 continuum-faint Lyman-α emitters from deep VLT/MUSE spectroscopy and directly measure their H α emission using stacked Spitzer/IRAC Ch. 1 photometry. Based on Hubble Space Telescope imaging, we determine that the average UV continuum magnitude is fainter than −16 (≈ 0.01 L⋆), implying a median Lyman-α equivalent width of 259 Å. By combining the H α measurement with the UV magnitude, we determine the ionizing photon production efficiency, ξion, a first for such faint galaxies. The measurement of log10 (ξion [Hz erg−1]) = 26.28 (+0.28−0.40) is in excess of literature measurements of both continuum- and emission line-selected samples, implying a more efficient production of ionizing photons in these lower luminosity, Lyman-α-selected systems. We conclude that this elevated efficiency can be explained by stellar populations with metallicities between 4 × 10−4 and 0.008, with light-weighted ages less than 3 Myr., Monthly Notices of the Royal Astronomical Society, 493 (4), ISSN:0035-8711, ISSN:1365-2966, ISSN:1365-8711

We present James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI) integral-field spectroscopy of the nearby merging, luminous infrared galaxy, NGC 7469. This galaxy hosts a Seyfert type-1.5 nucleus, a highly ionized outflow, and a bright, circumnuclear star-forming ring, making it an ideal target to study active galactic nucleus (AGN) feedback in the local universe. We take advantage of the high spatial/spectral resolution of JWST/MIRI to isolate the star-forming regions surrounding the central active nucleus and study the properties of the dust and warm molecular gas on ∼100 pc scales. The starburst ring exhibits prominent polycyclic aromatic hydrocarbon (PAH) emission, with grain sizes and ionization states varying by only ∼30%, and a total star formation rate of 10–30 M⊙ yr−1 derived from fine structure and recombination emission lines. Using pure rotational lines of H2 we detect 1.2 × 107 M⊙ of warm molecular gas at a temperature higher than 200 K in the ring. All PAH bands get significantly weaker toward the central source, where larger and possibly more ionized grains dominate the emission, likely the result of the ionizing radiation and/or the fast wind emerging from the AGN. The small grains and warm molecular gas in the bright regions of the ring however display properties consistent with normal star-forming regions. These observations highlight the power of JWST to probe the inner regions of dusty, rapidly evolving galaxies for signatures of feedback and inform models that seek to explain the coevolution of supermassive black holes and their hosts., This work is based on observations made with the NASA/ESA/CSA JWST. T.S.-Y.L. acknowledges funding support from NASA grant JWST-ERS-01328. The data were obtained from the Mikulski Archive for Space Telescopes at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-03127 for JWST. These observations are associated with program #1328 and can be accessed via doi:10.17909/0fe2-cf33. V.U. acknowledges funding support from NASA Astrophysics Data Analysis Program (ADAP) grant 80NSSC20K0450. The Flatiron Institute is supported by the Simons Foundation. H.I. and T.B. acknowledge support from JSPS KAKENHI grant No. JP21H01129 and the Ito Foundation for Promotion of Science. A.M.M. acknowledges support from the National Science Foundation under grant No. 2009416. A.S.E. and S.L. acknowledge support from NASA grant HST-GO15472 and HST-GO16914. Y.S. was funded in part by the NSF through the Grote Reer Fellowship Program administered by Associated Universities, Inc./National Radio Astronomy Observatory. S.A. gratefully acknowledges support from an ERC Advanced grant 789410, from the Swedish Research Council and from the Knut and Alice Wallenberg (KAW) Foundation. F.K. acknowledges support from the Spanish program Unidad de Excelencia María de Maeztu CEX2020-001058-M, financed by MCIN/AEI/10.13039/501100011033. K.I. acknowledges support by the Spanish MCIN under grant PID2019-105510GB-C33/AEI/10.13039/501100011033. F.M-S. acknowledges support from NASA through ADAP award 80NSSC19K1096. Finally, this research has made use of the NASA/IPAC Extragalactic Database (NED), which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration.

Recent work has shown that UV-luminous reionization-era galaxies often exhibit strong Lyman-alpha emission despite being situated at redshifts where the IGM is thought to be substantially neutral. It has been argued that this enhanced Ly α transmission reflects the presence of massive galaxies in overdense regions which power large ionized bubbles. An alternative explanation is that massive galaxies shift more of their Ly α profile to large velocities (relative to the systemic redshift) where the IGM damping wing absorption is reduced. Such a mass-dependent trend is seen at lower redshifts, but whether one exists at z ∼ 7 remains unclear owing to the small number of existing systemic redshift measurements in the reionization era. This is now changing with the emergence of [C ii]-based redshifts from ALMA. Here, we report MMT/Binospec Ly α spectroscopy of eight UV-bright (MUV ∼ −22) galaxies at z ≃ 7 selected from the ALMA REBELS survey. We detect Ly α in four of eight galaxies and use the [C ii] systemic redshifts to investigate the Ly α velocity profiles. The Ly α lines are significantly redshifted from systemic (average velocity offset = 223 km s–1) and broad (FWHM ≈ 300–650 km s−1), with two sources showing emission extending to ≈750 km s−1. We find that the broadest Ly α profiles are associated with the largest [C ii] line widths, suggesting a potential link between the Ly α FWHM and the dynamical mass. Since Ly α photons at high velocities transmit efficiently through the z = 7 IGM, our data suggest that velocity profiles play a significant role in boosting the Ly α visibility of the most UV-luminous reionization-era galaxies.

We present specific star-formation rates for 40 UV-bright galaxies at $z\sim7-8$ observed as part of the Reionization Era Bright Emission Line Survey (REBELS) ALMA large program. The sSFRs are derived using improved measures of SFR and stellar masses, made possible by measurements of far-infrared (FIR) continuum emission and [CII]-based spectroscopic redshifts. For each source in the sample, we derive stellar masses from SED fitting and total SFRs from calibrations of the UV and FIR emission. The median sSFR is $18_{-5}^{+7}$ Gyr$^{-1}$, significantly larger than literature measurements lacking constraints in the FIR. The increase in sSFR reflects the larger obscured SFRs we derive from the dust continuum relative to that implied by the UV+optical SED. We suggest that such differences may reflect spatial variations in dust across these luminous galaxies, with the component dominating the FIR distinct from that dominating the UV. We demonstrate that the inferred stellar masses (and hence sSFRs) are strongly-dependent on the assumed star formation history (SFH) in reionization-era galaxies. When large sSFR galaxies are modeled with non-parametric SFHs, the derived stellar masses can increase by an order of magnitude relative to constant star formation models, owing to the presence of a significant old stellar population that is outshined by the recent burst. The [CII] line widths in the largest sSFR systems are often very broad, suggesting dynamical masses that are easily able to accommodate the dominant old stellar population suggested by non-parametric models. Regardless of these systematic uncertainties in the derived parameters, we find that the sSFR increases rapidly toward higher redshifts for massive galaxies ($9.6, 18 pages, 10 figures, submitted to MNRAS

Cosmic dust is an essential component shaping both the evolution of galaxies and their observational signatures. How quickly dust builds up in the early Universe remains an open question that requires deep observations at (sub-)millimeter wavelengths to resolve. Here we use Atacama Large Millimeter Array observations of 45 galaxies from the Reionization Era Bright Emission Line Survey (REBELS) and its pilot programs, designed to target [CII] and dust emission in UV-selected galaxies at $z\sim7$, to investigate the dust content of high-redshift galaxies through a stacking analysis. We find that the typical fraction of obscured star formation $f_\mathrm{obs} = \mathrm{SFR}_\mathrm{IR} / \mathrm{SFR}_\mathrm{UV + IR}$ depends on stellar mass, similar to what is observed at lower redshift, and ranges from $f_\mathrm{obs} \approx 0.3 - 0.6$ for galaxies with $\log_{10}\left(M_\star / M_\odot\right) = 9.4 - 10.4$. We further adopt the $z\sim7$ stellar mass function from the literature to extract the obscured cosmic star formation rate density (SFRD) from the REBELS survey. Our results suggest only a modest decrease in the SFRD between $3\lesssim z \lesssim 7$, with dust-obscured star formation still contributing $\sim30\%$ at $z\sim7$. While we extensively discuss potential caveats, our analysis highlights the continued importance of dust-obscured star formation even well into the epoch of reionization., Comment: MNRAS Accepted (01 Nov 2022)

Ricci, C., et al., The merger of two or more galaxies can enhance the inflow of material from galactic scales into the close environments of active galactic nuclei (AGNs), obscuring and feeding the supermassive black hole (SMBH). Both recent simulations and observations of AGN in mergers have confirmed that mergers are related to strong nuclear obscuration. However, it is still unclear how AGN obscuration evolves in the last phases of the merger process. We study a sample of 60 luminous and ultra-luminous IR galaxies (U/LIRGs) from the GOALS sample observed by NuSTAR. We find that the fraction of AGNs that are Compton thick (CT; NH ≥ 1024 cm-2) peaks at 74+14 −19 per cent at a late merger stage, prior to coalescence, when the nuclei have projected separations (dsep) of 0.4-6 kpc. A similar peak is also observed in the median NH [(1.6 ± 0.5) × 1024 cm−2]. The vast majority (85+7 −9 per cent) of the AGNs in the final merger stages (dsep ≤ 10 kpc) are heavily obscured (NH ≥ 1023 cm−2), and the median NH of the accreting SMBHs in our sample is systematically higher than that of local hard X-ray-selected AGN, regardless of the merger stage. This implies that these objects have very obscured nuclear environments, with the NH ≥ 1023 cm−2 gas almost completely covering the AGN in late mergers. CT AGNs tend to have systematically higher absorption-corrected X-ray luminosities than less obscured sources. This could either be due to an evolutionary effect, with more obscured sources accreting more rapidly because they have more gas available in their surroundings, or to a selection bias. The latter scenario would imply that we are still missing a large fraction of heavily obscured, lower luminosity (L2−10 ≤ 1043 erg s−1) AGNs in U/LIRGs., LCH was supported by the National Key Research and Development Program of China (2016YFA0400702) and the National Science Foundation of China (11721303 and 11991052). CR acknowledges support from the Fondecyt Iniciacion grant 11190831. ET acknowledges support from CATA-Basal AFB-170002, FONDECYT Regular grant 1190818, ANID Anillo ACT172033, and Millennium Nucleus NCN19_058 (TITANs). FEB acknowledges support from ANID – Millennium Science Initiative Program – ICN12_009, CATA-Basal – AFB-170002, and FONDECYT Regular – 1190818 and 1200495. SA gratefully acknowledges funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No 789410). VU acknowledges funding support from NASA Astrophysics Data Analysis Program (ADAP) Grant 80NSSC20K0450. AMM acknowledges support from the National Science Foundation under grant number 2009416. KI acknowledges support by the Spanish MICINN under grant Proyecto/AEI/10.13039/501100011033 and ‘Unit of excellence María de Maeztu 2020-2023’ awarded to ICCUB (CEX2019-000918-M). PA acknowledges financial support from ANID Millennium Nucleus NCN19-058 (TITANS) and the Max Planck Society through a Partner Group. HI acknowledges support from JSPS KAKENHI Grant Number JP19K23462.

We study the physical properties of the ionized gas in local disks using the sample of 38 nearby $\sim10^{8.5-11.2}$M$_\odot$ Star-Forming Main Sequence (SFMS) galaxies observed so far as part of the MUSE Atlas of Disks (MAD). Specifically, we use all strong emission lines in the MUSE wavelength range 4650-9300 \AA\ to investigate the resolved ionized gas properties on $\sim$100 pc scales. This spatial resolution enables us to disentangle HII regions from the Diffuse Ionized Gas (DIG) in the computation of gas metallicities and Star Formation Rates (SFRs) of star forming regions. The gas metallicities generally decrease with radius. The metallicity of the HII regions is on average $\sim$0.1 dex higher than that of the DIG, but the metallicity radial gradient in both components is similar. The mean metallicities within the inner galaxy cores correlate with the total stellar mass of the galaxies. On our, Comment: 21 pages + appendices. Accepted for publication in MNRAS

The mid-infrared spectrometer and camera transit spectrometer (MISC-T) is one of the three baseline instruments for Origins Space Telescope (Origins) and provides the capability to assess the habitability of nearby exoplanets and search for signs of life. MISC-T employs a densified pupil optical design, and HgCdTe and Si:As detector arrays. This optical design allows the instrument to be relatively insensitive to minor line-of-sight pointing drifts and telescope aberrations, and the detectors do not require a sub-Kelvin refrigerator. MISC-T has three science spectral channels that share the same field-of-view by means of beam splitters, and all channels are operated simultaneously to cover the full spectral range from 2.8 to 20 μm at once with exquisite stability and precision (

We report the discovery of a 10 comoving megaparsec (cMpc)-scale structure traced by massive submillimeter galaxies (SMGs) at z ∼ 4.6. These galaxies are selected from an emission line search of ALMA Band 7 observations targeting 184 luminous submillimeter sources (S 850μm ≥ 6.2 mJy) across 1.6 degrees2 in the COSMOS field. We identify four [C ii] emitting SMGs and two probable [C ii] emitting SMG candidates at z = 4.60–4.64 with velocity-integrated signal-to-noise ratio of S/N > 8. Four of the six emitters are near-infrared blank SMGs. After excluding one SMG whose emission line is falling at the edge of the spectral window, all galaxies show clear velocity gradients along the major axes that are consistent with rotating gas disks. The estimated rotation velocities of the disks are 330–550 km s−1 and the inferred host dark-matter halo masses are ∼2–8 × 1012 M ⊙. From their estimated halo masses and [C ii] luminosity function, we suggest that these galaxies have a high (50%–100%) duty cycle and high (∼0.1) baryon conversion efficiency (SFR relative to baryon accretion rate), and that they contribute ≃2% to the total star formation rate density at z = 4.6. These SMGs are concentrated within just 0.3% of the full survey volume, suggesting they are strongly clustered. The extent of this structure and the individual halo masses suggest that these SMGs will likely evolve into members of a ∼1015 M ⊙ cluster at z = 0. This survey reveals a synchronized dusty starburst in massive halos at z > 4, which could be driven by mergers or fed by smooth gas accretion.

Star formation in galaxies is regulated by the heating and cooling in the interstellar medium. In particular, the processing of molecular gas into stars will depend strongly on the ratio of gas heating to gas cooling in the neutral gas around sites of recent star-formation. In this work, we combine mid-infrared (mid-IR) observations of Polycyclic Aromatic Hydrocarbons (PAHs), the dominant heating mechanism of gas in the interstellar medium (ISM), with [C II], [O I], and [Si II] fine-structure emission, the strongest cooling channels in dense, neutral gas. The ratio of IR cooling line emission to PAH emission measures the photoelectric efficiency, a property of the ISM which dictates how much energy carried by ultraviolet photons gets transferred into the gas. We find that star-forming, IR luminous galaxies in the Great Observatories All-Sky LIRG Survey (GOALS) with high IR surface densities have low photoelectric efficiencies. These systems also have, on average, higher ratios of radiation field strength to gas densities, and larger average dust grain size distributions. The data support a scenario in which the most compact galaxies have more young star-forming regions per unit area, which exhibit less efficient gas heating. These conditions may be more common at high-z, and may help explain the higher star-formation rates at cosmic noon. We make predictions on how this can be investigated with JWST., 17 pages, 8 figures, 2 tables, accepted for publication in ApJ

We study the molecular gas content of 24 star-forming galaxies at $z=3-4$, with a median stellar mass of $10^{9.1}$ M$_{\odot}$, from the MUSE Hubble Ultra Deep Field (HUDF) Survey. Selected by their Lyman-alpha-emission and H-band magnitude, the galaxies show an average EW $\approx 20$ angstrom, below the typical selection threshold for Lyman Alpha Emitters (EW $> 25$ angstrom), and a rest-frame UV spectrum similar to Lyman Break Galaxies. We use rest-frame optical spectroscopy from KMOS and MOSFIRE, and the UV features observed with MUSE, to determine the systemic redshifts, which are offset from Lyman alpha by 346 km s$^{-1}$, with a 100 to 600 km s$^{-1}$ range. Stacking CO(4-3) and [CI](1-0) (and higher-$J$ CO lines) from the ALMA Spectroscopic Survey of the HUDF (ASPECS), we determine $3\sigma$ upper limits on the line luminosities of $4.0\times10^{8}$ K km s$^{-1}$pc$^{2}$ and $5.6\times10^{8}$ K km s$^{-1}$pc$^{2}$, respectively (for a 300 km s$^{-1}$ linewidth). Stacking the 1.2 mm and 3 mm dust continuum flux densities, we find a $3\sigma$ upper limits of 9 $\mu$Jy and $1.2$ $\mu$Jy, respectively. The inferred gas fractions, under the assumption of a 'Galactic' CO-to-H$_{2}$ conversion factor and gas-to-dust ratio, are in tension with previously determined scaling relations. This implies a substantially higher $\alpha_{\rm CO} \ge 10$ and $\delta_{\rm GDR} \ge 1200$, consistent with the sub-solar metallicity estimated for these galaxies ($12 + \log(O/H) \approx 7.8 \pm 0.2$). The low metallicity of $z \ge 3$ star-forming galaxies may thus make it very challenging to unveil their cold gas through CO or dust emission, warranting further exploration of alternative tracers, such as [CII]., Comment: 26 pages, 12 figures, 3 tables. Accepted for publication in ApJ

SMI (SPICA Mid-infrared Instrument) is one of the three focal-plane science instruments for SPICA. SMI is the Japanese-led instrument proposed and managed by a university consortium. SMI covers the wavelength range from 10 to 36 μm with four separate channels: the low-resolution (R = 60 – 160) spectroscopy function for 17 – 36 μm, the broad-band (R = 5) imaging function at 34 μm, the mid-resolution (R = 1400 – 2600) spectroscopy function for 18 – 36 μm, and the high-resolution (R = 29000) spectroscopy function for 10 – 18 μm. In this presentation, we will show the latest design and specifications of SMI as a result of feasibility studies.

In this paper, we extend the source detection in the GOODS-ALMA field (69 arcmin$^2$, rms sensitivity $\sigma$ $\simeq$ 0.18 mJy.beam$^{-1}$), to deeper levels than presented in Franco et al. (2018). Using positional information at 3.6 and 4.5 $\mu$m (from Spitzer-IRAC), we explore the presence of galaxies detected at 1.1 mm with ALMA below our original blind detection limit of 4.8-$\sigma$ at which the number of spurious sources starts to dominate over that of real sources. In this Supplementary Catalog, we find a total of 16 galaxies, including 2 galaxies with no counterpart in HST images (also known as optically-dark galaxies) down to a 5$\sigma$ limiting depth of H = 28.2 AB (HST/WFC3 F160W). This brings the total sample of GOODS-ALMA 1.1 mm sources to 35 galaxies. Galaxies in the new sample cover a wider dynamic range in redshift ($z$ = 0.65 - 4.73), are on average twice as large (1.3 vs 0.65 kpc) and and have lower stellar mass (M$_{\star}^{\rm SC}$ = 7.6$\times$10$^{10}$M$_\odot$ vs M$_{\star}^{\rm MC}$ = 1.2$\times$10$^{11}$M$_\odot$). Although exhibiting larger physical sizes, these galaxies have still far-infrared sizes significantly more compact than inferred from their optical emission. We show that the astrometry of the HST image does not only suffer from a global astrometric shift, as already discussed in previous papers, but also from local shifts. These distortions were artificially introduced in the process of building the mosaic of the GOODS-South HST image. By comparing the positions of almost 400 galaxies detected by HST, Pan-STARRS and ALMA, we create a distortion map which can be used to correct for these astrometric issues., Comment: 18 pages, 13 figures. Accepted for publication in A&A

In this paper we study the properties of the six optically dark galaxies detected in the 69 arcmin^2 GOODS-ALMA 1.1mm continuum survey. While none of them are listed in the deepest H-band based CANDELS catalog in the GOODS-South field down to H=28.16AB, we were able to de-blend two of them from their bright neighbor and measure an $H$-band flux for them. We note that AGS4 and AGS15 have H=25.23, 27.11AB respectively. Their extreme proximity (0.50", 0.27") to a bright optical source and their extreme faintness prevented them from being included in the H-band catalog. We present the spectroscopic scan follow-up of five of the six sources with ALMA band 4. All are detected in the 2mm continuum with signal-to-noise ratios higher than eight. One emission line is detected in AGS4 (��_{obs} =151.44GHz with a S/N=8.58) and AGS17 (��_{obs} =154.78GHz with a S/N=10.23), which we interpret in both cases as being due to the CO(6-5) line at z^{AGS4}_{spec}=3.556 and z^{AGS4}_{spec}=3.467, respectively. These redshifts match both the probability distribution of the photometric redshifts derived from the UV to near-infrared spectral energy distributions (SEDs) and the far-infrared SEDs for typical dust temperatures of galaxies at these redshifts. We present evidence that nearly 70% (4/6 of galaxies) of the optically dark galaxies belong to the same overdensity of galaxies at z~3.5. overdensity The most massive one, AGS24 (M_{\star} = 10^{11.32^{+0.02}_{-0.19}} M_{\odot}), is the most massive galaxy without an active galactic nucleus (AGN) at $z$\,>\,3 in the GOODS-ALMA field. It falls in the very center of the peak of the galaxy surface density, which suggests that the surrounding overdensity is a proto-cluster in the process of virialization and that AGS24 is the candidate progenitor of the future brightest cluster galaxy (BCG)., 17 pages, 13 figures, 3 tables, accepted by A&A

We report molecular gas mass estimates obtained from a stacking analysis of CO line emission in the ALMA Spectroscopic Survey (ASPECS) using the spectroscopic redshifts from the optical integral field spectroscopic survey by the Multi Unit Spectroscopic Explorer (MUSE) of the Hubble Ultra Deep Field (HUDF). Stacking was performed on subsets of the sample of galaxies classified by their stellar mass and position relative to the main-sequence relation (on, above, below). Among all the CO emission lines, from CO(2-1) to CO(6-5), with redshifts accessible via the ASPECS Band 3 and the MUSE data, CO(2-1) provides the strongest constraints on the molecular gas content. We detect CO(2-1) emission in galaxies down to stellar masses of $\mathrm{log}({M}_{* }/{M}_{\odot })=10.0$ . Below this stellar mass, we present a new constraint on the molecular gas content of $z\sim 1.5$ main-sequence galaxies by stacking based on the MUSE detections. We find that the molecular gas mass of main-sequence galaxies continuously decreases with stellar mass down to $\mathrm{log}({M}_{* }/{M}_{\odot })\approx 9.0$ . Assuming a metallicity-based CO-to-H2 conversion factor, the molecular gas-to-stellar mass ratio from $\mathrm{log}({M}_{* }/{M}_{\odot })\sim 9.0$ to ˜10.0 does not seem to decrease as fast as for $\mathrm{log}({M}_{* }/{M}_{\odot })\gt 10.0$ , which is in line with simulations and studies at lower redshift. The inferred molecular gas density $\rho ({{\rm{H}}}_{2})=(0.49\pm 0.09)\times {10}^{8}\,{{\rm{M}}}_{\odot }\,{\mathrm{Mpc}}^{-3}$ of MUSE-selected galaxies at $z\sim 1.5$ is comparable with the one derived in the HUDF with a different CO selection. Using the MUSE data we recover most of the CO emission in our deep ALMA observations through stacking, demonstrating the synergy between volumetric surveys obtained at different wave bands.

We make use of sensitive (9.3 microJy/beam RMS) 1.2mm-continuum observations from the ASPECS ALMA large program of the Hubble Ultra Deep Field (HUDF) to probe dust-enshrouded star formation from 1362 Lyman-break galaxies spanning the redshift range z=1.5-10 (to ~7-28 Msolar/yr at 4 sigma over the entire range). We find that the fraction of ALMA-detected galaxies in our z=1.5-10 samples increases steeply with stellar mass, with the detection fraction rising from 0% at 10^9 Msolar to 85(-18)(+9)% at >10^{10} Msolar. Moreover, stacking all 1253 low-mass (10^{9.5} M_{solar} and a SMC-like relation at lower masses. Using stellar-mass and beta measurements for z~2 galaxies over CANDELS, we derive a new empirical relation between beta and stellar mass and then use this correlation to show that our IRX-beta and IRX-stellar mass relations are consistent with each other. We then use these constraints to express the infrared excess as a bivariate function of beta and stellar mass. Finally, we present updated estimates of star-formation rate density determinations at z>3, leveraging current improvements in the measured infrared excess and recent probes of ultra-luminous far-IR galaxies at z>2., 32 pages, 22 figures, 11 pages, accepted for publication in the Astrophysical Journal, key results shown in Figures 7, 8, 11, and 20, updated to match version in press

We present the results from our analysis of NuSTAR data of the luminous infrared galaxy Mrk 266, which contains two nuclei, south-western (SW) and north-eastern (NE), which were resolved in previous Chandra imaging. Combining this with the Chandra data, we intepret the hard X-ray spectrum obtained from a NuSTAR observation to result from a steeply rising flux from a Compton-thick active galactic nuclei (AGN) in the SW nucleus which is very faint in the Chandra band, confirming the previous claim. This hard X-ray component is dominated by reflection, and its intrinsic 2-10 keV luminosity is likely to be ∼1 × 1043 erg s-1. Although it is bright in soft X-ray, only a moderately absorbed NE nucleus has a 2-10 keV luminosity of 4 × 1041 erg s-1, placing it in the low-luminosity AGN class. These results have implications for understanding the detectability and duty cycles of emission from dual AGN in heavily obscured mergers., This research made use of data obtained from NuS-TAR, XMM-Newton and Chandra X-ray Observatory, software packages of HEASoft and R (R Core Team 2017), and the NASA/IPAC Extragalactic Databases (NED), which is funded by the National Aeronatutics and Space Administration and operated by the California Institute of Technology. KI acknowledges support by the Spanish MICINN under grant PID2019-105510GB-C33. T. D-S. acknowledges support from the CASSACA and CONICYT fund CAS-CONICYT Call 2018. C. R. acknowledges support from the Fondecyt Iniciacion grant 11190831. ASE, GCP and KI acknowledge NASA Astrophysics Data Analysis Program (ADAP) Grant 80NSSC20K0450 (PI: U). ASE was also supported by NSF grant AST 1816838.

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., 16 pages, 7 figures, 4 tables, published in Nature

The Lya emitter (LAE) fraction, X_LAE, is a potentially powerful probe of the evolution of the intergalactic neutral hydrogen gas fraction. However, uncertainties in the measurement of X_LAE are still debated. Thanks to deep data obtained with MUSE, we can measure the evolution of X_LAE homogeneously over a wide redshift range of z~3-6 for UV-faint galaxies (down to M_1500~-17.75). This is significantly fainter than in former studies, and allows us to probe the bulk of the population of high-z star-forming galaxies. We construct a UV-complete photo-redshift sample following UV luminosity functions and measure the Lya emission with MUSE using the second data release from the MUSE HUDF Survey. We derive the redshift evolution of X_LAE for M_1500 in [-21.75;-17.75] for the first time with a equivalent width range EW(Lya)>=65 A and find low values of X_ LAE~50 A at z~3-4, in contrast with previous work. The differences in X_LAE mainly arise from selection biases for Lyman Break Galaxies (LBGs) in the literature: UV-faint LBGs are more easily selected if they have strong Lya emission, hence X_LAE is biased towards higher values. Our results suggest either a lower increase of X_LAE towards z~6 than previously suggested, or even a turnover of X_LAE at z~5.5, which may be the signature of a late or patchy reionization process. We compared our results with predictions from a cosmological galaxy evolution model. We find that a model with a bursty star formation (SF) can reproduce our observed X_LAE much better than models where SF is a smooth function of time., Comment: 25 pages, 17 figures, 5 tables, accepted for publication in A&A

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

Galaxy interactions and mergers are thought to play an important role in the evolution of galaxies. Studies in the nearby universe show a higher fraction of active galactic nuclei (AGNs) in interacting and merging galaxies than in their isolated counterparts, indicating that such interactions are important contributors to black hole growth. To investigate the evolution of this role at higher redshifts, we have compiled the largest known sample of major spectroscopic galaxy pairs (2381 with ΔV < 5000 km s-1) at 0.5 < z < 3.0 from observations in the COSMOS and CANDELS surveys. We identify X-ray and IR AGNs among this kinematic pair sample, a visually identified sample of mergers and interactions, and a mass-, redshift-, and environment-matched control sample for each in order to calculate AGN fractions and the level of AGN enhancement as a function of relative velocity, redshift, and X-ray luminosity. While we see a slight increase in AGN fraction with decreasing projected separation, overall, we find no significant enhancement relative to the control sample at any separation. In the closest projected separation bin (< 25 kpc, ΔV < 1000 km s-1), we find enhancements of a factor of ${0.94}_{-0.16}^{+0.21}$ and ${1.00}_{-0.31}^{+0.58}$ for X-ray and IR-selected AGNs, respectively. While we conclude that galaxy interactions do not significantly enhance AGN activity on average over 0.5 < z < 3.0 at these separations, given the errors and the small sample size at the closest projected separations, our results would be consistent with the presence of low-level AGN enhancement....

23 pages, 11 figures, 5 tables. Accepted for publication in the Astrophysical Journal; International audience; We present a CO and atomic fine-structure line luminosity function analysis using the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS). ASPECS consists of two spatially-overlapping mosaics that cover the entire ALMA 3mm and 1.2mm bands. We combine the results of a line candidate search of the 1.2mm data cube with those previously obtained from the 3mm cube. Our analysis shows that $\sim$80% of the line flux observed at 3mm arises from CO(2-1) or CO(3-2) emitters at $z$=1-3 (`cosmic noon'). At 1.2mm, more than half of the line flux arises from intermediate-J CO transitions ($J_{\rm up}$=3-6); $\sim12$% from neutral carbon lines; and $< 1$% from singly-ionized carbon, [CII]. This implies that future [CII] intensity mapping surveys in the epoch of reionization will need to account for a highly significant CO foreground. The CO luminosity functions probed at 1.2mm show a decrease in the number density at a given line luminosity (in units of $L'$) at increasing $J_{\rm up}$ and redshift. Comparisons between the CO luminosity functions for different CO transitions at a fixed redshift reveal sub-thermal conditions on average in galaxies up to $z\sim 4$. In addition, the comparison of the CO luminosity functions for the same transition at different redshifts reveals that the evolution is not driven by excitation. The cosmic density of molecular gas in galaxies, $\rho_{\rm H2}$, shows a redshift evolution with an increase from high redshift up to $z\sim1.5$ followed by a factor $\sim 6$ drop down to the present day. This is in qualitative agreement with the evolution of the cosmic star-formation rate density, suggesting that the molecular gas depletion time is approximately constant with redshift, after averaging over the star-forming galaxy population.

International audience; We compare the extent of the dust, molecular gas and stars in three star-forming galaxies, at $z= 1.4, 1.6$ and $2.7$, selected from the Hubble Ultra Deep Field based on their bright CO and dust-continuum emission as well as their large rest-frame optical sizes. The galaxies have high stellar masses, $\mathrm{M}_*>10^{11}\mathrm{M}_\odot$, and reside on, or slightly below, the main sequence of star-forming galaxies at their respective redshifts. We probe the dust and molecular gas using subarcsecond Atacama Large Millimeter/submillimeter Array observations of the 1.3 mm continuum and CO line emission, respectively, and probe the stellar distribution using \emph{Hubble Space Telescope} observations at 1.6 \textmu m. We find that for all three galaxies the CO emission appears $\gtrsim 30\%$ more compact than the stellar emission. For the $z= 1.4$ and $2.7$ galaxies, the dust emission is also more compact, by $\gtrsim 50\%$, than the stellar emission, whereas for the $z=1.6$ galaxy, the dust and stellar emission have similar spatial extents. This similar spatial extent is consistent with observations of local disk galaxies. However, most high redshift observations show more compact dust emission, likely due to the ubiquity of central starbursts at high redshift and the limited sensitivity of many of these observations. Using the CO emission line, we also investigate the kinematics of the cold interstellar medium in the galaxies, and find that all three have kinematics consistent with a rotation-dominated disk.

We analyze the interstellar medium properties of a sample of sixteen bright CO line emitting galaxies identified in the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS) Large Program. This CO$-$selected galaxy sample is complemented by a couple of additional CO line emitters in the UDF that are identified based on their MUSE optical spectroscopic redshifts. The ASPECS CO$-$selected galaxies cover a larger range of star-formation rates and stellar masses compared to literature CO emitting galaxies at $z>1$ for which scaling relations have been established previously. Most of ASPECS CO-selected galaxies follow these established relations in terms of gas depletion timescales and gas fractions as a function of redshift, as well as the star-formation rate-stellar mass relation (`galaxy main sequence'). However, we find that $\sim30\%$ of the galaxies (5 out of 16) are offset from the galaxy main sequence at their respective redshift, with $\sim12\%$ (2 out of 16) falling below this relationship. Some CO-rich galaxies exhibit low star-formation rates, and yet show substantial molecular gas reservoirs, yielding long gas depletion timescales. Capitalizing on the well-defined cosmic volume probed by our observations, we measure the contribution of galaxies above, below, and on the galaxy main sequence to the total cosmic molecular gas density at different lookback times. We conclude that main sequence galaxies are the largest contributor to the molecular gas density at any redshift probed by our observations (z$\sim$1$-$3). The respective contribution by starburst galaxies above the main sequence decreases from z$\sim$2.5 to z$\sim$1, whereas we find tentative evidence for an increased contribution to the cosmic molecular gas density from the passive galaxies below the main sequence., Accepted for publication in ApJ. Only minor differences respect to previous version

International audience; The Atacama Large Millimeter/submillimeter Array (ALMA) SPECtroscopic Survey in the Hubble Ultra Deep Field (HUDF) is an ALMA large program that obtained a frequency scan in the 3 mm band to detect emission lines from the molecular gas in distant galaxies. Here we present our search strategy for emission lines and continuum sources in the HUDF. We compare several line search algorithms used in the literature, and critically account for the line widths of the emission line candidates when assessing significance. We identify 16 emission lines at high fidelity in our search. Comparing these sources to multiwavelength data we find that all sources have optical/ infrared counterparts. Our search also recovers candidates of lower significance that can be used statistically to derive, e.g., the CO luminosity function. We apply the same detection algorithm to obtain a sample of six 3 mm continuum sources. All of these are also detected in the 1.2 mm continuum with optical/near-infrared counterparts. We use the continuum sources to compute 3 mm number counts in the sub-millijansky regime, and find them to be higher by an order of magnitude than expected for synchrotron-dominated sources. However, the number counts are consistent with those derived at shorter wavelengths (0.85-1.3 mm) once extrapolating to 3 mm with a dust emissivity index of β=1.5, dust temperature of 35 K, and an average redshift of z=2.5. These results represent the best constraints to date on the faint end of the 3 mm number counts.

We discuss the nature and physical properties of gas-mass selected galaxies in the ALMA spectroscopic survey (ASPECS) of the Hubble Ultra Deep Field (HUDF). We capitalize on the deep optical integral-field spectroscopy from the MUSE HUDF Survey and multi-wavelength data to uniquely associate all 16 line-emitters, detected in the ALMA data without preselection, with rotational transitions of carbon monoxide (CO). We identify ten as CO(2-1) at $1 < z < 2$, five as CO(3-2) at $2 < z < 3$ and one as CO(4-3) at $z = 3.6$. Using the MUSE data as a prior, we identify two additional CO(2-1)-emitters, increasing the total sample size to 18. We infer metallicities consistent with (super-)solar for the CO-detected galaxies at $z \le 1.5$, motivating our choice of a Galactic conversion factor between CO luminosity and molecular gas mass for these galaxies. Using deep Chandra imaging of the HUDF, we determine an X-ray AGN fraction of 20% and 60% among the CO-emitters at $z \sim 1.4$ and $z \sim 2.6$, respectively. Being a CO-flux limited survey, ASPECS-LP detects molecular gas in galaxies on, above and below the main sequence (MS) at $z \sim 1.4$. For stellar masses $\ge 10^{10} (10^{10.5})$ M$_{\odot}$, we detect about 40% (50%) of all galaxies in the HUDF at $1 < z < 2$ ($2 < z < 3$). The combination of ALMA and MUSE integral-field spectroscopy thus enables an unprecedented view on MS galaxies during the peak of galaxy formation., 31 pages, 20 figures, 5 tables. Resubmitted to ApJ after addressing the first round of comments by the referee

[Abbreviated] The amount of integral field spectrograph (IFS) data has grown considerable over the last few decades. The demand for tools to analyze such data is therefore bigger now than ever. We present TDOSE; a flexible Python tool for Three Dimensional Optimal Spectral Extraction from IFS data cubes. TDOSE works on any three-dimensional data cube and bases the spectral extractions on morphological reference image models. In each wavelength layer of the IFS data cube, TDOSE simultaneously optimizes all sources in the morphological model to minimize the difference between the scaled model components and the IFS data. The flux optimization produces individual data cubes containing the scaled three-dimensional source models. This allows for efficient de-blending of flux in both the spatial and spectral dimensions of the IFS data cubes, and extraction of the corresponding one-dimensional spectra. We present an example of how the three-dimensional source models generated by TDOSE can be used to improve two-dimensional maps of physical parameters. By extracting TDOSE spectra of $\sim$150 [OII] emitters from the MUSE-Wide survey we show that the median increase in line flux is $\sim$5% when using multi-component models as opposed to single-component models. However, the increase in recovered line emission in individual cases can be as much as 50%. Comparing the TDOSE model-based extractions of the MUSE-Wide [OII] emitters with aperture spectra, the TDOSE spectra provides a median flux (S/N) increase of 9% (14%). Hence, TDOSE spectra optimizes the S/N while still being able to recover the total emitted flux. TDOSE version 3.0 presented in this paper is available at https://github.com/kasperschmidt/TDOSE., Comment: Accepted by A&A, 24 pages and 11 figures, code available at https://github.com/kasperschmidt/TDOSE

By using the deepest available mid and far infrared surveys in the CANDELS, GOODS and COSMOS fields we study the evolution of the Main Sequence (MS) of star forming galaxies (SFGs) from z~0 to` ~2.5 at stellar masses larger than 10^{10} M_{\odot}. The MS slope and scatter are consistent with a re-scaled version of the local relation and distribution, shifted at higher values of SFR according to ~(1+z)^{3.2}. The relation exhibits a bending at the high mass end and a slightly increasing scatter as a function of the stellar mass. We show that the previously reported evolution of the MS slope, in the considered mass and redshift range, is due to a selection effect. The distribution of galaxies in the MS region at fixed stellar mass is well represented by a single log-normal distribution at all redshifts and masses, with starburst galaxies (SBs) occupying the tail at high SFR., Comment: Accepted for publication by MNRAS

It is still a challenge to assess the merger fraction of galaxies at different cosmic epochs in order to probe the evolution of their mass assembly. Using the Illustris cosmological simulations, we investigate the relation between the separation of galaxies in a pair, both in velocity and projected spatial separation space, and the probability that these interacting galaxies will merge in the future. From this analysis, we propose a new set of criteria to select close pairs of galaxies along with a new corrective term to be applied to the computation of the galaxy merger fraction. We then probe the evolution of the major and minor merger fraction using the latest MUSE deep observations over the HUDF, HDFS, COSMOS-Gr30 and Abell 2744 regions. From a parent sample of 2483 galaxies with spectroscopic redshifts, we identify 366 close pairs spread over a large range of redshifts ($0.2, Comment: 26 pages, 15 figures. Accepted for publication in A&A

We present a detailed census of galaxies in and around PC217.96+32.3, a spectroscopically confirmed Coma analog at z=3.78. Diverse galaxy types identified in the field include Lya emitters (LAEs), massive star-forming galaxies, and ultra-massive galaxies (log (Mstar/Msun)>= 11) which may have already halted their star formation. The sky distribution of the star-forming galaxies suggests the presence of a significant overdensity (delta_g=8+/-2), which is spatially offset from the previously confirmed members by 3-4 Mpc to the west. Candidate quiescent and post-starburst galaxies are also found in large excess (a factor of ~ 8-15 higher surface density than the field) although their redshifts are less certain. We estimate that the total enclosed mass traced by star-forming galaxy candidates is roughly comparable to that of PC217.96+32.3 traced by the LAEs. We speculate that the true extent of PC217.96+32.3 may be larger than previously known, a half of which is missed by our LAE selection. Alternatively, the newly discovered overdensity may belong to another Coma progenitor not associated with PC217.96+32.3. Expectations from theory suggest that both scenarios are equally unlikely (, Comment: 21 pages, 11 figures, accepted for publication in Astrophysical Journal

By using a set of different SFR indicators, including WISE mid-infrared and Halpha emission, we study the slope of the Main Sequence (MS) of local star forming galaxies at stellar masses larger than 10^{10} M_{\odot}. The slope of the relation strongly depends on the SFR indicator used. In all cases, the local MS shows a bending at high stellar masses with respect to the slope obtained in the low mass regime. While the distribution of galaxies in the upper envelope of the MS is consistent with a log-normal distribution, the lower envelope shows an excess of galaxies, which increases as a function of the stellar mass but varies as a function of the SFR indicator used. The scatter of the best log-normal distribution increases with stellar mass from ~0.3 dex at 10^{10} M_{\odot} to ~0.45 at 10^{11} M_{\odot}. The MS high-mass end is dominated by central galaxies of group sized halos with a red bulge and a disk redder than the lower mass counterparts. We argue that the MS bending in this region is due to two processes: i) the formation of a bulge component as a consequence of the increased merger activity in groups, and ii) the cold gas starvation induced by the hot halo environment, which cuts off the gas inflow onto the disk. Similarly, the increase of the MS scatter at high stellar masses would be explained by the larger spread of star formation histories of central group and cluster galaxies with respect to lower mass systems., 16 pages with 14 figures, accepted for publication on MNRAS