Showing 1–50 of 110 results for author: Spilker, J

Robust Nitrogen and Oxygen Abundances of Haro 3 from Optical and Infrared Emission

Authors: Yuguang Chen, Tucker Jones, Ryan L. Sanders, Dario Fadda, Jessica Sutter, Robert Minchin, Nikolaus Z. Prusinski, Sunny Rhoades, Keerthi Vasan GC, Charles C. Steidel, Erin Huntzinger, Paige Kelly, Danielle A. Berg, Fabio Bresolin, Rodrigo Herrera-Camus, Ryan J. Rickards Vaught, Guido Roberts-Borsani, Peter Senchyna, Justin S. Spilker, Daniel P. Stark, Benjamin Weiner, D. Christopher Martin, Mateusz Matuszewski, Rosalie C. McGurk, James D. Neill

Abstract: Accurate chemical compositions of star-forming regions are a critical diagnostic tool to characterize the star formation history and gas flows which regulate galaxy formation. However, the abundance discrepancy factor (ADF) between measurements from the "direct" optical electron temperature ($T_e$) method and from the recombination lines (RL) represents $\sim0.2$ dex systematic uncertainty in oxyg… ▽ More Accurate chemical compositions of star-forming regions are a critical diagnostic tool to characterize the star formation history and gas flows which regulate galaxy formation. However, the abundance discrepancy factor (ADF) between measurements from the "direct" optical electron temperature ($T_e$) method and from the recombination lines (RL) represents $\sim0.2$ dex systematic uncertainty in oxygen abundance. The degree of uncertainty for other elements is unknown. We conduct a comprehensive analysis of O$^{++}$ and N$^+$ ion abundances using optical and far-infrared spectra of a star-forming region within the nearby dwarf galaxy Haro 3, which exhibits a typical ADF. Assuming homogeneous conditions, the far-IR emission indicates an O abundance which is higher than the $T_e$ method and consistent with the RL value, as would be expected from temperature fluctuations, whereas the N abundance is too large to be explained by temperature fluctuations. Instead a component of highly obscured gas is likely required to explain the high far-IR to optical flux ratios. Accounting for this obscured component reduces both the IR-based metallicities and the inferred magnitude of temperature fluctuations, such that they cannot fully explain the ADF in Haro 3. Additionally, we find potential issues when predicting the RL fluxes from current atomic data. Our findings underscore the critical importance of resolving the cause of abundance discrepancies and understanding the biases between different metallicity methods. This work represents a promising methodology, and we identify further approaches to address the current dominant uncertainties. △ Less

Submitted 28 May, 2024; originally announced May 2024.

Comments: 18 pages, 7 figures, and 2 tables. Submitted to ApJ

Single Aperture Large Telescope for Universe Studies (SALTUS): Science Overview

Authors: Gordon Chin, Carrie M. Anderson, Jennifer Bergner, Nicolas Biver, Gordon L. Bjoraker, Thibault Cavalie, Michael DiSanti, Jian-Rong Gao, Paul Hartogh, Leon K. Harding, Qing Hu, Daewook Kim, Craig Kulesa, Gert de Lange, David T. Leisawitz, Rebecca C. Levy, Arthur Lichtenberger, Daniel P. Marronh, Joan Najita, Trent Newswander, George H. Rieke, Dimitra Rigopoulou, Peter Roefsema, Nathan X. Roth, Kamber Schwarz , et al. (11 additional authors not shown)

Abstract: The SALTUS Probe mission will provide a powerful far-infrared (far-IR) pointed space observatory to explore our cosmic origins and the possibility of life elsewhere. The observatory employs an innovative deployable 14-m aperture, with a sunshield that will radiatively cool the off-axis primary to <45K. This cooled primary reflector works in tandem with cryogenic coherent and incoherent instruments… ▽ More The SALTUS Probe mission will provide a powerful far-infrared (far-IR) pointed space observatory to explore our cosmic origins and the possibility of life elsewhere. The observatory employs an innovative deployable 14-m aperture, with a sunshield that will radiatively cool the off-axis primary to <45K. This cooled primary reflector works in tandem with cryogenic coherent and incoherent instruments that span the 34 to 660 micron far-IR range at both high and moderate spectral resolutions. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: 49 pages, 10 figures, 3 tables, submitted to SPIE JATIS

Milky Way and Nearby Galaxies Science with the Single Aperture Large Telescope for Universe Studies (\textit{SALTUS}) Space Observatory

Authors: Rebecca C. Levy, Alexander Tielens, Justin Spilker, Daniel P. Marrone, Desika Narayanan, Christopher K. Walker

Abstract: This paper presents an overview of the Milky Way and nearby galaxies science case for the \textit{Single Aperture Large Telescope for Universe Studies} (SALTUS) far-infrared NASA probe-class mission concept. SALTUS offers enormous gains in spatial resolution and spectral sensitivity over previous far-IR missions, thanks to its cold ($<$40~K) 14-m primary mirror. Key Milky Way and nearby galaxies s… ▽ More This paper presents an overview of the Milky Way and nearby galaxies science case for the \textit{Single Aperture Large Telescope for Universe Studies} (SALTUS) far-infrared NASA probe-class mission concept. SALTUS offers enormous gains in spatial resolution and spectral sensitivity over previous far-IR missions, thanks to its cold ($<$40~K) 14-m primary mirror. Key Milky Way and nearby galaxies science goals for SALTUS focus on understanding the role of star formation in feedback in the Local Universe. In addition to this science case, SALTUS would open a new window to to of Galactic and extragalactic communities in the 2030s, enable fundamentally new questions to be answered, and be a far-IR analog to the near- and mid-IR capabilities of JWST. This paper summarizes the Milky Way and nearby galaxies science case and plans for notional observing programs in both guaranteed and guest (open) time. △ Less

Submitted 25 April, 2024; originally announced April 2024.

Comments: Submitted to JATIS as part of a special issue on SALTUS

High-Redshift Extragalactic Science with the Single Aperture Large Telescope for Universe Studies (SALTUS) Space Observatory

Authors: Justin Spilker, Rebecca C. Levy, Daniel Marrone, Stacey Alberts, Scott C. Chapman, Mark Dickinson, Eiichi Egami, Ryan Endsley, Desika Narayanan, George Rieke, Antony A. Stark, Alexander Tielens, Christopher K. Walker

Abstract: This paper presents an overview of the high-redshift extragalactic science case for the Single Aperture Large Telescope for Universe Studies (SALTUS) far-infrared NASA probe-class mission concept. Enabled by its 14m primary reflector, SALTUS offers enormous gains in spatial resolution and spectral sensitivity over previous far-IR missions. SALTUS would be a versatile observatory capable of respond… ▽ More This paper presents an overview of the high-redshift extragalactic science case for the Single Aperture Large Telescope for Universe Studies (SALTUS) far-infrared NASA probe-class mission concept. Enabled by its 14m primary reflector, SALTUS offers enormous gains in spatial resolution and spectral sensitivity over previous far-IR missions. SALTUS would be a versatile observatory capable of responding to the scientific needs of the extragalactic community in the 2030s, and a natural follow-on to the near- and mid-IR capabilities of JWST. Key early-universe science goals for SALTUS focus on understanding the role of galactic feedback processes in regulating galaxy growth across cosmic time, and charting the rise of metals and dust from the early universe to the present. This paper summarizes these science cases and the performance metrics most relevant for high-redshift observations. △ Less

Submitted 25 April, 2024; originally announced April 2024.

Comments: Submitted to JATIS as part of a special issue on SALTUS. Key performance metrics in Figs. 1 and 2

Medium Bands, Mega Science: a JWST/NIRCam Medium-Band Imaging Survey of Abell 2744

Authors: Katherine A. Suess, John R. Weaver, Sedona H. Price, Richard Pan, Bingjie Wang, Rachel Bezanson, Gabriel Brammer, Sam E. Cutler, Ivo Labbe, Joel Leja, Christina C. Williams, Katherine E. Whitaker, Pratika Dayal, Anna de Graaff, Robert Feldmann, Marijn Franx, Yoshinobu Fudamoto, Seiji Fujimoto, Lukas J. Furtak, Andy D. Goulding, Jenny E. Greene, Gourav Khullar, Vasily Kokorev, Mariska Kriek, Brian Lorenz , et al. (17 additional authors not shown)

Abstract: In this paper, we describe the "Medium Bands, Mega Science" JWST Cycle 2 survey (JWST-GO-4111) and demonstrate the power of these data to reveal both the spatially-integrated and spatially-resolved properties of galaxies from the local universe to the era of cosmic dawn. Executed in November 2023, MegaScience obtained ~30 arcmin^2 of deep multiband NIRCam imaging centered on the z~0.3 Abell 2744 c… ▽ More In this paper, we describe the "Medium Bands, Mega Science" JWST Cycle 2 survey (JWST-GO-4111) and demonstrate the power of these data to reveal both the spatially-integrated and spatially-resolved properties of galaxies from the local universe to the era of cosmic dawn. Executed in November 2023, MegaScience obtained ~30 arcmin^2 of deep multiband NIRCam imaging centered on the z~0.3 Abell 2744 cluster, including eleven medium-band filters and the two shortest-wavelength broad-band filters, F070W and F090W. Together, MegaScience and the UNCOVER Cycle 1 treasury program provide a complete set of deep (~28-30 mag) images in all NIRCam medium- and broad-band filters. This unique dataset allows us to precisely constrain photometric redshifts, map stellar populations and dust attenuation for large samples of distant galaxies, and examine the connection between galaxy structures and formation histories. MegaScience also includes ~17 arcmin^2 of NIRISS parallel imaging in two broad-band and four medium-band filters from 0.9-4.8um, expanding the footprint where robust spectral energy distribution (SED) fitting is possible. We provide example SEDs and multi-band cutouts at a variety of redshifts, and use a catalog of JWST spectroscopic redshifts to show that MegaScience improves both the scatter and catastrophic outlier rate of photometric redshifts by factors of 2-3. Additionally, we demonstrate the spatially-resolved science enabled by MegaScience by presenting maps of the [OIII] line emission and continuum emission in three spectroscopically-confirmed z>6 galaxies. We show that line emission in reionization-era galaxies can be clumpy, extended, and spatially offset from continuum emission, implying that galaxy assembly histories are complex even at these early epochs. We publicly release fully reduced mosaics and photometric catalogs for both the NIRCam primary and NIRISS parallel fields. △ Less

Submitted 19 April, 2024; originally announced April 2024.

Comments: 21 pages, 10 figures. Fully reduced imaging, photometric catalogs, and photometric redshift fits publicly available at https://jwst-uncover.github.io/megascience/

The ALMA-CRISTAL survey: Extended [CII] emission in an interacting galaxy system at z ~ 5.5

Authors: A. Posses, M. Aravena, J. González-López, N. M. Förster Schreiber, D. Liu, L. Lee, M. Solimano, T. Díaz-Santos, R. J. Assef, L. Barcos-Muñoz, S. Bovino, R. A. A. Bowler, G. Calistro Rivera, E. da Cunha, R. L. Davies, M. Killi, I. De Looze, A. Ferrara, D. B. Fisher, R. Herrera-Camus, R. Ikeda, T. Lambert, J. Li, D. Lutz, I. Mitsuhashi , et al. (9 additional authors not shown)

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

Submitted 5 March, 2024; originally announced March 2024.

Comments: Submitted to A&A - comments are welcome! - 19 pages, 13 figures

Reply to: Effects of density and temperature variations on the metallicity of Mrk 71

Authors: Yuguang Chen, Tucker Jones, Ryan Sanders, Dario Fadda, Jessica Sutter, Robert Minchin, Erin Huntzinger, Peter Senchyna, Daniel Stark, Justin Spilker, Benjamin Weiner, Guido Roberts-Borsani

Abstract: In Chen et al., 2023 (C23; arXiv:2304.09898), we introduced a new method to directly measure temperature fluctuations and applied it to a nearby dwarf galaxy, Mrk 71, finding a temperature fluctuation parameter $t^2 = 0.008\pm 0.043$. This result is lower by $\sim 2σ$ than the value required to explain the abundance discrepancy (AD) in this object. In the Matters Arising article submitted by Mende… ▽ More In Chen et al., 2023 (C23; arXiv:2304.09898), we introduced a new method to directly measure temperature fluctuations and applied it to a nearby dwarf galaxy, Mrk 71, finding a temperature fluctuation parameter $t^2 = 0.008\pm 0.043$. This result is lower by $\sim 2σ$ than the value required to explain the abundance discrepancy (AD) in this object. In the Matters Arising article submitted by Mendez-Delgado et al. (arXiv:2310.01197), the authors claim that using the same data presented in C23 in a different way, it is possible to conclude that the measurements are consistent with a larger $t^2 \simeq 0.1$ inferred indirectly from recombination lines (RLs). However, this requires a higher density such that the infrared [O III] 52 $μ$m and [O III] 88 $μ$m lines -- which form the basis of the direct measurement method -- are mutually inconsistent. Moreover, to reach agreement between the direct $t^2$ measurement and the larger $t^2$ value inferred from RLs requires systematically varying four parameters by $\sim 1σ$ from their best-determined values, which collectively amount to a $\sim2σ$ difference, consistent with the significance ($\sim 2 σ$) originally reported in C23. Therefore, we conclude that the results of C23 hold, and that the combined optical and infrared [O III] data disfavour $t^2 \simeq 0.1$ at the $\approx2σ$ level in Mrk 71. Future work is nonetheless warranted to better understand the AD associated with both optical and infrared emission line analysis. △ Less

Submitted 5 February, 2024; originally announced February 2024.

Comments: Response requested by the editorial board. Peer-reviewed version published in Nature Astronomy

TEMPLATES: Direct Abundance Constraints for Two Lensed Lyman-Break Galaxies

Authors: Brian Welch, Grace M. Olivier, Taylor A. Hutchison, Jane R. Rigby, Danielle A. Berg, Manuel Aravena, Matthew B. Bayliss, Jack E. Birkin, Scott C. Chapman, Håkon Dahle, Gourav Khullar, Keunho J. Kim, Guillaume Mahler, Matthew A. Malkan, Desika Narayanan, Kedar A. Phadke, Keren Sharon, J. D. T. Smith, Manuel Solimano, Justin S. Spilker, Joaquin D. Viera, David Vizgan

Abstract: Using integrated spectra for two gravitationally lensed galaxies from the JWST TEMPLATES Early Release Science program, we analyze faint auroral lines, which provide direct measurements of the gas-phase chemical abundance. For the brighter galaxy, SGAS1723$+$34 ($z = 1.3293$), we detect the [OIII]$\lambda4363$, [SIII]$\lambda6312$, and [OII]$λλ$7320,7330 auroral emission lines, and set an upper li… ▽ More Using integrated spectra for two gravitationally lensed galaxies from the JWST TEMPLATES Early Release Science program, we analyze faint auroral lines, which provide direct measurements of the gas-phase chemical abundance. For the brighter galaxy, SGAS1723$+$34 ($z = 1.3293$), we detect the [OIII]$\lambda4363$, [SIII]$\lambda6312$, and [OII]$λλ$7320,7330 auroral emission lines, and set an upper limit for the [NII]$\lambda5755$ line. For the second galaxy, SGAS1226$+$21 ($z = 2.925$), we do not detect any auroral lines, and report upper limits. With these measurements and upper limits, we constrain the electron temperatures in different ionization zones within both of these galaxies. For SGAS1723$+$34, where auroral lines are detected, we calculate direct oxygen and nitrogen abundances, finding an N/O ratio consistent with observations of nearby ($z\sim 0$) galaxies. These observations highlight the potent combination of JWST and gravitational lensing to measure faint emission lines in individual distant galaxies and to directly study the chemical abundance patterns in those galaxies. △ Less

Submitted 23 January, 2024; originally announced January 2024.

Comments: Submitted to the Astrophysical Journal's Focus Issue on the TEMPLATES JWST Early Release Science Program. 15 pages, 6 figures. Comments welcome

The ALMA-CRISTAL survey. Discovery of a 15 kpc-long gas plume in a $z=4.54$ Lyman-$α$ blob

Authors: M. Solimano, J. González-López, M. Aravena, R. Herrera-Camus, I. De Looze, N. M. Förster Schreiber, J. Spilker, K. Tadaki, R. J. Assef, L. Barcos-Muñoz, R. L. Davies, T. Díaz-Santos, A. Ferrara, D. B. Fisher, L. Guaita, R. Ikeda, E. J. Johnston, D. Lutz, I. Mitsuhashi, C. Moya-Sierralta, M. Relaño, T. Naab, A. C. Posses, K. Telikova, H. Übler , et al. (2 additional authors not shown)

Abstract: Massive star-forming galaxies in the high-redshift universe host large reservoirs of cold gas in their circumgalactic medium (CGM). Traditionally, these reservoirs have been linked to diffuse H I Lyman-$α$ (Ly$α)$ emission extending beyond $\approx 10$ kpc scales. In recent years, millimeter/submillimeter observations are starting to identify even colder gas in the CGM through molecular and/or ato… ▽ More Massive star-forming galaxies in the high-redshift universe host large reservoirs of cold gas in their circumgalactic medium (CGM). Traditionally, these reservoirs have been linked to diffuse H I Lyman-$α$ (Ly$α)$ emission extending beyond $\approx 10$ kpc scales. In recent years, millimeter/submillimeter observations are starting to identify even colder gas in the CGM through molecular and/or atomic tracers such as the [C II] $158\,μ$m transition. In this context, we study the well-known J1000+0234 system at $z=4.54$ that hosts a massive dusty star-forming galaxy (DSFG), a UV-bright companion, and a Ly$α$ blob. We combine new ALMA [C II] line observations taken by the CRISTAL survey with data from previous programs targeting the J1000+0234 system, and achieve a deep view into a DSFG and its rich environment at a 0.2" resolution. We identify an elongated [C II]-emitting structure with a projected size of 15 kpc stemming from the bright DSFG at the center of the field, with no clear counterpart at any other wavelength. The plume is oriented $\approx 40^{\circ}$ away from the minor axis of the DSFG, and shows significant spatial variation of its spectral parameters. In particular, the [C II] emission shifts from 180 km/s to 400 km/s between the bottom and top of the plume, relative to the DSFG's systemic velocity. At the same time, the line width starts at 400-600 km/s but narrows down to 190 km/s at top end of the plume. We discuss four possible scenarios to interpret the [C II] plume: a conical outflow, a cold accretion stream, ram pressure stripping, and gravitational interactions. While we cannot strongly rule out any of these with the available data, we disfavor the ram pressure stripping scenario due to the requirement of special hydrodynamic conditions. △ Less

Submitted 9 January, 2024; originally announced January 2024.

Comments: 17 pages (14 main text, 2 for references and 1 appendix page), 7 figures and 4 tables. Submitted to A&A

JWST Early Release Science Program TEMPLATES: Targeting Extremely Magnified Panchromatic Lensed Arcs and their Extended Star formation

Authors: Jane R. Rigby, Joaquin D. Vieira, Kedar A. Phadke, Taylor A. Hutchison, Brian Welch, Jared Cathey, Justin S. Spilker, Anthony H. Gonzalez, Prasanna Adhikari, M. Aravena, Matthew B. Bayliss, Jack E. Birkin, Emmy Bursk, Scott C. Chapman, Håkon Dahle, Lauren A. Elicker, Travis C. Fischer, Michael K. Florian, Michael D. Gladders, Christopher C. Hayward, Rose Hewald, Lily A. Kettler, Gourav Khullar, Seonwoo Kim, David R. Law , et al. (11 additional authors not shown)

Abstract: This paper gives an overview of TEMPLATES, a JWST Early Release Science program that targeted four extremely bright, gravitationally lensed galaxies: two extremely dusty, two with low attenuation, as templates for galaxy evolution studies with JWST. TEMPLATES obtains a common set of spectral diagnostics for these 1.3 < z < 4.2 galaxies, in particular H alpha, Paschen alpha, and the rest-frame opti… ▽ More This paper gives an overview of TEMPLATES, a JWST Early Release Science program that targeted four extremely bright, gravitationally lensed galaxies: two extremely dusty, two with low attenuation, as templates for galaxy evolution studies with JWST. TEMPLATES obtains a common set of spectral diagnostics for these 1.3 < z < 4.2 galaxies, in particular H alpha, Paschen alpha, and the rest-frame optical and near-infrared continua. In addition, two of the four targets have JWST coverage of [O III] 5007 Angstrom and H beta; the other two targets have have JWST coverage of PAH 3.3 micron and complementary ALMA data covering the [C II] 158 micron emission line. The science goals of TEMPLATES are to demonstrate attenuation-robust diagnostics of star formation, map the distribution of star formation, compare the young and old stellar populations, and measure the physical conditions of star formation and their spatial variation across the galaxies. In addition, TEMPLATES has technical goals to establish best practices for the Integral Field Units (IFU) within the NIRSpec and MIRI instruments, both in terms of observing strategy and in terms of data reduction. The paper describes TEMPLATES's observing program, scientific and technical goals, data reduction methods, and deliverables, including high-level data products and data reduction cookbooks. △ Less

Submitted 16 December, 2023; originally announced December 2023.

Comments: 28 pages, including 8 figures and 3 tables. Submitted to The Astrophysical Journal's Focus Issue on the TEMPLATES JWST Early Release Science Program, https://iopscience.iop.org/collections/apj-231204-01-templates

The ALMA-CRISTAL survey: Widespread dust-obscured star formation in typical star-forming galaxies at z=4-6

Authors: Ikki Mitsuhashi, Ken-ichi Tadaki, Ryota Ikeda, Rodrigo Herrera-Camus, Manuel Aravena, Ilse De Looze, Natascha M. Förster Schreiber, Jorge González-López, Justin Spilker, Roberto J. Assef, Rychard Bouwens, Loreto Barcos-Munoz, Jack Birkin, Rebecca A. A. Bowler, Gabriela Calistro Rivera, Rebecca Davies, Elisabete Da Cunha, Tanio Díaz-Santos, Andrea Ferrara, Deanne Fisher, Lilian L. Lee, Juno Li, Dieter Lutz, Monica Relaño, Thorsten Naab , et al. (7 additional authors not shown)

Abstract: We present the morphological parameters and global properties of dust-obscured star formation in typical star-forming galaxies at z=4-6. Among 26 galaxies composed of 20 galaxies observed by the Cycle-8 ALMA Large Program, CRISTAL, and six galaxies from archival data, we have individually detected rest-frame 158$μ$m dust continuum emission from 19 galaxies, nine of which are reported for the first… ▽ More We present the morphological parameters and global properties of dust-obscured star formation in typical star-forming galaxies at z=4-6. Among 26 galaxies composed of 20 galaxies observed by the Cycle-8 ALMA Large Program, CRISTAL, and six galaxies from archival data, we have individually detected rest-frame 158$μ$m dust continuum emission from 19 galaxies, nine of which are reported for the first time. The derived far-infrared luminosities are in the range $\log_{10} L_{\rm IR}\,[L_{\odot}]=$10.9-12.4, an order of magnitude lower than previously detected massive dusty star-forming galaxies (DSFGs). The average relationship between the fraction of dust-obscured star formation ($f_{\rm obs}$) and the stellar mass is consistent with previous results at z=4-6 in a mass range of $\log_{10} M_{\ast}\,[M_{\odot}]\sim$9.5-11.0 and show potential evolution from z=6-9. The individual $f_{\rm obs}$ exhibits a significant diversity, and it shows a correlation with the spatial offset between the dust and the UV continuum, suggesting the inhomogeneous dust reddening may cause the source-to-source scatter in $f_{\rm obs}$. The effective radii of the dust emission are on average $\sim$1.5 kpc and are $\sim2$ times more extended than the rest-frame UV. The infrared surface densities of these galaxies ($Σ_{\rm IR}\sim2.0\times10^{10}\,L_{\odot}\,{\rm kpc}^{-2}$) are one order of magnitude lower than those of DSFGs that host compact central starbursts. On the basis of the comparable contribution of dust-obscured and dust-unobscured star formation along with their similar spatial extent, we suggest that typical star-forming galaxies at z=4-6 form stars throughout the entirety of their disks. △ Less

Submitted 29 November, 2023; originally announced November 2023.

SuperCAM CO(3-2) APEX survey at 6 pc resolution in the Small Magellanic Clouds

Authors: H. P. Saldaño, M. Rubio, A. D. Bolatto, K. Sandstrom, B. J. Swift, C. Verdugo, K. Jameson, C. K. Walker, C. Kulesa, J. Spilker, P. Bergman, G. A. Salazar

Abstract: We present the CO(3-2) APEX survey at 6 pc resolution of the bar of the SMC. We aboard the CO analysis in the SMC-Bar comparing the CO(3-2) survey with that of the CO(2-1) of similar resolution. We study the CO(3-2)-to-CO(2-1) ratio (R32) that is very sensitive to the environment properties (e.g., star-forming regions). We analyzed the correlation of this ratio with observational quantities that t… ▽ More We present the CO(3-2) APEX survey at 6 pc resolution of the bar of the SMC. We aboard the CO analysis in the SMC-Bar comparing the CO(3-2) survey with that of the CO(2-1) of similar resolution. We study the CO(3-2)-to-CO(2-1) ratio (R32) that is very sensitive to the environment properties (e.g., star-forming regions). We analyzed the correlation of this ratio with observational quantities that trace the star formation as the local CO emission, the Spitzer color [70/160], and the total IR surface brightness measured from the Spitzer and Herschel bands. For the identification of the CO(3-2) clouds, we used the CPROPS algorithm, which allowed us to measure the physical properties of the clouds. We analyzed the scaling relationships of such physical properties. We obtained an R32 of 0.65 as a median value for the SMC, with a standard deviation of 0.3. We found that R32 varies from region to region, depending on the star formation activity. In regions dominated by HII and photo-dissociated regions (e.g., N22, N66), R32 tends to be higher than the median values. Meanwhile, lower values were found toward quiescent clouds. We also found that R32 correlates positively with the IR color [70/160] and the total IR surface brightness. This finding indicates that R32 increases with environmental properties like the dust temperature, the total gas density, and the radiation field. We have identified 225 molecular clouds with sizes R > 1.5 pc and signal-to-noise (S/N) ratio > 3 and only 17 well-resolved CO(3-2) clouds increasing the S/N ratio to > 5. These 17 clouds follow consistent scaling relationships to the inner Milky Way clouds but with some departure. The CO(3-2) tends to be less turbulent and less luminous than the inner Milky Way clouds of similar size. Finally, we estimated a median virial-based CO-to-H2 conversion factor of 12.6_{-7}^{+10} Msun/(K km s^{-1} pc^{2}) for the total sample. △ Less

Submitted 31 October, 2023; originally announced November 2023.

Journal ref: A&A 687, A26 (2024)

Recovering simulated planet and disk signals using SCALES aperture masking

Authors: Mackenzie Lach, Steph Sallum, Ravinder Banyal, Natalie Batalha, Geoff Blake, Tim Brandt, Zackery Briesemeister, Aditi Desai, Josh Eisner, Wen-fai Fong, Tom Greene, Mitsuhiko Honda, Isabel Kain, Charlie Kilpatrick, Katherine de Kleer, Michael Liu, Bruce Macintosh, Raquel Martinez, Dimitri Mawet, Brittany Miles, Caroline Morley, Imke de Pater, Diana Powell, Patrick Sheehan, Andrew Skemer , et al. (7 additional authors not shown)

Abstract: The Slicer Combined with Array of Lenslets for Exoplanet Spectroscopy (SCALES) instrument is a lenslet-based integral field spectrograph that will operate at 2 to 5 microns, imaging and characterizing colder (and thus older) planets than current high-contrast instruments. Its spatial resolution for distant science targets and/or close-in disks and companions could be improved via interferometric t… ▽ More The Slicer Combined with Array of Lenslets for Exoplanet Spectroscopy (SCALES) instrument is a lenslet-based integral field spectrograph that will operate at 2 to 5 microns, imaging and characterizing colder (and thus older) planets than current high-contrast instruments. Its spatial resolution for distant science targets and/or close-in disks and companions could be improved via interferometric techniques such as sparse aperture masking. We introduce a nascent Python package, NRM-artist, that we use to design several SCALES masks to be non-redundant and to have uniform coverage in Fourier space. We generate high-fidelity mock SCALES data using the scalessim package for SCALES' low spectral resolution modes across its 2 to 5 micron bandpass. We include realistic noise from astrophysical and instrument sources, including Keck adaptive optics and Poisson noise. We inject planet and disk signals into the mock datasets and subsequently recover them to test the performance of SCALES sparse aperture masking and to determine the sensitivity of various mask designs to different science signals. △ Less

Submitted 19 October, 2023; originally announced October 2023.

Journal ref: SPIE Proceedings Volume 12680, Techniques and Instrumentation for Detection of Exoplanets XI; 1268024 (2023)

Simulating medium-spectral-resolution exoplanet characterization with SCALES angular/reference differential imaging

Authors: Aditi Desai, Stephanie E. Sallum, Ravinder Banyal, Natalie Batalha, Natasha Batalha, Geoff Blake, Tim Brandt, Zack Briesemeister, Katherine de Kleer, Imke de Pater, Josh Eisner, Wen-fai Fong, Tom Greene, Mitsuhiko Honda, Isabel Kain, Charlie Kilpatrick, Mackenzie Lach, Mike Liu, Bruce Macintosh, Raquel A. Martinez, Dimitri Mawet, Brittany Miles, Caroline Morley, Diana Powell, Patrick Sheehan , et al. (8 additional authors not shown)

Abstract: SCALES (Slicer Combined with Array of Lenslets for Exoplanet Spectroscopy) is a 2 - 5 micron high-contrast lenslet-based integral field spectrograph (IFS) designed to characterize exoplanets and their atmospheres. The SCALES medium-spectral-resolution mode uses a lenslet subarray with a 0.34 x 0.36 arcsecond field of view which allows for exoplanet characterization at increased spectral resolution… ▽ More SCALES (Slicer Combined with Array of Lenslets for Exoplanet Spectroscopy) is a 2 - 5 micron high-contrast lenslet-based integral field spectrograph (IFS) designed to characterize exoplanets and their atmospheres. The SCALES medium-spectral-resolution mode uses a lenslet subarray with a 0.34 x 0.36 arcsecond field of view which allows for exoplanet characterization at increased spectral resolution. We explore the sensitivity limitations of this mode by simulating planet detections in the presence of realistic noise sources. We use the SCALES simulator scalessim to generate high-fidelity mock observations of planets that include speckle noise from their host stars, as well as other atmospheric and instrumental noise effects. We employ both angular and reference differential imaging as methods of disentangling speckle noise from the injected planet signals. These simulations allow us to assess the feasibility of speckle deconvolution for SCALES medium resolution data, and to test whether one approach outperforms another based on planet angular separations and contrasts. △ Less

Submitted 18 October, 2023; originally announced October 2023.

The Slicer Combined with Array of Lenslets for Exoplanet Spectroscopy (SCALES): driving science cases and expected outcomes

Authors: Steph Sallum, Andrew Skemer, Deno Stelter, Ravinder Banyal, Natalie Batalha, Natasha Batalha, Geoff Blake, Tim Brandt, Zack Briesemeister, Katherine de Kleer, Imke de Pater, Aditi Desai, Josh Eisner, Wen-fai Fong, Tom Greene, Mitsuhiko Honda, Rebecca Jensen-Clem, Isabel Kain, Charlie Kilpatrick, Renate Kupke, Mackenzie Lach, Michael C. Liu, Bruce Macintosh, Raquel A. Martinez, Dimitri Mawet , et al. (12 additional authors not shown)

Abstract: The Slicer Combined with Array of Lenslets for Exoplanet Spectroscopy (SCALES) is a $2-5~μ$m, high-contrast integral field spectrograph (IFS) currently being built for Keck Observatory. With both low ($R\lesssim250$) and medium ($R\sim3500-7000$) spectral resolution IFS modes, SCALES will detect and characterize significantly colder exoplanets than those accessible with near-infrared ($\sim1-2~μ$m… ▽ More The Slicer Combined with Array of Lenslets for Exoplanet Spectroscopy (SCALES) is a $2-5~μ$m, high-contrast integral field spectrograph (IFS) currently being built for Keck Observatory. With both low ($R\lesssim250$) and medium ($R\sim3500-7000$) spectral resolution IFS modes, SCALES will detect and characterize significantly colder exoplanets than those accessible with near-infrared ($\sim1-2~μ$m) high-contrast spectrographs. This will lead to new progress in exoplanet atmospheric studies, including detailed characterization of benchmark systems that will advance the state of the art of atmospheric modeling. SCALES' unique modes, while designed specifically for direct exoplanet characterization, will enable a broader range of novel (exo)planetary observations as well as galactic and extragalactic studies. Here we present the science cases that drive the design of SCALES. We describe an end-to-end instrument simulator that we use to track requirements, and show simulations of expected science yields for each driving science case. We conclude with a discussion of preparations for early science when the instrument sees first light in $\sim2025$. △ Less

Submitted 10 October, 2023; originally announced October 2023.

Comments: 10 pages, 16 figures, submitted to Proceedings of the SPIE

JWST's TEMPLATES for Star Formation: The First Resolved Gas-Phase Metallicity Maps of Dust-Obscured Star-Forming Galaxies at $z$ $\sim$ 4

Authors: Jack E. Birkin, Taylor A. Hutchison, Brian Welch, Justin S. Spilker, Manuel Aravena, Matthew B. Bayliss, Jared Cathey, Scott C. Chapman, Anthony H. Gonzalez, Gayathri Gururajan, Christopher C. Hayward, Gourav Khullar, Keunho J. Kim, Guillaume Mahler, Matthew A. Malkan, Desika Narayanan, Grace M. Olivier, Kedar A. Phadke, Cassie Reuter, Jane R. Rigby, J. D. T. Smith, Manuel Solimano, Nikolaus Sulzenauer, Joaquin D. Vieira, David Vizgan , et al. (1 additional authors not shown)

Abstract: We present the first spatially resolved maps of gas-phase metallicity for two dust-obscured star-forming galaxies (DSFGs) at $z\sim$ 4, from the JWST TEMPLATES Early Release Science program, derived from NIRSpec integral field unit spectroscopy of the H$α$ and [NII] emission lines. Empirical optical line calibrations are used to determine that the sources are globally enriched to near-solar levels… ▽ More We present the first spatially resolved maps of gas-phase metallicity for two dust-obscured star-forming galaxies (DSFGs) at $z\sim$ 4, from the JWST TEMPLATES Early Release Science program, derived from NIRSpec integral field unit spectroscopy of the H$α$ and [NII] emission lines. Empirical optical line calibrations are used to determine that the sources are globally enriched to near-solar levels. While one source shows elevated [NII]/H$α$ ratios and broad H$α$ emission consistent with the presence of an AGN in a $\gtrsim$1kpc region, we argue that both systems have already undergone significant metal enrichment as a result of their extremely high star formation rates. Utilizing ALMA rest-frame 380$μ$m continuum and [CI]($^3$P$_2$-$^3$P$_1$) line maps we compare the spatial variation of the metallicity and gas-to-dust ratio in the two galaxies, finding the two properties to be anticorrelated on highly resolved spatial scales, consistent with various literature studies of $z\sim$ 0 galaxies. The data are indicative of the enormous potential of JWST to probe the enrichment of the interstellar medium on $\sim$kpc scales in extremely dust-obscured systems at $z\sim$ 4 and beyond. △ Less

Submitted 28 August, 2023; v1 submitted 19 July, 2023; originally announced July 2023.

Comments: 14 pages, 6 figures, accepted for publication in ApJ

TEMPLATES: Characterization of a Merger in the Dusty Lensing SPT0418-47 System

Authors: Jared Cathey, Anthony H. Gonzalez, Sidney Lower, Kedar A. Phadke, Justin Spilker, Manuel Aravena, Jack E. Birkin, Simon Birrer, Scott Chapman, Håkon Dahle, Cristopher C. Hayward, Yashar Hezaveh, Ryley Hill, Taylor A. Hutchison, Guillaume Mahler, Daniel P. Marrone, Desika Narayanan, Alexander Navarre, Cassie Reuter, Jane R. Rigby, Keren Sharon, Manuel Solimano, Nikolaus Sulzenauer, Joaquin Vieira, David Vizgan

Abstract: We present JWST and ALMA results for the lensing system SPT0418-47, which includes a strongly-lensed, dusty star-forming galaxy at redshift z=4.225 and an associated multiply-imaged companion. JWST NIRCam and MIRI imaging observations presented in this paper were acquired as part of the Early Release Science program Targeting Extremely Magnified Panchromatic Lensed Arcs and Their Extended Star For… ▽ More We present JWST and ALMA results for the lensing system SPT0418-47, which includes a strongly-lensed, dusty star-forming galaxy at redshift z=4.225 and an associated multiply-imaged companion. JWST NIRCam and MIRI imaging observations presented in this paper were acquired as part of the Early Release Science program Targeting Extremely Magnified Panchromatic Lensed Arcs and Their Extended Star Formation (TEMPLATES). This data set provides robust, mutiwavelength detection of stellar light in both the main (SPT0418A) and companion (SPT0418B) galaxies, while the ALMA detection of [C II] emission confirms that SPT0418B lies at the same redshift as SPT0418A. From a source plane reconstruction, we infer that the projected physical separation of the two galaxies is $4.42\pm 0.05$ kpc. We derive total magnifications of $μ=29.5\pm1.2$ and $μ=4.2\pm 0.9$ for SPT0418A and SPT0418B, respectively. We use both CIGALE and PROSPECTOR to derive stellar masses. The stellar mass ratio of SPT0418A and SPT0418B is approximately 4 to 1 ($4.5\pm 1.0$ for CIGALE and $4.2^{+1.9}_{-1.6}$ for PROSPECTOR). We also see evidence of extended structure associated with SPT0418A in the lensing reconstruction that is suggestive of a tidal feature. Interestingly, the star formation rates and stellar masses of both galaxies are consistent with the main sequence of star-forming galaxies at this epoch, indicating that this ongoing interaction has not noticeably elevated the star formation levels. △ Less

Submitted 19 July, 2023; originally announced July 2023.

Comments: Submitted to AAS Journals,14 pages, 7 figures

Observations of neutral carbon in 29 high-z lensed dusty star forming galaxies and the comparison of gas mass tracers

Authors: G. Gururajan, M. Béthermin, N. Sulzenauer, P. Theulé, J. S. Spilker, M. Aravena, S. C. Chapman, A. Gonzalez, T. R. Greve, D. Narayanan, C. Reuter, J. D. Vieira, A. Weiss

Abstract: The nature and evolution of high-redshift dusty star-forming galaxies (high-z DSFGs) remain an open question. Their massive gas reservoirs play an important role in driving the intense star-formation rates hosted in these galaxies. We aim to estimate the molecular gas content of high-z DSFGs by using various gas mass tracers such as the [CI], CO, [CII] emission lines and the dust content. These tr… ▽ More The nature and evolution of high-redshift dusty star-forming galaxies (high-z DSFGs) remain an open question. Their massive gas reservoirs play an important role in driving the intense star-formation rates hosted in these galaxies. We aim to estimate the molecular gas content of high-z DSFGs by using various gas mass tracers such as the [CI], CO, [CII] emission lines and the dust content. These tracers need to be well calibrated as they are all limited by uncertainties on factors such as aCO, XCI, aCII and GDR, thereby affecting the determination of the gas mass accurately. The main goal of our work is to check the consistency between the gas mass tracers and cross-calibrate the uncertain factors. We observe the two [CI] line transitions for 29 SPT-SMGs with the ALMA-ACA. Additionally, we also present new APEX observations of [CII] line for 9 of these galaxies. We find a nearly linear relation between the infrared luminosity and [CI] luminosity if we fit the starbursts and main-sequence galaxies separately. We measure a median [CI]-derived excitation temperature of 34.5+/-2.1 K. We probe the properties of the interstellar medium (ISM) such as density and radiation field intensity using [CI] to mid- or high-J CO lines and [CI] to infrared luminosity ratio, and find similar values to the SMG populations in literature. Finally, the gas masses estimated from [CI], CO, dust, and [CII] do not exhibit any significant trend with the infrared luminosity or the dust temperature. We provide the various cross-calibrations between these tracers. Our study confirms that [CI] is a suitable tracer of the molecular gas content, and shows an overall agreement between all the classical gas tracers used at high redshift. However, their absolute calibration and thus the gas depletion timescale measurements remain uncertain. △ Less

Submitted 5 June, 2023; originally announced June 2023.

Comments: Accepted for publication in A&A, 25 pages, 11 figures, 6 tables

Journal ref: A&A 676, A89 (2023)

Spatial variations in aromatic hydrocarbon emission in a dust-rich galaxy

Authors: Justin S. Spilker, Kedar A. Phadke, Manuel Aravena, Melanie Archipley, Matthew B. Bayliss, Jack E. Birkin, Matthieu Bethermin, James Burgoyne, Jared Cathey, Scott C. Chapman, Hakon Dahle, Anthony H. Gonzalez, Gayathri Gururajan, Christopher C. Hayward, Yashar D. Hezaveh, Ryley Hill, Taylor A. Hutchison, Keunho J. Kim, Seonwoo Kim, David Law, Ronan Legin, Matthew A. Malkan, Daniel P. Marrone, Eric J. Murphy, Desika Narayanan , et al. (13 additional authors not shown)

Abstract: Dust grains absorb half of the radiation emitted by stars throughout the history of the universe, re-emitting this energy at infrared wavelengths. Polycyclic aromatic hydrocarbons (PAHs) are large organic molecules that trace millimeter-size dust grains and regulate the cooling of the interstellar gas within galaxies. Observations of PAH features in very distant galaxies have been difficult due to… ▽ More Dust grains absorb half of the radiation emitted by stars throughout the history of the universe, re-emitting this energy at infrared wavelengths. Polycyclic aromatic hydrocarbons (PAHs) are large organic molecules that trace millimeter-size dust grains and regulate the cooling of the interstellar gas within galaxies. Observations of PAH features in very distant galaxies have been difficult due to the limited sensitivity and wavelength coverage of previous infrared telescopes. Here we present JWST observations that detect the 3.3um PAH feature in a galaxy observed less than 1.5 billion years after the Big Bang. The high equivalent width of the PAH feature indicates that star formation, rather than black hole accretion, dominates the infrared emission throughout the galaxy. The light from PAH molecules, large dust grains, and stars and hot dust are spatially distinct from one another, leading to order-of-magnitude variations in the PAH equivalent width and the ratio of PAH to total infrared luminosity across the galaxy. The spatial variations we observe suggest either a physical offset between the PAHs and large dust grains or wide variations in the local ultraviolet radiation field. Our observations demonstrate that differences in the emission from PAH molecules and large dust grains are a complex result of localized processes within early galaxies. △ Less

Submitted 5 June, 2023; originally announced June 2023.

Comments: Published in Nature 5 June 2023 at https://www.nature.com/articles/s41586-023-05998-6. MIRI MRS reduction notebook is available at https://github.com/jwst-templates

ALMA 1.1mm Observations of a Conservative Sample of High Redshift Massive Quiescent Galaxies in SHELA

Authors: Katherine Chworowsky, Steven L. Finkelstein, Justin S. Spilker, Gene C. K. Leung, Micaela B. Bagley, Caitlin M. Casey, Caryl Gronwall, Shardha Jogee, Rebecca L. Larson, Casey Papovich, Rachel S. Somerville, Matthew Stevans, Isak G. B. Wold, L. Y. Aaron Yung

Abstract: We present a sample of 30 massive (log$(M_{\ast}/M_\odot) >11$) $z=3-5$ quiescent galaxies selected from the \textit{Spitzer-}HETDEX Exploratory Large Area (SHELA) Survey and observed at 1.1mm with Atacama Large Millimeter/submillimeter Array (ALMA) Band 6 observations. These ALMA observations would detect even modest levels of dust-obscured star-formation, on order of… ▽ More We present a sample of 30 massive (log$(M_{\ast}/M_\odot) >11$) $z=3-5$ quiescent galaxies selected from the \textit{Spitzer-}HETDEX Exploratory Large Area (SHELA) Survey and observed at 1.1mm with Atacama Large Millimeter/submillimeter Array (ALMA) Band 6 observations. These ALMA observations would detect even modest levels of dust-obscured star-formation, on order of $\sim 20 \ M_\odot \textrm{yr}^{-1}$ at $z\sim4$ at a $1σ$ level, allowing us to quantify the amount of contamination from dusty star-forming sources in our quiescent sample. Starting with a parent sample of candidate massive quiescent galaxies from the Stevans et al. 2021 v1 SHELA catalog, we use the Bayesian \textsc{Bagpipes} spectral energy distribution fitting code to derive robust stellar masses ($M_*$) and star-formation rates (SFRs) for these sources, and select a conservative sample of 36 candidate massive ($M_* > 10^{11}M_\odot$) quiescent galaxies, with specific SFRs at $>2σ$ below the star-forming main sequence at $z\sim4$. Based on ALMA imaging, six of these candidate quiescent galaxies have the presence of significant dust-obscured star-formation, thus were removed from our final sample. This implies a $\sim 17\%$ contamination rate from dusty star-forming galaxies with our selection criteria using the v1 SHELA catalog. This conservatively-selected quiescent galaxy sample at $z=3-5$ will provide excellent targets for future observations to better constrain how massive galaxies can both grow and shut-down their star-formation in a relatively short time period. △ Less

Submitted 10 May, 2023; originally announced May 2023.

Comments: 16 pages, 13 figures

Accurate Oxygen Abundance of Interstellar Gas in Mrk 71 from Optical and Infrared Spectra

Authors: Yuguang Chen, Tucker Jones, Ryan Sanders, Dario Fadda, Jessica Sutter, Robert Minchin, Erin Huntzinger, Peter Senchyna, Daniel Stark, Justin Spilker, Benjamin Weiner, Guido Roberts-Borsani

Abstract: The heavy element content ("metallicity") of the Universe is a record of the total star formation history. Gas-phase metallicity in galaxies, as well as its evolution with time, is of particular interest as a tracer of accretion and outflow processes. However, metallicities from the widely-used electron temperature ($T_e$) method are typically ~2x lower than the values based on the recombination l… ▽ More The heavy element content ("metallicity") of the Universe is a record of the total star formation history. Gas-phase metallicity in galaxies, as well as its evolution with time, is of particular interest as a tracer of accretion and outflow processes. However, metallicities from the widely-used electron temperature ($T_e$) method are typically ~2x lower than the values based on the recombination line method. This "abundance discrepancy factor" (ADF) is well known and is commonly ascribed to bias due to temperature fluctuations. We present a measurement of oxygen abundance in the nearby (3.4 Mpc) system, Mrk 71, using a combination of optical and far-IR emission lines to measure and correct for temperature fluctuation effects. Our far-IR result is inconsistent ($> 2 σ$ significance) with the metallicity from recombination lines and instead indicates little to no bias in the standard $T_e$ method, ruling out the long-standing hypothesis that the ADF is explained by temperature fluctuations for this object. Our results provide a framework to accurately measure metallicity across cosmic history, including with recent data reaching within the first billion years with JWST and the Atacama Large Millimeter Array (ALMA). △ Less

Submitted 19 April, 2023; originally announced April 2023.

Comments: 26 pages, 7 figures. Accepted for publication

Brightest Cluster Galaxy Formation in the z=4.3 Protocluster SPT2349-56: Discovery of a Radio-Loud AGN

Authors: Scott C. Chapman, Ryley Hill, Manuel Aravena, Melanie Archipley, Arif Babul, James Burgoyne, Rebecca E. A. Canning, Carlos De Breuck, Anthony H. Gonzalez, Christopher C. Hayward, Seon Woo Kim, Matt Malkan, Dan P. Marrone, Vincent McIntyre, Eric Murphy, Emily Pass, Ryan W. Perry, Kedar A. Phadke, Douglas Rennehan, Cassie Reuter, Kaja M. Rotermund, Douglas Scott, Nick Seymour, Manuel Solimano, Justin Spilker , et al. (7 additional authors not shown)

Abstract: We have observed the z=4.3 protocluster SPT2349-56 with ATCA with the aim of detecting radio-loud active galactic nuclei (AGN) amongst the ~30 submillimeter galaxies identified in the structure. We detect the central complex of SMGs at 2.2\,GHz with a luminosity of L_2.2=(4.42pm0.56)x10^{25} W/Hz. The ASKAP also detects the source at 888 MHz, constraining the radio spectral index to alpha=-1.6pm0.… ▽ More We have observed the z=4.3 protocluster SPT2349-56 with ATCA with the aim of detecting radio-loud active galactic nuclei (AGN) amongst the ~30 submillimeter galaxies identified in the structure. We detect the central complex of SMGs at 2.2\,GHz with a luminosity of L_2.2=(4.42pm0.56)x10^{25} W/Hz. The ASKAP also detects the source at 888 MHz, constraining the radio spectral index to alpha=-1.6pm0.3, consistent with ATCA non-detections at 5.5 and 9GHz, and implying L_1.4(rest)=(2.4pm0.3)x10^{26}W/Hz. This radio luminosity is about 100 times higher than expected from star formation, assuming the usual FIR-radio correlation, which is a clear indication of an AGN driven by a forming brightest cluster galaxy (BCG). None of the SMGs in SPT2349-56 show signs of AGN in any other diagnostics available to us (notably 12CO out to J=16, OH163um, CII/IR, and optical spectra), highlighting the radio continuum as a powerful probe of obscured AGN in high-z protoclusters. No other significant radio detections are found amongst the cluster members, consistent with the FIR-radio correlation. We compare these results to field samples of radio sources and SMGs, along with the 22 SPT-SMG gravitational lenses also observed in the ATCA program, as well as powerful radio galaxies at high redshifts. Our results allow us to better understand the effects of this gas-rich, overdense environment on early supermassive black hole (SMBH) growth and cluster feedback. We estimate that (3.3pm0.7)x10^{38} W of power are injected into the growing ICM by the radio-loud AGN, whose energy over 100Myr is comparable to the binding energy of the gas mass of the central halo. The AGN power is also comparable to the instantaneous energy injection from supernova feedback from the 23 catalogued SMGs in the core region of 120kpc projected radius. The SPT2349-56 radio-loud AGN may be providing strong feedback on a nascent ICM. △ Less

Submitted 4 January, 2023; v1 submitted 3 January, 2023; originally announced January 2023.

Comments: 31 pages, submitted to ApJ, Dec17,2022

Extended Lyman-$α$ emission towards the SPT2349-56 protocluster at $z=4.3$

Authors: Yordanka Apostolovski, Manuel Aravena, Timo Anguita, Matthieu Bethermin, James Burgoyne, Scott Chapman, Carlos De Breuck, Anthony Gonzalez, Max Gronke, Lucia Guaita, Yashar Hezaveh, Ryley Hill, Sreevani Jarugula, Evelyn Johnston, Matt Malkan, Desika Narayanan, Cassie Reuter, Manuel Solimano, Justin Spilker, Nikolaus Sulzenauer, Joaquin Vieira, David Vizgan, Axel Weiß

Abstract: Context. Deep spectroscopic surveys with the Atacama Large Millimeter/submillimeter Array (ALMA) have revealed that some of the brightest infrared sources in the sky correspond to concentrations of dusty star-forming galaxies (DSFG) at high redshift. Among these, the SPT2349-56 protocluster system at z = 4.304 is amongst the most extreme examples due to its high source density and integrated star… ▽ More Context. Deep spectroscopic surveys with the Atacama Large Millimeter/submillimeter Array (ALMA) have revealed that some of the brightest infrared sources in the sky correspond to concentrations of dusty star-forming galaxies (DSFG) at high redshift. Among these, the SPT2349-56 protocluster system at z = 4.304 is amongst the most extreme examples due to its high source density and integrated star formation rate. Aims. We conducted a deep Lyman-$α$ line emission survey around SPT2349-56 using the Multi-Unit Spectroscopic Explorer (MUSE) at Very Large Telescope (VLT) in order to characterize this uniquely dense environment. Methods. Taking advantage of the deep three-dimensional nature of this survey, we performed a sensitive search for Lyman-$α$ emitters (LAEs) toward the core and northern extension of the protocluster, which correspond to the brightest infrared regions in this field. Using a smoothed narrowband image extracted from the MUSE datacube around the protocluster redshift, we searched for possible extended structures. Results. We identify only three LAEs at z = 4.3 in this field, in concordance with expectations for blank-fields, and an extended Lyman-$α$ structure spatially associated with core of the protocluster. All the previously-identified DSFGs in this field are undetected in Lyman-$α$ emission, consistent with the conspicuous dust obscuration in these systems. We find an extended Lyman-$α$ structure, about $60 \times 60$ kpc$^{2}$ in size, and located 56 kpc west of the protocluster core. Three DSFGs coincide spatially with the location of this structure. We conclude that either the three co-spatial DSFGs or the protocluster core itself are feeding ionizing photons to the Lyman-$α$ structure. △ Less

Submitted 3 January, 2023; originally announced January 2023.

Comments: 14 pages, 9 figures, 2 tables, submitted to A&A

Merger Signatures are Common, but not Universal, In Massive, Recently-Quenched Galaxies at z~0.7

Authors: Margaret Verrico, David J. Setton, Rachel Bezanson, Jenny E. Greene, Katherine A. Suess, Andy D. Goulding, Justin S. Spilker, Mariska Kriek, Robert Feldmann, Desika Narayanan, Vincenzo Donofrio, Gourav Khullar

Abstract: We present visual classifications of merger-induced tidal disturbances in 143 $\rm{M}_* \sim 10^{11}\rm{M}_\odot$ post-starburst galaxies at z$\sim$0.7 identified in the SQuIGG$\vec{L}$E Sample. This sample spectroscopically selects galaxies from the Sloan Digital Sky Survey that have stopped their primary epoch of star formation within the past $\sim$500 Myrs. Visual classifications are performed… ▽ More We present visual classifications of merger-induced tidal disturbances in 143 $\rm{M}_* \sim 10^{11}\rm{M}_\odot$ post-starburst galaxies at z$\sim$0.7 identified in the SQuIGG$\vec{L}$E Sample. This sample spectroscopically selects galaxies from the Sloan Digital Sky Survey that have stopped their primary epoch of star formation within the past $\sim$500 Myrs. Visual classifications are performed on Hyper Suprime Cam (HSC) i-band imaging. We compare to a control sample of mass- and redshift-matched star-forming and quiescent galaxies from the Large Early Galaxy Census and find that post-starburst galaxies are more likely to be classified as disturbed than either category. This corresponds to a factor of $3.6^{+2.9}_{-1.3}$ times the disturbance rate of older quiescent galaxies and $2.1^{+1.9}_{-.73}$ times the disturbance rate of star-forming galaxies. Assuming tidal features persist for $\lesssim500$ Myr, this suggests merging is coincident with quenching in a significant fraction of these post-starbursts. Galaxies with tidal disturbances are younger on average than undisturbed post-starburst galaxies in our sample, suggesting tidal features from a major merger may have faded over time. This may be exacerbated by the fact that, on average, the undisturbed subset is fainter, rendering low surface brightness tidal features harder to identify. However, the presence of ten young ($\lesssim150$ Myr since quenching) undisturbed galaxies suggests that major mergers are not the only fast physical mechanism that shut down the primary epoch of star formation in massive galaxies at intermediate redshift. △ Less

Submitted 17 March, 2023; v1 submitted 29 November, 2022; originally announced November 2022.

Comments: 13 pages, 8 figures, accepted to ApJ

The Rest-Frame Submillimeter Spectrum of High Redshift, Dusty, Star-Forming Galaxies from the SPT-SZ Survey

Authors: C. Reuter, J. S. Spilker, J. D. Vieira, D. P. Marrone, A. Weiss, M. Aravena, M. A. Archipley, S. C. Chapman, A. Gonzalez, T. R. Greve, C. C. Hayward, R. Hill, S. Jarugula, S. Kim, M. Malkan, K. A. Phadke, A. A. Stark, N. Sulzenauer, D. Vizgan

Abstract: We present the average rest-frame spectrum of the final catalog of dusty star-forming galaxies (DSFGs) selected from the South Pole Telescope SZ survey (SPT-SZ) and measured with Band 3 of the Atacama Large Millimeter/submillimeter Array (ALMA). This work builds on the previous average rest-frame spectrum, given in Spilker et al. (2014) for the first 22 sources, and is comprised of a total of 78 s… ▽ More We present the average rest-frame spectrum of the final catalog of dusty star-forming galaxies (DSFGs) selected from the South Pole Telescope SZ survey (SPT-SZ) and measured with Band 3 of the Atacama Large Millimeter/submillimeter Array (ALMA). This work builds on the previous average rest-frame spectrum, given in Spilker et al. (2014) for the first 22 sources, and is comprised of a total of 78 sources, normalized by their respective apparent dust masses. The spectrum spans $1.9$$<$z$<$$6.9$ and covers rest-frame frequencies of 240$-$800 GHz. Combining this data with low-J CO observations from the Australia Telescope Compact Array (ATCA), we detect multiple bright line features from $^{12}$CO, $[$CI$]$, and H$_2$O, as well as fainter molecular transitions from $^{13}$CO, HCN, HCO$^+$, HNC, CN, H$_2$O$^+$, and CH. We use these detections, along with limits from other molecules, to characterize the typical properties of the interstellar medium (ISM) for these high redshift DSFGs. We are able to divide the large sample into subsets in order to explore how the average spectrum changes with various galaxy properties, such as effective dust temperature. We find that systems with hotter dust temperatures exhibit differences in the bright $^{12}$CO emission lines, and contain either warmer and more excited dense gas tracers, or larger dense gas reservoirs. These observations will serve as a reference point to studies of the ISM in distant luminous DSFGs (L$_{\mathrm{IR}}$$>$$10^{12}$L$_\odot$), and will inform studies of chemical evolution before the peak epoch of star formation at $z=2-3$. △ Less

Submitted 3 January, 2023; v1 submitted 20 October, 2022; originally announced October 2022.

Comments: 26 pages, 10 figures, accepted for publication in ApJ

Does the lockstep growth between black holes and bulges create their mass relation?

Authors: Guang Yang, W. N. Brandt, David M. Alexander, Médéric Boquien, Qingling Ni, Casey Papovich, Justin S. Spilker, Fabio Vito, Jonelle L. Walsh, Chengpeng Zhang

Abstract: Recent studies have revealed a strong relation between sample-averaged black-hole (BH) accretion rate (BHAR) and star formation rate (SFR) among bulge-dominated galaxies, i.e., "lockstep" BH-bulge growth, in the distant universe. This relation might be closely related to the BH-bulge mass correlation observed in the local universe. To understand further BH-bulge coevolution, we present ALMA CO(2-1… ▽ More Recent studies have revealed a strong relation between sample-averaged black-hole (BH) accretion rate (BHAR) and star formation rate (SFR) among bulge-dominated galaxies, i.e., "lockstep" BH-bulge growth, in the distant universe. This relation might be closely related to the BH-bulge mass correlation observed in the local universe. To understand further BH-bulge coevolution, we present ALMA CO(2-1) or CO(3-2) observations of 7 star-forming bulge-dominated galaxies at z=0.5-2.5. Using the ALMA data, we detect significant ($>3σ$) CO emission from 4 objects. For our sample of 7 galaxies, we measure (or constrain with upper limits) their CO line fluxes and estimate molecular gas masses ($M_{gas}$). We also estimate their stellar masses ($M_{star}$) and SFRs by modelling their spectral energy distributions (SEDs). Using these physical properties, we derive the gas-depletion timescales ($t_{dep} = M_{gas}/SFR$) and compare them with the bulge/BH growth timescales ($t_{grow} = M_{star}/SFR \sim M_{BH}/BHAR$). Our sample generally has $t_{dep}$ shorter than $t_{grow}$ by a median factor of $\gtrsim 4$, indicating that the cold gas will be depleted before significant bulge/BH growth takes place. This result suggests that the BH-bulge lockstep growth is mainly responsible for maintaining their mass relation, not creating it. We note that our sample is small and limited to $z<2.5$; JWST and ALMA will be able to probe to higher redshifts in the near future. △ Less

Submitted 18 October, 2022; originally announced October 2022.

Comments: 15 pages, 7 figures, and 3 tables. Accepted by ApJ

Molecular Gas Reservoirs in Massive Quiescent Galaxies at $\mathrm{z\sim0.7}$ Linked to Late Time Star Formation

Authors: Charity Woodrum, Christina C. Williams, Marcia Rieke, Joel Leja, Benjamin D. Johnson, Rachel Bezanson, Robert Kennicutt, Justin Spilker, Sandro Tacchella

Abstract: We explore how the presence of detectable molecular gas depends on the inferred star formation histories (SFHs) in 8 massive, quiescent galaxies at $\mathrm{z\sim0.7}$. Half of the sample have clear detections of molecular gas, traced by CO(2-1). We find that the molecular gas content is unrelated to the rate of star formation decline prior to the most recent 1 Gyr, suggesting that the gas reservo… ▽ More We explore how the presence of detectable molecular gas depends on the inferred star formation histories (SFHs) in 8 massive, quiescent galaxies at $\mathrm{z\sim0.7}$. Half of the sample have clear detections of molecular gas, traced by CO(2-1). We find that the molecular gas content is unrelated to the rate of star formation decline prior to the most recent 1 Gyr, suggesting that the gas reservoirs are not leftover from their primary star formation epoch. However, the recent SFHs of CO-detected galaxies demonstrate evidence for secondary bursts of star formation in their last Gyr. The fraction of stellar mass formed in these secondary bursts ranges from $\mathrm{f_{burst}\approx0.3-6\%}$, and ended between $\mathrm{t_{end\mbox{-}burst}\approx0-330~Myr}$ ago. The CO-detected galaxies form a higher fraction of mass in the last Gyr ($\mathrm{f_{M_{1Gyr}}=2.6\pm1.8\%}$) compared to the CO-undetected galaxies ($\mathrm{f_{M_{1Gyr}}=0.2\pm0.1\%}$). The galaxies with gas reservoirs have enhanced late-time star formation, highlighting this as a contributing factor to the observed heterogeneity in the gas reservoirs in high-redshift quiescent galaxies. We find that the amount of gas and star formation driven by these secondary bursts are inconsistent with that expected from dry minor mergers, and instead are likely driven by recently-accreted gas i.e., gas-rich minor mergers. This conclusion would not have been made based on $\mathrm{SFR_{UV+IR}}$ measurements alone, highlighting the power of detailed SFH modeling in the interpretation of gas reservoirs. Larger samples are needed to understand the frequency of low-level rejuvenation among quiescent galaxies at intermediate redshifts, and to what extent this drives the diversity of molecular gas reservoirs. △ Less

Submitted 4 January, 2023; v1 submitted 7 October, 2022; originally announced October 2022.

Journal ref: ApJ, 940, 29 (2022)

Star Formation Suppresion by Tidal Removal of Cold Molecular Gas from an Intermediate-Redshift Massive Post-Starburst Galaxy

Authors: Justin S. Spilker, Katherine A. Suess, David J. Setton, Rachel Bezanson, Robert Feldmann, Jenny E. Greene, Mariska Kriek, Sidney Lower, Desika Narayanan, Margaret Verrico

Abstract: Observations and simulations have demonstrated that star formation in galaxies must be actively suppressed to prevent the formation of over-massive galaxies. Galactic outflows driven by stellar feedback or supermassive black hole accretion are often invoked to regulate the amount of cold molecular gas available for future star formation, but may not be the only relevant quenching processes in all… ▽ More Observations and simulations have demonstrated that star formation in galaxies must be actively suppressed to prevent the formation of over-massive galaxies. Galactic outflows driven by stellar feedback or supermassive black hole accretion are often invoked to regulate the amount of cold molecular gas available for future star formation, but may not be the only relevant quenching processes in all galaxies. We present the discovery of vast molecular tidal features extending up to 64 kpc outside of a massive z=0.646 post-starburst galaxy that recently concluded its primary star-forming episode. The tidal tails contain (1.2 +/- 0.1)x10^10 Msun of molecular gas, 47 +/- 5 % of the total cold gas reservoir of the system. Both the scale and magnitude of the molecular tidal features are unprecedented compared to all known nearby or high-redshift merging systems. We infer that the cold gas was stripped from the host galaxies during the merger, which is most likely responsible for triggering the initial burst phase and the subsequent suppression of star formation. While only a single example, this result shows that galaxy mergers can regulate the cold gas contents in distant galaxies by directly removing a large fraction of the molecular gas fuel, and plausibly suppress star formation directly, a qualitatively different physical mechanism than feedback-driven outflows. △ Less

Submitted 29 August, 2022; originally announced August 2022.

Comments: Published in ApJL, 8 pages, 5 figures

Recovering the star formation histories of recently-quenched galaxies: the impact of model and prior choices

Authors: Katherine A. Suess, Joel Leja, Benjamin D. Johnson, Rachel Bezanson, Jenny E. Greene, Mariska Kriek, Sidney Lower, Desika Narayanan, David J. Setton, Justin S. Spilker

Abstract: Accurate models of the star formation histories (SFHs) of recently-quenched galaxies can provide constraints on when and how galaxies shut down their star formation. The recent development of "non-parametric" SFH models promises the flexibility required to make these measurements. However, model and prior choices significantly affect derived SFHs, particularly for post-starburst galaxies (PSBs) wh… ▽ More Accurate models of the star formation histories (SFHs) of recently-quenched galaxies can provide constraints on when and how galaxies shut down their star formation. The recent development of "non-parametric" SFH models promises the flexibility required to make these measurements. However, model and prior choices significantly affect derived SFHs, particularly for post-starburst galaxies (PSBs) which have sharp changes in their recent SFH. In this paper, we create mock PSBs, then use the Prospector SED fitting software to test how well four different SFH models recover key properties. We find that a two-component parametric model performs well for our simple mock galaxies, but is sensitive to model mismatches. The fixed- and flexible-bin non-parametric models included in Prospector are able to rapidly quench a major burst of star formation, but systematically underestimate the post-burst age by up to 200 Myr. We develop a custom SFH model that allows for additional flexibility in the recent SFH. Our flexible non-parametric model is able to constrain post-burst ages with no significant offset and just ~90 Myr of scatter. Our results suggest that while standard non-parametric models are able to recover first-order quantities of the SFH (mass, SFR, average age), accurately recovering higher-order quantities (burst fraction, quenching time) requires careful consideration of model flexibility. These mock recovery tests are a critical part of future SFH studies. Finally, we show that our new, public SFH model is able to accurately recover the properties of mock star-forming and quiescent galaxies and is suitable for broader use in the SED fitting community. △ Less

Submitted 6 July, 2022; originally announced July 2022.

Comments: 24 pages, 10 figures. Posted after referee report; comments welcome. Key results in Figures 8 & 9. New SFH model publicly available as "continuity_psb_sfh" in Prospector: https://github.com/bd-j/prospector

Probing cold gas in a massive, compact star-forming galaxy at z=6

Authors: Jorge A. Zavala, Caitlin M. Casey, Justin Spilker, Ken-ichi Tadaki, Akiyoshi Tsujita, Jaclyn Champagne, Daisuke Iono, Kotaro Kohno, Sinclaire Manning, Alfredo Montana

Abstract: Observations of low order CO transitions represent the most direct way to study galaxies' cold molecular gas, the fuel of star formation. Here we present the first detection of CO(2-1) in a galaxy lying on the main-sequence of star-forming galaxies at z>6. Our target, G09-83808 at z=6.03, has a short depletion time-scale of T_dep~50Myr and a relatively low gas fraction of M_gas/M_star=0.30 that co… ▽ More Observations of low order CO transitions represent the most direct way to study galaxies' cold molecular gas, the fuel of star formation. Here we present the first detection of CO(2-1) in a galaxy lying on the main-sequence of star-forming galaxies at z>6. Our target, G09-83808 at z=6.03, has a short depletion time-scale of T_dep~50Myr and a relatively low gas fraction of M_gas/M_star=0.30 that contrasts with those measured for lower redshift main-sequence galaxies. We conclude that this galaxy is undergoing a starburst episode with a high star formation efficiency that might be the result of gas compression within its compact rotating disk. Its starburst-like nature is further supported by its high star formation rate surface density, thus favoring the use of the Kennicutt-Schmidt relation as a more precise diagnostic diagram for starbursts. Without further significant gas accretion, this galaxy would become a compact, massive quiescent galaxy at z~5.5. In addition, we find that the calibration for estimating ISM masses from dust continuum emission satisfactorily reproduces the gas mass derived from the CO(2-1) transition (within a factor of ~2). This is in line with previous studies claiming a small redshift evolution in the gas-to-dust ratio of massive, metal-rich galaxies. In the absence of gravitational amplification, this detection would have required of order ~1000h of observing time. The detection of cold molecular gas in unlensed star-forming galaxies at high redshifts is thus prohibitive with current facilities and requires a ten-fold improvement in sensitivity, such as that envisaged for the ngVLA. △ Less

Submitted 6 June, 2022; originally announced June 2022.

Comments: Accepted for publication in The Astrophysical Journal

Searching Far and Long I: Pilot ALMA 2mm Follow-up of Bright Dusty Galaxies as a Redshift Filter

Authors: Olivia R. Cooper, Caitlin M. Casey, Jorge A. Zavala, Jaclyn B. Champagne, Elisabete da Cunha, Arianna S. Long, Justin S. Spilker, Johannes Staguhn

Abstract: A complete census of dusty star-forming galaxies (DSFGs) at early epochs is necessary to constrain the obscured contribution to the cosmic star formation rate density (CSFRD), however DSFGs beyond $z \sim 4$ are both rare and hard to identify from photometric data alone due to degeneracies in submillimeter photometry with redshift. Here, we present a pilot study obtaining follow-up Atacama Large M… ▽ More A complete census of dusty star-forming galaxies (DSFGs) at early epochs is necessary to constrain the obscured contribution to the cosmic star formation rate density (CSFRD), however DSFGs beyond $z \sim 4$ are both rare and hard to identify from photometric data alone due to degeneracies in submillimeter photometry with redshift. Here, we present a pilot study obtaining follow-up Atacama Large Millimeter Array (ALMA) $2\,$mm observations of a complete sample of 39 $850\,\rmμm$-bright dusty galaxies in the SSA22 field. Empirical modeling suggests $2\,$mm imaging of existing samples of DSFGs selected at $850\,\rmμm - 1\,$mm can quickly and easily isolate the "needle in a haystack" DSFGs that sit at $z>4$ or beyond. Combining archival submillimeter imaging with our measured ALMA $2\,$mm photometry ($1σ\sim 0.08\,$mJy$\,$beam$^{-1}$ rms), we characterize the galaxies' IR SEDs and use them to constrain redshifts. With available redshift constraints fit via the combination of six submillimeter bands, we identify 6/39 high-$z$ candidates each with $>50\%$ likelihood to sit at $z > 4$, and find a positive correlation between redshift and $2\,$mm flux density. Specifically, our models suggest the addition of $2\,$mm to a moderately constrained IR SED will improve the accuracy of a millimeter-derived redshift from $Δz/(1+z) = 0.3$ to $Δz/(1+z) = 0.2$. Our IR SED characterizations provide evidence for relatively high emissivity spectral indices ($\langle β\rangle = 2.4\pm0.3$) in the sample. We measure that especially bright ($S_{850\rmμm}>5.55\,$mJy) DSFGs contribute $\sim10$% to the cosmic-averaged CSFRD from $2<z<5$, confirming findings from previous work with similar samples. △ Less

Submitted 28 March, 2022; originally announced March 2022.

Comments: 22 pages, 7 figures, accepted for publication in ApJ

Chaotic and Clumpy Galaxy Formation in an Extremely Massive Reionization-Era Halo

Authors: Justin S. Spilker, Christopher C. Hayward, Daniel P. Marrone, Manuel Aravena, Matthieu Bethermin, James Burgoyne, Scott C. Chapman, Thomas R. Greve, Gayathri Gururajan, Yashar D. Hezaveh, Ryley Hill, Katrina C. Litke, Christopher C. Lovell, Matthew A. Malkan, Eric J. Murphy, Desika Narayanan, Kedar A. Phadke, Cassie Reuter, Antony A. Stark, Nikolaus Sulzenauer, Joaquin D. Vieira, David Vizgan, Axel Weiss

Abstract: The SPT0311-58 system at z=6.900 is an extremely massive structure within the reionization epoch, and offers a chance to understand the formation of galaxies in an extreme peak in the primordial density field. We present 70mas Atacama Large Millimeter/submillimeter Array observations of the dust continuum and CII 158um emission in the central pair of galaxies and reach physical resolution ~100-350… ▽ More The SPT0311-58 system at z=6.900 is an extremely massive structure within the reionization epoch, and offers a chance to understand the formation of galaxies in an extreme peak in the primordial density field. We present 70mas Atacama Large Millimeter/submillimeter Array observations of the dust continuum and CII 158um emission in the central pair of galaxies and reach physical resolution ~100-350pc, among the most detailed views of any reionization-era system to date. The observations resolve the source into at least a dozen kiloparsec-size clumps. The global kinematics and high turbulent velocity dispersion within the galaxies present a striking contrast to recent claims of dynamically cold thin-disk kinematics in some dusty galaxies just 800Myr later at z~4. We speculate that both gravitational interactions and fragmentation from massive parent disks have likely played a role in the overall dynamics and formation of clumps in the system. Each clump individually is comparable in mass to other 6<z<8 galaxies identified in rest-UV/optical deep field surveys, but with star formation rates elevated by ~3-5x. Internally, the clumps themselves bear close resemblance to greatly scaled-up versions of virialized cloud-scale structures identified in low-redshift galaxies. Our observations are qualitatively similar to the chaotic and clumpy assembly within massive halos seen in simulations of high-redshift galaxies. △ Less

Submitted 28 March, 2022; originally announced March 2022.

Comments: Accepted for publication in ApJ Letters, 9 pages, 5 figures. Interactive versions of Fig. 2 and data products are available at https://github.com/spt-smg/publicdata

The Compact Structures of Massive $z\sim0.7$ Post-Starburst Galaxies in the SQuIGG$\vec{L}$E Sample

Authors: David J. Setton, Margaret Verrico, Rachel Bezanson, Jenny E. Greene, Katherine A. Suess, Andy D. Goulding, Justin S. Spilker, Mariska Kriek, Robert Feldmann, Desika Narayanan, Khalil Hall-Hooper, Erin Kado-Fong

Abstract: We present structural measurements of 145 spectroscopically selected intermediate-redshift (z$\sim$0.7), massive ($M_\star \sim 10^{11} \ M_\odot$) post-starburst galaxies from the SQuIGG$\vec{L}$E Sample measured using wide-depth Hyper Suprime-Cam i-band imaging. This deep imaging allows us to probe the sizes and structures of these galaxies, which we compare to a control sample of star forming a… ▽ More We present structural measurements of 145 spectroscopically selected intermediate-redshift (z$\sim$0.7), massive ($M_\star \sim 10^{11} \ M_\odot$) post-starburst galaxies from the SQuIGG$\vec{L}$E Sample measured using wide-depth Hyper Suprime-Cam i-band imaging. This deep imaging allows us to probe the sizes and structures of these galaxies, which we compare to a control sample of star forming and quiescent galaxies drawn from the LEGA-C Survey. We find that post-starburst galaxies systematically lie $\sim0.1$ dex below the quiescent mass-size (half-light radius) relation, with a scatter of $\sim0.2$ dex. This finding is bolstered by non-parametric measures, such as the Gini coefficient and the concentration, which also reveal these galaxies to have more compact light profiles than both quiescent and star-forming populations at similar mass and redshift. The sizes of post-starburst galaxies show either negative or no correlation with the time since quenching, such that more recently quenched galaxies are larger or similarly sized. This empirical finding disfavors the formation of post-starburst galaxies via a purely central burst of star formation that simultaneously shrinks the galaxy and shuts off star formation. We show that the central densities of post-starburst and quiescent galaxies at this epoch are very similar, in contrast with their effective radii. The structural properties of z$\sim$0.7 post-starburst galaxies match those of quiescent galaxies that formed in the early universe, suggesting that rapid quenching in the present epoch is driven by a similar mechanism to the one at high redshift. △ Less

Submitted 13 June, 2022; v1 submitted 16 March, 2022; originally announced March 2022.

Comments: 28 pages, 19 figures, 2 appendices, updated to accepted ApJ version

REQUIEM-2D: A diversity of formation pathways in a sample of spatially-resolved massive quiescent galaxies at z~2

Authors: Mohammad Akhshik, Katherine E. Whitaker, Joel Leja, Johan Richard, Justin S. Spilker, Mimi Song, Gabriel Brammer, Rachel Bezanson, Harald Ebeling, Anna R. Gallazzi, Guillaume Mahler, Lamiya A. Mowla, Erica J. Nelson, Camilla Pacifici, Keren Sharon, Sune Toft, Christina C. Williams, Lillian Wright, Johannes Zabl

Abstract: REQUIEM-2D (REsolving QUIEscent Magnified galaxies with 2D grism spectroscopy) is comprised of a sample of 8 massive ($\log M_*/M_\odot > 10.6$) strongly lensed quiescent galaxies at $z\sim2$. REQUIEM-2D combines the natural magnification from strong gravitational lensing with the high spatial-resolution grism spectroscopy of \emph{Hubble Space Telescope} through a spectrophotometric fit to study… ▽ More REQUIEM-2D (REsolving QUIEscent Magnified galaxies with 2D grism spectroscopy) is comprised of a sample of 8 massive ($\log M_*/M_\odot > 10.6$) strongly lensed quiescent galaxies at $z\sim2$. REQUIEM-2D combines the natural magnification from strong gravitational lensing with the high spatial-resolution grism spectroscopy of \emph{Hubble Space Telescope} through a spectrophotometric fit to study spatially resolved stellar populations. We show that quiescent galaxies in the REQUIEM-2D survey have diverse formation histories manifesting as a gradient in stellar ages, including examples of (1) a younger central region supporting outside-in formation, (2) flat age gradients that show evidence for both spatially-uniform early formation or inside-out quenching, and (3) regions at a fixed radial distance having different ages (such asymmetries cannot be recovered when averaging stellar population measurements azimuthally). The typical dust attenuation curve for the REQUIEM-2D galaxies is constrained to be steeper than Calzetti's law in the UV and generally consistent with $A_V<1$. Combined together and accounting for the different physical radial distances and formation time-scales, we find that the REQUIEM-2D galaxies that formed earlier in the universe exhibit slow and uniform growth in their inner core, whereas the galaxies that formed later have rapid inner growth in their inner core with younger ages relative to the outskirts. These results challenge the currently accepted paradigm of how massive quiescent galaxies form, where the earliest galaxies are thought to form most rapidly. Significantly larger samples close to the epoch of formation with similar data quality and higher spectral resolution are required to validate this finding. △ Less

Submitted 9 March, 2022; originally announced March 2022.

Comments: Submitted to ApJ

Multi-Phase ISM in the z = 5.7 Hyperluminous Starburst SPT0346-52

Authors: Katrina C. Litke, Daniel P. Marrone, Manuel Aravena, Matthieu Bethermin, Scott C. Chapman, Chenxing Dong, Christopher C. Hayward, Ryley Hill, Sreevani Jarugula, Matthew A. Malkan, Desika Narayanan, Cassie A. Reuter, Justin S. Spilker, Nikolaus Sulzenauer, Joaquin D. Vieira, Axel Weiss

Abstract: SPT0346-52 (z=5.7) is the most intensely star-forming galaxy discovered by the South Pole Telescope, with Sigma_SFR ~ 4200 Msol yr^-1 kpc^-2. In this paper, we expand on previous spatially-resolved studies, using ALMA observations of dust continuum, [NII]205 micron, [CII]158 micron, [OI]146 micron, and undetected [NII]122 micron and [OI]63 micron emission to study the multi-phase interstellar medi… ▽ More SPT0346-52 (z=5.7) is the most intensely star-forming galaxy discovered by the South Pole Telescope, with Sigma_SFR ~ 4200 Msol yr^-1 kpc^-2. In this paper, we expand on previous spatially-resolved studies, using ALMA observations of dust continuum, [NII]205 micron, [CII]158 micron, [OI]146 micron, and undetected [NII]122 micron and [OI]63 micron emission to study the multi-phase interstellar medium (ISM) in SPT0346-52. We use pixelated, visibility-based lens modeling to reconstruct the source-plane emission. We also model the source-plane emission using the photoionization code CLOUDY and find a supersolar metallicity system. We calculate T_dust = 48.3 K and lambda_peak = 80 micron, and see line deficits in all five lines. The ionized gas is less dense than comparable galaxies, with n_e < 32 cm^-3, while ~20% of the [CII]158 emission originates from the ionized phase of the ISM. We also calculate the masses of several phases of the ISM. We find that molecular gas dominates the mass of the ISM in SPT0346-52, with the molecular gas mass ~4x higher than the neutral atomic gas mass and ~100x higher than the ionized gas mass. △ Less

Submitted 25 March, 2022; v1 submitted 24 February, 2022; originally announced February 2022.

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

ALMA measures molecular gas reservoirs comparable to field galaxies in a low-mass galaxy cluster at z=1.3

Authors: Christina C. Williams, Stacey Alberts, Justin S. Spilker, Allison G. Noble, Mauro Stefanon, Christopher N. A. Willmer, Rachel Bezanson, Desika Narayanan, Katherine E. Whitaker

Abstract: We report the serendipitous discovery of an overdensity of CO emitters in an X-ray-identified cluster (Log$_{10}$M$_{\rm halo}/M_{\odot}\sim13.6$ at z=1.3188) using ALMA. We present spectroscopic confirmation of 6 new cluster members exhibiting CO(2-1) emission, adding to 2 existing optical/IR spectroscopic members undetected in CO. This is the lowest mass cluster to date at z>1 with molecular gas… ▽ More We report the serendipitous discovery of an overdensity of CO emitters in an X-ray-identified cluster (Log$_{10}$M$_{\rm halo}/M_{\odot}\sim13.6$ at z=1.3188) using ALMA. We present spectroscopic confirmation of 6 new cluster members exhibiting CO(2-1) emission, adding to 2 existing optical/IR spectroscopic members undetected in CO. This is the lowest mass cluster to date at z>1 with molecular gas measurements, bridging the observational gap between galaxies in the more extreme, well-studied clusters (Log$_{10}$~M$_{\rm halo}/M_{\odot}\gtrsim14$) and those in group or field environments at cosmic noon. The CO sources are concentrated on the sky (within ~1-arcmin diameter) and phase space analysis indicates the gas resides in galaxies already within the cluster environment. We find that CO sources sit in similar phase space as CO-rich galaxies in more massive clusters at similar redshifts (have similar accretion histories) while maintaining field-like molecular gas reservoirs, compared to scaling relations. This work presents the deepest CO survey to date in a galaxy cluster at z>1, uncovering gas reservoirs down to M$_{\rm H_{2}}>1.6\times10^{10}$M$_{\odot}$ (5$σ$ at 50% primary beam). Our deep limits rule out the presence of gas content in excess of the field scaling relations; however, combined with literature CO detections, cluster gas fractions in general appear systematically high, on the upper envelope or above the field. This study is the first demonstration that low mass clusters at z~1-2 can host overdensities of CO emitters with surviving gas reservoirs, in line with the prediction that quenching is delayed after first infall while galaxies consume the gas bound to the disk. △ Less

Submitted 4 January, 2022; originally announced January 2022.

Comments: Submitted to AAS Journals on December 15, 2021. Comments welcome

SQuIGGLE: Studying Quenching in Intermediate-z Galaxies -- Gas, AnguLar Momentum, and Evolution

Authors: Katherine A. Suess, Mariska Kriek, Rachel Bezanson, Jenny E. Greene, David Setton, Justin S. Spilker, Robert Feldmann, Andy D. Goulding, Benjamin D. Johnson, Joel Leja, Desika Narayanan, Khalil Hall-Hooper, Qiana Hunt, Sidney Lower, Margaret Verrico

Abstract: We describe the SQuIGGLE survey of intermediate-redshift post-starburst galaxies. We leverage the large sky coverage of the SDSS to select ~1300 recently-quenched galaxies at 0.5<z<~0.9 based on their unique spectral shapes. These bright, intermediate-redshift galaxies are ideal laboratories to study the physics responsible for the rapid quenching of star formation: they are distant enough to be u… ▽ More We describe the SQuIGGLE survey of intermediate-redshift post-starburst galaxies. We leverage the large sky coverage of the SDSS to select ~1300 recently-quenched galaxies at 0.5<z<~0.9 based on their unique spectral shapes. These bright, intermediate-redshift galaxies are ideal laboratories to study the physics responsible for the rapid quenching of star formation: they are distant enough to be useful analogs for high-redshift quenching galaxies, but low enough redshift that multi-wavelength follow-up observations are feasible with modest telescope investments. We use the Prospector code to infer the stellar population properties and non-parametric star formation histories of all galaxies in the sample. We find that SQuIGGLE galaxies are both very massive (M* ~ 10^11.25 Msun) and quenched, with inferred star formation rates <~1Msun/yr, more than an order of magnitude below the star-forming main sequence. The best-fit star formation histories confirm that these galaxies recently quenched a major burst of star formation: >75% of SQuIGGLE galaxies formed at least a quarter of their total stellar mass in the recent burst, which ended just ~200Myr before observation. We find that SQuIGGLE galaxies are on average younger and more burst-dominated than most other z<~1 post-starburst samples. This large sample of bright post-starburst galaxies at intermediate redshift opens a wide range of studies into the quenching process. In particular, the full SQuIGGLE survey will investigate the molecular gas reservoirs, morphologies, kinematics, resolved stellar populations, AGN incidence, and infrared properties of this unique sample of galaxies in order to place definitive constraints on the quenching process. △ Less

Submitted 29 November, 2021; originally announced November 2021.

Comments: 23 pages, 16 figures, accepted to ApJ

Now you see it, now you don't: Star formation truncation precedes the loss of molecular gas by ~100 Myr in massive post-starburst galaxies at z~0.6

Authors: Rachel Bezanson, Justin S. Spilker, Katherine A. Suess, David J. Setton, Robert Feldmann, Jenny E. Greene, Mariska Kriek, Desika Narayanan, Margaret Verrico

Abstract: We use ALMA observations of CO(2-1) in 13 massive ($M_{\star}\gtrsim 10^{11} M_{\odot}$) post-starburst galaxies at $z\sim0.6$ to constrain the molecular gas content in galaxies shortly after they quench their major star-forming episode. The post-starburst galaxies in this study are selected from the Sloan Digital Sky Survey spectroscopic samples (DR14) based on their spectral shapes, as part of t… ▽ More We use ALMA observations of CO(2-1) in 13 massive ($M_{\star}\gtrsim 10^{11} M_{\odot}$) post-starburst galaxies at $z\sim0.6$ to constrain the molecular gas content in galaxies shortly after they quench their major star-forming episode. The post-starburst galaxies in this study are selected from the Sloan Digital Sky Survey spectroscopic samples (DR14) based on their spectral shapes, as part of the SQuIGGLE program. Early results showed that two post-starburst galaxies host large H$_2$ reservoirs despite their low inferred star formation rates. Here we expand this analysis to a larger statistical sample of 13 galaxies. Six of the primary targets (45%) are detected, with $M_{H_2}\gtrsim10^9 M_{\odot}$. Given their high stellar masses, this mass limit corresponds to an average gas fraction of $\langle f_{H_2} \equiv M_{H_2}/M_{\star} \rangle \sim7\%$, or ${\sim}14\%$ using lower stellar masses estimates derived from analytic, exponentially declining star formation histories. The gas fraction correlates with the $D_n4000$ spectral index, suggesting that the cold gas reservoirs decrease with time since burst, as found in local K+A galaxies. Star formation histories derived from flexible stellar population synthesis modeling support this empirical finding: galaxies that quenched $\lesssim 150$ Myr prior to observation host detectable CO(2-1) emission, while older post-starburst galaxies are undetected. The large $\mathrm{H_2}$ reservoirs and low star formation rates in the sample imply that the quenching of star formation precedes the disappearance of the cold gas reservoirs. However, within the following 100-200 Myrs, the SQuIGGLE galaxies require the additional and efficient heating or removal of cold gas to bring their low star formation rates in line with standard $\mathrm{H_2}$ scaling relations. △ Less

Submitted 29 November, 2021; originally announced November 2021.

Comments: 20 pages, 13 figures, accepted for publication in ApJ

High Molecular-Gas to Dust Mass Ratios Predicted in Most Quiescent Galaxies

Authors: Katherine E. Whitaker, Desika Narayanan, Christina C. Williams, Qi Li, Justin S. Spilker, Romeel Davé, Mohammad Akhshik, Hollis B. Akins, Rachel Bezanson, Neal Katz, Joel Leja, Georgios E. Magdis, Lamiya Mowla, Erica J. Nelson, Alexandra Pope, George C. Privon, Sune Toft, Francesco Valentino

Abstract: Observations of cold molecular gas reservoirs are critical for understanding the shutdown of star formation in massive galaxies. While dust continuum is an efficient and affordable tracer, this method relies upon the assumption of a "normal" molecular-gas to dust mass ratio, $δ_{\mathrm{GDR}}$, typically of order one hundred. Recent null detections of quiescent galaxies in deep dust continuum obse… ▽ More Observations of cold molecular gas reservoirs are critical for understanding the shutdown of star formation in massive galaxies. While dust continuum is an efficient and affordable tracer, this method relies upon the assumption of a "normal" molecular-gas to dust mass ratio, $δ_{\mathrm{GDR}}$, typically of order one hundred. Recent null detections of quiescent galaxies in deep dust continuum observations support a picture where the cold gas and dust has been rapidly depleted or expelled. In this work, we present another viable explanation: a significant fraction of galaxies with low star formation per unit stellar mass are predicted to have extreme $δ_{\mathrm{GDR}}$ ratios. We show that simulated massive quiescent galaxies at $0 < z < 3$ in the \textsc{simba} cosmological simulations have $δ_{\mathrm{GDR}}$ values that extend $>$4 orders of magnitude. The dust in most simulated quiescent galaxies is destroyed significantly more rapidly than the molecular gas depletes, and cannot be replenished. The transition from star-forming to quiescent halts dust formation via star formation processes, with dust subsequently destroyed by supernova shocks and thermal sputtering of dust grains embedded in hot plasma. After this point, the dust growth rate in the models is not sufficient to overcome the loss of $>$3 orders of magnitude in dust mass to return to normal values of $δ_{\mathrm{GDR}}$ despite having high metallicity. Our results indicate that it is not straight forward to use a single observational indicator to robustly pre-select exotic versus normal ratios. These simulations make strong predictions that can be tested with millimeter facilities. △ Less

Submitted 10 November, 2021; v1 submitted 9 November, 2021; originally announced November 2021.

Comments: 9 pages, 5 figures, version accepted for publication in the Astrophysical Journal Letters

The ramp-up of interstellar medium enrichment at z>4

Authors: M. Franco, K. E. K. Coppin, J. E. Geach, C. Kobayashi, S. C. Chapman, C. Yang, E. González-Alfonso, J. S. Spilker, A. Cooray, M. J. Michałowski

Abstract: Fluorine is one of the most interesting elements for nuclear and stellar astrophysics. Fluorine abundance was first measured for stars other than the Sun in 1992, then for a handful metal-poor stars, which are likely to have formed in the early Universe. The main production sites of fluorine are under debate and include asymptotic giant branch (AGB) stars, $ν$-process in core-collapse supernovae,… ▽ More Fluorine is one of the most interesting elements for nuclear and stellar astrophysics. Fluorine abundance was first measured for stars other than the Sun in 1992, then for a handful metal-poor stars, which are likely to have formed in the early Universe. The main production sites of fluorine are under debate and include asymptotic giant branch (AGB) stars, $ν$-process in core-collapse supernovae, and Wolf-Rayet (WR) stars. Due to the difference in the mass and lifetime of progenitor stars, high redshift observations of fluorine can help constrain the mechanism of fluorine production in massive galaxies. Here, we report the detection of HF (S/N = 8) in absorption in a gravitationally lensed dusty star-forming galaxy at redshift z=4.4 with $N_{\rm HF}$/$N_{\rm{H_2}}$ as high as $\sim2\times10^{-9}$, indicating a very quick ramp-up of the chemical enrichment in this high-z galaxy. At z=4.4, AGB stars of a few solar masses are very unlikely to dominate the enrichment. Instead, we show that WR stars are required to produce the observed fluorine abundance at this time, with other production mechanisms becoming important at later times. These observations therefore provide an insight into the underlying processes driving the `ramp-up' phase of chemical enrichment alongside rapid stellar mass assembly in a young massive galaxy. △ Less

Submitted 4 November, 2021; originally announced November 2021.

Comments: 21 pages, 4 figures, published in Nature Astronomy

Characterization of Two 2mm-detected Optically-Obscured Dusty Star-Forming Galaxies

Authors: Sinclaire M. Manning, Caitlin M. Casey, Jorge A. Zavala, Georgios E. Magdis, Patrick M. Drew, Jaclyn B. Champagne, Manuel Aravena, Matthieu Béthermin, David L. Clements, Steven L. Finkelstein, Seiji Fujimoto, Christopher C. Hayward, Jacqueline A. Hodge, Olivier Ilbert, Jeyhan S. Kartaltepe, Kirsten K. Knudsen, Anton M. Koekemoer, Allison W. S. Man, David B. Sanders, Kartik Sheth, Justin S. Spilker, Johannes Staguhn, Margherita Talia, Ezequiel Treister, Min S. Yun

Abstract: The 2mm Mapping Obscuration to Reionization with ALMA (MORA) Survey was designed to detect high redshift ($z\gtrsim4$), massive, dusty star-forming galaxies (DSFGs). Here we present two, likely high redshift sources, identified in the survey whose physical characteristics are consistent with a class of optical/near-infrared (OIR) invisible DSFGs found elsewhere in the literature. We first perform… ▽ More The 2mm Mapping Obscuration to Reionization with ALMA (MORA) Survey was designed to detect high redshift ($z\gtrsim4$), massive, dusty star-forming galaxies (DSFGs). Here we present two, likely high redshift sources, identified in the survey whose physical characteristics are consistent with a class of optical/near-infrared (OIR) invisible DSFGs found elsewhere in the literature. We first perform a rigorous analysis of all available photometric data to fit spectral energy distributions and estimate redshifts before deriving physical properties based on our findings. Our results suggest the two galaxies, called MORA-5 and MORA-9, represent two extremes of the "OIR-dark" class of DSFGs. MORA-5 ($z_{\rm phot}=4.3^{+1.5}_{-1.3}$) is a significantly more active starburst with a star-formation rate of 830$^{+340}_{-190}$M$_\odot$yr$^{-1}$ compared to MORA-9 ($z_{\rm phot}=4.3^{+1.3}_{-1.0}$) whose star-formation rate is a modest 200$^{+250}_{-60}$M$_\odot$yr$^{-1}$. Based on the stellar masses (M$_{\star}\approx10^{10-11}$M$_\odot$), space density ($n\sim(5\pm2)\times10^{-6}$Mpc$^{-3}$, which incorporates two other spectroscopically confirmed OIR-dark DSFGs in the MORA sample at $z=4.6$ and $z=5.9$), and gas depletion timescales ($<1$Gyr) of these sources, we find evidence supporting the theory that OIR-dark DSFGs are the progenitors of recently discovered $3<z<4$ massive quiescent galaxies. △ Less

Submitted 3 November, 2021; originally announced November 2021.

Comments: 18 pages, 7 figures, accepted for publication in ApJ

Mapping Obscuration to Reionization with ALMA (MORA): 2mm Efficiently Selects the Highest-Redshift Obscured Galaxies

Authors: Caitlin M. Casey, Jorge A. Zavala, Sinclaire M. Manning, Manuel Aravena, Matthieu Béthermin, Karina I. Caputi, Jaclyn B. Champagne, David L. Clements, Patrick Drew, Steven L. Finkelstein, Seiji Fujimoto, Christopher C. Hayward, Anton M. Koekemoer, Vasily Kokorev, Claudia del P. Lagos, Arianna S. Long, Georgios E. Magdis, Allison W. S. Man, Ikki Mitsuhashi, Gergö Popping, Justin Spilker, Johannes Staguhn, Margherita Talia, Sune Toft, Ezequiel Treister , et al. (2 additional authors not shown)

Abstract: We present the characteristics of 2mm-selected sources from the largest Atacama Large Millimeter and submillimeter Array (ALMA) blank-field contiguous survey conducted to-date, the Mapping Obscuration to Reionization with ALMA (MORA) survey covering 184arcmin$^2$ at 2mm. Twelve of the thirteen detections above 5$σ$ are attributed to emission from galaxies, eleven of which are dominated by cold dus… ▽ More We present the characteristics of 2mm-selected sources from the largest Atacama Large Millimeter and submillimeter Array (ALMA) blank-field contiguous survey conducted to-date, the Mapping Obscuration to Reionization with ALMA (MORA) survey covering 184arcmin$^2$ at 2mm. Twelve of the thirteen detections above 5$σ$ are attributed to emission from galaxies, eleven of which are dominated by cold dust emission. These sources have a median redshift of $\langle z_{\rm 2mm}\rangle=3.6^{+0.4}_{-0.3}$ primarily based on optical/near-infrared (OIR) photometric redshifts with some spectroscopic redshifts, with 77$\pm$11% of sources at $z>3$ and 38$\pm$12% of sources at $z>4$. This implies that 2mm selection is an efficient method for identifying the highest redshift dusty star-forming galaxies (DSFGs). Lower redshift DSFGs ($z<3$) are far more numerous than those at $z>3$ yet likely to drop out at 2mm. MORA shows that DSFGs with star-formation rates in excess of 300M$_\odot$ yr$^{-1}$ and relative rarity of $\sim$10$^{-5}$ Mpc$^{-3}$ contribute $\sim$30% to the integrated star-formation rate density between $3<z<6$. The volume density of 2mm-selected DSFGs is consistent with predictions from some cosmological simulations and is similar to the volume density of their hypothesized descendants: massive, quiescent galaxies at $z>2$. Analysis of MORA sources' spectral energy distributions hint at steeper empirically-measured dust emissivity indices than typical literature studies, with $\langleβ\rangle=2.2^{+0.5}_{-0.4}$. The MORA survey represents an important step in taking census of obscured star-formation in the Universe's first few billion years, but larger area 2mm surveys are needed to more fully characterize this rare population and push to the detection of the Universe's first dusty galaxies. △ Less

Submitted 13 October, 2021; originally announced October 2021.

Comments: 42 pages, 14 figures, accepted for publication in ApJ

Quenching of star formation from a lack of inflowing gas to galaxies

Authors: Katherine E. Whitaker, Christina C. Williams, Lamiya Mowla, Justin S. Spilker, Sune Toft, Desika Narayanan, Alexandra Pope, Georgios E. Magdis, Pieter G. van Dokkum, Mohammad Akhshik, Rachel Bezanson, Gabriel B. Brammer, Joel Leja, Allison Man, Erica J. Nelson, Johan Richard, Camilla Pacifici, Keren Sharon, Francesco Valentino

Abstract: Star formation in half of massive galaxies was quenched by the time the Universe was three billion years old. Very low amounts of molecular gas appear responsible for this, at least in some cases, though morphological gas stabilization, shock heating, or activity associated with accretion onto a central supermassive black hole is invoked in other cases. Recent studies of quenching by gas depletion… ▽ More Star formation in half of massive galaxies was quenched by the time the Universe was three billion years old. Very low amounts of molecular gas appear responsible for this, at least in some cases, though morphological gas stabilization, shock heating, or activity associated with accretion onto a central supermassive black hole is invoked in other cases. Recent studies of quenching by gas depletion have been based upon upper limits that are insufficiently sensitive to determine this robustly, or stacked emission with its problems of averaging. Here we report 1.3mm observations of dust emission from six strongly lensed galaxies where star formation has been quenched, with magnifications of up to a factor of 30. Four of the six galaxies are undetected in dust emission, with an estimated upper limit on the dust mass of 0.0001 times the stellar mass, and by proxy (assuming a Milky Way molecular gas-to-dust ratio) 0.01 times the stellar mass in molecular gas. This is two orders of magnitude less molecular gas per unit stellar mass than seen in star forming galaxies at similar redshifts. It remains difficult to extrapolate from these small samples, but these observations establish that gas depletion is responsible for a cessation of star formation in some fraction of high-redshift galaxies. △ Less

Submitted 21 September, 2021; originally announced September 2021.

Comments: 17 pages, 3 figures. Authors' version. Published online by Nature on September 22, 2021

Rapid build-up of the stellar content in the protocluster core SPT2349$-$56 at $z\,{=}\,4.3$

Authors: Ryley Hill, Scott Chapman, Kedar A. Phadke, Manuel Aravena, Melanie Archipley, Matthew L. N. Ashby, Matthieu Bethermin, Rebecca E. A. Canning, Anthony Gonzalez, Thomas R. Greve, Gayathri Gururajan, Christopher C. Hayward, Yashar Hezaveh, Sreevani Jarugula, Duncan MacIntyre, Daniel P. Marrone, Tim Miller, Douglas Rennehan, Cassie Reuter, Kaja Rotermund, Douglas Scott, Justin Spilker, Joaquin D. Vieira, George Wang, Axel Weiss

Abstract: The protocluster SPT2349$-$56 at $z\,{=}\,4.3$ contains one of the most actively star-forming cores known, yet constraints on the total stellar mass of this system are highly uncertain. We have therefore carried out deep optical and infrared observations of this system, probing rest-frame ultraviolet to infrared wavelengths. Using the positions of the spectroscopically-confirmed protocluster membe… ▽ More The protocluster SPT2349$-$56 at $z\,{=}\,4.3$ contains one of the most actively star-forming cores known, yet constraints on the total stellar mass of this system are highly uncertain. We have therefore carried out deep optical and infrared observations of this system, probing rest-frame ultraviolet to infrared wavelengths. Using the positions of the spectroscopically-confirmed protocluster members, we identify counterparts and perform detailed source deblending, allowing us to fit spectral energy distributions in order to estimate stellar masses. We show that the galaxies in SPT2349$-$56 have stellar masses proportional to their high star-formation rates, consistent with other protocluster galaxies and field submillimetre galaxies (SMGs) around redshift 4. The galaxies in SPT2349$-$56 have on average lower molecular gas-to-stellar mass fractions and depletion timescales than field SMGs, although with considerable scatter. We construct the stellar-mass function for SPT2349$-$56 and compare it to the stellar-mass function of $z\,{=}\,1$ galaxy clusters, finding consistent shapes between the two. We measure rest-frame galaxy ultraviolet half-light radii from our HST-F160W imaging, finding that on average the galaxies in our sample are similar in size to typical star-forming galaxies at these redshifts. However, the brightest HST-detected galaxy in our sample, found near the luminosity-weighted centre of the protocluster core, remains unresolved at this wavelength. Hydrodynamical simulations predict that the core galaxies will quickly merge into a brightest cluster galaxy, thus our observations provide a direct view of the early formation mechanisms of this class of object. △ Less

Submitted 7 April, 2022; v1 submitted 9 September, 2021; originally announced September 2021.

Comments: Submitted to MNRAS Sept. 3, 2021 Accepted by MNRAS Dec, 4, 2021

Journal ref: MNRAS 512(3), 4352, 2022

High resolution spectral imaging of CO(7-6), [CI](2-1) and continuum of three high-z lensed dusty star-forming galaxies using ALMA

Authors: G. Gururajan, M. Béthermin, P. Theulé, J. S. Spilker, M. Aravena, M. A. Archipley, S. C. Chapman, C. DeBreuck, A. Gonzalez, C. C. Hayward, Y. Hezaveh, R. Hill, S. Jarugula, K. C. Litke, M. Malkan, D. Marrone, D. Narayanan, K. A. Phadke, C. Reuter, J. Vieira, D. Vizgan, A. Weiß

Abstract: High-redshift dusty star-forming galaxies with very high star formation rates (500 -- 3000 M$_{\odot}$ yr$^{-1}$) are key to understanding the formation of the most extreme galaxies in the early Universe. Characterising the gas reservoir of these systems can reveal the driving factor behind the high star formation. Using molecular gas tracers like high-J CO lines, neutral carbon lines and the dust… ▽ More High-redshift dusty star-forming galaxies with very high star formation rates (500 -- 3000 M$_{\odot}$ yr$^{-1}$) are key to understanding the formation of the most extreme galaxies in the early Universe. Characterising the gas reservoir of these systems can reveal the driving factor behind the high star formation. Using molecular gas tracers like high-J CO lines, neutral carbon lines and the dust continuum, we can estimate the gas density and radiation field intensity in their interstellar media. In this paper, we present high resolution ($\sim$0.4$^{\prime\prime}$) observations of CO(7-6), [CI](2-1) and dust continuum of 3 lensed galaxies from the SPT-SMG sample at $z\sim$ 3 with the Atacama Large Millimeter/submillimeter Array. Our sources have high intrinsic star-formation rates ($>$850 M$_{\odot}$yr$^{-1}$) and rather short depletion timescales ($<$100 Myr). Based on the L$_{[\rm CI](2-1)}$/L$_{\rm CO(7-6)}$ and L$_{[\rm CI](2-1)}$/L$_{\rm IR}$ ratios, our galaxy sample has similar radiation field intensities and gas densities compared to other submillimetre galaxies. We perform visibility-based lens modelling on these objects to reconstruct the kinematics in the source plane. We find that the cold gas masses of the sources are compatible with simple dynamical mass estimates using ULIRG-like values of the CO-H$_2$ conversion factor $α_{\rm CO}$ but not Milky Way-like values. We find diverse source kinematics in our sample: SPT0103-45 and SPT2147-50 are likely rotating disks while SPT2357-51 is possibly a major merger. The analysis presented in the paper could be extended to a larger sample to determine better statistics of morphologies and interstellar medium properties of high-$z$ dusty star-forming galaxies. △ Less

Submitted 15 March, 2022; v1 submitted 8 September, 2021; originally announced September 2021.

Comments: 29 pages, 26 figures, Accepted for publication in A&A. Figure 16 and the related conclusions have been revised

Journal ref: A&A 663, A22 (2022)

Molecular Line Observations in Two Dusty Star-Forming Galaxies at z = 6.9

Authors: Sreevani Jarugula, Joaquin D. Vieira, Axel Weiß, Justin S. Spilker, Manuel Aravena, Melanie Archipley, Matthieu Béthermin, Scott C. Chapman, Chenxing Dong, Thomas R. Greve, Kevin Harrington, Christopher C. Hayward, Yashar Hezaveh, Ryley Hill, Katrina C. Litke, Matthew A. Malkan, Daniel P. Marrone, Desika Narayanan, Kedar A. Phadke, Cassie Reuter, Kaja M. Rotermund

Abstract: SPT0311-58 is the most massive infrared luminous system discovered so far during the Epoch of Reionization (EoR). In this paper, we present a detailed analysis of the molecular interstellar medium at z = 6.9, through high-resolution observations of the CO(6-5), CO(7-6), CO(10-9), [CI](2-1), and p-H2O(211-202) lines and dust continuum emission with the Atacama Large Millimeter/submillimeter Array (… ▽ More SPT0311-58 is the most massive infrared luminous system discovered so far during the Epoch of Reionization (EoR). In this paper, we present a detailed analysis of the molecular interstellar medium at z = 6.9, through high-resolution observations of the CO(6-5), CO(7-6), CO(10-9), [CI](2-1), and p-H2O(211-202) lines and dust continuum emission with the Atacama Large Millimeter/submillimeter Array (ALMA). The system consists of a pair of intensely star-forming gravitationally lensed galaxies (labelled West and East). The intrinsic far-infrared luminosity is (16 $\pm$ 4)$\times\rm 10^{12} \ \rm L_{\odot}$ in West and (27 $\pm$ 4)$\times\rm 10^{11} \ \rm L_{\odot}$ in East. We model the dust, CO, and [CI] using non-local thermodynamic equilibrium radiative transfer models and estimate the intrinsic gas mass to be (5.4 $\pm$ 3.4)$\times\rm 10^{11} \ \rm M_{\odot}$ in West and (3.1 $\pm$ 2.7)$\times\rm 10^{10} \ \rm M_{\odot}$ in East. We find that the CO spectral line energy distribution in West and East are typical of high-redshift sub-millimeter galaxies (SMGs). The CO-to-H2 conversion factor ($α_{CO}$) and the gas depletion time scales estimated from the model are consistent with the high-redshift SMGs in the literature within the uncertainties. We find no evidence of evolution of depletion time with redshift in SMGs at z > 3. This is the most detailed study of molecular gas content of a galaxy in the EoR to-date, with the most distant detection of H2O in a galaxy without any evidence for active galactic nuclei in the literature. △ Less

Submitted 25 August, 2021; originally announced August 2021.

Comments: Accepted for publication in The Astrophysical Journal. 34 pages, 17 figures

ALMA Observations of the Sub-kpc Structure of the Host Galaxy of a z= 6.5 Lensed Quasar: A Rotationally-Supported Hyper-Starburst System at the Epoch of Reionization

Authors: Minghao Yue, Jinyi Yang, Xiaohui Fan, Feige Wang, Justin Spilker, Iskren Y. Georgiev, Charles R. Keeton, Katrina C. Litke, Daniel P. Marrone, Fabian Walter, Ran Wang, Xue-Bing Wu, Bram P. Venemans, Ann Zabludoff

Abstract: We report ALMA observations of the dust continuum and {\cii} emission of the host galaxy of J0439+1634, a gravitationally lensed quasar at $z=6.5$. Gravitational lensing boosts the source-plane resolution to $\sim0\farcs15$ $(\sim0.8\text{ kpc})$. The lensing model derived from the ALMA data is consistent with the fiducial model in \citet{fan19} based on {\it HST} imaging. The host galaxy of J0439… ▽ More We report ALMA observations of the dust continuum and {\cii} emission of the host galaxy of J0439+1634, a gravitationally lensed quasar at $z=6.5$. Gravitational lensing boosts the source-plane resolution to $\sim0\farcs15$ $(\sim0.8\text{ kpc})$. The lensing model derived from the ALMA data is consistent with the fiducial model in \citet{fan19} based on {\it HST} imaging. The host galaxy of J0439+1634 can be well-fitted by a Sérsic profile consistent with an exponential disk, both in the far-infrared (FIR) continuum and the {\cii} emission. The overall magnification is $4.53\pm0.05$ for the continuum and $3.44\pm0.05$ for the {\cii} line. The host galaxy of J0439+1634 is a compact ultra-luminous infrared galaxy, with a total star formation rate (SFR) of $1.56\times10^{3}M_\odot/\text{year}$ after correcting for lensing and an effective radius of $0.74$ kpc. The resolved regions in J0439+1634 follow the ``{\cii} deficit," where the {\cii}-to-FIR ratio decreases with FIR surface brightness. The reconstructed velocity field of J0439+1634 appears to be rotation-like. The maximum line-of-sight rotation velocity of 130 km/s at a radius of 2 kpc. However, our data cannot be fit by an axisymmetric thin rotating disk, and the inclination of the rotation axis, $i$, remains unconstrained. We estimate the dynamical mass of the host galaxy to be $7.9\sin^{-2}(i)\times10^{9}M_\odot$. J0439+1634 is likely to have a high gas-mass fraction and an oversized SMBH compared to local relations. The SFR of J0439+1634 reaches the maximum possible values, and the SFR surface density is close to the highest value seen in any star-forming galaxy currently known in the universe. △ Less

Submitted 15 June, 2021; v1 submitted 14 June, 2021; originally announced June 2021.

Comments: 14 pages, 7 figures. Accepted by ApJ

Comprehensive Gas Characterization of a $z= 2.5$ Protocluster: A Cluster Core Caught in the Beginning of Virialization?

Authors: Jaclyn B. Champagne, Caitlin M. Casey, Jorge A. Zavala, Asantha Cooray, Helmut Dannerbauer, Andrew Fabian, Christopher C. Hayward, Arianna S. Long, Justin S. Spilker

Abstract: In order to connect galaxy clusters to their progenitor protoclusters, we must constrain the star formation histories within their member galaxies and the timescale of virial collapse. In this paper we characterize the complex star-forming properties of a $z=2.5$ protocluster in the COSMOS field using ALMA dust continuum and new VLA CO(1-0) observations of two filaments associated with the structu… ▽ More In order to connect galaxy clusters to their progenitor protoclusters, we must constrain the star formation histories within their member galaxies and the timescale of virial collapse. In this paper we characterize the complex star-forming properties of a $z=2.5$ protocluster in the COSMOS field using ALMA dust continuum and new VLA CO(1-0) observations of two filaments associated with the structure, sometimes referred to as the "Hyperion" protocluster. We focus in particular on the protocluster "core" which has previously been suggested as the highest redshift bona fide galaxy cluster traced by extended X-ray emission in a stacked Chandra/XMM image. We re-analyze this data and refute these claims, finding that at least 40 $\pm$ 17% of extended X-ray sources of similar luminosity and size at this redshift arise instead from Inverse Compton scattering off recently extinguished radio galaxies rather than intracluster medium. Using ancillary COSMOS data, we also constrain the SEDs of the two filaments' eight constituent galaxies from the rest-frame UV to radio. We do not find evidence for enhanced star formation efficiency in the core and conclude that the constituent galaxies are already massive (M$_{\star} \approx 10^{11} M_{\odot}$), with molecular gas reservoirs $>10^{10} M_{\odot}$ that will be depleted within 200-400 Myr. Finally, we calculate the halo mass of the nested core at $z=2.5$ and conclude that it will collapse into a cluster of 2-9 $\times 10^{14} M_{\odot}$, comparable to the size of the Coma cluster at $z=0$ and accounting for at least 50% of the total estimated halo mass of the extended "Hyperion" structure. △ Less

Submitted 5 April, 2021; originally announced April 2021.

Comments: 30 pages, 12 figures, 1 appendix. Accepted for publication in ApJ

Early Science with the Large Millimeter Telescope: Constraining the Gas Fraction of a Compact Quiescent Galaxy at z=1.883

Authors: Joyce N. Caliendo, Katherine E. Whitaker, Mohammad Akhshik, Grant Wilson, Christina C. Williams, Justin S. Spilker, Guillaume Mahler, Alexandra Pope, Keren Sharon, Emmaly Aguilar, Rachel Bezanson, Miguel Chavez Dagastino, Arturo I. Gómez-Ruiz, Alfredo Montaña, Sune Toft, Miguel Velazquez De La Rosa, Milagros Zeballos

Abstract: We present constraints on the dust continuum flux and inferred gas content of a gravitationally lensed massive quiescent galaxy at $z$=1.883$\pm$0.001 using AzTEC 1.1mm imaging with the Large Millimeter Telescope. MRG-S0851 appears to be a prototypical massive compact quiescent galaxy, but has evidence that it experienced a centrally concentrated rejuvenation event in the last 100 Myr (see Akhshik… ▽ More We present constraints on the dust continuum flux and inferred gas content of a gravitationally lensed massive quiescent galaxy at $z$=1.883$\pm$0.001 using AzTEC 1.1mm imaging with the Large Millimeter Telescope. MRG-S0851 appears to be a prototypical massive compact quiescent galaxy, but has evidence that it experienced a centrally concentrated rejuvenation event in the last 100 Myr (see Akhshik et al. 2020). This galaxy is undetected in the AzTEC image but we calculate an upper limit on the millimeter flux and use this to estimate the H$_2$ mass limit via an empirically calibrated relation that assumes a constant molecular gas-to-dust ratio of 150. We constrain the 3$σ$ upper limit of the H$_2$ fraction from the dust continuum in MRG-S0851 to be ${M_{H_2}/M_{\star}}$ $\leq$ 6.8%. MRG-S0851 has a low gas fraction limit with a moderately low sSFR owing to the recent rejuvenation episode, which together results in a relatively short depletion time of $<$0.6 Gyr if no further H$_2$ gas is accreted. Empirical and analytical models both predict that we should have detected molecular gas in MRG-S0851, especially given the rejuvenation episode; this suggests that cold gas and/or dust is rapidly depleted in at least some early quiescent galaxies. △ Less

Submitted 4 March, 2021; originally announced March 2021.

Comments: 8 pages, 4 figures, accepted for publication in Astrophysical Journal Letters in January 2021 (in press)

The Physical Drivers of the Luminosity-Weighted Dust Temperatures in High-Redshift Galaxies

Authors: Anne D. Burnham, Caitlin M. Casey, Jorge A. Zavala, Sinclaire M. Manning, Justin S. Spilker, Scott C. Chapman, Chian-Chou Chen, Asantha Cooray, David B. Sanders, Nick Z. Scoville

Abstract: The underlying distribution of galaxies' dust SEDs (i.e., their spectra re-radiated by dust from rest-frame $\sim$3$μ$m-3mm) remains relatively unconstrained due to a dearth of FIR/(sub)mm data for large samples of galaxies. It has been claimed in the literature that a galaxy's dust temperature -- observed as the wavelength where the dust SED peaks ($λ_{peak}$) -- is traced most closely by its spe… ▽ More The underlying distribution of galaxies' dust SEDs (i.e., their spectra re-radiated by dust from rest-frame $\sim$3$μ$m-3mm) remains relatively unconstrained due to a dearth of FIR/(sub)mm data for large samples of galaxies. It has been claimed in the literature that a galaxy's dust temperature -- observed as the wavelength where the dust SED peaks ($λ_{peak}$) -- is traced most closely by its specific star-formation rate (sSFR) or parameterized 'distance' to the SFR-M$_\star$ relation (the galaxy 'main sequence'). We present 0.24" resolved 870$μ$m ALMA dust continuum observations of seven $z=1.4-4.6$ dusty star-forming galaxies (DSFGs) chosen to have a large range of well-constrained luminosity-weighted dust temperatures. We also draw on similar resolution dust continuum maps from a sample of ALESS submillimeter galaxies from Hodge et al. (2016). We constrain the physical scales over which the dust radiates and compare those measurements to characteristics of the integrated SED. We confirm significant correlations of $λ_{peak}$ with both L$_{IR}$ (or SFR) and $Σ_{\rm IR}$ ($\propto$SFR surface density). We investigate the correlation between $\log_{10}$($λ_{peak}$) and $\log_{10}$($Σ_{\rm IR}$) and find the relation to hold as would be expected from the Stefan-Boltzmann Law, or the effective size of an equivalent blackbody. The correlations of $λ_{peak}$ with sSFR and distance from the SFR-M$_\star$ relation are less significant than those for $Σ_{\rm IR}$ or L$_{IR}$; therefore, we conclude that the more fundamental tracer of galaxies' luminosity-weighted integrated dust temperatures are indeed their star-formation surface densities in line with local Universe results, which relate closely to the underlying geometry of dust in the ISM. △ Less

Submitted 11 February, 2021; originally announced February 2021.

Comments: 18 pages, 7 figures, accepted for publication in ApJ