Showing 1–50 of 524 results for author: Van der Werf, P

MICONIC: JWST/MIRI MRS observations of the nuclear and circumnuclear regions of Mrk231

Authors: A. Alonso-Herrero, L. Hermosa Muñoz, A. Labiano, P. Guillard, V. A. Buiten, D. Dicken, P. van der Werf, J. Álvarez-Márquez, T. Böker, L. Colina, A. Eckart, M. García-Marín, O. C. Jones, L. Pantoni, P. G. Pérez-González, D. Rouan, M. J. Ward, M. Baes, G. Östlin, P. Royer, G. S. Wright, M. Güdel, Th. Henning, P. -O. Lagage, E. F. van Dishoeck

Abstract: We present JWST/MIRI MRS spatially resolved $\sim 5-28\,μ$m observations of the central ~4-8kpc of the ultraluminous infrared galaxy and broad absorption line quasar Mrk231. These are part of the Mid-Infrared Characterization of Nearby Iconic galaxy Centers (MICONIC) program of the MIRI European Consortium guaranteed time observations. No high excitation lines (i.e., [MgV] at 5.61$μ$m or [NeV] at… ▽ More We present JWST/MIRI MRS spatially resolved $\sim 5-28\,μ$m observations of the central ~4-8kpc of the ultraluminous infrared galaxy and broad absorption line quasar Mrk231. These are part of the Mid-Infrared Characterization of Nearby Iconic galaxy Centers (MICONIC) program of the MIRI European Consortium guaranteed time observations. No high excitation lines (i.e., [MgV] at 5.61$μ$m or [NeV] at 14.32$μ$m) typically associated with the presence of an active galactic nucleus (AGN) are detected in the nuclear region of Mrk231. This is likely due to the intrinsically X-ray weak nature of its quasar. Some intermediate ionization potential lines, for instance, [ArIII] at 8.99$μ$m and [SIV] at 10.51$μ$m, are not detected either, even though they are clearly observed in a star-forming region ~920pc south-east of the AGN. Thus, the strong nuclear mid-infrared (mid-IR) continuum is also in part hampering the detection of faint lines in the nuclear region. The nuclear [NeIII]/[NeII]line ratio is consistent with values observed in star-forming galaxies. Moreover, we resolve for the first time the nuclear starburst in the mid-IR low-excitation line emission (size of ~400pc, FWHM). Several pieces of evidence also indicate that it is partly obscured even at these wavelengths. At the AGN position, the ionized and warm molecular gas emission lines have modest widths (W_80~300km/s). There are, however, weak blueshifted wings reaching velocities v_02~-400km/s in [NeII]. The nuclear starburst is at the center of a large (~8kpc), massive rotating disk with widely-spread, low velocity outflows. Given the high star formation rate of Mrk231, we speculate that part of the nuclear outflows and the large-scale non-circular motions observed in the mid-IR are driven by its powerful nuclear starburst. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: Accepted for publication in Astronomy & Astrophysics

MIDIS. Near-infrared rest-frame morphology of massive galaxies at $3<z<5.5$ in the Hubble eXtreme Deep Field

Authors: L. Costantin, S. Gillman, L. A. Boogaard, P. G. Pérez-González, E. Iani, P. Rinaldi, J. Melinder, A. Crespo Gómez, L. Colina, T. R. Greve, G. Östlin, G. Wright, A. Alonso-Herrero, J. Álvarez-Márquez, M. Annunziatella, A. Bik., K. I. Caputi, D. Dicken, A. Eckart, J. Hjorth, O. Ilbert, I. Jermann, A. Labiano, D. Langeroodi, F. Peißker , et al. (7 additional authors not shown)

Abstract: Thanks to decades of observations using HST, the structure of galaxies at redshift $z>2$ has been widely studied in the rest-frame ultraviolet regime, which traces recent star formation from young stellar populations. But, we still have little information about the spatial distribution of the older, more evolved, stellar populations, constrained by the rest-frame infrared portion of galaxies' spec… ▽ More Thanks to decades of observations using HST, the structure of galaxies at redshift $z>2$ has been widely studied in the rest-frame ultraviolet regime, which traces recent star formation from young stellar populations. But, we still have little information about the spatial distribution of the older, more evolved, stellar populations, constrained by the rest-frame infrared portion of galaxies' spectral energy distribution. We present the morphological characterization of a sample of 21 massive galaxies ($\log(M_{\star}/M_{\odot})>9.5$) at redshift $3<z<5.5$. These galaxies are observed as part of the GTO program MIDIS with the Mid-Infrared Instrument (MIRI) onboard JWST. The deep MIRI 5.6~$μ$m imaging allows us to characterize for the first time the rest-frame near-infrared structure of galaxies beyond cosmic noon, at higher redshifts than possible with NIRCam, tracing their older stellar populations. We derive the galaxies' non-parametric morphology and model the galaxies' light distribution with a Sérsic component. We find that at $z>3$ massive galaxies show a smooth distribution of their rest-infrared light, strongly supporting the increasing number of regular disk galaxies already in place at early epochs. On the contrary, the ultraviolet structure obtained from HST observations is generally more irregular, catching the most recent episodes of star formation. Importantly, we find a segregation of morphologies across cosmic time, having massive galaxies at redshift $z>4$ later-type morphologies compared to $z\sim3$ galaxies. These findings suggest a transition phase in galaxy assembly and central mass build up already taking place at $z\sim3-4$. MIRI provides unique information about the structure of the mature stellar population of high-redshift galaxies, unveiling that massive galaxies beyond cosmic noon are prevalently compact disk galaxies with smooth mass distribution. △ Less

Submitted 28 June, 2024; originally announced July 2024.

Comments: Submitted to A&A. 7 pages, 4 figures. Comments welcome

MIDIS: MIRI uncovers Virgil, an extended source at $z\simeq 6.6$ with the photometric properties of Little Red Dots

Authors: Edoardo Iani, Pierluigi Rinaldi, Karina I. Caputi, Marianna Annunziatella, Danial Langeroodi, Jens Melinder, Pablo G. Pérez-González, Javier Álvarez-M��rquez, Leindert A. Boogaard, Sarah E. I. Bosman, Luca Costantin, Thibaud Moutard, Luis Colina, Göran Östlin, Thomas R. Greve, Gillian Wright, Almudena Alonso-Herrero, Arjan Bik, Steven Gillman, Alejandro Crespo Gómez, Jens Hjorth, Alvaro Labiano, John P. Pye, Tuomo V. Tikkanen, Paul P. van der Werf

Abstract: We present Virgil, a MIRI extremely red object (MERO) detected with the F1000W filter as part of the MIRI Deep Imaging Survey (MIDIS) observations of the Hubble Ultra Deep Field (HUDF). Virgil is a Lyman-$α$ emitter (LAE) at $z_{spec} = 6.6312\pm 0.0019$ (from VLT/MUSE) with a rest-frame UV-to-optical spectral energy distribution (SED) typical of LAEs at similar redshifts. However, MIRI observatio… ▽ More We present Virgil, a MIRI extremely red object (MERO) detected with the F1000W filter as part of the MIRI Deep Imaging Survey (MIDIS) observations of the Hubble Ultra Deep Field (HUDF). Virgil is a Lyman-$α$ emitter (LAE) at $z_{spec} = 6.6312\pm 0.0019$ (from VLT/MUSE) with a rest-frame UV-to-optical spectral energy distribution (SED) typical of LAEs at similar redshifts. However, MIRI observations reveal an unexpected extremely red color at rest-frame near-infrared wavelengths, $\rm F444W - F1000W = 2.33 \pm 0.06$. Such steep rise in the near-infrared, completely missed without MIRI imaging, is poorly reproduced by models including only stellar populations and hints towards the presence of an Active Galactic Nucleus (AGN). Interestingly, the overall SED shape of Virgil resembles that of the recently discovered population of Little Red Dots (LRDs) but does not meet their compactness criterion: at rest-frame UV-optical wavelengths Virgil's morphology follows a 2D-Sérsic profile with average index $n = 0.93^{+0.85}_{-0.31}$ and $r_e = 0.43$~pkpc. Only at MIRI wavelengths Virgil is unresolved due to the coarser PSF. We also estimate a bolometric luminosity $L_{\rm bol} = (8.4-11.1)\times 10^{44}\rm~erg~s^{-1}$ and a supermassive black hole mass $M_{\rm BH} = (4-7)\times 10^7\rm ~ M_\odot$ in agreement with recently reported values for LRDs. This discovery demonstrates the crucial importance of deep MIRI surveys to find AGN amongst high-$z$ galaxies that otherwise would be completely missed and raises the question of how common Virgil-like objects could be in the early Universe. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: 17 pages, 10 figures, 3 tables. Submitted to ApJ

Imaging of I Zw 18 by JWST: II. Spatially resolved star formation history

Authors: Giacomo Bortolini, Göran Östlin, Nolan Habel, Alec S. Hirschauer, Olivia C. Jones, Kay Justtanont, Margaret Meixner, Martha L. Boyer, Joris A. D. L. Blommaert, Nicolas Crouzet, Lenkić, Conor Nally, Beth A. Sargent, Paul van der Werf, Manuel Güdel, Thomas Henning, Pierre O. Lagage

Abstract: The blue compact dwarf galaxy I Zw 18 is one of the most metal-poor ($Z \sim 3% Z_{\sun}$) star-forming galaxies in the local Universe. Its evolutionary status has sparked debate within the astronomical community. We aim to investigate the stellar populations of I Zw 18 in the near-IR using JWST/NIRCam's high spatial resolution and sensitivity. Additionally, we aim to derive the galaxy's spatially… ▽ More The blue compact dwarf galaxy I Zw 18 is one of the most metal-poor ($Z \sim 3% Z_{\sun}$) star-forming galaxies in the local Universe. Its evolutionary status has sparked debate within the astronomical community. We aim to investigate the stellar populations of I Zw 18 in the near-IR using JWST/NIRCam's high spatial resolution and sensitivity. Additionally, we aim to derive the galaxy's spatially resolved star formation history (SFH) over the last 1 Gyr and provide constraints for older epochs. We used DOLPHOT to measure positions and fluxes of point sources in the F115W and F200W filters' images of I Zw 18. To derive I Zw 18's SFH, we applied the color-magnitude diagram (CMD) fitting technique SFERA 2.0, using two independent sets of stellar models. Our analysis reveals three main stellar populations: one younger than $\sim30$ Myr, mainly in the northwest star-forming (SF) region; an intermediate-age population ($\sim 100 - 800$ Myr) in the southeast SF region; and a red and faint population linked to the underlying halo, older than 1 Gyr and possibly as old as 13.8 Gyr. The main body of the galaxy shows a very low star formation rate (SFR) of $\sim 10^{-4} M_{\odot} \text{yr}^{-1}$ between 1 and 13.8 Gyr ago. In the last billion years, I Zw 18 shows increasing SF, with strong bursts around $\sim10$ and $\sim100$ Myr ago. Component C mirrors the main body's evolution but with lower SFRs. Our findings confirm that I Zw 18 contains stars of all ages, indicating it is not a young galaxy but has an old stellar halo, similar to other BCDs. The low SF activity over the past billion years supports the "slow cooking" dwarf scenario, explaining its low metal content. Currently, the galaxy is undergoing its strongest SF episode ($\sim 0.6 M_{\odot} \text{yr}^{-1}$) mainly in the northwest region, likely due to a recent gravitational interaction with Component C. △ Less

Submitted 25 June, 2024; originally announced June 2024.

Comments: 14 pages, 11 figures, accepted for publications in Astronomy & Astrophysics (section "4. Extragalactic astronomy")

The Structure of Massive Star-Forming Galaxies from JWST and ALMA: Dusty, High Redshift Disk Galaxies

Authors: Steven Gillman, Ian Smail, Bitten Gullberg, A. M. Swinbank, Aswin P. Vijayan, Minju Lee, Gabe Brammer, U. Dudzevičiūtė, Thomas R. Greve, Omar Almaini, Malte Brinch, Scott C. Chapman, Chian-Chou Chen, Soh Ikarashi, Yuichi Matsuda, Wei-Hao Wang, Fabian Walter, Paul P. van der Werf

Abstract: We present an analysis of the JWST NIRCam and MIRI morphological properties of 80 massive ($\log_{10}(M_\ast[M_{\odot}])$=11.2$\pm$0.1) dusty star-forming galaxies at $z$$=$2.7$^{+1.2}_{-0.7}$, identified as sub-millimetre galaxies (SMGs) by ALMA, that have been observed as part of the JWST PRIMER project. To compare the structure of these massive, active galaxies to more typical star-forming gala… ▽ More We present an analysis of the JWST NIRCam and MIRI morphological properties of 80 massive ($\log_{10}(M_\ast[M_{\odot}])$=11.2$\pm$0.1) dusty star-forming galaxies at $z$$=$2.7$^{+1.2}_{-0.7}$, identified as sub-millimetre galaxies (SMGs) by ALMA, that have been observed as part of the JWST PRIMER project. To compare the structure of these massive, active galaxies to more typical star-forming galaxies, we define a sample of 850 field galaxies with matched redshifts and specific star formation rates. We identify 20$\pm$5% of the SMGs as candidate late-stage major mergers, a further 40$\pm$10% as potential minor mergers and 40$\pm$10% which have comparatively undisturbed disk-like morphologies, with no obvious massive neighbours on $\lesssim$20-30kpc (projected) scales. These rates are comparable to those for the field sample and indicate that the majority of the sub-millimetre-detected galaxies are not late-stage major mergers, but have interaction rates similar to the less-active population at $z$$\sim$2-3. Through a multi-wavelength morphological analysis, we establish that SMGs have comparable near-infrared sizes to the less active population, but exhibit lower Sérsic indices, consistent with bulge-less disks and have more structured morphologies at 2$μ$m relative to 4$μ$m. We find evidence for dust reddening as the origin of the morphological differences between the populations, identifying a strong correlation between the F200W$-$F444W pixel colour and the 870$μ$m surface brightness. We conclude that SMGs and less active galaxies at the same epochs share a common disk-like structure, but the weaker bulge components of the SMGs results in a lower dynamical stability. Consequently, instabilities triggered either secularly or by minor external perturbations result in higher levels of activity (and dust content) in SMGs compared to typical star-forming galaxies. [Abridged] △ Less

Submitted 5 June, 2024; originally announced June 2024.

Comments: 29 pages, 13 figures including appendices. Submitted to A&A. Comments welcome

An ALCHEMI inspection of sulphur-bearing species towards the central molecular zone of NGC 253

Authors: M. Bouvier, S. Viti, E. Behrens, J. Butterworth, K. -Y. Huang, J. G. Mangum, N. Harada, S. Martín, V. M. Rivilla, S. Muller, K. Sakamoto, Y. Yoshimura, K. Tanaka, K. Nakanishi, R. Herrero-Illana, L. Colzi, M. D. Gorski, C. Henkel, P. K. Humire, D. S. Meier, P. P. van der Werf, Y. T. Yan

Abstract: Sulphur-bearing species are detected in various environments within Galactic star-forming regions and are particularly abundant in the gas phase of outflows and shocks, and photo-dissociation regions. In this work, we aim to investigate the nature of the emission from the most common sulphur-bearing species observable at millimetre wavelengths towards the nuclear starburst of the galaxy NGC 253. W… ▽ More Sulphur-bearing species are detected in various environments within Galactic star-forming regions and are particularly abundant in the gas phase of outflows and shocks, and photo-dissociation regions. In this work, we aim to investigate the nature of the emission from the most common sulphur-bearing species observable at millimetre wavelengths towards the nuclear starburst of the galaxy NGC 253. We intend to understand which type of regions are probed by sulphur-bearing species and which process(es) dominate(s) the release of sulphur into the gas phase. We used the high-angular resolution (1.6" or 27 pc) observations from the ALCHEMI ALMA Large Program to image several sulphur-bearing species towards the central molecular zone (CMZ) of NGC 253. We performed local thermodynamic equilibrium (LTE) and non-LTE large velocity gradient (LVG) analyses to derive the physical conditions of the gas in which S-bearing species are emitted, and their abundance ratios across the CMZ. Finally, we compared our results with previous ALCHEMI studies and a few selected Galactic environments. We found that not all sulphur-bearing species trace the same type of gas: strong evidence indicates that H2S and part of the emission of OCS, H2CS, and SO, are tracing shocks whilst part of SO and CS emission rather trace the dense molecular gas. For some species, such as CCS and SO2, we could not firmly conclude on their origin of emission. The present analysis indicates that the emission from most sulphur-bearing species throughout the CMZ is likely dominated by shocks associated with ongoing star formation. In the inner part of the CMZ where the presence of super star clusters was previously indicated, we could not distinguish between shocks or thermal evaporation as the main process releasing the S-bearing species. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: 44 pages, 20 figures, Accepted for publication in A&A

REBELS-25: Discovery of a dynamically cold disc galaxy at z = 7.31

Authors: Lucie E. Rowland, Jacqueline Hodge, Rychard Bouwens, Pavel Mancera Piña, Alexander Hygate, Hiddo Algera, Manuel Aravena, Rebecca Bowler, Elisabete da Cunha, Pratika Dayal, Andrea Ferrara, Thomas Herard-Demanche, Hanae Inami, Ivana van Leeuwen, Ilse de Looze, Pascal Oesch, Andrea Pallottini, Siân Phillips, Matus Rybak, Sander Schouws, Renske Smit, Laura Sommovigo, Mauro Stefanon, Paul van der Werf

Abstract: We present high resolution ($\sim0.14$" = 710 pc) ALMA [CII] 158$μ$m and dust continuum follow-up observations of REBELS-25, a [CII]-luminous ($L_{\mathrm{[CII]}}=(1.7\pm0.2)\times 10^9 \mathrm{L_{\odot}}$) galaxy at redshift $z=7.3065\pm0.0001$. These high resolution, high signal-to-noise observations allow us to study the sub-kpc morphology and kinematics of this massive (… ▽ More We present high resolution ($\sim0.14$" = 710 pc) ALMA [CII] 158$μ$m and dust continuum follow-up observations of REBELS-25, a [CII]-luminous ($L_{\mathrm{[CII]}}=(1.7\pm0.2)\times 10^9 \mathrm{L_{\odot}}$) galaxy at redshift $z=7.3065\pm0.0001$. These high resolution, high signal-to-noise observations allow us to study the sub-kpc morphology and kinematics of this massive ($M_* = 8^{+4}_{-2} \times 10^9 \mathrm{M_{\odot}}$) star-forming (SFR$_{\mathrm{UV+IR}} = 199^{+101}_{-63} \mathrm{M_{\odot}} \mathrm{yr}^{-1}$) galaxy in the Epoch of Reionisation. By modelling the kinematics with $^{\mathrm{3D}}$BAROLO, we find it has a low velocity dispersion ($\barσ = 33 \pm 9$ km s$^{-1}$) and a high ratio of ordered-to-random motion ($V_{\mathrm{rot, ~max}}/\barσ = 11 ^{+8}_{-4}$), indicating that REBELS-25 is a dynamically cold disc. Additionally, we find that the [CII] distribution is well fit by a near-exponential disc model, with a Sérsic index, $n$, of $1.3 \pm 0.2$, and we see tentative evidence of more complex non-axisymmetric structures suggestive of a bar in the [CII] and dust continuum emission. By comparing to other high spatial resolution cold gas kinematic studies, we find that dynamically cold discs seem to be more common in the high redshift Universe than expected based on prevailing galaxy formation theories, which typically predict more turbulent and dispersion-dominated galaxies in the early Universe as an outcome of merger activity, gas accretion and more intense feedback. This higher degree of rotational support seems instead to be consistent with recent cosmological simulations that have highlighted the contrast between cold and warm ionised gas tracers, particularly for massive galaxies. We therefore show that dynamically settled disc galaxies can form as early as 700 Myr after the Big Bang. △ Less

Submitted 9 May, 2024; originally announced May 2024.

Comments: Submitted to MNRAS

JWST Observations of Starbursts: Cold Clouds and Plumes Launching in the M82 Outflow

Authors: Deanne B. Fisher, Alberto D. Bolatto, John Chisholm, Drummond Fielding, Rebecca C. Levy, Elizabeth Tarantino, Martha L. Boyer, Serena A. Cronin, Laura A. Lopez, J. D. Smith, Danielle A. Berg, Sebastian Lopez, Sylvain Veilleux, Paul P. van der Werf, Torsten Böker, Leindert A. Boogaard, Laura Lenkić, Simon C. O. Glover, Vicente Villanueva, Divakara Mayya, Thomas S. -Y. Lai, Daniel A. Dale, Kimberly L. Emig, Fabian Walter, Monica Relaño , et al. (6 additional authors not shown)

Abstract: In this paper we study the filamentary substructure of 3.3 $μ$m PAH emission from JWST/NIRCam observations in the base of the M82 star-burst driven wind. We identify plume-like substructure within the PAH emission with widths of $\sim$50 pc. Several of the plumes extend to the edge of the field-of-view, and thus are at least 200-300 pc in length. In this region of the outflow, the vast majority (… ▽ More In this paper we study the filamentary substructure of 3.3 $μ$m PAH emission from JWST/NIRCam observations in the base of the M82 star-burst driven wind. We identify plume-like substructure within the PAH emission with widths of $\sim$50 pc. Several of the plumes extend to the edge of the field-of-view, and thus are at least 200-300 pc in length. In this region of the outflow, the vast majority ($\sim$70\%) of PAH emission is associated with the plumes. We show that those structures contain smaller scale "clouds" with widths that are $\sim$5-15 pc, and they are morphologically similar to the results of "cloud-crushing" simulations. We estimate the cloud-crushing time-scales of $\sim$0.5-3 Myr, depending on assumptions. We show this time scale is consistent with a picture in which these observed PAH clouds survived break-out from the disk rather than being destroyed by the hot wind. The PAH emission in both the midplane and the outflow is shown to tightly correlate with that of Pa$α$ emission (from HST/NICMOS data), at the scale of both plumes and clouds, though the ratio of PAH-to-Pa$α$ increases at further distances from the midplane. Finally, we show that the outflow PAH emission is suppressed in regions of the M82 wind that are bright in X-ray emission. Overall, our results are broadly consistent with a picture in which cold gas in galactic outflows is launched via hierarchically structured plumes, and those small scale clouds are more likely to survive the wind environment when collected into the larger plume structure. △ Less

Submitted 6 May, 2024; originally announced May 2024.

Comments: Submitted to MNRAS

A dusty proto-cluster surrounding the binary galaxy HerBS-70 at $z = 2.3$

Authors: Tom J. L. C. Bakx, S. Berta, H. Dannerbauer, P. Cox, K. M. Butler, M. Hagimoto, D. H. Hughes, D. A. Riechers, P. P. van der Werf, C. Yang, A. J. Baker, A. Beelen, G. J. Bendo, E. Borsato, V. Buat, A. R. Cooray, L. Dunne, S. Dye, S. Eales, R. Gavazzi, A. I. Harris, D. Ismail, R. J. Ivison, B. Jones, M. Krips , et al. (16 additional authors not shown)

Abstract: We report on deep SCUBA-2 observations at 850$μ$m and NOEMA spectroscopic measurements at 2 mm of the environment surrounding the luminous, massive ($M_{*} \approx 2 \times 10^{11}$ M$_{\odot}$) Herschel-selected source HerBS-70. This source was revealed by previous NOEMA observations to be a binary system of dusty star-forming galaxies at $z= 2.3$, with the East component (HerBS-70E) hosting an A… ▽ More We report on deep SCUBA-2 observations at 850$μ$m and NOEMA spectroscopic measurements at 2 mm of the environment surrounding the luminous, massive ($M_{*} \approx 2 \times 10^{11}$ M$_{\odot}$) Herschel-selected source HerBS-70. This source was revealed by previous NOEMA observations to be a binary system of dusty star-forming galaxies at $z= 2.3$, with the East component (HerBS-70E) hosting an Active Galactic Nucleus (AGN). The SCUBA-2 observations detected, in addition to the binary system, twenty-one sources at $> 3.5 σ$ over an area of $\sim 25$ square comoving Mpc with a sensitivity of $σ_{850} = 0.75$ mJy. The surface density of continuum sources around HerBS-70 is three times higher than for field galaxies. The NOEMA spectroscopic measurements confirm the protocluster membership of three of the nine brightest sources through their CO(4 - 3) line emission, yielding a volume density 36 times higher than for field galaxies. All five confirmed sub-mm galaxies in the HerBS-70 system have relatively short gas depletion times ($80 - 500$ Myr), indicating the onset of quenching for this protocluster core due to the depletion of gas. The dark matter halo mass of the HerBS-70 system is estimated around $5 \times{} 10^{13}$ M$_{\odot}$, with a projected current-day mass of $10^{15}$ M$_{\odot}$, similar to the local Virgo and Coma clusters. These observations support the claim that DSFGs, in particular the ones with observed multiplicity, can trace cosmic overdensities. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Comments: 19 pages, 13 figures, accepted for publication in MNRAS

Physical Properties of the Southwest Outflow Streamer in the Starburst Galaxy NGC 253 with ALCHEMI

Authors: Min Bao, Nanase Harada, Kotaro Kohno, Yuki Yoshimura, Fumi Egusa, Yuri Nishimura, Kunihiko Tanaka, Kouichiro Nakanishi, Sergio Martín, Jeffrey G. Mangum, Kazushi Sakamoto, Sébastien Muller, Mathilde Bouvier, Laura Colzi, Kimberly L. Emig, David S. Meier, Christian Henkel, Pedro Humire, Ko-Yun Huang, Víctor M. Rivilla, Paul van der Werf, Serena Viti

Abstract: The physical properties of galactic molecular outflows are important as they could constrain outflow formation mechanisms. We study the properties of the southwest (SW) outflow streamer including gas kinematics, optical depth, dense gas fraction, and shock strength in the central molecular zone of the starburst galaxy NGC 253. We image the molecular emission at a spatial resolution of $\sim$27 pc… ▽ More The physical properties of galactic molecular outflows are important as they could constrain outflow formation mechanisms. We study the properties of the southwest (SW) outflow streamer including gas kinematics, optical depth, dense gas fraction, and shock strength in the central molecular zone of the starburst galaxy NGC 253. We image the molecular emission at a spatial resolution of $\sim$27 pc based on data from the ALCHEMI program. We trace the kinematics of molecular gas with CO(1-0) line. We constrain the optical depth of CO emission with CO/$^{13}$CO(1-0) ratio, the dense gas fraction with HCN/CO(1-0) ratio, as well as the shock strength with SiO(2-1)/$^{13}$CO(1-0) ratio. The CO/$^{13}$CO(1-0) integrated intensity ratio is $\sim$21 in the SW streamer region, which approximates the C/$^{13}$C isotopic abundance ratio. The higher integrated intensity ratio compared to the disk can be attributed to the optically thinner environment for CO(1-0) emission inside the SW streamer. The HCN/CO(1-0) and SiO(2-1)/$^{13}$CO(1-0) integrated intensity ratios both approach $\sim$0.2 in three giant molecular clouds (GMCs) at the base of the outflow streamers, which implies the higher dense gas fraction and enhanced strength of fast shocks in those GMCs than in the disk. The contours of those two integrated intensity ratios are extended towards the directions of outflow streamers, which connects the enhanced dense gas fraction and shock strength with molecular outflow. Moreover, the molecular gas with enhanced dense gas fraction and shock strength located at the base of the SW streamer shares the same velocity with the outflow. These phenomena suggest that the star formation inside the GMCs can trigger the shocks and further drive the molecular outflow. △ Less

Submitted 6 April, 2024; originally announced April 2024.

Comments: Accepted for publication in A&A, 14 pages, 11 figures

A spectacular galactic scale magnetohydrodynamic powered wind in ESO 320-G030

Authors: M. D. Gorski, S. Aalto, S. König, C. F. Wethers, C. Yang, S. Muller, K. Onishi, M. Sato, N. Falstad, Jeffrey G. Mangum, S. T. Linden, F. Combes, S. Martín, M. Imanishi, Keiichi Wada, L. Barcos-Muñoz, F. Stanley, S. García-Burillo, P. P. van der Werf, A. S. Evans, C. Henkel, S. Viti, N. Harada, T. Díaz-Santos, J. S. Gallagher , et al. (1 additional authors not shown)

Abstract: How galaxies regulate nuclear growth through gas accretion by supermassive black holes (SMBHs) is one of the most fundamental questions in galaxy evolution. One potential way to regulate nuclear growth is through a galactic wind that removes gas from the nucleus. It is unclear whether galactic winds are powered by jets, mechanical winds, radiation, or via magnetohydrodynamic (MHD) processes. Compa… ▽ More How galaxies regulate nuclear growth through gas accretion by supermassive black holes (SMBHs) is one of the most fundamental questions in galaxy evolution. One potential way to regulate nuclear growth is through a galactic wind that removes gas from the nucleus. It is unclear whether galactic winds are powered by jets, mechanical winds, radiation, or via magnetohydrodynamic (MHD) processes. Compact obscured nuclei (CONs) represent a significant phase of galactic nuclear growth. These galaxies hide growing SMBHs or unusual starbursts in their very opaque, extremely compact (r $<$ 100 pc) centres. They are found in approximately 30 % of the luminous and ultra-luminous infrared galaxy (LIRG and ULIRG) population. Here, we present high-resolution ALMA observations ($\sim$30 mas, $\sim$5 pc) of ground-state and vibrationally excited HCN towards ESO 320-G030 (IRAS 11506-3851). ESO 320-G030 is an isolated luminous infrared galaxy known to host a compact obscured nucleus and a kiloparsec-scale molecular wind. Our analysis of these high-resolution observations excludes the possibility of a starburst-driven wind, a mechanically or energy driven active galactic nucleus (AGN) wind, and exposes a molecular MDH wind. These results imply that the nuclear evolution of galaxies and the growth of SMBHs are similar to the growth of hot cores or protostars where gravitational collapse of the nuclear torus drives a MHD wind. These results mean galaxies are capable, in part, of regulating the evolution of their nuclei without feedback. △ Less

Submitted 25 March, 2024; originally announced March 2024.

Comments: 9 pages, 10 figures

Journal ref: A&A 684, L11 (2024)

GOALS-JWST: The Warm Molecular Outflows of the Merging Starburst Galaxy NGC 3256

Authors: Thomas Bohn, Hanae Inami, Aditya Togi, Lee Armus, Thomas S. -Y. Lai, Loreto Barcos-Munoz, Yiqing Song, Sean T. Linden, Jason Surace, Marina Bianchin, Vivian U, Aaron S. Evans, Torsten Böker, Matthew A. Malkan, Kirsten L. Larson, Sabrina Stierwalt, Victorine A. Buiten, Vassilis Charmandaris, Tanio Diaz-Santos, Justin H. Howell, George C. Privon, Claudio Ricci, Paul P. van der Werf, Susanne Aalto, Christopher C. Hayward , et al. (4 additional authors not shown)

Abstract: We present James Webb Space Telescope (JWST) Integral Field Spectrograph observations of NGC 3256, a local infrared-luminous late-stage merging system with two nuclei about 1 kpc apart, both of which have evidence of cold molecular outflows. Using JWST NIRSpec and MIRI datasets, we investigate this morphologically complex system on spatial scales of $<$100 pc, where we focus on the warm molecular… ▽ More We present James Webb Space Telescope (JWST) Integral Field Spectrograph observations of NGC 3256, a local infrared-luminous late-stage merging system with two nuclei about 1 kpc apart, both of which have evidence of cold molecular outflows. Using JWST NIRSpec and MIRI datasets, we investigate this morphologically complex system on spatial scales of $<$100 pc, where we focus on the warm molecular H$_2$ gas surrounding the nuclei. We detect collimated outflowing warm H$_2$ gas originating from the southern nucleus, though we do not find significant outflowing warm H$_2$ gas surrounding the northern nucleus. Within the observed region, the maximum intrinsic velocities of the outflow reach up to $\sim$1,000 km s$^{-1}$, and extend out to a distance of 0.7 kpc. Based on H$_2$ S(7)/S(1) line ratios, we find a larger fraction of warmer gas near the S nucleus, which decreases with increasing distance from the nucleus, signifying the S nucleus as a primary source of H$_2$ heating. The gas mass of the warm H$_2$ outflow component is estimated to be $M\rm{_{warm,out}}$ = 8.9$\times$10$^5\;M_{\odot}$, as much as 4$\%$ of the cold H$_2$ mass as estimated using ALMA CO data. The outflow time scale is about $7\times10^5$ yr, resulting in a mass outflow rate of $\dot{M}\rm{_{warm,out}}$ = 1.3 M$_{\odot}$ yr$^{-1}$ and kinetic power of $P\rm{_{warm,out}}\;\sim\;2\times10^{41}$ erg s$^{-1}$. Lastly, the regions where the outflowing gas reside show high [FeII]/Pa$β$ and H$_2$/Br$γ$ line ratios, indicating enhanced mechanical heating caused by the outflows. At the same time, the 3.3 $μ$m and 6.2 $μ$m Polycyclic Aromatic Hydrocarbon fluxes in these regions are not significantly suppressed compared to those outside the outflows, suggesting the outflows have no clear negative feedback effect on the local star formation. △ Less

Submitted 21 March, 2024; originally announced March 2024.

Comments: 16 pages, 8 figures. Submitted to ApJ

Imaging of I Zw 18 by JWST. I. Strategy and First Results of Dusty Stellar Populations

Authors: Alec S. Hirschauer, Nicolas Crouzet, Nolan Habel, Laura Lenkić, Conor Nally, Olivia C. Jones, Giacomo Bortolini, Martha L. Boyer, Kay Justtanont Margaret Meixner, Göran Östlin, Gillian S. Wright, Ruyman Azzollini, Joris A. D. L. Blommaert, Bernhard Brandl, Leen Decin, Omnarayani Nayak, Pierre Royer, B. A. Sargent, Paul van der Werf

Abstract: We present a JWST imaging survey of I Zw 18, the archetypal extremely metal-poor, star-forming (SF), blue compact dwarf galaxy. With an oxygen abundance of only $\sim$3% $Z_{\odot}$, it is among the lowest-metallicity systems known in the local Universe, and is, therefore, an excellent accessible analog for the galactic building blocks which existed at early epochs of ionization and star formation… ▽ More We present a JWST imaging survey of I Zw 18, the archetypal extremely metal-poor, star-forming (SF), blue compact dwarf galaxy. With an oxygen abundance of only $\sim$3% $Z_{\odot}$, it is among the lowest-metallicity systems known in the local Universe, and is, therefore, an excellent accessible analog for the galactic building blocks which existed at early epochs of ionization and star formation. These JWST data provide a comprehensive infrared (IR) view of I Zw 18 with eight filters utilizing both Near Infrared Camera (F115W, F200W, F356W, and F444W) and Mid-Infrared Instrument (F770W, F1000W, F1500W, and F1800W) photometry, which we have used to identify key stellar populations that are bright in the near- and mid-IR. These data allow for a better understanding of the origins of dust and dust-production mechanisms in metal-poor environments by characterizing the population of massive, evolved stars in the red supergiant (RSG) and asymptotic giant branch (AGB) phases. In addition, it enables the identification of the brightest dust-enshrouded young stellar objects (YSOs), which provide insight into the formation of massive stars at extremely low metallicities typical of the very early Universe. This paper provides an overview of the observational strategy and data processing, and presents first science results, including identifications of dusty AGB, RSG, and bright YSO candidates. These first results assess the scientific quality of JWST data and provide a guide for obtaining and interpreting future observations of the dusty and evolved stars inhabiting compact dwarf SF galaxies in the local Universe. △ Less

Submitted 26 June, 2024; v1 submitted 11 March, 2024; originally announced March 2024.

Comments: Accepted for publication in The Astronomical Journal; 25 pages, 14 figures, 4 tables

JWST/MIRI unveils the stellar component of the GN20 dusty galaxy overdensity at $z$=4.05

Authors: A. Crespo Gómez, L. Colina, J. Álvarez-Márquez, A. Bik, L. Boogaard, G. Östlin, F. Peißker, F. Walter, A. Labiano, P. G. Pérez-González, T. R. Greve, G. Wright, A. Alonso-Herrero, K. I. Caputi, L. Costantin, A. Eckart, M. García-Marín, S. Gillman, J. Hjorth, E. Iani, D. Langeroodi, J. P. Pye, P. Rinaldi, T. Tikkanen, P. van der Werf , et al. (2 additional authors not shown)

Abstract: Despite the importance of the dusty star-forming galaxies (DSFGs) at $z$>2 for understanding the galaxy evolution in the early Universe, their stellar distributions traced by the near-IR emission were spatially unresolved until the arrival of the JWST. In this work we present, for the first time, a spatially-resolved morphological analysis of the rest-frame near-IR (~1.1-3.5$μ$m) emission in DSFGs… ▽ More Despite the importance of the dusty star-forming galaxies (DSFGs) at $z$>2 for understanding the galaxy evolution in the early Universe, their stellar distributions traced by the near-IR emission were spatially unresolved until the arrival of the JWST. In this work we present, for the first time, a spatially-resolved morphological analysis of the rest-frame near-IR (~1.1-3.5$μ$m) emission in DSFGs traced with the JWST/MIRI. In particular, we study the mature stellar component for the three DSFGs and a Lyman-break galaxy (LBG) present in an overdensity at $z$=4.05. Moreover, we use MIRI images along with UV to (sub)-mm ancillary photometric data to model their SEDs and extract their main physical properties. The sub-arcsec resolution MIRI images have revealed that the stellar component present a wide range of morphologies, from disc-like to compact and clump-dominated structures. These near-IR structures contrast with their UV emission, which is usually diffuse and off-centered. The SED fitting analysis shows that GN20 dominates the total SFR with a value ~2500 $M_\odot$yr$^{-1}$ while GN20.2b has the highest stellar mass in the sample ($M_*$~2$\times$10$^{11}$ $M_\odot$). The two DSFGs classified as LTGs (GN20 and GN20.2a) show high specific SFR (sSFR>30 Gyr$^{-1}$) placing them above the star-forming main sequence (SFMS) at z~4 by >0.5 dex while the ETG (i.e.,GN20.2b) is compatible with the high-mass end of the main sequence. When comparing with other DSFGs in overdensities at $z$~2-7 we observe that our objects present similar SFRs, depletion times and projected separations. Nevertheless, the effective radii computed for GN20 and GN20.2a are up to two times larger than those of isolated galaxies observed in CEERS and ALMA-HUDF at similar redshifts. We interpret this difference in size as an effect of rapid growth induced by the dense environment. △ Less

Submitted 26 June, 2024; v1 submitted 28 February, 2024; originally announced February 2024.

Comments: 15 pages, 6 figures. Submitted to A&A

A NIRCam-dark galaxy detected with the MIRI/F1000W filter in the MIDIS/JADES Hubble Ultra Deep Field

Authors: Pablo G. Pérez-González, Pierluigi Rinaldi, Karina I. Caputi, Javier Álvarez-Márquez, Marianna Annunziatella, Danial Langeroodi, Thibaud Moutard, Leindert Boogaard, Edoardo Iani, Jens Melinder, Luca Costantin, Goran Östlin, Luis Colina, Thomas R. Greve, Gillian Wright, Almudena Alonso-Herrero, Arjan Bik, Sarah E. I. Bosman, Alejandro Crespo Gómez, Daniel Dicken, Andreas Eckart, Macarena García-Marín, Steven Gillman, Manuel Güdel, Thomas Henning , et al. (10 additional authors not shown)

Abstract: We report the discovery of Cerberus, an extremely red object detected with the MIRI Deep Imaging Survey (MIDIS) observations in the F1000W filter of the Hubble Ultra Deep Field. The object is detected at $S/N\sim6$, with $\mathrm{F1000W}\sim27$ mag, and undetected in the NIRCam data gathered by the JWST Advanced Deep Extragalactic Survey, JADES, fainter than the 30.0-30.5 mag $5σ$ detection limits… ▽ More We report the discovery of Cerberus, an extremely red object detected with the MIRI Deep Imaging Survey (MIDIS) observations in the F1000W filter of the Hubble Ultra Deep Field. The object is detected at $S/N\sim6$, with $\mathrm{F1000W}\sim27$ mag, and undetected in the NIRCam data gathered by the JWST Advanced Deep Extragalactic Survey, JADES, fainter than the 30.0-30.5 mag $5σ$ detection limits in individual bands, as well as in the MIDIS F560W ultra-deep data ($\sim$29 mag, $5σ$). Analyzing the spectral energy distribution built with low-$S/N$ ($<5$) measurements in individual optical-to-mid-infrared filters and higher $S/N$ ($\gtrsim5$) in stacked NIRCam data, we discuss the possible nature of this red NIRCam-dark source using a battery of codes. We discard the possibility of Cerberus being a Solar System body based on the $<0.016$" proper motion in the 1-year apart JADES and MIDIS observations. A sub-stellar Galactic nature is deemed unlikely, given that the Cerberus' relatively flat NIRCam-to-NIRCam and very red NIRCam-to-MIRI flux ratios are not consistent with any brown dwarf model. The extragalactic nature of Cerberus offers 3 possibilities: (1) A $z\sim0.4$ galaxy with strong emission from polycyclic aromatic hydrocarbons; the very low inferred stellar mass, $\mathrm{M}_\star=10^{5-6}$ M$_\odot$, makes this possibility highly improbable. (2) A dusty galaxy at $z\sim4$ with an inferred stellar mass $\mathrm{M}_\star\sim10^{8}$ M$_\odot$. (3) A galaxy with observational properties similar to those of the reddest little red dots discovered around $z\sim7$, but Cerberus lying at $z\sim15$, with the rest-frame optical dominated by emission from a dusty torus or a dusty starburst. △ Less

Submitted 26 May, 2024; v1 submitted 26 February, 2024; originally announced February 2024.

Comments: Accepted for publication in ApJL, results unchanged, manuscript shortened to publish as a letter

CON-quest II. Spatially and spectrally resolved HCN/HCO+ line ratios in local luminous and ultraluminous infrared galaxies

Authors: Y. Nishimura, S. Aalto, M. D. Gorski, S. König, K. Onishi, C. Wethers, C. Yang, L. Barcos-Muñoz, F. Combes, T. Díaz-Santos, J. S. Gallagher, S. García-Burillo, E. González-Alfonso, T. R. Greve, N. Harada, C. Henkel, M. Imanishi, K. Kohno, S. T. Linden, J. G. Mangum, S. Martín, S. Muller, G. C. Privon, C. Ricci, F. Stanley , et al. (2 additional authors not shown)

Abstract: Nuclear regions of ultraluminous and luminous infrared galaxies (U/LIRGs) are powered by starbursts and/or active galactic nuclei (AGNs). These regions are often obscured by extremely high columns of gas and dust. Molecular lines in the submillimeter windows have the potential to determine the physical conditions of these compact obscured nuclei (CONs). We aim to reveal the distributions of HCN an… ▽ More Nuclear regions of ultraluminous and luminous infrared galaxies (U/LIRGs) are powered by starbursts and/or active galactic nuclei (AGNs). These regions are often obscured by extremely high columns of gas and dust. Molecular lines in the submillimeter windows have the potential to determine the physical conditions of these compact obscured nuclei (CONs). We aim to reveal the distributions of HCN and HCO$^+$ emission in local U/LIRGs and investigate whether and how they are related to galaxy properties. Using ALMA, we have conducted sensitive observations of the HCN J=3--2 and HCO$^+$ J=3--2 lines toward 23 U/LIRGs in the local Universe (z < 0.07) with a spatial resolution of ~0.3" (~50--400 pc). We detected both HCN and HCO$^+$ in 21 galaxies, only HCN in one galaxy, and neither in one galaxy. The global HCN/HCO$^+$ line ratios, averaged over scales of ~0.5--4 kpc, range from 0.4 to 2.3, with an unweighted mean of 1.1. These line ratios appear to have no systematic trend with bolometric AGN luminosity or star formation rate. The line ratio varies with position and velocity within each galaxy, with an average interquartile range of 0.38 on a spaxel-by-spaxel basis. In eight out of ten galaxies known to have outflows and/or inflows, we found spatially and kinematically symmetric structures of high line ratios. These structures appear as a collimated bicone in two galaxies and as a thin spherical shell in six galaxies. Non-LTE analysis suggests that the high HCN/HCO$^+$ line ratio in outflows is predominantly influenced by the abundance ratio. Chemical model calculations indicate that the enhancement of HCN abundance in outflows is likely due to high-temperature chemistry triggered by shock heating. These results imply that the HCN/HCO$^+$ line ratio can aid in identifying the outflow geometry when the shock velocity of the outflows is sufficiently high to heat the gas. △ Less

Submitted 25 April, 2024; v1 submitted 23 February, 2024; originally announced February 2024.

Comments: 52 pages, 35 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 686, A48 (2024)

Molecular isotopologue measurements toward super star clusters and the relation to their ages in NGC253 with ALCHEMI

Authors: J. Butterworth, S. Viti, P. P. Van der Werf, J. G. Mangum, S. Martín, N. Harada, K. L. Emig, S. Muller, K. Sakamoto, Y. Yoshimura, K. Tanaka, R. Herrero-Illana, L. Colzi, V. M. Rivilla, K. Y. Huang, M. Bouvier, E. Behrens, C. Henkel, Y. T. Yan, D. S. Meier, D. Zhou

Abstract: Determining the evolution of the CNO isotopes in the interstellar medium (ISM) of starburst galaxies can yield important constraints on the ages of superstar clusters (SSCs), or on other aspects and contributing factors of their evolution. Due to the time-dependent nature of the abundances of isotopes within the ISM as they are supplied from processes such as nucleosynthesis or chemical fractionat… ▽ More Determining the evolution of the CNO isotopes in the interstellar medium (ISM) of starburst galaxies can yield important constraints on the ages of superstar clusters (SSCs), or on other aspects and contributing factors of their evolution. Due to the time-dependent nature of the abundances of isotopes within the ISM as they are supplied from processes such as nucleosynthesis or chemical fractionation, this provides the possible opportunity to probe the ability of isotopes ratios to trace the ages of high star forming regions, such as SSCs. The goal of this study is to investigate whether the isotopic variations in SSC regions within NGC253 are correlated with their different ages as derived from stellar population modelling. We have measured abundance ratios of CO, HCN and HCO$^+$ isotopologues in six regions containing SSCs within NGC253 using high spatial resolution (1.6",$\sim 28$pc) data from the ALCHEMI (ALma Comprehensive High-resolution Extragalactic Molecular Inventory) ALMA Large program. We have then analysed these ratios using RADEX radiative transfer modelling, with the parameter space sampled using the nested sampling Monte Carlo algorithm MLFriends. These abundance ratios were then compared to ages predicted in each region via the fitting of observed star formation tracers (such as Br$γ$) to starburst stellar population evolution models. We do not find any significant trend with age for the CO and HCN isotopologue ratios on the timescales for the ages of the SSC* regions observed. The driving factors of these ratios within SSCs could be the Initial Mass Function as well as possibly fractionation effects. To further probe these effects in SSCs over time a larger sample of SSCs must be observed spanning a larger age range. △ Less

Submitted 16 February, 2024; originally announced February 2024.

Comments: 44 pages, 43 Figures, Accepted for Publication to A&A

Journal ref: A&A 686, A31 (2024)

JWST Observations of Starbursts: Polycyclic Aromatic Hydrocarbon Emission at the Base of the M 82 Galactic Wind

Authors: Alberto D. Bolatto, Rebecca C. Levy, Elizabeth Tarantino, Martha L. Boyer, Deanne B. Fisher, Adam K. Leroy, Serena A. Cronin, Ralf S. Klessen, J. D. Smith, Dannielle A. Berg, Torsten Boeker, Leindert A. Boogaard, Eve C. Ostriker, Todd A. Thompson, Juergen Ott, Laura Lenkic, Laura A. Lopez, Daniel A. Dale, Sylvain Veilleux, Paul P. van der Werf, Simon C. O. Glover, Karin M. Sandstrom, Evan D. Skillman, John Chisholm, Vicente Villanueva , et al. (15 additional authors not shown)

Abstract: We present new observations of the central 1 kpc of the M 82 starburst obtained with the James Webb Space Telescope (JWST) near-infrared camera (NIRCam) instrument at a resolution ~0.05"-0.1" (~1-2 pc). The data comprises images in three mostly continuum filters (F140M, F250M, and F360M), and filters that contain [FeII] (F164N), H2 v=1-0 (F212N), and the 3.3 um PAH feature (F335M). We find promine… ▽ More We present new observations of the central 1 kpc of the M 82 starburst obtained with the James Webb Space Telescope (JWST) near-infrared camera (NIRCam) instrument at a resolution ~0.05"-0.1" (~1-2 pc). The data comprises images in three mostly continuum filters (F140M, F250M, and F360M), and filters that contain [FeII] (F164N), H2 v=1-0 (F212N), and the 3.3 um PAH feature (F335M). We find prominent plumes of PAH emission extending outward from the central starburst region, together with a network of complex filamentary substructure and edge-brightened bubble-like features. The structure of the PAH emission closely resembles that of the ionized gas, as revealed in Paschen alpha and free-free radio emission. We discuss the origin of the structure, and suggest the PAHs are embedded in a combination of neutral, molecular, and photoionized gas. △ Less

Submitted 21 April, 2024; v1 submitted 29 January, 2024; originally announced January 2024.

Comments: Submitted to The Astrophysical Journal

At the end of cosmic noon: Short gas depletion times in unobscured quasars at $z \sim$ 1

Authors: M. Frias Castillo, M. Rybak, J. Hodge, P. van der Werf, L. J. Abbo, F. J. Ballieux, S. Ward, C. Harrison, G. Calistro Rivera, J. P. McKean, H. R. Stacey

Abstract: Unobscured quasars (QSOs) are predicted to be the final stage in the evolutionary sequence from gas-rich mergers to gas-depleted, quenched galaxies. Studies of this population, however, find a high incidence of far-infrared-luminous sources -suggesting significant dust-obscured star formation-but direct observations of the cold molecular gas fuelling this star formation are still necessary. We pre… ▽ More Unobscured quasars (QSOs) are predicted to be the final stage in the evolutionary sequence from gas-rich mergers to gas-depleted, quenched galaxies. Studies of this population, however, find a high incidence of far-infrared-luminous sources -suggesting significant dust-obscured star formation-but direct observations of the cold molecular gas fuelling this star formation are still necessary. We present a NOEMA study of CO(2-1) emission, tracing the cold molecular gas, in ten lensed z=1-1.5 unobscured QSOs. We detected CO(2-1) in seven of our targets, four of which also show continuum emission (λ_rest = 1.3mm). After subtracting the foreground galaxy contribution to the photometry, spectral energy distribution fitting yielded stellar masses of 10^9-11 M_\odot, with star formation rates of 25-160 M_\odot yr^-1 for the host galaxies. These QSOs have lower $L'_\mathrm{CO}$ than star-forming galaxies with the same L_IR, and show depletion times spanning a large range (50-900 Myr), but with a median of just 90 Myr. We find molecular gas masses in the range 2-40 x 10^9(alpha_CO/4) M_\odot, which suggest gas fractions above ~50% for most of the targets. Despite the presence of an unobscured QSO, the host galaxies are able to retain significant amounts of cold gas. However, with a median depletion time of ~90 Myr, the intense burst of star formation taking place in these targets will quickly deplete their molecular gas reservoirs in the absence of gas replenishment, resulting in a quiescent host galaxy. The non-detected QSOs are three of the four radio-loud QSOs in the sample, and their properties indicate that they are likely already transitioning into quiescence. Recent cosmological simulations tend to overestimate the depletion times expected for these z~1 QSO-host galaxies, which is likely linked to their difficulty producing starbursts across the general high-redshift galaxy population. △ Less

Submitted 12 January, 2024; v1 submitted 9 January, 2024; originally announced January 2024.

Comments: 20 pages. Accepted for publication in A&A

The ALCHEMI atlas: principal component analysis reveals starburst evolution in NGC 253

Authors: Nanase Harada, David S. Meier, Sergio Martín, Sebastien Muller, Kazushi Sakamoto, Toshiki Saito, Mark D. Gorski, Christian Henkel, Kunihiko Tanaka, Jeffrey G. Mangum, Susanne Aalto, Rebeca Aladro, Mathilde Bouvier, Laura Colzi, Kimberly L. Emig, Rubén Herrero-Illana, Ko-Yun Huang, Kotaro Kohno, Sabine König, Kouichiro Nakanishi, Yuri Nishimura, Shuro Takano, Víctor M. Rivilla, Serena Viti, Yoshimasa Watanabe , et al. (2 additional authors not shown)

Abstract: Molecular lines are powerful diagnostics of the physical and chemical properties of the interstellar medium (ISM). These ISM properties, which affect future star formation, are expected to differ in starburst galaxies from those of more quiescent galaxies. We investigate the ISM properties in the central molecular zone of the nearby starburst galaxy NGC 253 using the ultra-wide millimeter spectral… ▽ More Molecular lines are powerful diagnostics of the physical and chemical properties of the interstellar medium (ISM). These ISM properties, which affect future star formation, are expected to differ in starburst galaxies from those of more quiescent galaxies. We investigate the ISM properties in the central molecular zone of the nearby starburst galaxy NGC 253 using the ultra-wide millimeter spectral scan survey from the ALMA Large Program ALCHEMI. We present an atlas of velocity-integrated images at a 1".6 resolution of 148 unblended transitions from 44 species, including the first extragalactic detection of HCNH$^+$ and the first interferometric images of C$_3$H$^+$, NO, HCS$^+$. We conduct a principal component analysis (PCA) on these images to extract correlated chemical species and to identify key groups of diagnostic transitions. To the best of our knowledge, our dataset is currently the largest astronomical set of molecular lines to which PCA has been applied. The PCA can categorize transitions coming from different physical components in NGC 253 such as i) young starburst tracers characterized by high-excitation transitions of HC$_3$N and complex organic molecules (COMs) versus tracers of on-going star formation (radio recombination lines) and high-excitation transitions of CCH and CN tracing PDRs, ii) tracers of cloud-collision-induced shocks (low-excitation transitions of CH$_3$OH, HNCO, HOCO$^+$, and OCS) versus shocks from star-formation-induced outflows (high-excitation transitions of SiO), as well as iii) outflows showing emission from HOC$^+$, CCH, H$_3$O$^+$, CO isotopologues, HCN, HCO$^+$, CS, and CN. Our findings show these intensities vary with galactic dynamics, star formation activities, and stellar feedback. △ Less

Submitted 4 January, 2024; originally announced January 2024.

Comments: 65 pages, 39 figures. Accepted for publication in ApJS

The kinematics of massive high-redshift dusty star-forming galaxies

Authors: A. Amvrosiadis, J. L. Wardlow, J. E. Birkin, I. Smail, A. M. Swinbank, J. Nightingale, F. Bertoldi, W. N. Brandt, C. M. Casey, S. C. Chapman, C. -C. Chen, P. Cox, E. da Cunha, H. Dannerbauer, U. Dudzevičiūtė, B. Gullberg, J. A. Hodge, K. K. Knudsen, K. Menten, F. Walter, P. van der Werf

Abstract: 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) emissi… ▽ More 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$μ$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 $α_{\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}\,-\,σ$ plane and their relative space densities. △ Less

Submitted 14 December, 2023; originally announced December 2023.

Clumpy star formation and an obscured nuclear starburst in the luminous dusty z=4 galaxy GN20 seen by MIRI/JWST

Authors: A. Bik, J. Álvarez-Márquez, L. Colina, A. Crespo Gómez, F. Peissker, F. Walter, L. A. Boogaard, G. Östlin, T. R. Greve, G. Wright, A. Alonso-Herrero, K. I. Caputi, L. Costantin, A. Eckart, S. Gillman, J. Hjorth, E. Iani, I. Jermann, A. Labiano, D. Langeroodi, J. Melinder, P. G. Pérez-González, J. P. Pye, P. Rinaldi, T. Tikkanen , et al. (6 additional authors not shown)

Abstract: Dusty star-forming galaxies emit most of their light at far-IR to mm wavelengths as their star formation is highly obscured. Far-IR and mm observations have revealed their dust, neutral and molecular gas properties. The sensitivity of JWST at rest-frame optical and near-infrared wavelengths now allows the study of the stellar and ionized gas content. We investigate the spatially resolved distribut… ▽ More Dusty star-forming galaxies emit most of their light at far-IR to mm wavelengths as their star formation is highly obscured. Far-IR and mm observations have revealed their dust, neutral and molecular gas properties. The sensitivity of JWST at rest-frame optical and near-infrared wavelengths now allows the study of the stellar and ionized gas content. We investigate the spatially resolved distribution and kinematics of the ionized gas in GN20, a dusty star forming galaxy at $z$=4.0548. We present deep MIRI/MRS integral field spectroscopy of the near-infrared rest-frame emission of GN20. We detect spatially resolved \paa, out to a radius of 6 kpc, distributed in a clumpy morphology. The star formation rate derived from \paa\ (144 $\pm$ 9 \msunperyear) is only 7.7 $\pm 0.5 $\% of the infrared star formation rate (1860 $\pm$ 90 \msunperyear). We attribute this to very high extinction (A$_V$ = 17.2 $\pm$ 0.4 mag, or A$_{V,mixed}$ = 44 $\pm$ 3 mag), especially in the nucleus of GN20, where only faint \paa\ is detected, suggesting a deeply buried starburst. We identify four, spatially unresolved, clumps in the \paa\ emission. Based on the double peaked \paa\ profile we find that each clump consist of at least two sub-clumps. We find mass upper limits consistent with them being formed in a gravitationally unstable gaseous disk. The UV bright region of GN20 does not have any detected \paa\ emission, suggesting an age of more than 10 Myrs for this region of the galaxy. From the rotation profile of \paa\ we conclude that the gas kinematics are rotationally dominated and the $v_{rot}/σ_{m} = 3.8 \pm 1.4$ is similar to low-redshift LIRGs. We speculate that the clumps seen in GN20 could contribute to building up the inner disk and bulge of GN20. △ Less

Submitted 2 March, 2024; v1 submitted 5 December, 2023; originally announced December 2023.

Comments: 13 pages, 7 figures, accepted for publication in A&A

Journal ref: A&A 686, A3 (2024)

GOALS-JWST: Mid-Infrared Molecular Gas Excitation Probes the Local Conditions of Nuclear Star Clusters and the AGN in the LIRG VV 114

Authors: Victorine A. Buiten, Paul P. van der Werf, Serena Viti, Lee Armus, Andrew G. Barr, Loreto Barcos-Muñoz, Aaron S. Evans, Hanae Inami, Sean T. Linden, George C. Privon, Yiqing Song, Jeffrey A. Rich, Susanne Aalto, Philip N. Appleton, Torsten Böker, Vassilis Charmandaris, Tanio Diaz-Santos, Christopher C. Hayward, Thomas S. -Y. Lai, Anne M. Medling, Claudio Ricci, Vivian U

Abstract: The enormous increase in mid-IR sensitivity and spatial and spectral resolution provided by the JWST spectrographs enables, for the first time, detailed extragalactic studies of molecular vibrational bands. This opens an entirely new window for the study of the molecular interstellar medium in luminous infrared galaxies (LIRGs). We present a detailed analysis of rovibrational bands of gas-phase CO… ▽ More The enormous increase in mid-IR sensitivity and spatial and spectral resolution provided by the JWST spectrographs enables, for the first time, detailed extragalactic studies of molecular vibrational bands. This opens an entirely new window for the study of the molecular interstellar medium in luminous infrared galaxies (LIRGs). We present a detailed analysis of rovibrational bands of gas-phase CO, H$_2$O, C$_2$H$_2$ and HCN towards the heavily-obscured eastern nucleus of the LIRG VV 114, as observed by NIRSpec and MIRI MRS. Spectra extracted from apertures of 130 pc in radius show a clear dichotomy between the obscured AGN and two intense starburst regions. We detect the 2.3 $μ$m CO bandheads, characteristic of cool stellar atmospheres, in the star-forming regions, but not towards the AGN. Surprisingly, at 4.7 $\mathrmμ$m we find highly-excited CO ($T_\mathrm{ex} \approx 700-800$ K out to at least rotational level $J = 27$) towards the star-forming regions, but only cooler gas ($T_\mathrm{ex} \approx 200$ K) towards the AGN. We conclude that only mid-infrared pumping through the rovibrational lines can account for the equilibrium conditions found for CO and H$_2$O in the deeply-embedded starbursts. Here the CO bands probe regions with an intense local radiation field inside dusty young massive star clusters or near the most massive young stars. The lack of high-excitation molecular gas towards the AGN is attributed to geometric dilution of the intense radiation from the bright point source. An overview of the relevant excitation and radiative transfer physics is provided in an appendix. △ Less

Submitted 8 March, 2024; v1 submitted 4 December, 2023; originally announced December 2023.

Comments: 28 pages, 17 figures. Accepted for publication in ApJ. This version includes small revisions following the referee report

Metal and dust evolution in ALMA REBELS galaxies: insights for future JWST observations

Authors: Marco Palla, Ilse De Looze, Monica Relaño, Stefan van der Giessen, Pratika Dayal, Andrea Ferrara, Raffaella Schneider, Luca Graziani, Hiddo S. B. Algera, Manuel Aravena, Rebecca A. A. Bowler, Alexander P. S. Hygate, Hanae Inami, Ivana van Leeuwen, Rychard Bouwens, Jacqueline Hodge, Renske Smit, Mauro Stefanon, Paul van der Werf

Abstract: ALMA observations revealed the presence of significant amounts of dust in the first Gyr of Cosmic time. However, the metal and dust buildup picture remains very uncertain due to the lack of constraints on metallicity. JWST has started to reveal the metal content of high-redshift targets, which may lead to firmer constraints on high-redshift dusty galaxies evolution. In this work, we use detailed c… ▽ More ALMA observations revealed the presence of significant amounts of dust in the first Gyr of Cosmic time. However, the metal and dust buildup picture remains very uncertain due to the lack of constraints on metallicity. JWST has started to reveal the metal content of high-redshift targets, which may lead to firmer constraints on high-redshift dusty galaxies evolution. In this work, we use detailed chemical and dust evolution models to explore the evolution of galaxies within the ALMA REBELS survey, testing different metallicity scenarios that could be inferred from JWST observations. In the models, we track the buildup of stellar mass by using non-parametric SFHs for REBELS galaxies. Different scenarios for metal and dust evolution are simulated by allowing different prescriptions for gas flows and dust processes. The model outputs are compared with measured dust scaling relations, by employing metallicity-dependent calibrations for the gas mass based on the [CII]158micron line. Independently of the galaxies metal content, we found no need for extreme dust prescriptions to explain the dust masses revealed by ALMA. However, different levels of metal enrichment will lead to different dominant dust production mechanisms, with stardust production dominant over other ISM dust processes only in the metal-poor case. This points out how metallicity measurements from JWST will significantly improve our understanding of the dust buildup in high-redshift galaxies. We also show that models struggle to reproduce observables such as dust-to-gas and dust-to-stellar ratios simultaneously, possibly indicating an overestimation of the gas mass through current calibrations, especially at high metallicities. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: 16 pages + appendices, 9 Figures, 1 Table. Resubmitted to MNRAS after moderate revision

MIDIS: The Relation between Strong (Hb+[OIII]) Emission, Star Formation and Burstiness Around the Epoch of Reionization

Authors: Karina I. Caputi, Pierluigi Rinaldi, Edoardo Iani, Pablo G. Pérez-González, Göran Ostlin, Luis Colina, Thomas R. Greve, Hans-Ulrik Nørgaard-Nielsen, Gillian S. Wright, Javier Alvarez-Márquez, Andreas Eckart, Jens Hjorth, Alvaro Labiano, Olivier Le Fèvre, Fabian Walter, Paul van der Werf, Leindert Boogaard, Luca Costantin, Alejandro Crespo-Gómez, Steven Gillman, Iris Jermann, Danial Langeroodi, Jens Melinder, Florian Peissker, Manuel Güdel , et al. (3 additional authors not shown)

Abstract: We investigate the properties of strong (Hb+[OIII]) emitters before and after the end of the Epoch of Reionization from z=8 to z=5.5. We make use of ultra-deep JWST/NIRCam imaging in the Parallel Field of the MIRI Deep Imaging Survey (MIDIS) in the Hubble eXtreme Deep Field (P2-XDF), in order to select prominent (Hb+[OIII]) emitters (with rest EW_0 > 100 Angstroms) at z=5.5-7, based on their flux… ▽ More We investigate the properties of strong (Hb+[OIII]) emitters before and after the end of the Epoch of Reionization from z=8 to z=5.5. We make use of ultra-deep JWST/NIRCam imaging in the Parallel Field of the MIRI Deep Imaging Survey (MIDIS) in the Hubble eXtreme Deep Field (P2-XDF), in order to select prominent (Hb+[OIII]) emitters (with rest EW_0 > 100 Angstroms) at z=5.5-7, based on their flux density enhancement in the F356W band with respect to the spectral energy distribution continuum. We complement our selection with other (Hb+[OIII]) emitters from the literature at similar and higher (z=7-8) redshifts. We find (non-independent) anti-correlations between EW_0(Hb+[OIII]) and both galaxy stellar mass and age, in agreement with previous studies, and a positive correlation with specific star formation rate (sSFR). On the SFR-M* plane, the (Hb+[OIII]) emitters populate both the star-formation main sequence and the starburst region, which become indistinguishable at low stellar masses (log10(M*) < 7.5). We find tentative evidence for a non-monotonic relation between EW_0(Hb+[OIII]) and SFR, such that both parameters correlate with each other at SFR > 1 Msun/yr, while the correlation flattens out at lower SFRs. This suggests that low metallicities producing high EW_0(Hb+[OIII]) could be important at low SFR values. Interestingly, the properties of the strong emitters and other galaxies (33% and 67% of our z=5.5-7 sample, respectively) are similar, including, in many cases, high sSFR. Therefore, it is crucial to consider both emitters and non-emitters to obtain a complete picture of the cosmic star formation activity around the Epoch of Reionization. △ Less

Submitted 18 June, 2024; v1 submitted 21 November, 2023; originally announced November 2023.

Comments: 14 pages, 8 figures. Accepted for publication at the ApJ. Version v3 includes minor update to match journal version

Volume density structure of the NGC253 CMZ through ALCHEMI excitation analysis

Authors: Kunihiko Tanaka, Jeffrey G. Mangum, Serena Viti, Sergio Martin, Nanase Harada, Kazushi Sakamoto, Sebastien Muller, Yuki Yoshimura, Kouichiro Nakanishi, Ruben Herrero Illana, Kimberly L. Emig, S. Muhle, Hiroyuki Kaneko, Tomoka Tosaki, Erica Behrens, Victor M. Rivilla, Laura Colzi, Yuri Nishimura, P. K. Humire, Mathilde Bouvier, Ko-Yun Huang, Joshua Butterworth, David S. Meier, Paul P. van der Werf

Abstract: We present a spatially-resolved excitation analysis for the central molecular zone (CMZ) of the starburst galaxy NGC 253 using the data from the ALMA Large program ALCHEMI, whereby we explore parameters distinguishing NGC 253 from the quiescent Milky Way's Galactic Center (GC). Non-LTE analyses employing a hierarchical Bayesian framework are applied to Band 3-7 transitions from nine molecular spec… ▽ More We present a spatially-resolved excitation analysis for the central molecular zone (CMZ) of the starburst galaxy NGC 253 using the data from the ALMA Large program ALCHEMI, whereby we explore parameters distinguishing NGC 253 from the quiescent Milky Way's Galactic Center (GC). Non-LTE analyses employing a hierarchical Bayesian framework are applied to Band 3-7 transitions from nine molecular species to delineate the position-position-velocity distributions of column density ($N_\mathrm{H_2}$), volume density ($n_\mathrm{H_2}$), and temperature ($T_\mathrm{kin}$) at 27 pc resolution. Two distinct components are detected: a low-density component with $(n_\mathrm{H_2},\ T_\mathrm{kin})\sim(10^{3.3}\ \mathrm{cm}^{-3}, 85 K)$ and a high-density component with $(n_\mathrm{H_2},\ T_\mathrm{kin})\sim (10^{4.4}\ \mathrm{cm}^{-3}, 110\ \mathrm{K})$, separated at $n_\mathrm{H_2}\sim10^{3.8}\ \mathrm{cm}^{-3}$. NGC 253 has $\sim10$ times the high-density gas mass and $\sim3$ times the dense-gas mass fraction of the GC. These properties are consistent with their HCN/CO ratio but cannot alone explain the factor of $\sim30$ difference in their star formation efficiencies (SFEs), contradicting the dense-gas mass to star formation rate scaling law. The $n_\mathrm{H_2}$ histogram toward NGC 253 exhibits a shallow declining slope up to $n_\mathrm{H_2}\sim10^6\ \mathrm{cm}^{-3}$, while that of the GC steeply drops in $n_\mathrm{H_2}\gtrsim10^{4.5}\ \mathrm{cm}^{-3}$ and vanishes at $10^5\ \mathrm{cm}^{-3}$. Their dense-gas mass fraction ratio becomes consistent with their SFEs when the threshold $n_\mathrm{H_2}$ for the dense gas is taken at $\sim 10^{4.2\mbox{-}4.6}\ \mathrm{cm}^{-3}$. The rich abundance of gas above this density range in the NGC 253 CMZ, or its scarcity in the GC, is likely to be the critical difference characterizing the contrasting star formation in the centers of the two galaxies. △ Less

Submitted 20 November, 2023; originally announced November 2023.

Comments: 49 pages, 27 figures, 7 tables, accepted for publication in the Astrophysical Journal

Characterisation of Herschel-selected strong lens candidates through HST and sub-mm/mm observations

Authors: Edoardo Borsato, Lucia Marchetti, Mattia Negrello, Enrico Maria Corsini, David Wake, Aristeidis Amvrosiadis, Andrew Baker, Tom Bakx, Alexandre Beelen, Stefano Berta, David Clements, Asantha Cooray, Pierre Cox, Helmut Dannerbauer, Gianfranco de Zotti, Simon Dye, Stephen Eales, Andrea Enia, Duncan Farrah, Joaquin Gonzalez-Nuevo, David Hughes, Diana Ismail, Shuowen Jin, Andrea Lapi, Matthew Lehnert , et al. (12 additional authors not shown)

Abstract: We have carried out HST snapshot observations at 1.1 $μ$m of 281 candidate strongly lensed galaxies identified in the wide-area extragalactic surveys conducted with the Herschel space observatory. Our candidates comprise systems with flux densities at $500\,μ$m$ S_{500}\geq 80$ mJy. We model and subtract the surface brightness distribution for 130 systems, where we identify a candidate for the for… ▽ More We have carried out HST snapshot observations at 1.1 $μ$m of 281 candidate strongly lensed galaxies identified in the wide-area extragalactic surveys conducted with the Herschel space observatory. Our candidates comprise systems with flux densities at $500\,μ$m$ S_{500}\geq 80$ mJy. We model and subtract the surface brightness distribution for 130 systems, where we identify a candidate for the foreground lens candidate. After combining visual inspection, archival high-resolution observations, and lens subtraction, we divide the systems into different classes according to their lensing likelihood. We confirm 65 systems to be lensed. Of these, 30 are new discoveries. We successfully perform lens modelling and source reconstruction on 23 systems, where the foreground lenses are isolated galaxies and the background sources are detected in the HST images. All the systems are successfully modelled as a singular isothermal ellipsoid. The Einstein radii of the lenses and the magnifications of the background sources are consistent with previous studies. However, the background source circularised radii (between 0.34 kpc and 1.30 kpc) are $\sim$3 times smaller than the ones measured in the sub-mm/mm for a similarly selected and partially overlapping sample. We compare our lenses with those in the SLACS survey, confirming that our lens-independent selection is more effective at picking up fainter and diffuse galaxies and group lenses. This sample represents the first step towards characterising the near-IR properties and stellar masses of the gravitationally lensed dusty star-forming galaxies. △ Less

Submitted 2 November, 2023; originally announced November 2023.

Comments: 57 pages, 18 figures, 11 tables. Accepted for publication in MNRAS

The ALMA REBELS survey: obscured star formation in massive Lyman-break galaxies at z = 4-8 revealed by the IRX-$β$ and $M_{\star}$ relations

Authors: R. A. A. Bowler, H. Inami, L. Sommovigo, R. Smit, H. S. B. Algera, M. Aravena, L. Barrufet, R. Bouwens, E. da Cunha, F. Cullen, P. Dayal, I. de Looze, J. S. Dunlop, Y. Fudamoto, V. Mauerhofer, R. J. McLure, M. Stefanon, R. Schneider, A. Ferrara, L. Graziani, J. A. Hodge, T. Nanayakkara, M. Palla, S. Schouws, D. P. Stark , et al. (1 additional authors not shown)

Abstract: 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-… ▽ More 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 ($β\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 $β_{\rm d} = 2.0$), we find that the REBELS sources are in agreement with the local ''Calzetti-like'' starburst Infrared-excess (IRX)-$β$ 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-$β$ 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-$β$ 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. △ Less

Submitted 28 November, 2023; v1 submitted 29 September, 2023; originally announced September 2023.

Comments: 18 pages, 10 figures, 2 tables (plus 1 figure and 2 tables in the appendix). Updated to match MNRAS accepted version after minor corrections

The ALMA Reionization Era Bright Emission Line Survey (REBELS): The molecular gas content of galaxies at z~7

Authors: M. Aravena, K. E. Heintz, M. Dessauges-Zavadsky, P. A. Oesch, H. S. B. Algera, R. J. Bouwens, E. Da Cunha, P. Dayal, I. De Looze, A. Ferrara, Y. Fudamoto, V. Gonzalez, L. Graziani, H. Inami, A. Pallottini, R. Schneider, S. Schouws, L. Sommovigo, M. Topping, P. van der Werf, M. Palla

Abstract: 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$μ$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 Brigh… ▽ More 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$μ$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 ($μ_{\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 $μ_{\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$(ρ_{\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. △ Less

Submitted 29 September, 2023; v1 submitted 27 September, 2023; originally announced September 2023.

Comments: Accepted for publication in Astronomy & Astrophysics. 11 pages, 7 figures. Typo fixed in an author's name

MIDIS: Unveiling the Role of Strong Ha-emitters during the Epoch of Reionization with JWST

Authors: P. Rinaldi, K. I. Caputi, E. Iani, L. Costantin, S. Gillman, P. G. Perez-Gonzalez, G. Ostlin, L. Colina, T. R. Greve, H. U. Noorgard-Nielsen, G. S. Wright, J. Alvarez-Marquez, A. Eckart, M. Garcia-Marin, J. Hjorth, O. Ilbert, S. Kendrew, A. Labiano, O. Le Fevre, J. Pye, T. Tikkanen, F. Walter, P. van der Werf, M. Ward, M. Annunziatella , et al. (18 additional authors not shown)

Abstract: By using the ultra-deep \textit{JWST}/MIRI image at 5.6 $μm$ in the Hubble eXtreme Deep Field, we constrain the role of strong H$α$-emitters (HAEs) during Cosmic Reionization at $z\simeq7-8$. Our sample of HAEs is comprised of young ($<35\;\rm Myr$) galaxies, except for one single galaxy ($\approx 300\;\rm Myr$), with low stellar masses ($\lesssim 10^{9}\;\rm M_{\odot}$). These HAEs show a wide ra… ▽ More By using the ultra-deep \textit{JWST}/MIRI image at 5.6 $μm$ in the Hubble eXtreme Deep Field, we constrain the role of strong H$α$-emitters (HAEs) during Cosmic Reionization at $z\simeq7-8$. Our sample of HAEs is comprised of young ($<35\;\rm Myr$) galaxies, except for one single galaxy ($\approx 300\;\rm Myr$), with low stellar masses ($\lesssim 10^{9}\;\rm M_{\odot}$). These HAEs show a wide range of UV-$β$ slopes, with a median value of $β= -2.15\pm0.21$ which broadly correlates with stellar mass. We estimate the ionizing photon production efficiency ($ξ_{ion,0}$) of these sources (assuming $f_{esc,LyC} = 0\%$), which yields a median value $\rm log_{10}(ξ_{ion,0}/(Hz\;erg^{-1})) = 25.50^{+0.10}_{-0.12}$. We show that $ξ_{ion,0}$ positively correlates with EW$_{0}$(H$α$) and specific star formation rate (sSFR). Instead $ξ_{ion,0}$ weakly anti-correlates with stellar mass and $β$. Based on the $β$ values, we predict $f_{esc, LyC}=4\%^{+3}_{-2}$, which results in $\rm log_{10}(ξ_{ion}/(Hz\;erg^{-1})) = 25.55^{+0.11}_{-0.13}$. Considering this and related findings from the literature, we find a mild evolution of $ξ_{ion}$with redshift. Additionally, our results suggest that these HAEs require only modest escape fractions ($f_{esc, rel}$) of 6$-$15\% to reionize their surrounding intergalactic medium. By only considering the contribution of these HAEs, we estimated their total ionizing emissivity ($\dot{N}_{ion}$) as $\dot{N}_{ion} = 10^{50.53 \pm 0.45}; \text{s}^{-1}\text{Mpc}^{-3}$. When comparing their $\dot{N}_{ion}$ with "non-H$α$ emitter" galaxies across the same redshift range, we find that that strong, young, and low-mass emitters may have played an important role during Cosmic Reionization. △ Less

Submitted 13 June, 2024; v1 submitted 27 September, 2023; originally announced September 2023.

Comments: 21 pages, 11 Figures, 1 table. Accepted for publication in ApJ. Updated version 13/06/2024 -- found a mistake in the reported eq. 2 (now corrected)

Bright beacons? ALMA non-detection of a supposedly bright [OI] 63-um line in a redshift-6 dusty galaxy

Authors: M. Rybak, L. Lemsom, A. Lundgren, J. Zavala, J. A. Hodge, C. de Breuck, C. M. Casey, R. Decarli, K. Torstensson, J. L. Wardlow, P. P. van der Werf

Abstract: We report a non-detection of the [OI] 63-um emission line from the z = 6.03 galaxy G09.83808 using ALMA Band 9 observations, refuting the previously claimed detection with APEX by (Rybak et al. 2020); the new upper limit on the [OI] 63-um flux is almost 20-times lower. [OI] 63-um line could be a powerful tracer of neutral gas in the Epoch of Reionisation: yet our null result shows that detecting [… ▽ More We report a non-detection of the [OI] 63-um emission line from the z = 6.03 galaxy G09.83808 using ALMA Band 9 observations, refuting the previously claimed detection with APEX by (Rybak et al. 2020); the new upper limit on the [OI] 63-um flux is almost 20-times lower. [OI] 63-um line could be a powerful tracer of neutral gas in the Epoch of Reionisation: yet our null result shows that detecting [OI] 63-um from z$\geq$6 galaxies is more challenging than previously hypothesised. △ Less

Submitted 22 September, 2023; originally announced September 2023.

Journal ref: RNAAS, 7, 188 (2023)

MIDIS: JWST NIRCam and MIRI unveil the stellar population properties of Ly$α$-emitters and Lyman-Break galaxies at z ~ 3-7

Authors: Edoardo Iani, Karina I. Caputi, Pierluigi Rinaldi, Marianna Annunziatella, Leindert A. Boogaard, Göran Östlin, Luca Costantin, Steven Gillman, Pablo G. Pérez-González, Luis Colina, Gillian Wright, Almudena Alonso-Herrero, Javier Álvarez-Márquez, Arjan Bik, Sarah E. I. Bosman, Alejandro Crespo-Gómez, Andreas Eckart, Thomas R. Greve, Thomas K. Henning, Jens Hjorth, Iris Jermann, Alvaro Labiano, Danial Langeroodi, Jens Melinder, Thibaud Moutard , et al. (7 additional authors not shown)

Abstract: We study the stellar population properties of 182 spectroscopically-confirmed (MUSE/VLT) Lyman-$α$ emitters (LAEs) and 450 photometrically-selected Lyman-Break galaxies (LBGs) at z = 2.8 - 6.7 in the Hubble eXtreme Deep Field (XDF). Leveraging the combined power of HST and JWST NIRCam and MIRI observations, we analyse their rest-frame UV-through-near-IR spectral energy distributions (SEDs) with MI… ▽ More We study the stellar population properties of 182 spectroscopically-confirmed (MUSE/VLT) Lyman-$α$ emitters (LAEs) and 450 photometrically-selected Lyman-Break galaxies (LBGs) at z = 2.8 - 6.7 in the Hubble eXtreme Deep Field (XDF). Leveraging the combined power of HST and JWST NIRCam and MIRI observations, we analyse their rest-frame UV-through-near-IR spectral energy distributions (SEDs) with MIRI playing a crucial role in robustly assessing the LAE's stellar mass and ages. Our LAEs are low-mass objects (log$_{10}$(M$_\star$[M$_\odot$]) ~ 7.5), with little or no dust extinction (E(B - V) ~ 0.1) and a blue UV continuum slope ($β$ ~ -2.2). While 75% of our LAEs are young (< 100 Myr), the remaining 25% have significantly older stellar populations (> 100 Myr). These old LAEs are statistically more massive, less extinct and have lower specific star formation rate (sSFR) compared to young LAEs. Besides, they populate the M$_\star$ - SFR plane along the main-sequence (MS) of star-forming galaxies, while young LAEs populate the starburst region. The comparison between the LAEs properties to those of a stellar-mass matched sample of LBGs shows no statistical difference between these objects, except for the LBGs redder UV continuum slope and marginally larger E(B - V) values. Interestingly, 48% of the LBGs have ages < 10 Myr and are classified as starbursts, but lack detectable Ly$α$ emission. This is likely due to HI resonant scattering and/or selective dust extinction. Overall, we find that JWST observations are crucial in determining the properties of LAEs and shedding light on the properties and similarities between LAEs and LBGs. △ Less

Submitted 25 January, 2024; v1 submitted 15 September, 2023; originally announced September 2023.

Comments: 22 pages, 15 figures, 5 tables. Accepted for publication in APJ

Spatially-resolved H$α$ and ionizing photon production efficiency in the lensed galaxy MACS1149-JD1 at a redshift of 9.11

Authors: J. Álvarez-Márquez, L. Colina, A. Crespo Gómez, P. Rinaldi, J. Melinder, G. Östlin, M. Annunziatella, A. Labiano, A. Bik, S. Bosman, T. R. Greve, G. Wright, A. Alonso-Herrero, L. Boogaard, R. Azollini, K. I. Caputi, L. Costantin, A. Eckart, M. GarcÍa-MarÍn, S. Gillman, J. Hjorth, E. Iani, O. Ilbert, I. Jermann, D. Langeroodi , et al. (10 additional authors not shown)

Abstract: We present MIRI/JWST medium-resolution spectroscopy (MRS) and imaging (MIRIM) of the lensed galaxy MACS1149-JD1 at a redshift of $z$=9.1092$\pm$0.0002 (Universe age about 530 Myr). We detect, for the first time, spatially resolved H$α$ emission in a galaxy at a redshift above nine. The structure of the H$α$ emitting gas consists of two clumps, S and N. The total H$α$ luminosity implies an instanta… ▽ More We present MIRI/JWST medium-resolution spectroscopy (MRS) and imaging (MIRIM) of the lensed galaxy MACS1149-JD1 at a redshift of $z$=9.1092$\pm$0.0002 (Universe age about 530 Myr). We detect, for the first time, spatially resolved H$α$ emission in a galaxy at a redshift above nine. The structure of the H$α$ emitting gas consists of two clumps, S and N. The total H$α$ luminosity implies an instantaneous star-formation of 5.3$\pm$0.4 $M_{\odot}$ yr$^{-1}$ for solar metallicities. The ionizing photon production efficiency, $\log(ζ_\mathrm{ion})$, shows a spatially resolved structure with values of 25.55$\pm$0.03, 25.47$\pm$0.03, and 25.91$\pm$0.09 Hz erg$^{-1}$ for the integrated galaxy, and clumps S and N, respectively. The H$α$ rest-frame equivalent width, EW$_{0}$(H$α$), is 726$^{+660}_{-182}$ Ángstrom for the integrated galaxy, but presents extreme values of 531$^{+300}_{-96}$ Ángstrom and $\geq$1951 Ángstrom for clumps S and N, respectively. The spatially resolved ionizing photon production efficiency is within the range of values measured in galaxies at redshift above six, and well above the canonical value (25.2$\pm$0.1 Hz erg$^{-1}$). The extreme difference of EW$_{0}$(H$α$) for Clumps S and N indicates the presence of a recent (<5 Myrs) burst in clump N and a star formation over a larger period of time (e.g., $\sim$50 Myr) in clump S. Finally, clump S and N show very different H$α$ kinematics with velocity dispersions of 56$\pm$4 km s$^{-1}$ and 113$\pm$33 km s$^{-1}$, likely indicating the presence of outflows or increased turbulence in the clump N. The dynamical mass, $M_\mathrm{dyn}$= (2.4$\pm$0.5)$\times$10$^{9}$ $M_{\odot}$, is within the range measured with the spatially resolved [OIII]88$μ$m line. △ Less

Submitted 24 March, 2024; v1 submitted 12 September, 2023; originally announced September 2023.

Comments: 11 pages, 6 figures, published version in A&A

Journal ref: A&A 686, A85 (2024)

NOEMA observations of GN-z11: Constraining Neutral Interstellar Medium and Dust Formation in the Heart of Cosmic Reionization at $z=10.6$

Authors: Y. Fudamoto, P. A. Oesch, F. Walter, R. Decarli, C. L. Carilli, A. Ferrara, L. Barrufet, R. Bouwens, M. Dessauges-Zavadsky, E. J. Nelson, H. Dannerbauer, G. Illingworth, A. K. Inoue, R. Marques-Chaves, I. Pérez-Fournon, D. A. Riechers, D. Schaerer, R. Smit, Y. Sugahara, P. van der Werf

Abstract: We present results of dust continuum and [CII]$\,158\,{\rm μm}$ emission line observations of a remarkably UV-luminous ($M_{\rm UV}=-21.6$) galaxy at $z=10.603$: GN-z11. Using the Northern Extended Millimeter Array (NOEMA), observations have been carried out over multiple observing cycles. We achieved a high sensitivity resulting in a $λ_{\rm rest}=160\,{\rm μm}$ continuum $1\,σ$ depth of… ▽ More We present results of dust continuum and [CII]$\,158\,{\rm μm}$ emission line observations of a remarkably UV-luminous ($M_{\rm UV}=-21.6$) galaxy at $z=10.603$: GN-z11. Using the Northern Extended Millimeter Array (NOEMA), observations have been carried out over multiple observing cycles. We achieved a high sensitivity resulting in a $λ_{\rm rest}=160\,{\rm μm}$ continuum $1\,σ$ depth of $13.0\,\rm{μJy/beam}$ and a [CII] emission line $1\,σ$ sensitivity of $31\,\rm{mJy/beam\,km/s}$ using $50\,\rm{km/s}$ binning with a $\sim 2\,{\rm arcsec}$ synthesized beam. Neither dust continuum nor [CII]$\,158\,{\rm μm}$ line emission are detected at the expected frequency of $ν_{\rm [CII]} = 163.791\,\rm{GHz}$ and the sky location of GN-z11. The upper limits show that GN-z11 is neither luminous in $L_{\rm IR}$ nor $L_{\rm [CII]}$, with a dust mass $3\,σ$ limit of ${\rm log}(M_{\rm dust}/{\rm M_{\odot}}) < 6.5-6.9$ and with a [CII] based molecular gas mass $3\,σ$ limit of ${\rm log}(M_{\rm mol,[CII]}/{\rm M_{\odot}}) < 9.3$. Together with radiative transfer calculations, we also investigated the possible cause of the dust poor nature of the GN-z11 showed by the blue color in the UV continuum of GN-z11 ($β_{\rm UV}=-2.4$), and found that $\gtrsim3\times$ deeper observations are crucial to study dust production at very high-redshift. Nevertheless, our observations show the crucial role of deep mm/submm observations of very high redshift galaxies to constrain multiple phases in the interstellar medium. △ Less

Submitted 5 September, 2023; originally announced September 2023.

Comments: submitted to MNRAS, 7 pages, 4 figures, 1 table

MIDIS: JWST/MIRI reveals the Stellar Structure of ALMA-selected Galaxies in the Hubble-UDF at Cosmic Noon

Authors: Leindert A. Boogaard, Steven Gillman, Jens Melinder, Fabian Walter, Luis Colina, Göran Östlin, Karina I. Caputi, Edoardo Iani, Pablo Pérez-González, Paul van der Werf, Thomas R. Greve, Gillian Wright, Almudena Alonso-Herrero, Javier Álvarez-Márquez, Marianna Annunziatella, Arjan Bik, Sarah Bosman, Luca Costantin, Alejandro Crespo Gómez, Dan Dicken, Andreas Eckart, Jens Hjorth, Iris Jermann, Alvaro Labiano, Danial Langeroodi , et al. (8 additional authors not shown)

Abstract: We present deep James Webb Space Telescope (JWST)/MIRI F560W observations of a flux-limited, ALMA-selected sample of 28 galaxies at z=0.5-3.6 in the Hubble Ultra Deep Field (HUDF). The data from the MIRI Deep Imaging Survey (MIDIS) reveal the stellar structure of the HUDF galaxies at rest-wavelengths of >1 micron for the first time. We revise the stellar mass estimates using new JWST photometry an… ▽ More We present deep James Webb Space Telescope (JWST)/MIRI F560W observations of a flux-limited, ALMA-selected sample of 28 galaxies at z=0.5-3.6 in the Hubble Ultra Deep Field (HUDF). The data from the MIRI Deep Imaging Survey (MIDIS) reveal the stellar structure of the HUDF galaxies at rest-wavelengths of >1 micron for the first time. We revise the stellar mass estimates using new JWST photometry and find good agreement with pre-JWST analysis; the few discrepancies can be explained by blending issues in the earlier lower-resolution Spitzer data. At z~2.5, the resolved rest-frame near-infrared (1.6 micron) structure of the galaxies is significantly more smooth and centrally concentrated than seen by HST at rest-frame 450 nm (F160W), with effective radii of Re(F560W)=1-5 kpc and Sérsic indices mostly close to an exponential (disk-like) profile (n~1), up to n~5 (excluding AGN). We find an average size ratio of Re(F560W)/Re(F160W)~0.7 that decreases with stellar mass. The stellar structure of the ALMA-selected galaxies is indistinguishable from a HUDF reference sample of galaxies with comparable MIRI flux density. We supplement our analysis with custom-made, position-dependent, empirical PSF models for the F560W observations. The results imply that a smoother stellar structure is in place in massive gas-rich, star-forming galaxies at Cosmic Noon, despite a more clumpy rest-frame optical appearance, placing additional constraints on galaxy formation simulations. As a next step, matched-resolution, resolved ALMA observations will be crucial to further link the mass- and light-weighted galaxy structures to the dusty interstellar medium. △ Less

Submitted 26 April, 2024; v1 submitted 31 August, 2023; originally announced August 2023.

Comments: 19 pages, 12 figures, 1 table, Accepted for publication in ApJ

SUNRISE: The rich molecular inventory of high-redshift dusty galaxies revealed by broadband spectral line surveys

Authors: Chentao Yang, Alain Omont, Sergio Martín, Thomas G. Bisbas, Pierre Cox, Alexandre Beelen, Eduardo González-Alfonso, Raphaël Gavazzi, Susanne Aalto, Paola Andreani, Cecilia Ceccarelli, Yu Gao, Mark Gorski, Michel Guélin, Hai Fu, R. J. Ivison, Kirsten K. Knudsen, Matthew Lehnert, Hugo Messias, Sebastien Muller, Roberto Neri, Dominik Riechers, Paul van der Werf, Zhi-Yu Zhang

Abstract: Understanding the nature of high-$z$ dusty galaxies requires a comprehensive view of their ISM and molecular complexity. However, the molecular ISM at high-$z$ is commonly studied using only a few species beyond CO, limiting our understanding. In this paper, we present the results of deep 3 mm spectral line surveys using the NOEMA targeting two lensed dusty galaxies: APM 08279+5255 (APM), a quasar… ▽ More Understanding the nature of high-$z$ dusty galaxies requires a comprehensive view of their ISM and molecular complexity. However, the molecular ISM at high-$z$ is commonly studied using only a few species beyond CO, limiting our understanding. In this paper, we present the results of deep 3 mm spectral line surveys using the NOEMA targeting two lensed dusty galaxies: APM 08279+5255 (APM), a quasar at redshift $z=3.911$, and NCv1.143 (NC), a $z=3.565$ starburst galaxy. The spectral line surveys cover rest-frame frequencies from about 330-550 GHz. We report the detection of 38 and 25 emission lines in APM and NC, respectively. The spectra reveal the chemical richness and the complexity of the physical properties of the ISM. By comparing the spectra of the two sources and combining the gas excitation analysis, we find that the physical properties and the chemical imprints of the ISM are different between them: the molecular gas is more excited in APM, exhibiting higher molecular-gas temperatures and densities compared to NC; the chemical abundances in APM are akin to the values of local AGN, showing boosted relative abundances of the dense gas tracers that might be related to high-temperature chemistry and/or XDRs, while NC more closely resembles local starburst galaxies. The most significant differences are found in H2O, where the 448GHz H2O line is significantly brighter in APM, likely linked to the intense far-infrared radiation from the dust powered by AGN. Our astrochemical model suggests that at such high column densities, FUV radiation is less important in regulating the ISM, while CRs (X-rays/shocks) are the key players in shaping the abundance of the molecules and the initial conditions of star formation. Such deep spectral line surveys open a new window to study the physical and chemical properties of the ISM and the radiation field of galaxies in the early Universe. (abridged) △ Less

Submitted 22 October, 2023; v1 submitted 14 August, 2023; originally announced August 2023.

Comments: Accepted for publication in A&A. 39 pages (including 9 pages of Appendix). Updated after language editing

Journal ref: A&A 680, A95 (2023)

PRUSSIC II -- ALMA imaging of dense-gas tracers in SDP.81: Evidence for low mechanical heating and a sub-solar metallicity in a z=3.04 dusty galaxy

Authors: M. Rybak, J. van Marrewijk, J. A. Hodge, P. Andreani, G. Calistro Rivera, L. Graziani, J. P. McKean, S. Viti, P. P. van der Werf

Abstract: We present deep ALMA Band 3 observations of the HCN, HCO+, and HNC (4-3) emission in SDP.81, a well-studied z = 3.042 strongly lensed galaxy. These lines trace the high-density gas, which remains almost entirely unexplored in z$\geq$1 galaxies. Additionally, these dense-gas tracers are potentially powerful diagnostics of the mechanical heating of the interstellar medium. While the HCN(4-3) and HNC… ▽ More We present deep ALMA Band 3 observations of the HCN, HCO+, and HNC (4-3) emission in SDP.81, a well-studied z = 3.042 strongly lensed galaxy. These lines trace the high-density gas, which remains almost entirely unexplored in z$\geq$1 galaxies. Additionally, these dense-gas tracers are potentially powerful diagnostics of the mechanical heating of the interstellar medium. While the HCN(4-3) and HNC(4-3) lines are not detected, the HCO+(4-3) emission is clearly detected and resolved. This is the third detection of this line in a high-redshift star-forming galaxy. We find an unusually high HCO+/HCN intensity ratio of $\geq$2.2. Based on the photodissociation region modelling, the most likely explanation for the elevated HCO+/HCN ratio is that SDP.81 has low mechanical heating - less than 10% of the total energy budget - and a sub-solar metallicity, Z=0.5 Z$_\odot$. While such conditions might not be representative of the general population of high-redshift dusty galaxies, lower-than-solar metallicity might have a significant impact on gas masses inferred from CO observations. In addition, we report the detection of CO(0-1) absorption from the foreground lensing galaxy and CO(1-0) emission from a massive companion to the lensing galaxy, approximately 50 kpc to the southeast. △ Less

Submitted 25 September, 2023; v1 submitted 5 August, 2023; originally announced August 2023.

Comments: A&A accepted, corrected version

Journal ref: A&A 679, A119 (2023)

GOALS-JWST: Gas Dynamics and Excitation in NGC7469 revealed by NIRSpec

Authors: Marina Bianchin, Vivian U, Yiqing Song, Thomas S. -Y. Lai, Raymond P. Remigio, Loreto Barcos-Munoz, Tanio Diaz-Santos, Lee Armus, Hanae Inami, Kirsten L. Larson, Aaron S. Evans, Torsten Boker, Justin A. Kader, Sean T. Linden, Vassilis Charmandaris, Matthew A. Malkan, Jeff Rich, Thomas Bohn, Anne M. Medling, Sabrina Stierwalt, Joseph M. Mazzarella, David R. Law, George C. Privon, Susanne Aalto, Philip Appleton , et al. (14 additional authors not shown)

Abstract: We present new JWST-NIRSpec IFS data for the luminous infrared galaxy NGC7469: a nearby (70.6Mpc) active galaxy with a Sy 1.5 nucleus that drives a highly ionized gas outflow and a prominent nuclear star-forming ring. Using the superb sensitivity and high spatial resolution of the JWST instrument NIRSpec-IFS, we investigate the role of the Seyfert nucleus in the excitation and dynamics of the circ… ▽ More We present new JWST-NIRSpec IFS data for the luminous infrared galaxy NGC7469: a nearby (70.6Mpc) active galaxy with a Sy 1.5 nucleus that drives a highly ionized gas outflow and a prominent nuclear star-forming ring. Using the superb sensitivity and high spatial resolution of the JWST instrument NIRSpec-IFS, we investigate the role of the Seyfert nucleus in the excitation and dynamics of the circumnuclear gas. Our analysis focuses on the [Fe ii], H2, and hydrogen recombination lines that trace the radiation/shocked-excited molecular and ionized ISM around the AGN. We investigate the gas excitation through H2/Brγ and [Fe ii]/Pa\b{eta} emission line ratios and find that photoionization by the AGN dominates within the central 300 pc of the galaxy and together with a small region show ing signatures of shock-heated gas; these shock-heated regions are likely associated with a compact radio jet. In addition, the velocity field and velocity dispersion maps reveal complex gas kinematics. Rotation is the dominant feature, but we also identify non-circular motions consistent with gas inflows as traced by the velocity residuals and the spiral pattern in the Paα velocity dispersion map. The inflow is consistent with the mass outflow rate and two orders of magnitude higher than the AGN accretion rate. The compact nuclear radio jet has enough power to drive the highly ionized outflow. This scenario suggests that the inflow and outflow are in a self-regulating feeding-feedback process, with a contribution from the radio jet helping to drive the outflow. △ Less

Submitted 15 February, 2024; v1 submitted 31 July, 2023; originally announced August 2023.

Comments: 15 pages, 4 figures, 3 tables, accepted for publication in ApJ

z-GAL -- A NOEMA spectroscopic redshift survey of bright Herschel galaxies: [III] Physical properties

Authors: S. Berta, F. Stanley, D. Ismail, P. Cox, R. Neri, C. Yang, A. J. Young, S. Jin, H. Dannerbauer, T. J. Bakx, A. Beelen, A. Weiss, A. Nanni, A. Omont, P. van der Werf, M. Krips, A. J. Baker, G. Bendo, E. Borsato, V. Buat, K. M. Butler, N. Chartab, A. Cooray, S. Dye, S. Eales , et al. (13 additional authors not shown)

Abstract: The z-GAL survey observed 137 bright Herschel-selected targets with the IRAM NOrthern Extended Millimeter Array, with the aim to measure their redshift and study their properties. Several of them have been resolved into multiple sources. Consequently, robust spectroscopic redshifts have been measured for 165 individual galaxies in the range 0.8<z<6.5. In this paper we analyse the millimetre spectr… ▽ More The z-GAL survey observed 137 bright Herschel-selected targets with the IRAM NOrthern Extended Millimeter Array, with the aim to measure their redshift and study their properties. Several of them have been resolved into multiple sources. Consequently, robust spectroscopic redshifts have been measured for 165 individual galaxies in the range 0.8<z<6.5. In this paper we analyse the millimetre spectra of the z-GAL sources, using both their continuum and line emission to derive their physical properties. At least two spectral lines are detected for each source, including transitions of 12CO, [CI], and H2O. The observed 12CO line ratios and spectral line energy distributions of individual sources resemble those of local starbursts. In seven sources the para-H2O(2_11-2_02) transition is detected and follows the IR versus H2O luminosity relation of sub-millimetre galaxies. The molecular gas mass of the z-GAL sources is derived from their 12CO, [CI], and sub-millimetre dust continuum emission. The three tracers lead to consistent results, with the dust continuum showing the largest scatter when compared to 12CO. The gas-to-dust mass ratio of these sources was computed by combining the information derived from 12CO and the dust continuum and has a median value of 107, similar to star-forming galaxies of near-solar metallicity. The same combined analysis leads to depletion timescales in the range between 0.1 and 1.0 Gyr, which place the z-GAL sources between the `main sequence' of star formation and the locus of starbursts. Finally, we derived a first estimate of stellar masses - modulo possible gravitational magnification - by inverting known gas scaling relations: the z-GAL sample is confirmed to be mostly composed by starbursts, whereas ~25% of its members lie on the main sequence of star-forming galaxies (within +/- 0.5 dex). △ Less

Submitted 28 July, 2023; originally announced July 2023.

Comments: Accepted for publication on A&A; 26 pages; 12 figures

Journal ref: A&A 678, A28 (2023)

z-GAL -- A NOEMA spectroscopic redshift survey of bright Herschel galaxies: [II] Dust properties

Authors: D. Ismail, A. Beelen, V. Buat, S. Berta, P. Cox, F. Stanley, A. Young, S. Jin, R. Neri, T. Bakx, H. Dannerbauer, K. Butler, A. Cooray, A. Nanni, A. Omont, S. Serjeant, P. van der Werf, C. Vlahakis, A. Weiss, C. Yang, A. J. Baker, G. Bendo, E. Borsato, N. Chartab, S. Dye , et al. (12 additional authors not shown)

Abstract: (Abridged) We present the dust properties of 125 bright Herschel galaxies selected from the z-GAL survey. The large instantaneous bandwidth of NOEMA provides an exquisite sampling of the underlying dust continuum emission at 2 and 3 mm in the observed frame, with flux densities in at least four side bands for each source. Together with the available Herschel 250, 350, and 500 micron and SCUBA-2 85… ▽ More (Abridged) We present the dust properties of 125 bright Herschel galaxies selected from the z-GAL survey. The large instantaneous bandwidth of NOEMA provides an exquisite sampling of the underlying dust continuum emission at 2 and 3 mm in the observed frame, with flux densities in at least four side bands for each source. Together with the available Herschel 250, 350, and 500 micron and SCUBA-2 850 micron flux densities, the spectral energy distribution of each source can be analyzed from the far-infrared to the millimeter, with a fine sampling of the Rayleigh-Jeans tail. This wealth of data provides a solid basis to derive robust dust properties, in particular the dust emissivity index, beta, and the dust temperature, T(dust). In order to demonstrate our ability to constrain the dust properties, we used a flux-generated mock catalog and analyzed the results under the assumption of an optically thin and optically thick modified black body emission. For the z-GAL sources, we report a range of dust emissivities with beta ~ 1.5 - 3 estimated up to high precision with relative uncertainties that vary in the range 7% - 15%, and an average of 2.2 +/- 0.3. We find dust temperatures varying from 20 to 50 K with an average of T(dust) ~ 30 K for the optically thin case and ~38 K in the optically thick case. For all the sources, we estimate the dust masses and apparent infrared luminosities (based on the optically thin approach). An inverse correlation is found between T(dust) and beta, which is similar to what is seen in the local Universe. Finally, we report an increasing trend in the dust temperature as a function of redshift at a rate of 6.5 +/- 0.5 K/z for this 500 micron-selected sample. Based on this study, future prospects are outlined to further explore the evolution of dust temperature across cosmic time. △ Less

Submitted 28 July, 2023; originally announced July 2023.

Comments: Accepted for publication on A&A; 32 pages; 26 figures

Journal ref: A&A 678, A27 (2023)

z-GAL -- A NOEMA spectroscopic redshift survey of bright Herschel galaxies: [I] Overview

Authors: P. Cox, R. Neri, S. Berta, D. Ismail, F. Stanley, A. Young, S. Jin, T. Bakx, A. Beelen, H. Dannerbauer, M. Krips, M. Lehnert, A. Omont, D. A. Riechers, A. J. Baker, G. Bendo, E. Borsato, V. Buat, K. Butler, N. Chartab, A. Cooray, S. Dye, S. Eales, R. Gavazzi, D. Hughes , et al. (13 additional authors not shown)

Abstract: (Abridged) Using the IRAM NOEMA interferometer, we measures the redshifts of 126 bright galaxies detected in the Herschel H-ATLAS, HeLMS, and HerS surveys. We report reliable spectroscopic redshifts for a total of 124 of the Herschel-selected galaxies. The redshifts are estimated from scans of the 3 and 2-mm bands (and, in one case, the 1-mm band) and are based on the detection of at least two emi… ▽ More (Abridged) Using the IRAM NOEMA interferometer, we measures the redshifts of 126 bright galaxies detected in the Herschel H-ATLAS, HeLMS, and HerS surveys. We report reliable spectroscopic redshifts for a total of 124 of the Herschel-selected galaxies. The redshifts are estimated from scans of the 3 and 2-mm bands (and, in one case, the 1-mm band) and are based on the detection of at least two emission lines. Together with the Pilot Programme (Neri et al. 2020), including spectroscopic redshifts of 11 sources, our survey has derived precise redshifts for 135 bright Herschel-selected galaxies, making it the largest sample of high-z galaxies with robust redshifts to date. Most emission lines detected are from 12CO (mainly from J=2-1 to 5-4), with some sources seen in [CI] and H2O emission lines. The spectroscopic redshifts are in the range 0.8<z<6.55 with a median value of z=2.56 +/- 0.10. The line widths of the sources are large, with a mean value for the full width at half maximum Delta(V) of 590 +/- 25 km/s and with 35% of the sources having widths of 700 km/s < Delta(V) < 1800 km/s. Most of the sources are unresolved or barely resolved on scales of 2 to 3 arcsec (or linear sizes of 15-25 kpc, unlensed). Some fields reveal double or multiple sources and, in some cases, sources at different redshifts. Taking these sources into account, there are, in total, 165 individual sources with robust spectroscopic redshifts, including lensed galaxies, binary systems, and over-densities. We present an overview of the z-GAL survey and provide the observed properties of the emission lines, the derived spectroscopic redshifts, and an atlas of the entire sample. The data presented here will serve as a foundation for the other z-GAL papers in this series reporting on the dust emission, the molecular and atomic gas properties, and a detailed analysis of the nature of the sources. △ Less

Submitted 28 July, 2023; originally announced July 2023.

Comments: Accepted for publication on A&A; 63 pages

Journal ref: A&A 678, A26 (2023)

GOALS-JWST: Small neutral grains and enhanced 3.3 micron PAH emission in the Seyfert galaxy NGC 7469

Authors: Thomas S. -Y. Lai, Lee Armus, Marina Bianchin, Tanio Diaz-Santos, Sean T. Linden, George C. Privon, Hanae Inami, Vivian U, Thomas Bohn, Aaron S. Evans, Kirsten L. Larson, Brandon S. Hensley, J. -D. T. Smith, Matthew A. Malkan, Yiqing Song, Sabrina Stierwalt, Paul P. van der Werf, Jed McKinney, Susanne Aalto, Victorine A. Buiten, Jeff Rich, Vassilis Charmandaris, Philip Appleton, Loreto Barcos-Munoz, Torsten Boker , et al. (14 additional authors not shown)

Abstract: We present James Webb Space Telescope (JWST) Near Infrared Spectrograph (NIRSpec) integral-field spectroscopy of the nearby luminous infrared galaxy, NGC 7469. We take advantage of the high spatial/spectral resolution and wavelength coverage of JWST /NIRSpec to study the 3.3 um neutral polycyclic aromatic hydrocarbon (PAH) grain emission on ~60 pc scales. We find a clear change in the average grai… ▽ More We present James Webb Space Telescope (JWST) Near Infrared Spectrograph (NIRSpec) integral-field spectroscopy of the nearby luminous infrared galaxy, NGC 7469. We take advantage of the high spatial/spectral resolution and wavelength coverage of JWST /NIRSpec to study the 3.3 um neutral polycyclic aromatic hydrocarbon (PAH) grain emission on ~60 pc scales. We find a clear change in the average grain properties between the star-forming ring and the central AGN. Regions in the vicinity of the AGN, with [NeIII]/[NeII]>0.25, tend to have larger grain sizes and lower aliphatic-to-aromatic (3.4/3.3) ratios indicating that smaller grains are preferentially removed by photo-destruction in the vicinity of the AGN. We find an overall suppression of the total PAH emission relative to the ionized gas in the central 1 kpc region of the AGN in NGC 7469 compared to what has been observed with Spitzer on 3 kpc scales. However, the fractional 3.3 um to total PAH power is enhanced in the starburst ring, possibly due to a variety of physical effects on sub-kpc scales, including recurrent fluorescence of small grains or multiple photon absorption by large grains. Finally, the IFU data show that while the 3.3 um PAH-derived star formation rate (SFR) in the ring is 8% higher than that inferred from the [NeII] and [NeIII] emission lines, the integrated SFR derived from the 3.3 um feature would be underestimated by a factor of two due to the deficit of PAHs around the AGN, as might occur if a composite system like NGC 7469 were to be observed at high-redshift. △ Less

Submitted 27 July, 2023; originally announced July 2023.

Comments: 14 pages, 5 figures, 2 tables, Submitted to ApJL

First rest-frame infrared spectrum of a z>7 quasar: JWST/MRS observations of J1120+0641

Authors: Sarah E. I. Bosman, Javier Álvarez-Márquez, Luis Colina, Fabian Walter, Almudena Alonso-Herrero, Martin J. Ward, Göran Östlin, Thomas R. Greve, Gillian Wright, Arjan Bik, Leindert Boogaard, Karina I. Caputi, Luca Costantin, Andreas Eckart, Macarena García-Marín, Steven Gillman, Manuel Güdel, Thomas Henning, Jens Hjorth, Edoardo Iani, Olivier Ilbert, Iris Jermann, Alvaro Labiano, Pierre-Olivier Lagage, Danial Langeroodi , et al. (7 additional authors not shown)

Abstract: We present a JWST/MRS spectrum of the quasar J1120+0641 at z=7.0848, the first spectroscopic observation of a reionisation-era quasar in the rest-frame infrared ($0.6<λ<3.4μ$m). In the context of the mysterious fast assembly of the first supermassive black holes at z>7, our observations enable for the first time the detection of hot torus dust, the H$α$ emission line, and the Paschen-series broad… ▽ More We present a JWST/MRS spectrum of the quasar J1120+0641 at z=7.0848, the first spectroscopic observation of a reionisation-era quasar in the rest-frame infrared ($0.6<λ<3.4μ$m). In the context of the mysterious fast assembly of the first supermassive black holes at z>7, our observations enable for the first time the detection of hot torus dust, the H$α$ emission line, and the Paschen-series broad emission lines in a quasar at z>7. Hot torus dust is clearly detected as an upturn in the continuum emission at $λ_{\text{rest}}\simeq1.3μ$m, leading to a black-body temperature of $T=1413.5^{+5.7}_{-7.4}$K. Compared to similarly-luminous quasars at 0<z<6, the hot dust in J1120+0641 is somewhat elevated in temperature (top 1%). The temperature is more typical among 6<z<6.5 quasars (top 25%), leading us to postulate a weak evolution in the hot dust temperature at z>6 ($2σ$ significance). We measure the black hole mass of J1120+0641 based on the H$α$ Balmer line, $M_{\text{BH}}=1.52\pm0.17\cdot 10^9 M_\odot$, which is in good agreement with the previous rest-UV MgII black hole mass measurement. The black hole mass based on the Paschen-series lines is also consistent, indicating no significant extinction in the rest-frame UV measurement. The broad H$α$, Pa-$α$ and Pa-$β$ emission lines are consistent with an origin in a common broad-line region (BLR) with density log$N_H/\text{cm}^{-3}\geq 12$, ionisation parameter $-7<$log$U<-4$, and extinction E(B-V)$\lesssim 0.1$mag. These BLR parameters are consistent with similarly-bright quasars at 0<z<4. Overall, we find that both J1120+0641's hot dust torus and hydrogen BLR properties show no significant peculiarity when compared to luminous quasars down to z=0. The quasar accretion structures must have therefore assembled very quickly, as they appear fully "mature" less than 760 million years after the Big Bang. △ Less

Submitted 26 July, 2023; originally announced July 2023.

Comments: 19 pages, 10 figures; submitted

Neutral outflows in high-z QSOs

Authors: Kirsty M. Butler, Paul P. van der Werf, Alain Omont, Pierre Cox

Abstract: OH+ absorption is a powerful tracer of inflowing and outflowing gas in the predominantly atomic diffuse and turbulent halo surrounding galaxies. In this letter, we present observations of OH+(1_1-1_0), CO(9-8) and the underlying dust continuum in 5 strongly lensed z~2-4 QSOs, using ALMA to detect outflowing neutral gas. Blue-shifted OH+ absorption is detected in 3/5 QSOs and tentatively detected i… ▽ More OH+ absorption is a powerful tracer of inflowing and outflowing gas in the predominantly atomic diffuse and turbulent halo surrounding galaxies. In this letter, we present observations of OH+(1_1-1_0), CO(9-8) and the underlying dust continuum in 5 strongly lensed z~2-4 QSOs, using ALMA to detect outflowing neutral gas. Blue-shifted OH+ absorption is detected in 3/5 QSOs and tentatively detected in a 4th. Absorption at systemic velocities is also detected in one. OH+ emission is observed in 3/5 QSOs at systemic velocities and CO(9-8) is detected in all 5 QSOs at high S/N, providing information on the dense molecular gas within the host galaxy. We compare our sample to high-z far-infrared (FIR) luminous star-forming and active galaxies from the literature. We find no difference in OH+ absorption line properties between active and star-forming galaxies with both samples following the same optical depth-dust temperature relation, suggesting that these observables are driven by the same mechanism in both samples. Similarly, star-forming and active galaxies both follow the same OH+ emission-FIR relation. Obscured QSOs display broader (>800 km/s) emission than the unobscured QSOs and all but one of the high-z star-forming galaxies, likely caused by the warm molecular gas reservoir obscuring the accreting nucleus. Broader CO(9-8) emission (>500 km/s) is found in obscured versus unobscured QSOs, but overall cover a similar range in line widths as the star-forming galaxies and follow the CO(9-8)-FIR luminosity relation found in low-z galaxies. We find that outflows traced by OH+ are only detected in extreme star-forming galaxies (broad CO emission) and in both types of QSOs, which, in turn, display no red-shifted absorption. This suggests that diffuse neutral outflows in galaxy halos may be associated with the most energetic evolutionary phases leading up to and following the obscured QSO phase. △ Less

Submitted 6 May, 2023; originally announced May 2023.

Comments: 8 pages, 3 figures, 4 tables, accepted to A&A letters

Journal ref: A&A 674, L5 (2023)

Uncovering the stellar structure of the dusty star-forming galaxy GN20 at z=4.055 with MIRI/JWST

Authors: L. Colina, A. Crespo Gómez, J. Álvarez-Márquez, A. Bik, F. Walter, L. Boogaard, A. Labiano, F. Peissker, P. Pérez-González, G. Östlin, T. R. Greve, H. U. Nørgaard-Nielsen, G. Wright, A. Alonso-Herrero, R. Azollini, K. I. Caputi, D. Dicken, M. García-Marín, J. Hjorth, O. Ilbert, S. Kendrew, J. P. Pye, T. Tikkanen, P. van der Werf, L. Costantin , et al. (13 additional authors not shown)

Abstract: Luminous infrared galaxies at high redshifts ($z$>4) include extreme starbursts that build their stellar mass over short periods of time (>100 Myr). These galaxies are considered to be the progenitors of massive quiescent galaxies at intermediate redshifts ($z\sim$2) but their stellar structure and buildup is unknown. Here, we present the first spatially resolved near-infrared imaging of GN20, one… ▽ More Luminous infrared galaxies at high redshifts ($z$>4) include extreme starbursts that build their stellar mass over short periods of time (>100 Myr). These galaxies are considered to be the progenitors of massive quiescent galaxies at intermediate redshifts ($z\sim$2) but their stellar structure and buildup is unknown. Here, we present the first spatially resolved near-infrared imaging of GN20, one of the most luminous dusty star-forming galaxies known to date, observed at an epoch when the Universe was only 1.5 Gyr old. The 5.6$μ$m image taken with the JWST Mid-Infrared Instrument (MIRI/JWST) shows that GN20 is a very luminous galaxy (M$_\mathrm{1.1μm,AB}$=$-$25.01), with a stellar structure composed of a conspicuous central source and an extended envelope. The central source is an unresolved nucleus that carries 9% of the total flux. The nucleus is co-aligned with the peak of the cold dust emission, and offset by 3.9 kpc from the ultraviolet stellar emission. The diffuse stellar envelope is similar in size to the clumpy CO molecular gas distribution. The centroid of the stellar envelope is offset by 1 kpc from the unresolved nucleus, suggesting GN20 is involved in an interaction or merger event supported by its location as the brightest galaxy in a proto-cluster. The stellar size of GN20 is larger by a factor of about 3-5 than known spheroids, disks, and irregulars at $z\sim$4, while its size and low Sérsic index are similar to those measured in dusty, infrared luminous galaxies at $z\sim$2 of the same mass. GN20 has all the ingredients necessary for evolving into a massive spheroidal quiescent galaxy at intermediate $z$: it is a large, luminous galaxy at $z$=4.05 involved in a short and massive starburst centred in the stellar nucleus and extended over the entire galaxy, out to radii of 4 kpc, and likely induced by the interaction or merger with a member of the proto-cluster. △ Less

Submitted 26 April, 2023; originally announced April 2023.

Comments: 7 pages, 4 figures

Journal ref: A&A 673, L6 (2023)

The ALMA REBELS Survey: Discovery of a massive, highly star-forming and morphologically complex ULIRG at $z =7.31$

Authors: A. P. S. Hygate, J. A. Hodge, E. da Cunha, M. Rybak, S. Schouws, H. Inami, M. Stefanon, L. Graziani, R. Schneider, P. Dayal, R. J. Bouwens, R. Smit, R. A. A. Bowler, R. Endsley, V. Gonzalez, P. A. Oesch, D. P. Stark, H. S. B. Algera, M. Aravena, L. Barrufet, A. Ferrara, Y. Fudamoto, J. H. A, I. De Looze, T. Nanayakkara , et al. (5 additional authors not shown)

Abstract: We present Atacama Large Millimeter/Submillimeter Array (ALMA) [CII] and $\sim158$ $\rmμ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 sig… ▽ More We present Atacama Large Millimeter/Submillimeter Array (ALMA) [CII] and $\sim158$ $\rmμ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$. △ Less

Submitted 18 April, 2023; originally announced April 2023.

Comments: Accepted for publication in MNRAS. 21 pages, 8 figures

Reconstructing the shock history in the CMZ of NGC 253 with ALCHEMI

Authors: K. -Y. Huang, S. Viti, J. Holdship, J. G. Mangum, S. Martín, N. Harada, S. Muller, K. Sakamoto, K. Tanaka, Y. Yoshimura, R. Herrero-Illana, D. S. Meier, E. Behrens, P. P. van der Werf, C. Henkel, S. García-Burillo, V. M. Rivilla, K. L. Emig, L. Colzi, P. K. Humire, R. Aladro, M. Bouvier

Abstract: HNCO and SiO are well known shock tracers and have been observed in nearby galaxies, including the nearby (D=3.5 Mpc) starburst galaxy NGC 253. The simultaneous detection of these two species in regions where the star formation rate is high may be used to study the shock history of the gas. We perform a multi-line molecular study using these two shock tracers (SiO and HNCO) with the aim of charact… ▽ More HNCO and SiO are well known shock tracers and have been observed in nearby galaxies, including the nearby (D=3.5 Mpc) starburst galaxy NGC 253. The simultaneous detection of these two species in regions where the star formation rate is high may be used to study the shock history of the gas. We perform a multi-line molecular study using these two shock tracers (SiO and HNCO) with the aim of characterizing the gas properties. We also explore the possibility of reconstructing the shock history in NGC 253's Central Molecular Zone (CMZ). Six SiO transitions and eleven HNCO transitions were imaged at high resolution $1''.6$ (28 pc) with the Atacama Large Millimeter/submillimeter Array (ALMA) as part of the ALCHEMI Large Programme. Both non-LTE radiative transfer analysis and chemical modelling were performed in order to characterize the gas properties, and to investigate the chemical origin of the emission. The non-LTE radiative transfer analysis coupled with Bayesian inference shows clear evidence that the gas traced by SiO has different densities and temperatures than that traced by HNCO, with an indication that shocks are needed to produce both species. Chemical modelling further confirms such a scenario and suggests that fast and slow shocks are responsible for SiO and HNCO production, respectively, in most GMCs. We are also able to infer the physical characteristics of the shocks traced by SiO and HNCO for each GMC. Radiative transfer and chemical analysis of the SiO and HNCO in the CMZ of NGC 253 reveal a complex picture whereby most of the GMCs are subjected to shocks. We speculate on the possible shock scenarios responsible for the observed emission and provide potential history and timescales for each shock scenario. Higher spatial resolution observations of these two species are required in order to quantitatively differentiate between scenarios. △ Less

Submitted 22 March, 2023; originally announced March 2023.

Journal ref: A&A 675, A151 (2023)

The ALMA REBELS Survey: The First Infrared Luminosity Function Measurement at $\mathbf{z \sim 7}

Authors: L. Barrufet, P. A. Oesch, R. Bouwens, H. Inami, L. Sommovigo, H. Algera, E. da Cunha, M. Aravena, P. Dayal, A. Ferrara, Y. Fudamoto, V. Gonzalez, L. Graziani, A. Hygate, I. de Looze, T. Nanayakkara, A. Pallottini, R. Schneider, M. Stefanon, M. Topping, P. van Der Werf

Abstract: 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… ▽ More 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. △ Less

Submitted 20 March, 2023; originally announced March 2023.

Comments: 9 pages, 5 figures

Bright Extragalactic ALMA Redshift Survey (BEARS) III: Detailed study of emission lines from 71 Herschel targets

Authors: M. Hagimoto, T. J. L. C. Bakx, S. Serjeant, G. J. Bendo, S. A. Urquhart, S. Eales, K. C. Harrington, Y. Tamura, H. Umehata, S. Berta, A. R. Cooray, P. Cox, G. De Zotti, M. D. Lehnert, D. A. Riechers, D. Scott, P. Temi, P. P. van der Werf, C. Yang, A. Amvrosiadis, P. M. Andreani, A. J. Baker, A. Beelen, E. Borsato, V. Buat , et al. (33 additional authors not shown)

Abstract: We analyse the molecular and atomic emission lines of 71 bright Herschel-selected galaxies between redshifts 1.4 to 4.6 detected by the Atacama Large Millimetre/submillimetre Array. These lines include a total of 156 CO, [C I], and H2O emission lines. For 46 galaxies, we detect two transitions of CO lines, and for these galaxies we find gas properties similar to those of other dusty star-forming g… ▽ More We analyse the molecular and atomic emission lines of 71 bright Herschel-selected galaxies between redshifts 1.4 to 4.6 detected by the Atacama Large Millimetre/submillimetre Array. These lines include a total of 156 CO, [C I], and H2O emission lines. For 46 galaxies, we detect two transitions of CO lines, and for these galaxies we find gas properties similar to those of other dusty star-forming galaxy (DSFG) samples. A comparison to photo-dissociation models suggests that most of Herschel-selected galaxies have similar interstellar medium conditions as local infrared-luminous galaxies and high-redshift DSFGs, although with denser gas and more intense far-ultraviolet radiation fields than normal star-forming galaxies. The line luminosities agree with the luminosity scaling relations across five orders of magnitude, although the star-formation and gas surface density distributions (i.e., Schmidt-Kennicutt relation) suggest a different star-formation phase in our galaxies (and other DSFGs) compared to local and low-redshift gas-rich, normal star-forming systems. The gas-to-dust ratios of these galaxies are similar to Milky Way values, with no apparent redshift evolution. Four of 46 sources appear to have CO line ratios in excess of the expected maximum (thermalized) profile, suggesting a rare phase in the evolution of DSFGs. Finally, we create a deep stacked spectrum over a wide rest-frame frequency (220-890 GHz) that reveals faint transitions from HCN and CH, in line with previous stacking experiments. △ Less

Submitted 8 March, 2023; originally announced March 2023.

Comments: 30 pages, 17 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal. Comments are warmly welcomed

VLA Legacy Survey of Molecular Gas in Massive Star-forming Galaxies at High Redshift

Authors: Marta Frias Castillo, Jacqueline Hodge, Matus Rybak, Paul van der Werf, Ian Smail, Jack Birkin, Chian-Chou Chen, Scott Chapman, Ryley Hill, Claudia del P. Lagos, Cheng-Lin Liao, Elisabete da Cunha, Gabriela Calistro Rivera, Jianhang Chen, Eric Jimenez Andrade, Eric Murphy, Douglas Scott, Mark Swinbank, Fabian Walter, Rob Ivison, Helmut Dannerbauer

Abstract: We present initial results of an ongoing survey with the Karl G. Jansky Very Large Array targeting the CO($J$ = 1-0) transition in a sample of 30 submillimeter-selected, dusty star-forming galaxies at $z =$ 2-5 with existing mid--$J$ CO detections from ALMA and NOEMA, of which 17 have been fully observed. We detect CO(1-0) emission in 11 targets, along with three tentative ($\sim$1.5-2$σ$) detecti… ▽ More We present initial results of an ongoing survey with the Karl G. Jansky Very Large Array targeting the CO($J$ = 1-0) transition in a sample of 30 submillimeter-selected, dusty star-forming galaxies at $z =$ 2-5 with existing mid--$J$ CO detections from ALMA and NOEMA, of which 17 have been fully observed. We detect CO(1-0) emission in 11 targets, along with three tentative ($\sim$1.5-2$σ$) detections; three galaxies are undetected. Our results yield total molecular gas masses of 6-23$\times$10$^{10}$ ($α_\mathrm{CO}$/1) M$_\odot$, with gas mass fractions, $f_\mathrm{gas}$=$M_\mathrm{mol}$/($M_*$+$M_\mathrm{mol}$), of 0.1-0.8 and a median depletion time of (140$\pm$70) Myr. We find median CO excitation ratios of $r_{31}$ = 0.75$\pm$0.39 and $r_{41}$ = 0.63$\pm$0.44, with a significant scatter. We find no significant correlation between the excitation ratio and a number of key parameters such as redshift, CO(1-0) line width or $Σ_\mathrm{SFR}$. We only find a tentative positive correlation between $r_{41}$ and the star-forming efficiency, but we are limited by our small sample size. Finally, we compare our results to predictions from the SHARK semi-analytical model, finding a good agreement between the molecular gas masses, depletion times and gas fractions of our sources and their SHARK counterparts. Our results highlight the heterogeneous nature of the most massive star-forming galaxies at high-redshift, and the importance of CO(1--0) observations to robustly constrain their total molecular gas content and ISM properties. △ Less

Submitted 10 February, 2023; v1 submitted 7 February, 2023; originally announced February 2023.

Comments: 10 pages, 15 figures. Accepted for publication in ApJ