Showing 1–50 of 106 results for author: Popping, G

The ALMA-ALPAKA survey II. Evolution of turbulence in galaxy disks across cosmic time: difference between cold and warm gas

Authors: F. Rizzo, C. Bacchini, M. Kohandel, L. Di Mascolo, F. Fraternali, F. Roman-Oliveira, A. Zanella, G. Popping, F. Valentino, G. Magdis, K. Whitaker

Abstract: The gas in the interstellar medium (ISM) of galaxies is supersonically turbulent. Measurements of turbulence typically rely on cold gas emission lines for low-z galaxies and warm ionized gas observations for z>0 galaxies. Studies of warm gas kinematics at z>0 conclude that the turbulence strongly evolves as a function of redshift, due to the increasing impact of gas accretion and mergers in the ea… ▽ More The gas in the interstellar medium (ISM) of galaxies is supersonically turbulent. Measurements of turbulence typically rely on cold gas emission lines for low-z galaxies and warm ionized gas observations for z>0 galaxies. Studies of warm gas kinematics at z>0 conclude that the turbulence strongly evolves as a function of redshift, due to the increasing impact of gas accretion and mergers in the early Universe. However, recent findings suggest potential biases in turbulence measurements derived from ionized gas at high-z, impacting our understanding of turbulence origin, ISM physics and disk formation. We investigate the evolution of turbulence using velocity dispersion ($σ$) measurements from cold gas tracers (i.e., CO, [CI], [CII]) derived from a sample of 57 galaxy disks spanning the redshift range z=0-5. This sample consists of main-sequence and starburst galaxies with stellar masses $\gtrsim 10^{10} M_{\odot}$. The comparison with current H$α$ kinematic observations and existing models demonstrates that the velocity dispersion inferred from cold gas tracers differ by a factor of $\approx 3$ from those obtained using emission lines tracing warm gas. We show that stellar feedback is the main driver of turbulence measured from cold gas tracers. This is fundamentally different from the conclusions of studies based on warm gas, which had to consider additional turbulence drivers to explain the high values of $σ$. We present a model predicting the redshift evolution of turbulence in galaxy disks, attributing the increase of $σ$ with redshift to the higher energy injected by supernovae due to the elevated star-formation rate in high-z galaxies. This supernova-driven model suggests that turbulence is lower in galaxies with lower stellar mass compared to those with higher stellar mass. Additionally, it forecasts the evolution of $σ$ in Milky-Way like progenitors. △ Less

Submitted 13 July, 2024; v1 submitted 8 July, 2024; originally announced July 2024.

Comments: Accepted for publication in A&A. The abstract has been modified to comply with arXiv's character limit

ALMA reveals a dust-obscured galaxy merger at Cosmic Noon

Authors: I. Langan, G. Popping, M. Ginolfi, F. Gentile, F. Valentino, M. Kaasinen

Abstract: Galaxy mergers play a critical role in galaxy evolution - altering the size, morphology, dynamics and composition of galaxies. So far, galaxy mergers have mostly been identified through visual inspection of their rest-frame optical and NIR emission. But, dust can obscure this emission, resulting in the misclassification of mergers as single galaxies, and the incorrect interpretation of their baryo… ▽ More Galaxy mergers play a critical role in galaxy evolution - altering the size, morphology, dynamics and composition of galaxies. So far, galaxy mergers have mostly been identified through visual inspection of their rest-frame optical and NIR emission. But, dust can obscure this emission, resulting in the misclassification of mergers as single galaxies, and the incorrect interpretation of their baryonic properties. Having serendipitously discovered a dust-obscured galaxy merger at z = 1.17, we aim to determine the baryonic properties of the two merging galaxies, including the star formation rate, and stellar, molecular gas, and dust masses. Using Band 3 and 6 observations from the Atacama Large Millimeter and submillimeter Array (ALMA), and ancillary data, we study the morphology of this previously misclassified merger. We deblend the emission, derive the gas masses from CO observations, and model the spectral energy distributions, to determine the properties of each galaxy. Using the rare combination of ALMA CO(2-1), CO(5-4) and dust-continuum (rest-frame 520um) observations, we provide insights into the gas and dust content and ISM properties of each merger component. We find that only one of the two galaxies is highly dust-obscured, whereas both are massive (> 10^10.5 Msun), highly star-forming (SFR = 60-900Msun/yr), have a moderate-to-low depletion time (tdepl < 0.7Gyr) and high gas fraction ( fgas >= 1). These properties can be interpreted as the positive impact of the merger. With this serendipitous discovery, we highlight the power of (sub)millimeter observations to identify and characterise the individual components of obscured galaxy mergers. △ Less

Submitted 3 June, 2024; originally announced June 2024.

Comments: 10 pages, 8 figures, accepted for publication to A&A

Observational predictions for the survival of atomic hydrogen in simulated Fornax-like galaxy clusters

Authors: Avinash Chaturvedi, Stephanie Tonnesen, Greg L. Bryan, Gergö Popping, Michael Hilker, Paolo Serra, Shy Genel

Abstract: The presence of dense, neutral hydrogen clouds in the hot, diffuse intra-group and intra-cluster medium is an important clue to the physical processes controlling the survival of cold gas and sheds light on cosmological baryon flows in massive halos. Advances in numerical modeling and observational surveys means that theory and observational comparisons are now possible. In this paper, we use the… ▽ More The presence of dense, neutral hydrogen clouds in the hot, diffuse intra-group and intra-cluster medium is an important clue to the physical processes controlling the survival of cold gas and sheds light on cosmological baryon flows in massive halos. Advances in numerical modeling and observational surveys means that theory and observational comparisons are now possible. In this paper, we use the high-resolution TNG50 cosmological simulation to study the HI distribution in seven halos with masses similar to the Fornax galaxy cluster. Adopting observational sensitivities similar to the MeerKAT Fornax Survey (MFS), an ongoing HI survey that will probe to column densities of $10^{18}$ cm$^{-2}$, we find that Fornax-like TNG50 halos have an extended distribution of neutral hydrogen clouds. Within one virial radius, we predict the MFS will observe a total HI covering fraction around $\sim$ 12\% (mean value) for 10 kpc pixels and 6\% for 2 kpc pixels. If we restrict this to gas more than 10 half-mass radii from galaxies, the mean values only decrease mildly, to 10\% (4\%) for 10 (2) kpc pixels (albeit with significant halo-to-halo spread). Although there are large amounts of HI outside of galaxies, the gas seems to be associated with satellites, judging both by the visual inspection of projections and by comparison of the line of sight velocities of galaxies and intracluster HI. △ Less

Submitted 25 April, 2024; originally announced April 2024.

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

Atacama Large Aperture Submillimeter Telescope (AtLAST) Science: The hidden circumgalactic medium

Authors: Minju M. Lee, Alice Schimek, Claudia Cicone, Paola Andreani, Gergö Popping, Laura Sommovigo, Philip N. Appleton, Manuela Bischetti, Sebastiano Cantalupo, Chian-Chou Chen, Helmut Dannerbauer, Carlos De Breuck, Luca Di Mascolo, Bjorn H. C. Emonts, Evanthia Hatziminaoglou, Antonio Pensabene, Francesca Rizzo, Matus Rybak, Sijing Shen, Andreas Lundgren, Mark Booth, Pamela Klaassen, Tony Mroczkowski, Martin A. Cordiner, Doug Johnstone , et al. (7 additional authors not shown)

Abstract: Our knowledge of galaxy formation and evolution has incredibly progressed through multi-wavelength observational constraints of the interstellar medium (ISM) of galaxies at all cosmic epochs. However, little is known about the physical properties of the more diffuse and lower surface brightness reservoir of gas and dust that extends beyond ISM scales and fills dark matter haloes of galaxies up to… ▽ More Our knowledge of galaxy formation and evolution has incredibly progressed through multi-wavelength observational constraints of the interstellar medium (ISM) of galaxies at all cosmic epochs. However, little is known about the physical properties of the more diffuse and lower surface brightness reservoir of gas and dust that extends beyond ISM scales and fills dark matter haloes of galaxies up to their virial radii, the circumgalactic medium (CGM). New theoretical studies increasingly stress the relevance of the latter for understanding the feedback and feeding mechanisms that shape galaxies across cosmic times, whose cumulative effects leave clear imprints into the CGM. Recent studies are showing that a -- so far unconstrained -- fraction of the CGM mass may reside in the cold (T < 1e4 K) molecular and atomic phase, especially in high-redshift dense environments. These gas phases, together with the warmer ionised phase, can be studied in galaxies from z ~ 0 to z ~ 10 through bright far-infrared and sub-millimeter emission lines such as [C II] 158$μ$m, [O III] 88 $μ$m, [C I] 609$μ$m, [C I] 370$μ$m, and the rotational transitions of CO. Imaging such hidden cold CGM can lead to a breakthrough in galaxy evolution studies but requires a new facility with the specifications of the proposed Atacama Large Aperture Submillimeter Telescope (AtLAST). In this paper, we use theoretical and empirical arguments to motivate future ambitious CGM observations with AtLAST and describe the technical requirements needed for the telescope and its instrumentation to perform such science. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: Submitted to Open Research Europe as part of the AtLAST collection: https://open-research-europe.ec.europa.eu/collections/atlast/about

The ALPINE-ALMA [C II] survey: Characterisation of Spatial Offsets in Main-Sequence Galaxies at $z \sim$ 4-6

Authors: Meghana Killi, Michele Ginolfi, Gergö Popping, Darach Watson, Giovanni Zamorani, Brian C. Lemaux, Seiji Fujimoto, Andreas Faisst, Matthieu Bethermin, Michael Romano, Yoshinobu Fudamoto, Sandro Bardelli, Médéric Boquien, Stefano Carniani, Miroslava Dessauges-Zavadsky, Carlotta Gruppioni, Nimish Hathi, Eduardo Ibar, Gareth C. Jones, Anton M. Koekemoer, Ivanna Langan, Hugo Méndez-Hernández, Yuma Sugahara, Livia Vallini, Daniela Vergani

Abstract: Galaxy morphology is shaped by stellar activity, feedback, gas and dust properties, and interactions with surroundings, and can therefore provide insight into these processes. In this paper, we study the spatial offsets between stellar and interstellar medium emission in a sample of 54 main-sequence star-forming galaxies at $z\sim4-6$ observed with the Atacama Large Millimeter/submillimeter Array… ▽ More Galaxy morphology is shaped by stellar activity, feedback, gas and dust properties, and interactions with surroundings, and can therefore provide insight into these processes. In this paper, we study the spatial offsets between stellar and interstellar medium emission in a sample of 54 main-sequence star-forming galaxies at $z\sim4-6$ observed with the Atacama Large Millimeter/submillimeter Array (ALMA) and drawn from the ALMA Large Program to INvestigate C$^+$ at Early times (ALPINE). We find no significant spatial offset for the majority ($\sim$ 70 percent) of galaxies in the sample among any combination of [C II], far-infrared continuum, optical, and ultraviolet emission. However, a fraction of the sample ($\sim$ 30 percent) shows offsets larger than the median by more than 3$σ$ significance (compared to the uncertainty on the offsets), especially between [C II] and ultraviolet emission. We find that these significant offsets are of the order of $\sim$0.5-0.7 arcsec, corresponding to $\sim$3.5-4.5 kiloparsecs. The offsets could be caused by a complex dust geometry, strong feedback from stars and active galactic nuclei, large-scale gas inflow and outflow, or a combination of these phenomena. However, our current analysis does not definitively constrain the origin. Future, higher resolution ALMA and JWST observations may help resolve the ambiguity. Regardless, since there exist at least some galaxies that display such large offsets, galaxy models and spectral energy distribution fitting codes cannot assume co-spatial emission in all main-sequence galaxies, and must take into account that the observed emission across wavelengths may be spatially segregated. △ Less

Submitted 12 February, 2024; originally announced February 2024.

Comments: 20 pages, 8 figures. Submitted to MNRAS; second revision

Carbon envelopes around merging galaxies at z ~ 4.5

Authors: C. Di Cesare, M. Ginolfi, L. Graziani, R. Schneider, M. Romano, G. Popping

Abstract: Galaxies evolve through a dynamic exchange of material with their immediate surrounding environment, the circumgalactic medium (CGM). Understanding the physics of gas flows and the nature of the CGM is thus fundamental to studying galaxy evolution, especially at $4 \leq z \leq 6$ when galaxies rapidly assembled their masses and reached their chemical maturity. Galactic outflows are predicted to en… ▽ More Galaxies evolve through a dynamic exchange of material with their immediate surrounding environment, the circumgalactic medium (CGM). Understanding the physics of gas flows and the nature of the CGM is thus fundamental to studying galaxy evolution, especially at $4 \leq z \leq 6$ when galaxies rapidly assembled their masses and reached their chemical maturity. Galactic outflows are predicted to enrich the CGM with metals, although gas stripping in systems undergoing a major merger has also been suggested to play a role. In this work, we explore the metal enrichment of the medium around merging galaxies at $z\sim4.5$, observed by the ALMA Large Program to INvestigate [CII] at Early times (ALPINE) survey. To do so, we study the nature of the [CII]158 $μ$m emission in the CGM around these systems, using simulations to help disentangle the mechanisms contributing to the CGM metal pollution. By adopting an updated classification of major merger systems in the ALPINE survey, we select and analyse merging galaxies whose components can be spatially and/or spectrally resolved in a robust way. In this way, we can distinguish between the [CII] emission coming from the single components of the system and that coming from the system as a whole. We also make use of the dustyGadget cosmological simulation to select synthetic analogues of observed galaxies and guide the interpretation of the observational results. We find a large diffuse [CII] envelope (> 20 kpc) embedding all the merging systems, with around 50% of the total [CII] emission coming from the medium between the galaxies. Using predictions from dustyGadget we suggest that this emission has a two-fold nature: it is due to both dynamical interactions between the galaxies which result in tidal stripped gas and the presence of star-forming satellites (currently unresolved by ALMA) that enrich the medium with heavy elements. △ Less

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

Comments: 20 pages, 14 figure, 2 tables. Few updated references, comments welcome

$\texttt{slick}$: Modeling a Universe of Molecular Line Luminosities in Hydrodynamical Simulations

Authors: Karolina Garcia, Desika Narayanan, Gergö Popping, R. Anirudh, Sagan Sutherland, Melanie Kaasinen

Abstract: We present {\sc slick} (the Scalable Line Intensity Computation Kit), a software package that calculates realistic CO, [\ion{C}{1}], and [\ion{C}{2}] luminosities for clouds and galaxies formed in hydrodynamic simulations. Built on the radiative transfer code {\sc despotic}, {\sc slick} computes the thermal, radiative, and statistical equilibrium in concentric zones of model clouds, based on their… ▽ More We present {\sc slick} (the Scalable Line Intensity Computation Kit), a software package that calculates realistic CO, [\ion{C}{1}], and [\ion{C}{2}] luminosities for clouds and galaxies formed in hydrodynamic simulations. Built on the radiative transfer code {\sc despotic}, {\sc slick} computes the thermal, radiative, and statistical equilibrium in concentric zones of model clouds, based on their physical properties and individual environments. We validate our results applying {\sc slick} to the high-resolution run of the {\sc Simba} simulations, testing the derived luminosities against empirical and theoretical/analytic relations. To simulate the line emission from a universe of emitting clouds, we have incorporated random forest machine learning (ML) methods into our approach, allowing us to predict cosmologically evolving properties of CO, [\ion{C}{1}] and [\ion{C}{2}] emission from galaxies such as luminosity functions. We tested this model in 100,000 gas particles, and 2,500 galaxies, reaching an average accuracy of $\sim$99.8\% for all lines. Finally, we present the first model light cones created with realistic and ML-predicted CO, [\ion{C}{1}], and [\ion{C}{2}] luminosities in cosmological hydrodynamical simulations, from $z=0$ to $z=10$. △ Less

Submitted 2 November, 2023; originally announced November 2023.

Comments: 17 pages, 11 figures, comments are welcome

Theoretical strong line metallicity diagnostics for the JWST era

Authors: Prerak Garg, Desika Narayanan, Ryan L. Sanders, Romeel Davè, Gergö Popping, Alice E. Shapley, Daniel P. Stark, Jonathan R. Trump

Abstract: The ratios of strong rest-frame optical emission lines are the dominant indicator of metallicities in high-redshift galaxies. Since typical strong-line based metallicity indicators are calibrated on auroral lines at $z=0$, their applicability for galaxies in the distant Universe is unclear. In this paper, we make use of mock emission line data from cosmological simulations to investigate the calib… ▽ More The ratios of strong rest-frame optical emission lines are the dominant indicator of metallicities in high-redshift galaxies. Since typical strong-line based metallicity indicators are calibrated on auroral lines at $z=0$, their applicability for galaxies in the distant Universe is unclear. In this paper, we make use of mock emission line data from cosmological simulations to investigate the calibration of rest-frame optical emission lines as metallicity indicators at high redshift. Our model, which couples the SIMBA cosmological galaxy formation simulation with cloudy photoionization calculations, includes contributions from HII regions, post-AGB stars and Diffuse Ionized Gas (DIG). We find mild redshift evolution in the 12 indicators that we study, which implies that the dominant physical properties that evolve in our simulations do have a discernible impact on the metallicity calibrations at high redshifts. When comparing our calibrations with high redshift auroral line observations from James Webb Space Telescope we find a slight offset between our model results and the observations and find that a higher ionization parameter at high redshifts can be one of the possible explanations. We explore the physics that drives the shapes of strong-line metallicity relationships and propose calibrations for hitherto unexplored low-metallicity regimes. Finally, we study the contribution of DIG to total line fluxes. We find that the contribution of DIG increases with metallicity at z $\sim$ 0 for singly ionized oxygen and sulfur lines and can be as high as 70% making it crucial to include their contribution when modeling nebular emission. △ Less

Submitted 8 November, 2023; v1 submitted 12 October, 2023; originally announced October 2023.

Comments: 28 pages, 6 figures, typo corrected: in polynomial fits, x represents metallicity and y represents the line ratio

XLSSC 122 caught in the act of growing up: Spatially resolved SZ observations of a z=1.98 galaxy cluster

Authors: J. van Marrewijk, L. Di Mascolo, A. S. Gill, N. Battaglia, E. S. Battistelli, J. R. Bond, M. J. Devlin, P. Doze, J. Dunkley, K. Knowles, A. Hincks, J. P. Hughes, M. Hilton, K. Moodley, T. Mroczkowski, S. Naess, B. Partridge, G. Popping, C. Sifón, S. T. Staggs, E. J. Wollack

Abstract: How protoclusters evolved from sparse galaxy overdensities to mature galaxy clusters is still not well understood. In this context, detecting and characterizing the hot ICM at high redshifts (z~2) is key to understanding how the continuous accretion from and mergers along the filamentary large-scale structure impact the first phases of cluster formation. We study the dynamical state and morphology… ▽ More How protoclusters evolved from sparse galaxy overdensities to mature galaxy clusters is still not well understood. In this context, detecting and characterizing the hot ICM at high redshifts (z~2) is key to understanding how the continuous accretion from and mergers along the filamentary large-scale structure impact the first phases of cluster formation. We study the dynamical state and morphology of the z=1.98 galaxy cluster XLSSC 122 with high-resolution observations (~5") of the ICM through the SZ effect. Via Bayesian forward modeling, we map the ICM on scales from the virial radius down to the core of the cluster. To constrain such a broad range of spatial scales, we employ a new technique that jointly forward-models parametric descriptions of the pressure distribution to interferometric ACA and ALMA observations and multi-band imaging data from the 6-m, single-dish Atacama Cosmology Telescope. We detect the SZ effect with $11σ$ in the ALMA+ACA observations and find a flattened inner pressure profile that is consistent with a non-cool core classification with a significance of $>3σ$. In contrast to the previous works, we find better agreement between the SZ effect signal and the X-ray emission as well as the cluster member distribution. Further, XLSSC 122 exhibits an excess of SZ flux in the south of the cluster where no X-ray emission is detected. By reconstructing the interferometric observations and modeling in the uv-plane, we obtain a tentative detection of an infalling group or filamentary-like structure that is believed to boost and heat up the ICM while the density of the gas is low. In addition, we provide an improved SZ mass of $M_{500,\mathrm{c}} = 1.66^{+0.23}_{-0.20} \times 10^{14} \rm M_\odot$. Altogether, the observations indicate that we see XLSSC 122 in a dynamic phase of cluster formation while a large reservoir of gas is already thermalized. △ Less

Submitted 19 June, 2024; v1 submitted 9 October, 2023; originally announced October 2023.

Galaxy mergers in Subaru HSC-SSP: a deep representation learning approach for identification and the role of environment on merger incidence

Authors: Kiyoaki Christopher Omori, Connor Bottrell, Mike Walmsley, Hassen M. Yesuf, Andy D. Goulding, Xuheng Ding, Gergö Popping, John D. Silverman, Tsutomu T. Takeuchi, Yoshiki Toba

Abstract: We take a deep learning-based approach for galaxy merger identification in Subaru HSC-SSP, specifically through the use of deep representation learning and fine-tuning, with the aim of creating a pure and complete merger sample within the HSC-SSP survey. We can use this merger sample to conduct studies on how mergers affect galaxy evolution. We use Zoobot, a deep learning representation learning m… ▽ More We take a deep learning-based approach for galaxy merger identification in Subaru HSC-SSP, specifically through the use of deep representation learning and fine-tuning, with the aim of creating a pure and complete merger sample within the HSC-SSP survey. We can use this merger sample to conduct studies on how mergers affect galaxy evolution. We use Zoobot, a deep learning representation learning model pre-trained on citizen science votes on Galaxy Zoo DeCALS images. We fine-tune Zoobot for the purpose of merger classification of images of SDSS and GAMA galaxies in HSC-SSP PDR 3. Fine-tuning is done using 1200 synthetic HSC-SSP images of galaxies from the TNG simulation. We then find merger probabilities on observed HSC images using the fine-tuned model. Using our merger probabilities, we examine the relationship between merger activity and environment. We find that our fine-tuned model returns an accuracy on the synthetic validation data of 76%. This number is comparable to those of previous studies where convolutional neural networks were trained with simulation images, but with our work requiring a far smaller number of training samples. For our synthetic data, our model is able to achieve completeness and precision values of 80%. In addition, our model is able to correctly classify both mergers and non-mergers of diverse morphologies and structures, including those at various stages and mass ratios, while distinguishing between projections and merger pairs. For the relation between galaxy mergers and environment, we find two distinct trends. Using stellar mass overdensity estimates for TNG simulations and observations using SDSS and GAMA, we find that galaxies with higher merger scores favor lower density environments on scales of 0.5 to 8 h^-1 Mpc. However, below these scales in the simulations, we find that galaxies with higher merger scores favor higher density environments. △ Less

Submitted 27 September, 2023; originally announced September 2023.

Comments: 36 pages, 15 figures, accepted to Astronomy and Astrophysics

Faint [CI](1-0) emission in z $\sim$ 3.5 radio galaxies

Authors: S. Kolwa, C. De Breuck, J. Vernet, D. Wylezalek, W. Wang, G. Popping, A. W. S. Man, C. M. Harrison, P. Andreani

Abstract: We present Atacama Large Millimeter/sub-millimeter Array (ALMA) neutral carbon, [C I](1-0), line observations that probe molecular hydrogen gas (H$_2$) within seven radio galaxies at $z = 2.9 - 4.5$ surrounded by extended ($\gtrsim100$ kpc) Ly-$α$ nebulae. We extract [C I](1-0) emission from the radio-active galactic nuclei (AGN) host galaxies whose positions are set by near-infrared detections an… ▽ More We present Atacama Large Millimeter/sub-millimeter Array (ALMA) neutral carbon, [C I](1-0), line observations that probe molecular hydrogen gas (H$_2$) within seven radio galaxies at $z = 2.9 - 4.5$ surrounded by extended ($\gtrsim100$ kpc) Ly-$α$ nebulae. We extract [C I](1-0) emission from the radio-active galactic nuclei (AGN) host galaxies whose positions are set by near-infrared detections and radio detections of the cores. Additionally, we place constraints on the galaxies' systemic redshifts via He II $λ$1640 lines seen with the Multi-Unit Spectroscopic Explorer (MUSE). We detect faint [C I] emission in four out of seven sources. In two of these galaxies, we discover narrow line emission of full width at half maximum $\lesssim100$ km s$^{-1}$ which may trace emission from bright kpc-scale gas clouds within the ISM. In the other two [C I]-detected galaxies, line dispersions range from $\sim100 - 600$ km s$^{-1}$ and may be tracing the rotational component of the cold gas. Overall, the [C I] line luminosities correspond to H$_2$ masses of M$_{\rm H_2,[C I]} \simeq (0.5 - 3) \times 10^{10} M_\odot$ for the detections and M$_{H_2,[C I]} < 0.65 \times 10^{10} M_\odot$ for the [C I] non-detections in three out of seven galaxies within the sample. The molecular gas masses in our sample are relatively low in comparison to previously reported measures for similar galaxies which are M$_{H_2,[C I]} \simeq (3 - 4) \times 10^{10}.$ Our results imply that the observed faintness in carbon emission is representative of a decline in molecular gas supply from previous star-formation epochs and/or a displacement of molecular gas from the ISM due to jet-powered outflows. △ Less

Submitted 1 September, 2023; originally announced September 2023.

Comments: 16 pages, 4 figures and 5 tables. Accepted for publication in MNRAS

IllustrisTNG in the HSC-SSP: image data release and the major role of mini mergers as drivers of asymmetry and star formation

Authors: Connor Bottrell, Hassen M. Yesuf, Gergö Popping, Kiyoaki Christopher Omori, Shenli Tang, Xuheng Ding, Annalisa Pillepich, Dylan Nelson, Lukas Eisert, Hua Gao, Andy D. Goulding, Boris S. Kalita, Wentao Luo, Jenny E. Greene, Jingjing Shi, John D. Silverman

Abstract: At fixed galaxy stellar mass, there is a clear observational connection between structural asymmetry and offset from the star forming main sequence, $Δ$SFMS. Herein, we use the TNG50 simulation to investigate the relative roles of major mergers (stellar mass ratios $μ\geq0.25$), minor ($0.1 \leq μ< 0.25$), and mini mergers ($0.01 \leq μ< 0.1$) in driving this connection amongst star forming galaxi… ▽ More At fixed galaxy stellar mass, there is a clear observational connection between structural asymmetry and offset from the star forming main sequence, $Δ$SFMS. Herein, we use the TNG50 simulation to investigate the relative roles of major mergers (stellar mass ratios $μ\geq0.25$), minor ($0.1 \leq μ< 0.25$), and mini mergers ($0.01 \leq μ< 0.1$) in driving this connection amongst star forming galaxies (SFGs). We use dust radiative transfer post-processing with SKIRT to make a large, public collection of synthetic Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) images of simulated TNG galaxies over $0.1\leq z \leq 0.7$ with $\log M_{\star} / \mathrm{M}_{\odot}\geq9$ ($\sim750$k images). Using their instantaneous SFRs, known merger histories/forecasts, and HSC-SSP asymmetries, we show (1) that TNG50 SFGs qualitatively reproduce the observed trend between $Δ$SFMS and asymmetry and (2) a strikingly similar trend emerges between $Δ$SFMS and the time-to-coalescence for mini mergers. Controlling for redshift, stellar mass, environment, and gas fraction, we show that individual mini merger events yield small enhancements in SFRs and asymmetries that are sustained on long timescales (at least $\sim3$ Gyr after coalescence, on average) -- in contrast to major/minor merger remnants which peak at much greater amplitudes but are consistent with controls only $\sim1$ Gyr after coalescence. Integrating the boosts in SFRs and asymmetries driven by $μ\geq0.01$ mergers since $z=0.7$ in TNG50 SFGs, we show that mini mergers are responsible for (i) $55$ per cent of all merger-driven star formation and (ii) $70$ per cent of merger-driven asymmetric structure. Due to their relative frequency and prolonged boost timescales, mini mergers dominate over their minor and major counterparts in driving star formation and asymmetry in SFGs. △ Less

Submitted 7 October, 2023; v1 submitted 28 August, 2023; originally announced August 2023.

Comments: 31 pages; 19 figures; MNRAS accepted; Image data available via the TNG website: www.tng-project.org/bottrell23

Emulating Radiative Transfer with Artificial Neural Networks

Authors: Snigdaa S. Sethuram, Rachel K. Cochrane, Christopher C. Hayward, Viviana Acquaviva, Francisco Villaescusa-Navarro, Gergo Popping, John H. Wise

Abstract: Forward-modeling observables from galaxy simulations enables direct comparisons between theory and observations. To generate synthetic spectral energy distributions (SEDs) that include dust absorption, re-emission, and scattering, Monte Carlo radiative transfer is often used in post-processing on a galaxy-by-galaxy basis. However, this is computationally expensive, especially if one wants to make… ▽ More Forward-modeling observables from galaxy simulations enables direct comparisons between theory and observations. To generate synthetic spectral energy distributions (SEDs) that include dust absorption, re-emission, and scattering, Monte Carlo radiative transfer is often used in post-processing on a galaxy-by-galaxy basis. However, this is computationally expensive, especially if one wants to make predictions for suites of many cosmological simulations. To alleviate this computational burden, we have developed a radiative transfer emulator using an artificial neural network (ANN), ANNgelina, that can reliably predict SEDs of simulated galaxies using a small number of integrated properties of the simulated galaxies: star formation rate, stellar and dust masses, and mass-weighted metallicities of all star particles and of only star particles with age <10 Myr. Here, we present the methodology and quantify the accuracy of the predictions. We train the ANN on SEDs computed for galaxies from the IllustrisTNG project's TNG50 cosmological magnetohydrodynamical simulation. ANNgelina is able to predict the SEDs of TNG50 galaxies in the ultraviolet (UV) to millimetre regime with a typical median absolute error of ~7 per cent. The prediction error is the greatest in the UV, possibly due to the viewing-angle dependence being greatest in this wavelength regime. Our results demonstrate that our ANN-based emulator is a promising computationally inexpensive alternative for forward-modeling galaxy SEDs from cosmological simulations. △ Less

Submitted 25 August, 2023; originally announced August 2023.

The IR Compactness of Dusty Galaxies Set Star-formation and Dust Properties at z~0-2

Authors: Jed McKinney, Alexandra Pope, Allison Kirkpatrick, Lee Armus, Tanio Diaz-Santos, Carlos Gomez-Guijarro, Maximilien Franco, David Elbaz, Christopher C. Hayward, Hanae Inami, Gergo Popping, Mengyuan Xiao

Abstract: Surface densities of gas, dust and stars provide a window into the physics of star-formation that, until the advent of high-resolution far-infrared/sub-millimeter observations, has been historically difficult to assess amongst dusty galaxies. To study the link between infrared (IR) surface densities and dust properties, we leverage the Atacama Large Millimetre/Submillimetre Array (ALMA) archive to… ▽ More Surface densities of gas, dust and stars provide a window into the physics of star-formation that, until the advent of high-resolution far-infrared/sub-millimeter observations, has been historically difficult to assess amongst dusty galaxies. To study the link between infrared (IR) surface densities and dust properties, we leverage the Atacama Large Millimetre/Submillimetre Array (ALMA) archive to measure the extent of cold dust emission in 15 $z\sim2$ IR selected galaxies selected on the basis of having available mid-IR spectroscopy from Spitzer. We use the mid-IR spectra to constrain the relative balance between dust heating from star-formation and active galactic nuclei (AGN), and to measure emission from Polycylic Aromatic Hydrocarbons (PAHs) -- small dust grains that play a key role in the photoelectric heating of gas. In general, we find that dust-obscured star-formation at high IR surface densities exhibits similar properties at low- and high-redshift, namely: local luminous IR galaxies have comparable PAH luminosity to total dust mass ratios as high-$z$ galaxies, and star-formation at $z\sim0-2$ is more efficient at high IR surface densities despite the fact that our sample of high$-z$ galaxies are closer to the main-sequence than local luminous IR galaxies. High star-formation efficiencies are coincident with a decline in the PAH/IR luminosity ratio reminiscent of the deficit observed in far-infrared fine-structure lines. Changes in the gas and dust conditions arising from high star-formation surface densities might help drive the star-formation efficiency up. This could help explain high efficiencies needed to reconcile star-formation and gas volume densities in dusty galaxies at cosmic noon. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 18 pages, 10 figures, accepted to ApJ

ALMA High-Level Data Products: Submillimetre counterparts of SDSS quasars in the ALMA footprint

Authors: A. Wong, E. Hatziminaoglou, A. Borkar, G. Popping, I. Pérez-Fournon, F. Poidevin, F. Stoehr, H. Messias

Abstract: The Atacama Large Millimetre/submillimetre Array (ALMA) is the world's most advanced radio interferometric facility, producing science data with an average rate of about 1 TB per day. After a process of calibration, imaging and quality assurance, the scientific data are stored in the ALMA Science Archive (ASA), along with the corresponding raw data, making the ASA an invaluable resource for origin… ▽ More The Atacama Large Millimetre/submillimetre Array (ALMA) is the world's most advanced radio interferometric facility, producing science data with an average rate of about 1 TB per day. After a process of calibration, imaging and quality assurance, the scientific data are stored in the ALMA Science Archive (ASA), along with the corresponding raw data, making the ASA an invaluable resource for original astronomical research. Due to their complexity, each ALMA data set has the potential for scientific results that go well beyond the ideas behind the original proposal that led to each observation. For this reason, the European ALMA Regional Centre initiated the High-Level Data Products initiative to develop science-oriented data products derived from data sets publicly available in the ASA, that go beyond the formal ALMA deliverables. The first instance of this initiative is the creation of a catalogue of submillimetre (submm) detections of Sloan Digital Sky Survey (SDSS) quasars from the SDSS Data Release 14 that lie in the aggregate ALMA footprint observed since ALMA Cycle 0. The ALMA fluxes are extracted in an automatic fashion, using the ALMA Data Mining Toolkit. All extractions above a signal-to-noise cut of 3.5 are considered, they have been visually inspected and the reliable detections are presented in a catalogue of 376 entries, corresponding to 275 unique quasars. Interesting targets found in the process, i.e. lensed or jetted quasars as well as quasars with nearby submm counterparts are highlighted, to facilitate further studies or potential follow up observations. △ Less

Submitted 9 May, 2023; originally announced May 2023.

The ALMA-ALPAKA survey I: high-resolution CO and [CI] kinematics of star-forming galaxies at z = 0.5-3.5

Authors: F. Rizzo, F. Roman-Oliveira, F. Fraternali, D. Frickmann, F. Valentino, G. Brammer, A. Zanella, V. Kokorev, G. Popping, K. E. Whitaker, M. Kohandel, G. E. Magdis, L. Di Mascolo, R. Ikeda, S. Jin, S. Toft

Abstract: Spatially-resolved studies of the kinematics of galaxies provide crucial insights into their assembly and evolution, enabling to infer the properties of the dark matter halos, derive the impact of feedback on the ISM, characterize the outflow motions. To date, most of the kinematic studies at z=0.5-3.5 were obtained using emission lines tracing the warm, ionized gas. However, whether these provide… ▽ More Spatially-resolved studies of the kinematics of galaxies provide crucial insights into their assembly and evolution, enabling to infer the properties of the dark matter halos, derive the impact of feedback on the ISM, characterize the outflow motions. To date, most of the kinematic studies at z=0.5-3.5 were obtained using emission lines tracing the warm, ionized gas. However, whether these provide an exhaustive or only a partial view of the dynamics of galaxies and of the properties of the ISM is still debated. Complementary insights on the cold gas kinematics are therefore needed. We present ALPAKA, a project aimed at gathering high-resolution observations of CO and [CI] emission lines of star-forming galaxies at z=0.5-3.5 from the ALMA public archive. With 147 hours of total integration time, ALPAKA assembles ~0.25'' observations for 28 star-forming galaxies, the largest sample with spatially-resolved cold gas kinematics as traced by either CO or [CI] at z>0.5. By combining multi-wavelength ancillary data, we derive the stellar masses ($M_{\star}$) and star-formation rates (SFR) for our targets, finding values of $M_{\star}\gtrsim 10^{10}$ M$_{\odot}$ and SFR of 10-3000 M$_{\odot}$/yr. A large fraction of ALPAKA galaxies (19/28) lie in overdense regions (clusters, groups, and protoclusters). We exploit the ALMA data to infer their dynamical state and we find that 19/28 ALPAKA galaxies are rotating disks, 2 are interacting systems, while for the remaining 7 sources the classification is uncertain. The disks have velocity dispersion values that are typically larger in the innermost regions than in the outskirts, with a median value for the entire disk sample of 35$^{+11}_{-9}$ km/s. Despite the bias of our sample towards galaxies hosting very energetic mechanisms, the ALPAKA disks have high ratios of ordered-to-random motion ($V/σ$) with a median value of 9$^{+7}_{-2}$. △ Less

Submitted 13 August, 2023; v1 submitted 28 March, 2023; originally announced March 2023.

Comments: 36 pages, 23 figures, 5 tables; accepted for publication in A&A. The data and the outputs of the kinematic analysis will be made available at https://alpaka-survey.github.io/index.html

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

ALMA Lensing Cluster Survey: Deep 1.2 mm Number Counts and Infrared Luminosity Functions at $z\simeq1-8$

Authors: Seiji Fujimoto, Kotaro Kohno, Masami Ouchi, Masamune Oguri, Vasily Kokorev, Gabriel Brammer, Fengwu Sun, Jorge Gonzalez-Lopez, Franz E. Bauer, Gabriel B. Caminha, Bunyo Hatsukade, Johan Richard, Ian Smail, Akiyoshi Tsujita, Yoshihiro Ueda, Ryosuke Uematsu, Adi Zitrin, Dan Coe, Jean-Paul Kneib, Marc Postman, Keiichi Umetsu, Claudia del P. Lagos, Gergo Popping, Yiping Ao, Larry Bradley , et al. (18 additional authors not shown)

Abstract: We present a statistical study of 180 dust continuum sources identified in 33 massive cluster fields by the ALMA Lensing Cluster Survey (ALCS) over a total of 133 arcmin$^{2}$ area, homogeneously observed at 1.2 mm. ALCS enables us to detect extremely faint mm sources by lensing magnification, including near-infrared (NIR) dark objects showing no counterparts in existing {\it Hubble Space Telescop… ▽ More We present a statistical study of 180 dust continuum sources identified in 33 massive cluster fields by the ALMA Lensing Cluster Survey (ALCS) over a total of 133 arcmin$^{2}$ area, homogeneously observed at 1.2 mm. ALCS enables us to detect extremely faint mm sources by lensing magnification, including near-infrared (NIR) dark objects showing no counterparts in existing {\it Hubble Space Telescope} and {\it Spitzer} images. The dust continuum sources belong to a blind sample ($N=141$) with S/N $\gtrsim$ 5.0 (a purity of $>$ 0.99) or a secondary sample ($N=39$) with S/N= $4.0-5.0$ screened by priors. With the blind sample, we securely derive 1.2-mm number counts down to $\sim7$ $μ$Jy, and find that the total integrated 1.2mm flux is 20.7$^{+8.5}_{-6.5}$ Jy deg$^{-2}$, resolving $\simeq$ 80 % of the cosmic infrared background light. The resolved fraction varies by a factor of $0.6-1.1$ due to the completeness correction depending on the spatial size of the mm emission. We also derive infrared (IR) luminosity functions (LFs) at $z=0.6-7.5$ with the $1/V_{\rm max}$ method, finding the redshift evolution of IR LFs characterized by positive luminosity and negative density evolution. The total (=UV+IR) cosmic star-formation rate density (SFRD) at $z>4$ is estimated to be $161^{+25}_{-21}$ % of the established measurements, which were almost exclusively based on optical$-$NIR surveys. Although our general understanding of the cosmic SFRD is unlikely to change beyond a factor of 2, these results add to the weight of evidence for an additional ($\approx 60$ %) SFRD component contributed by the faint-mm population, including NIR dark objects. △ Less

Submitted 20 June, 2024; v1 submitted 2 March, 2023; originally announced March 2023.

Comments: ApJS in press, 34 pages, 8 Tables, and 18 Figures (+38 pages, 5 Tables, and 6 Figures in Appendix)

MusE GAs FLOw and Wind (MEGAFLOW) IX. The impact of gas flows on the relations between the mass, star formation rate and metallicity of galaxies

Authors: I. Langan, J. Zabl, N. F. Bouche, M. Ginolfi, G. Popping, I. Schroetter, M. Wendt, J. Schaye, L. Boogaard, J. Freundlich, J. Richard, J. Matthee, W. Mercier, T. Contini, Y. Guo, M. Cherrey

Abstract: We study the link between gas flow events and key galaxy scaling relations: the relations between star formation rate (SFR) and stellar mass (the main sequence, MS), gas metallicity and stellar mass (the mass-metallicity relation, MZR) and gas metallicity, stellar mass and SFR (the fundamental metallicity relation, FMR). Using all star-forming galaxies (SFGs) in the 22 MUSE fields of the MusE GAs… ▽ More We study the link between gas flow events and key galaxy scaling relations: the relations between star formation rate (SFR) and stellar mass (the main sequence, MS), gas metallicity and stellar mass (the mass-metallicity relation, MZR) and gas metallicity, stellar mass and SFR (the fundamental metallicity relation, FMR). Using all star-forming galaxies (SFGs) in the 22 MUSE fields of the MusE GAs FLOw and Wind (MEGAFLOW) survey, we derive the MS, MZR and FMR scaling relations for 385 SFGs with $M = 10^8 - 10^{11.5}$ $M_\odot$ at redshifts 0.35 < z < 0.85. Using the MUSE data and complementary X-Shooter spectra at 0.85 < z < 1.4, we determine the locations of 21 SFGs associated with inflowing or outflowing circumgalactic gas (i.e. with strong MgII absorption in background quasar spectra) relative to these scaling relations. Compared to a control sample of galaxies without gas flows (i.e., without MgII absorption within 70 kpc of the quasar), SFGs with inflow events (i.e., MgII absorption along the major axis) are preferentially located above the MS, while SFGs with ouflow events (i.e., MgII absorption along the minor axis) are preferentially more metal rich. Our observations support the scenario in which gas accretion increases the SFR while diluting the metal content and where circumgalactic outflows are found in more metal-rich galaxies. △ Less

Submitted 14 March, 2023; v1 submitted 25 January, 2023; originally announced January 2023.

Comments: 13 pages, 8 figures

A NOEMA molecular line scan of the Hubble Deep Field North: Improved constraints on the CO luminosity functions and cosmic density of molecular gas

Authors: Leindert A. Boogaard, Roberto Decarli, Fabian Walter, Axel Weiss, Gergö Popping, Roberto Neri, Manuel Aravena, Dominik Riechers, Richard S. Ellis, Chris Carilli, Pierre Cox, Jérôme Pety

Abstract: We present measurements of the CO luminosity functions (LFs) and the evolution of the cosmic molecular gas density out to z~6 based on an 8.5 arcmin^2 spectral scan survey at 3mm of the iconic Hubble Deep Field North (HDF-N) observed with the NOrthern Extended Millimeter Array (NOEMA). We use matched filtering to search for line emission from galaxies and determine their redshift probability distr… ▽ More We present measurements of the CO luminosity functions (LFs) and the evolution of the cosmic molecular gas density out to z~6 based on an 8.5 arcmin^2 spectral scan survey at 3mm of the iconic Hubble Deep Field North (HDF-N) observed with the NOrthern Extended Millimeter Array (NOEMA). We use matched filtering to search for line emission from galaxies and determine their redshift probability distributions exploiting the extensive multi-wavelength data for the HDF-N. We identify the 7 highest-fidelity sources as CO emitters at 1<z<6, including the well-known submillimeter galaxy HDF850.1 at z=5.18. Four high-fidelity 3mm continuum sources are all found to be radio galaxies at z<=1, plus HDF850.1. We constrain the CO LFs in the HDF-N out to z~6, including a first measurement of the CO(5-4) LF at <z>=5.0. The relatively large area and depth of the NOEMA HDF-N survey extends the existing luminosity functions at 1<z<4 above the knee, yielding a somewhat lower density by 0.15-0.4 dex at the overlap region for the CO(2-1) and CO(3-2) transitions, attributed to cosmic variance. We perform a joint analysis of the CO LFs in the HDF-N and Hubble Ultra Deep Field from ASPECS, finding that they can be well described by a single Schechter function. The evolution of the cosmic molecular gas density from a joint analysis is in good agreement with earlier determinations. This implies that the impact of cosmic field-to-field variance on the measurements is consistent with previous estimates, adding to the challenges for simulations that model galaxies from first principles. △ Less

Submitted 13 January, 2023; originally announced January 2023.

Comments: 19 pages, 13 figures, 7 tables. Accepted for publication in the Astrophysical Journal (ApJ)

A solar metallicity galaxy at $z >$ 7? Possible detection of the [N II] 122 $μ$m and [O III] 52 $μ$m lines

Authors: Meghana Killi, Darach Watson, Seiji Fujimoto, Hollis Akins, Kirsten Knudsen, Johan Richard, Yuichi Harikane, Dimitra Rigopoulou, Francesca Rizzo, Michele Ginolfi, Gergö Popping, Vasily Kokorev

Abstract: We present the first detection of the [N II] 122 $μ$m and [O III] 52 $μ$m lines for a reionisation-epoch galaxy. Based on these lines and previous [C II] 158 $μ$m and [O III] 88 $μ$m measurements, using two different radiative transfer models of the interstellar medium, we estimate an upper limit on electron density of $\lesssim$ 500 cm$^{-3}$ and a gas-phase metallicity… ▽ More We present the first detection of the [N II] 122 $μ$m and [O III] 52 $μ$m lines for a reionisation-epoch galaxy. Based on these lines and previous [C II] 158 $μ$m and [O III] 88 $μ$m measurements, using two different radiative transfer models of the interstellar medium, we estimate an upper limit on electron density of $\lesssim$ 500 cm$^{-3}$ and a gas-phase metallicity $Z/Z_\odot \sim 1.1 \pm 0.2$ for A1689-zD1, a gravitationally-lensed, dusty galaxy at $z$ = 7.133. Other measurements or indicators of metallicity so far in galaxy interstellar media at $z \gtrsim$ 6 are typically an order of magnitude lower than this. The unusually high metallicity makes A1689-zD1 inconsistent with the fundamental metallicity relation, although there is likely significant dust obscuration of the stellar mass, which may partly resolve the inconsistency. Given a solar metallicity, the dust-to-metals ratio is a factor of several lower than expected, hinting that galaxies beyond $z \sim$ 7 may have lower dust formation efficiency. Finally, the inferred nitrogen enrichment compared to oxygen, on which the metallicity measurement depends, indicates that star-formation in the system is older than about 250 Myr, pushing the beginnings of this galaxy to $z >$ 10. △ Less

Submitted 8 March, 2023; v1 submitted 2 November, 2022; originally announced November 2022.

Comments: 10 pages, 7 figures; Published MNRAS: 06 March 2023

ALMACAL IX: multi-band ALMA survey for dusty star-forming galaxies and the resolved fractions of the cosmic infrared background

Authors: Jianhang Chen, R. J. Ivison, Martin A. Zwaan, Ian Smail, Anne Klitsch, Céline Péroux, Gergö Popping, Andrew D. Biggs, Roland Szakacs, Aleksandra Hamanowicz, Claudia Lagos

Abstract: Wide, deep, blind continuum surveys at submillimetre/millimetre (submm/mm) wavelengths are required to provide a full inventory of the dusty, distant Universe. However, conducting such surveys to the necessary depth, with sub-arcsec angular resolution, is prohibitively time-consuming, even for the most advanced submm/mm telescopes. Here, we report the most recent results from the ALMACAL project,… ▽ More Wide, deep, blind continuum surveys at submillimetre/millimetre (submm/mm) wavelengths are required to provide a full inventory of the dusty, distant Universe. However, conducting such surveys to the necessary depth, with sub-arcsec angular resolution, is prohibitively time-consuming, even for the most advanced submm/mm telescopes. Here, we report the most recent results from the ALMACAL project, which exploits the 'free' calibration data from the Atacama Large Millimetre/submillimetre Array (ALMA) to map the lines of sight towards and beyond the ALMA calibrators. ALMACAL has now covered 1,001 calibrators, with a total sky coverage around 0.3 deg2, distributed across the sky accessible from the Atacama desert, and has accumulated more than 1,000h of integration. The depth reached by combining multiple visits to each field makes ALMACAL capable of searching for faint, dusty, star-forming galaxies (DSFGs), with detections at multiple frequencies to constrain the emission mechanism. Based on the most up-to-date ALMACAL database, we report the detection of 186 DSFGs with flux densities down to S870um ~ 0.2mJy, comparable with existing ALMA large surveys but less susceptible to cosmic variance. We report the number counts at five wavelengths between 870um and 3mm, in ALMA bands 3, 4, 5, 6 and 7, providing a benchmark for models of galaxy formation and evolution. By integrating the observed number counts and the best-fitting functions, we also present the resolved fraction of the cosmic infrared background (CIB) and the CIB spectral shape. Combining existing surveys, ALMA has currently resolved about half of the CIB in the submm/mm regime. △ Less

Submitted 28 November, 2022; v1 submitted 17 October, 2022; originally announced October 2022.

Comments: 22 pages, 16 figures, accepted for publication in MNRAS

Journal ref: 2023MNRAS.518.1378C

Semi-analytic forecasts for Roman -- the beginning of a new era of deep-wide galaxy surveys

Authors: L. Y. Aaron Yung, Rachel S. Somerville, Steven L. Finkelstein, Peter Behroozi, Romeel Davé, Henry C. Ferguson, Jonathan P. Gardner, Gergö Popping, Sangeeta Malhotra, Casey Papovich, James E. Rhoads, Micaela B. Bagley, Michaela Hirschmann, Anton M. Koekemoer

Abstract: The Nancy Grace Roman Space Telescope, NASA's next flagship observatory, will redefine deep-field galaxy survey with a field of view two orders of magnitude larger than Hubble and an angular resolution of matching quality. These future deep-wide galaxy surveys necessitate new simulations to forecast their scientific output and to optimise survey strategies. In this work, we present five realizatio… ▽ More The Nancy Grace Roman Space Telescope, NASA's next flagship observatory, will redefine deep-field galaxy survey with a field of view two orders of magnitude larger than Hubble and an angular resolution of matching quality. These future deep-wide galaxy surveys necessitate new simulations to forecast their scientific output and to optimise survey strategies. In this work, we present five realizations of 2-deg^2 lightcones, containing a total of >25 million simulated galaxies with -16 < MUV < -25 spanning z ~ 0 to 10. This dataset enables a new set of experiments with the impacts of survey size on the derived galaxy formation and cosmological constraints. The intrinsic and observable galaxy properties are predicted using a well-established, physics-based semi-analytic modelling approach. We provide forecasts for number density, cosmic SFR, field-to-field variance, and angular two-point correlation functions, and demonstrate how the future wide-field surveys will be able to improve these measurements relative to current generation surveys. We also present a comparison between these lightcones and others that have been constructed with empirical models. The mock lightcones are designed to facilitate the exploration of multi-instrument synergies and connecting with current generation instruments and legacy surveys. In addition to Roman, we also provide photometry for a number of other instruments on upcoming facilities, including Euclid and Rubin, as well as the instruments that are part of many legacy surveys. Full object catalogues and data tables for the results presented in this work are made available through a web-based, interactive portal https://www.simonsfoundation.org/semi-analytic-forecasts. △ Less

Submitted 5 December, 2022; v1 submitted 10 October, 2022; originally announced October 2022.

Comments: 25 pages, 16 figures. Accepted for publication in MNRAS

To see or not to see a $z\sim13$ galaxy, that is the question

Authors: Melanie Kaasinen, Joshiwa van Marrewijk, Gergö Popping, Michele Ginolfi, Luca Di Mascolo, Tony Mroczkowski, Alice Concas, Claudia Di Cesare, Meghana Killi, Ivanna Langan

Abstract: Determining when the first galaxies formed remains an outstanding goal of modern observational astronomy. Theory and current stellar population models imply that the first galaxies formed at least at $z = 14-15$. But to date, only one galaxy at $z > 13$ (GS-z13-0) has been spectroscopically confirmed.. The galaxy `HD1' was recently proposed to be a z=13.27 galaxy based on its potential Lyman break… ▽ More Determining when the first galaxies formed remains an outstanding goal of modern observational astronomy. Theory and current stellar population models imply that the first galaxies formed at least at $z = 14-15$. But to date, only one galaxy at $z > 13$ (GS-z13-0) has been spectroscopically confirmed.. The galaxy `HD1' was recently proposed to be a z=13.27 galaxy based on its potential Lyman break and tentative [O III] 88 μm detection with ALMA. We hereby aim to test this scenario with new ALMA Band 4, DDT observations of what would be the [C II] 158 μm emission, if HD1 is at z$\sim$13.27. We carefully analyse the new ALMA Band 4 observations and re-analyse the existing ALMA Band 6 data on the source to determine the proposed redshift. We find a tentative $4σ$ feature in the Band 4 data that is spatially offset by 1.7" and spectrally offset by 190 km s-1 from the previously-reported $3.8σ$ `[O III] 88 μm' feature. Through various statistical tests, we demonstrate that these tentative features are fully consistent with both being random noise features. We conclude that we are more likely to be recovering noise features than both [O III] and [C II] emission from a source at $z\sim 13.27$. Although we find no credible evidence of a $z\sim 13.27$ galaxy, we cannot entirely rule out this scenario. Non-detections are also possible for a $z\sim 13$ source with a low interstellar gas-phase metallicity or ionisation parameter and/or high gas density. Determining where and exactly what type of galaxy HD1 is, will now likely require JWST/NIRSpec spectroscopy. △ Less

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

Comments: Accepted by A&A on 16/01/2023, 13 pages, 9 figures,

Investigating the [C$\,{\rm \scriptsize II}$]-to-H$\,{\rm \scriptsize I}$ conversion factor and the H$\,{\rm \scriptsize I}$ gas budget of galaxies at $z\approx 6$ with hydrodynamical simulations

Authors: David Vizgan, Kasper E. Heintz, Thomas R. Greve, Desika Narayanan, Romeel Davé, Karen P. Olsen, Gergö Popping, Darach Watson

Abstract: One of the most fundamental baryonic matter components of galaxies is the neutral atomic hydrogen (H$\,{\rm \scriptsize I}$). At low redshifts, this component can be traced directly through the 21-cm transition, but to infer H$\,{\rm \scriptsize I}$ gas content of the most distant galaxies, a viable tracer is needed. We here investigate the fidelity of the fine structure transition of the (… ▽ More One of the most fundamental baryonic matter components of galaxies is the neutral atomic hydrogen (H$\,{\rm \scriptsize I}$). At low redshifts, this component can be traced directly through the 21-cm transition, but to infer H$\,{\rm \scriptsize I}$ gas content of the most distant galaxies, a viable tracer is needed. We here investigate the fidelity of the fine structure transition of the ($^2P_{3/2} - ^2P_{1/3}$) transition of singly-ionized carbon [C$\,{\rm \scriptsize II}$] at $158\,μ$m as a proxy for H$\,{\rm \scriptsize I}$ in a set simulated galaxies at $z\approx 6$, following the work by Heintz et al. (2021). We select 11,125 star-forming galaxies from the SIMBA simulations, with far-infrared line emissions post-processed and modeled within the SIGAME framework. We find a strong connection between [C$\,{\rm \scriptsize II}$] and H$\,{\rm \scriptsize I}$, with the relation between this [C$\,{\rm \scriptsize II}$]-to-H$\,{\rm \scriptsize I}$ relation ($β_{\rm [C\,{\rm \scriptsize II}]}$) being anti-correlated with the gas-phase metallicity of the simulated galaxies. We further use these simulations to make predictions for the total baryonic matter content of galaxies at $z\approx 6$, and specifically the HI gas mass fraction. We find mean values of $M_{\rm HI}/M_\star = 1.4$, and $M_{\rm HI}/M_{\rm bar,tot} = 0.45$. These results provide strong evidence for H$\,{\rm \scriptsize I}$ being the dominant baryonic matter component by mass in galaxies at $z\approx 6$. △ Less

Submitted 7 October, 2022; originally announced October 2022.

Comments: 7 pages, 3 figures. Accepted for publication by ApJL

An upper limit on [OIII] 88 $μ$m and 1.2 mm continuum emission from a JWST $z \approx 12-13$ galaxy candidate with ALMA

Authors: Gergö Popping

Abstract: A number of new $z>11$ galaxy candidates have recently been identified based on public James Webb Space Telescope (JWST) NIRCam observations. Spectroscopic confirmation of these candidates is necessary to robustly measure their redshift and put them in the context of our understanding of the buildup of galaxies in the early Universe. GLASS-z13 is one of these candidates, with a reported photometri… ▽ More A number of new $z>11$ galaxy candidates have recently been identified based on public James Webb Space Telescope (JWST) NIRCam observations. Spectroscopic confirmation of these candidates is necessary to robustly measure their redshift and put them in the context of our understanding of the buildup of galaxies in the early Universe. GLASS-z13 is one of these candidates, with a reported photometric redshift $z>11.9$. I present publicly available Atacama Large Millimeter/submillimeter Array (ALMA) band 6 Director's Discretionary Time observations (project 2021.A.00020.S; PI T. Bakx), taken to acquire a spectroscopic redshift for GLASS-z13 by searching for [OIII] line emission in the redshift range $z=11.9-13.5$. No [OIII] emission is detected in integrated spectra extracted within an aperture around GLASS-z13, nor when using an automated line finding algorithm (applying different uv-weighting strategies for the imaging). 1.2 mm continuum emission associated to GLASS-z13 is not detected either. If GLASS-z13 is at z$\approx$12-13, this implies a 3-$σ$ upper limit on the [OIII] and rest-frame $\sim$90 $μ$m continuum emission of $\sim1\times10^8\,\rm{L}_\odot$ and 10.8 $μ$Jy, respectively. The non-detection of [OIII] and continuum emission does not necessarily imply that GLASS-z13 is not at $z\approx12-13$. It can also be explained by a low metallicity ($\sim 0.2\,\rm{Z}_\odot$ or lower) and/or high-density (at least 100 $\rm{cm}^{-3}$) interstellar medium. This work demonstrates the synergy between ALMA and JWST to study the properties of the first galaxies, although JWST/NIRSpec spectroscopy will be necessary to confirm or reject the high photometric-redshift of GLASS-z13. △ Less

Submitted 15 December, 2022; v1 submitted 27 August, 2022; originally announced August 2022.

Comments: Version accepted for publication in A&A

Journal ref: A&A 669, L8 (2023)

The strongest cool core in REXCESS: Missing X-ray cavities in RXC J2014.8-2430

Authors: Tony Mroczkowski, Megan Donahue, Joshiwa van Marrewijk, Tracy E. Clarke, Aaron Hoffer, Huib Intema, Luca Di Mascolo, Gergö Popping, Gabriel W. Pratt, Ming Sun, Mark Voit

Abstract: We present a multiwavelength study of RXC J2014.8-2430, the most extreme cool-core cluster in the Representative $XMM-Newton$ Cluster Structure Survey (REXCESS), using $Chandra$ X-ray, Southern Astrophysical Research (SOAR) Telescope, Atacama Large Millimeter/submillimeter Array (ALMA), Very Large Array (VLA), and Giant Metrewave Radio Telescope (GMRT) observations. While feedback from an active g… ▽ More We present a multiwavelength study of RXC J2014.8-2430, the most extreme cool-core cluster in the Representative $XMM-Newton$ Cluster Structure Survey (REXCESS), using $Chandra$ X-ray, Southern Astrophysical Research (SOAR) Telescope, Atacama Large Millimeter/submillimeter Array (ALMA), Very Large Array (VLA), and Giant Metrewave Radio Telescope (GMRT) observations. While feedback from an active galactic nucleus (AGN) is thought to be the dominant mechanism by which a cooling flow is suppressed, the $Chandra$ imaging observations surprisingly do not reveal the bi-lateral X-ray cavities expected in the intracluster medium (ICM) of an extreme cool core hosting a powerful radio source. We discuss the limits on the presence of any radio bubbles associated with any undetected X-ray cavities. We place upper limits on any significant X-ray AGN in the brightest cluster galaxy, and show that the X-ray peak is offset from the central radio source, which exhibits a steep low frequency radio spectrum indicative of electron ageing. The SOAR data reveal an extended, luminous emission line source. From our narrowband H$α$ imaging of the BCG, the central H$α$ peak is coincident with the radio observations, yet offset from the X-ray peak, consistent with sloshing found previously in this cluster. ALMA observations reveal a large reservoir of molecular gas that traces the extended H$α$ emission. We conclude either that the radio source and its cavities in the X-ray gas are nearly aligned along the line of sight, or that ram pressure induced by sloshing has significantly displaced the cool molecular gas feeding it, perhaps preempting the AGN feedback cycle. We argue that the sloshing near the core is likely subsonic, as expected, given the co-location of the H$α$, CO(1-0), radio continuum, and stellar emission peaks and their proximity to the intact cool core seen in X-ray. △ Less

Submitted 6 October, 2023; v1 submitted 27 July, 2022; originally announced July 2022.

Comments: 20 pages (including the appendix), 15 figures. A&A accepted on 18 July 2022

Journal ref: A&A 665, A48 (2022)

Semi-analytic forecasts for JWST -- VI. Simulated lightcones and galaxy clustering predictions

Authors: L. Y. Aaron Yung, Rachel S. Somerville, Henry C. Ferguson, Steven L. Finkelstein, Jonathan P. Gardner, Romeel Davé, Micaela B. Bagley, Gergö Popping, Peter Behroozi

Abstract: In anticipation of the new era of high-redshift exploration marked by the commissioning of the James Webb Space Telescope (JWST), we present two sets of galaxy catalogues that are designed to aid the planning and interpretation of observing programs. We provide a set of 40 wide-field lightcones with footprints spanning approximately ~ 1,000 sq. arcmin, containing galaxies up to z = 10, and a new s… ▽ More In anticipation of the new era of high-redshift exploration marked by the commissioning of the James Webb Space Telescope (JWST), we present two sets of galaxy catalogues that are designed to aid the planning and interpretation of observing programs. We provide a set of 40 wide-field lightcones with footprints spanning approximately ~ 1,000 sq. arcmin, containing galaxies up to z = 10, and a new set of 8 ultra-deep lightcones with 132 sq. arcmin footprints, containing galaxies up to z ~ 12 down to the magnitudes expected to be reached in the deepest JWST surveys. These mock lightcones are extracted from dissipationless N-body simulations and populated with galaxies using the well-established, computationally efficient Santa Cruz semi-analytic model for galaxy formation. We provide a wide range of predicted physical properties, and simulated photometry from NIRCam and many other instruments. We explore the predicted counts and luminosity functions and angular two-point correlation functions for galaxies in these simulated lightcones. We also explore the predicted field-to-field variance using multiple lightcone realizations. We find that these lightcones reproduce the available measurements of observed clustering from 0.2 < z < 7.5 very well. We provide predictions for galaxy clustering at high redshift that may be obtained from future JWST observations. All of the lightcones presented here are made available through a web-based, interactive data release portal https://flathub.flatironinstitute.org/group/sam-forecasts. △ Less

Submitted 27 July, 2022; v1 submitted 27 June, 2022; originally announced June 2022.

Comments: Final paper of the Semi-analytic forecasts for JWST series. 24 pages, 14 figures, many years of exciting science to come, accepted for publication in MNRAS

The dust-to-gas mass ratio of luminous galaxies as a function of their metallicity at cosmic noon

Authors: Gergö Popping, Irene Shivaei, Ryan L. Sanders, Tucker Jones, Alexandra Pope, Naveen A. Reddy, Alice E. Shapley, Alison L. Coil, Mariska Kriek

Abstract: We aim to quantify the relation between the dust-to-gas mass ratio (DTG) and gas-phase metallicity of $z=$2.1-2.5 luminous galaxies and contrast this high-redshift relation against analogous constraints at z$=$0. We present a sample of ten star-forming main-sequence galaxies in the redshift range $2.1<z<2.5$ with rest-optical emission-line information available from the MOSDEF survey and with ALMA… ▽ More We aim to quantify the relation between the dust-to-gas mass ratio (DTG) and gas-phase metallicity of $z=$2.1-2.5 luminous galaxies and contrast this high-redshift relation against analogous constraints at z$=$0. We present a sample of ten star-forming main-sequence galaxies in the redshift range $2.1<z<2.5$ with rest-optical emission-line information available from the MOSDEF survey and with ALMA 1.2 millimetre and CO J$=$3-2 follow-up observations. The galaxies have stellar masses ranging from $10^{10.3}$ to $10^{10.6}\,\rm{M}_\odot$ and cover a range in star-formation rate from 35 to 145 $\rm{M}_\odot\,\rm{yr}^{-1}$. We calculated the gas-phase oxygen abundance of these galaxies from rest-optical nebular emission lines (8.4 < $12 + \log{(\rm{O/H})} < 8.8$, corresponding to 0.5 - 1.25 Z$_\odot$). We estimated the dust and H$_2$ masses of the galaxies (using a metallicity-dependent CO-to-H$_2$ conversion factor) from the 1.2~mm and CO J$=$3-2 observations, respectively, from which we estimated a DTG. We find that the galaxies in this sample follow the trends already observed between CO line luminosity and dust-continuum luminosity from $z=0$ to $z=3$, extending such trends to fainter galaxies at $2.1<z<2.5$ than observed to date. We find no second-order metallicity dependence in the CO - dust-continuum luminosity relation for the galaxies presented in this work. The DTGs of main-sequence galaxies at $2.1<z<2.5$ are consistent with an increase in the DTG with gas-phase metallicity. The metallicity dependence of the DTG is driven by the metallicity dependence of the CO-to-H$_2$ conversion factor. Galaxies at z$=$2.1-2.5 are furthermore consistent with the DTG-metallicity relation found at z$=$0 (i.e. with no significant evolution), providing relevant constraints for galaxy formation models. △ Less

Submitted 19 December, 2022; v1 submitted 18 April, 2022; originally announced April 2022.

Comments: Version accepted for publication in A&A

Journal ref: A&A 670, A138 (2023)

CO Emission, Molecular Gas, and Metallicity in Main-Sequence Star-Forming Galaxies at $z\sim2.3$

Authors: Ryan L. Sanders, Alice E. Shapley, Tucker Jones, Irene Shivaei, Gergö Popping, Naveen A. Reddy, Romeel Davé, Sedona H. Price, Bahram Mobasher, Mariska Kriek, Alison L. Coil, Brian Siana

Abstract: We present observations of CO(3-2) in 13 main-sequence $z=2.0-2.5$ star-forming galaxies at $\log(M_*/M_{\odot})=10.2-10.6$ that span a wide range in metallicity (O/H) based on rest-optical spectroscopy. We find that CO(3-2)/SFR decreases with decreasing metallicity, implying that the CO luminosity per unit gas mass is lower in low-metallicity galaxies at $z\sim2$. We constrain the CO-to-H$_2$ con… ▽ More We present observations of CO(3-2) in 13 main-sequence $z=2.0-2.5$ star-forming galaxies at $\log(M_*/M_{\odot})=10.2-10.6$ that span a wide range in metallicity (O/H) based on rest-optical spectroscopy. We find that CO(3-2)/SFR decreases with decreasing metallicity, implying that the CO luminosity per unit gas mass is lower in low-metallicity galaxies at $z\sim2$. We constrain the CO-to-H$_2$ conversion factor ($α_{\text{CO}}$) and find that $α_{\text{CO}}$ inversely correlates with metallicity at $z\sim2$. We derive molecular gas masses ($M_{\text{mol}}$) and characterize the relations among $M_*$, SFR, $M_{\text{mol}}$, and metallicity. At $z\sim2$, $M_{\text{mol}}$ increases and molecular gas fraction ($M_{\text{mol}}$/$M_*$) decrease with increasing $M_*$, with a significant secondary dependence on SFR. Galaxies at $z\sim2$ lie on a near-linear molecular KS law that is well-described by a constant depletion time of 700 Myr. We find that the scatter about the mean SFR-$M_*$, O/H-$M_*$, and $M_{\text{mol}}$-$M_*$ relations is correlated such that, at fixed $M_*$, $z\sim2$ galaxies with larger $M_{\text{mol}}$ have higher SFR and lower O/H. We thus confirm the existence of a fundamental metallicity relation at $z\sim2$ where O/H is inversely correlated with both SFR and $M_{\text{mol}}$ at fixed $M_*$. These results suggest that the scatter of the $z\sim2$ star-forming main sequence, mass-metallicity relation, and $M_{\text{mol}}$-$M_*$ relation are primarily driven by stochastic variations in gas inflow rates. We place constraints on the mass loading of galactic outflows and perform a metal budget analysis, finding that massive $z\sim2$ star-forming galaxies retain only 30% of metals produced, implying that a large mass of metals resides in the circumgalactic medium. △ Less

Submitted 14 April, 2022; originally announced April 2022.

Comments: 35 pages, 17 figures, submitted to ApJ

The UV 2175Å Attenuation Bump and its Correlation with PAH Emission at z~2

Authors: Irene Shivaei, Leindert Boogaard, Tanio Díaz-Santos, Andrew Battisti, Jarle Brinchmann, Elisabete da Cunha, Michael Maseda, Jorryt Matthee, Ana Monreal-Ibero, Themiya Nanayakkara, Gergö Popping, Alba Vidal-García, Peter M. Weilbacher

Abstract: The UV bump is a broad absorption feature centered at 2175Å that is seen in the attenuation/extinction curve of some galaxies, but its origin is not well known. Here, we use a sample of 86 star-forming galaxies at z=1.7-2.7 with deep rest-frame UV spectroscopy from the MUSE HUDF Survey to study the connection between the strength of the observed UV 2175Å bump and the Spitzer/MIPS 24 micron photome… ▽ More The UV bump is a broad absorption feature centered at 2175Å that is seen in the attenuation/extinction curve of some galaxies, but its origin is not well known. Here, we use a sample of 86 star-forming galaxies at z=1.7-2.7 with deep rest-frame UV spectroscopy from the MUSE HUDF Survey to study the connection between the strength of the observed UV 2175Å bump and the Spitzer/MIPS 24 micron photometry, which at the redshift range of our sample probes mid-IR polycyclic aromatic hydrocarbon (PAH) emission at ~6-8 micron. The sample has robust spectroscopic redshifts and consists of typical main-sequence galaxies with a wide range in stellar mass (log(Mstar/Msun) ~ 8.5-10.7) and star formation rates (SFRs; SFR ~ 1-100 Msun/yr). Galaxies with MIPS detections have strong UV bumps, except for those with mass-weighted ages younger than ~150 Myr. We find that the UV bump amplitude does not change with SFR at fixed stellar mass but increases with mass at fixed SFR. The UV bump amplitude and the PAH strength (defined as mid-IR emission normalized by SFR) are highly correlated and both also correlate strongly with stellar mass. We interpret these correlations as the result of the mass-metallicity relationship, such that at low metallicities PAH emission is weak due to a lower abundance of PAH molecules. The weak or complete absence of the 2175Å bump feature on top of the underlying smooth attenuation curve at low mass/metallicities is then expected if the PAH carriers are the main source of the additional UV absorption. △ Less

Submitted 16 May, 2022; v1 submitted 17 March, 2022; originally announced March 2022.

Comments: 9 pages, 7 figures, accepted for publication in mnras

Tracing Molecular Gas in z $\simeq$ 6 Galaxies with [C${\rm \scriptsize II}$]

Authors: David Vizgan, Thomas R. Greve, Karen P. Olsen, Anita Zanella, Desika Narayanan, Romeel Davè, Georgios E. Magdis, Gergö Popping, Francesco Valentino, Kasper E. Heintz

Abstract: We investigate the fine-structure [C${\rm \scriptsize II}$] line at $158\,μ$m as a molecular gas tracer by analyzing the relationship between molecular gas mass ($M_{\rm mol}$) and [C${\rm \scriptsize II}$] line luminosity ($L_{\rm [CII]}$) in 11,125 $z\simeq 6$ star-forming, main sequence galaxies from the SIMBA simulations, with line emission modeled by SÍGAME. Though most ($\sim 50-100\,\%$) of… ▽ More We investigate the fine-structure [C${\rm \scriptsize II}$] line at $158\,μ$m as a molecular gas tracer by analyzing the relationship between molecular gas mass ($M_{\rm mol}$) and [C${\rm \scriptsize II}$] line luminosity ($L_{\rm [CII]}$) in 11,125 $z\simeq 6$ star-forming, main sequence galaxies from the SIMBA simulations, with line emission modeled by SÍGAME. Though most ($\sim 50-100\,\%$) of the gas mass in our simulations is ionized, the bulk ($> 50\,\%$) of the [C${\rm \scriptsize II}$] emission comes from the molecular phase. We find a sub-linear (slope $0.78\pm 0.01$) $\log L_{\rm [CII]}-\log M_{\rm mol}$ relation, in contrast with the linear relation derived from observational samples of more massive, metal-rich galaxies at $z \lesssim 6$. We derive a median [C${\rm \scriptsize II}$]-to-$M_{\rm mol}$ conversion factor of $α_{\rm [CII]} \simeq 18\,{\rm M_{\rm \odot}/L_{\rm \odot}}$. This is lower than the average value of $\simeq 30\,{\rm M_{\rm \odot}/L_{\rm \odot}}$ derived from observations, which we attribute to lower gas-phase metallicities in our simulations. Thus, a lower, luminosity-dependent, conversion factor must be applied when inferring molecular gas masses from [C${\rm \scriptsize II}$] observations of low-mass galaxies. For our simulations, [C${\rm \scriptsize II}$] is a better tracer of the molecular gas than CO $J=1-0$, especially at the lowest metallicities, where much of the gas is 'CO-dark'. We find that $L_{\rm [CII]}$ is more tightly correlated with $M_{\rm mol}$ than with star-formation rate (${\rm SFR}$), and both the $\log L_{\rm [CII]}-\log M_{\rm mol}$ and $\log L_{\rm [CII]}-\log {\rm SFR}$ relations arise from the Kennicutt-Schmidt relation. Our findings suggest that $L_{\rm [CII]}$ is a promising tracer of the molecular gas at the earliest cosmic epochs. △ Less

Submitted 10 March, 2022; originally announced March 2022.

Comments: 13 pages, 9 figures. Accepted for publication in ApJ

Journal ref: 2022 ApJ 929 92

Observed cosmic evolution of galaxy dust properties with metallicity and tensions with models

Authors: Gergö Popping, Céline Péroux

Abstract: The dust abundance of the interstellar medium plays an important role in galaxy physics, the chemical evolution of matter and the absorption and re-emission of stellar light. The last years have seen a surge in observational and theoretical studies constraining the dust-abundance of galaxies up to $z\sim5$. In this work we gather the latest observational measurements (with a focus on absorption st… ▽ More The dust abundance of the interstellar medium plays an important role in galaxy physics, the chemical evolution of matter and the absorption and re-emission of stellar light. The last years have seen a surge in observational and theoretical studies constraining the dust-abundance of galaxies up to $z\sim5$. In this work we gather the latest observational measurements (with a focus on absorption studies covering metallicities in the range $6.8 < 12 + \log{(O/H)}<9$) and theoretical predictions (from six different galaxy formation models) for the dust-to-gas (DTG) and dust-to-metal (DTM) ratio of galaxies. The observed trend between DTG and DTM and gas-phase metallicity can be described by a linear relation and shows no evolution from $0<z<5$. Importantly, the fit to the DTG-metallicity relation provides a refined tool for robust dust-based gas mass estimates inferred from millimeter dust-continuum observations. The lack of evolution in the observed relations are indicative of a quickly reached balance (already when the Universe was 1.2 Gyr old) between the formation and destruction of dust and a constant timescale for star-formation at fixed metallicities over cosmic time. None of the models is able to reproduce the observed trends over the entire range in metallicity and redshift probed. The comparison between models and simulations furthermore rules out some of the current implementations for the growth and destruction of dust in galaxy formation models and places tight constraints on the predicted timescale for star-formation. △ Less

Submitted 17 March, 2022; v1 submitted 7 March, 2022; originally announced March 2022.

Comments: Accepted for publication in MNRAS. Some corrections added

IR SED and Dust Masses of Sub-solar Metallicity Galaxies at z~2.3

Authors: Irene Shivaei, Gergö Popping, George Rieke, Naveen Reddy, Alexandra Pope, Robert Kennicutt, Bahram Mobasher, Alison Coil, Yoshinobu Fudamoto, Mariska Kriek, Jianwei Lyu, Pascal Oesch, Ryan Sanders, Alice Shapley, Brian Siana

Abstract: We present results from ALMA 1.2mm continuum observations of a sample of 27 star-forming galaxies at z=2.1-2.5 from the MOSFIRE Deep Evolution Field (MOSDEF) survey. These galaxies have gas-phase metallicity and star-formation rate measurements from Hb, [OIII], Ha, and [NII]. Using stacks of Spitzer, Herschel, and ALMA photometry (rest-frame ~ 8-400$μ$m), we examine the IR SED of high-redshift sub… ▽ More We present results from ALMA 1.2mm continuum observations of a sample of 27 star-forming galaxies at z=2.1-2.5 from the MOSFIRE Deep Evolution Field (MOSDEF) survey. These galaxies have gas-phase metallicity and star-formation rate measurements from Hb, [OIII], Ha, and [NII]. Using stacks of Spitzer, Herschel, and ALMA photometry (rest-frame ~ 8-400$μ$m), we examine the IR SED of high-redshift subsolar metallicity (~0.5 $Z_{\odot}$) LIRGs. We find that the data agree well with an average SED template of higher luminosity local low-metallicity dwarf galaxies (reduced $χ^2$ of 1.8). When compared with the commonly used templates for solar-metallicity local galaxies or high-redshift LIRGs and ULIRGs, even in the most favorable case (with reduced $χ^2$ of 2.8), the templates are rejected at >98% confidence level. The broader and hotter IR SED of both the local dwarfs and high-redshift subsolar metallicity galaxies may result from different grain properties, a clumpy dust geometry, or a harder/more intense ionizing radiation field that heats the dust to higher temperatures. The obscured SFR indicated by the FIR emission of the subsolar metallicity galaxies is only ~ 60% of the total SFR, which is considerably lower than that of the local LIRGs with ~ 96-97% obscured fractions. Due to the evolving IR SED shape, the local LIRG templates fit to mid-IR data can overestimate the Rayleigh-Jeans tail measurements at z~2 by a factor of 2-20, and these templates underestimate IR luminosities if fit to the observed ALMA fluxes by >0.4dex. At a given stellar mass or metallicity, dust masses at z~2.3 are an order of magnitude higher than those at z~0. Given the predicted molecular gas mass fractions, the observed z~2.3 dust-to-stellar mass ratios suggest lower dust-to-molecular gas masses than in local galaxies at the same metallicity. △ Less

Submitted 21 February, 2022; v1 submitted 11 January, 2022; originally announced January 2022.

Comments: Accepted for publication in ApJ. 29 pages, 2 appendices, 10 figures, 3 tables

Journal ref: The Astrophysical Journal 2022, Volume 928, Issue 1, id.68, 26 pp

The BPT Diagram in Cosmological Galaxy Formation Simulations: Understanding the Physics Driving Offsets at High-Redshift

Authors: Prerak Garg, Desika Narayanan, Nell Byler, Ryan L. Sanders, Alice E. Shapley, Allison L. Strom, Romeel Davé, Michaela Hirschmann, Christopher C. Lovell, Justin Otter, Gergö Popping, George C. Privon

Abstract: The Baldwin, Philips, & Terlevich diagram of [O III]/H$β$ vs. [N II]/H$α$ (hereafter N2-BPT) has long been used as a tool for classifying galaxies based on the dominant source of ionizing radiation. Recent observations have demonstrated that galaxies at $z\sim2$ reside offset from local galaxies in the N2-BPT space. In this paper, we conduct a series of controlled numerical experiments to understa… ▽ More The Baldwin, Philips, & Terlevich diagram of [O III]/H$β$ vs. [N II]/H$α$ (hereafter N2-BPT) has long been used as a tool for classifying galaxies based on the dominant source of ionizing radiation. Recent observations have demonstrated that galaxies at $z\sim2$ reside offset from local galaxies in the N2-BPT space. In this paper, we conduct a series of controlled numerical experiments to understand the potential physical processes driving this offset. We model nebular line emission in a large sample of galaxies, taken from the SIMBA cosmological hydrodynamic galaxy formation simulation, using the CLOUDY photoionization code to compute the nebular line luminosities from H II regions. We find that the observed shift toward higher [O III]/H$β$ and [N II]/H$α$ values at high redshift arises from sample selection: when we consider only the most massive galaxies $M_* \sim 10^{10-11} M_\odot$, the offset naturally appears, due to their high metallicities. We predict that deeper observations that probe lower-mass galaxies will reveal galaxies that lie on a locus comparable to $z\sim 0$ observations. Even when accounting for sample selection effects, we find that there is a subtle mismatch between simulations and observations. To resolve this discrepancy, we investigate the impact of varying ionization parameters, H II region densities, gas-phase abundance patterns, and increasing radiation field hardness on N2-BPT diagrams. We find that either decreasing the ionization parameter or increasing the N/O ratio of galaxies at fixed O/H can move galaxies along a self-similar arc in N2-BPT space that is occupied by high-redshift galaxies. △ Less

Submitted 10 January, 2022; originally announced January 2022.

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

Exploring the High-Redshift Universe with ALMA

Authors: Evanthia Hatziminaoglou, Gergö Popping, Martin Zwaan

Abstract: The properties of the interstellar medium (ISM) of the highest-redshift galaxies and quasars provide important indications of the complex interplay between the accretion of baryons onto galaxies, the physics that drives the build-up of stars out of this gas, the subsequent chemical evolution and feedback processes and the reionisation of the Universe. The Atacama Large Millimeter/submillimeter Arr… ▽ More The properties of the interstellar medium (ISM) of the highest-redshift galaxies and quasars provide important indications of the complex interplay between the accretion of baryons onto galaxies, the physics that drives the build-up of stars out of this gas, the subsequent chemical evolution and feedback processes and the reionisation of the Universe. The Atacama Large Millimeter/submillimeter Array (ALMA) continues to play a pivotal role in the characterisation of the ISM of high-redshift galaxies. Observations of the dust continuum emission, atomic fine-structure and molecular lines arising from high-redshift galaxies are now carried out routinely, providing ever more constraints on the theoretical models of galaxy formation and evolution in the early Universe. The European Astronomical Society's EAS 2021 symposium dedicated to the exploration of the high-redshift Universe with ALMA provided a forum for the observational and theoretical high-redshift ALMA communities to exchange their views and recent results in this rapidly evolving field. The article summarises the exciting results that were presented at the meeting. △ Less

Submitted 14 December, 2021; originally announced December 2021.

Comments: 4 pages, 3 figures, EAS 2021 symposium "Exploring the High-Redshift Universe with ALMA"

Journal ref: The ESO Messenger, 2021, Vol. 185

Revisiting the [C II]$_{158 μ\text{m}}$ line-intensity mapping power spectrum from the EoR using non-uniform line-luminosity scatter

Authors: Chandra Shekhar Murmu, Karen P. Olsen, Thomas R. Greve, Suman Majumdar, Kanan K. Datta, Bryan R. Scott, T. K. Daisy Leung, Romeel Dave, Gergo Popping, Raul Ortega Ochoa, David Vizgan, Desika Narayanan

Abstract: Detecting the line-intensity mapping (LIM) signal from the galaxies of the Epoch of Reionization is an emerging tool to constrain their role in reionization. Ongoing and upcoming experiments target the signal fluctuations across the sky to reveal statistical and astrophysical properties of these galaxies via signal statistics, e.g., the power spectrum. Here, we revisit the [C II]… ▽ More Detecting the line-intensity mapping (LIM) signal from the galaxies of the Epoch of Reionization is an emerging tool to constrain their role in reionization. Ongoing and upcoming experiments target the signal fluctuations across the sky to reveal statistical and astrophysical properties of these galaxies via signal statistics, e.g., the power spectrum. Here, we revisit the [C II]$_{158 μ\text{m}}$ LIM power spectrum under non-uniform line-luminosity scatter, which has a halo-mass variation of statistical properties. Line-luminosity scatter from a cosmological hydrodynamic and radiative transfer simulation of galaxies at $z=6$ is considered in this study. We test the robustness of different model frameworks that interpret the impact of the line-luminosity scatter on the signal statistics. We use a simple power-law model to fit the scatter and demonstrate that the mean luminosity-halo mass correlation fit cannot preserve the mean intensity of the LIM signal (hence the clustering power spectrum) under non-uniform scatter. In our case, the mean intensity changes by $\sim 48$ per cent compared to the mean correlation fit in contrast to the general case with semi-analytic scatter. However, we find that the prediction for the mean intensity from the most-probable fit can be modelled robustly, considering the generalized and more realistic non-uniform scatter. We also explore the possibility of diminishing luminosity bias under non-uniform scatter, affecting the clustering power spectrum, although this phenomenon might not be statistically significant. Therefore, we should adopt appropriate approaches that can consistently interpret the LIM power spectrum from observations. △ Less

Submitted 9 November, 2022; v1 submitted 20 October, 2021; originally announced October 2021.

Comments: 9 pages, 8 figures, comments are welcome, accepted in MNRAS

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

Semi-analytic forecasts for JWST -- V. AGN luminosity functions and helium reionization at z = 2-7

Authors: L. Y. Aaron Yung, Rachel S. Somerville, Steven L. Finkelstein, Michaela Hirschmann, Romeel Davé, Gergö Popping, Jonathan P. Gardner, Aparna Venkatesan

Abstract: Active galactic nuclei (AGN) forming in the early universe are thought to be the primary source of hard ionizing photons contributing to the reionization of intergalactic helium. However, the number density and spectral properties of high-redshift AGN remain largely unconstrained. In this work, we make use of physically-informed models calibrated with a wide variety of available observations to pr… ▽ More Active galactic nuclei (AGN) forming in the early universe are thought to be the primary source of hard ionizing photons contributing to the reionization of intergalactic helium. However, the number density and spectral properties of high-redshift AGN remain largely unconstrained. In this work, we make use of physically-informed models calibrated with a wide variety of available observations to provide estimates for the role of AGN throughout the Epoch of Reionization. We present AGN luminosity functions in various bands between z = 2 to 7 predicted by the well-established Santa Cruz semi-analytic model, which includes modelling of black hole accretion and AGN feedback. We then combine the predicted AGN populations with a physical spectral model for self-consistent estimates of ionizing photon production rates, which depend on the mass and accretion rate of the accreting supermassive black hole. We then couple the predicted comoving ionizing emissivity with an analytic model to compute the subsequent reionization history of intergalactic helium and hydrogen. This work demonstrates the potential of coupling physically motivated analytic or semi-analytic techniques to capture multi-scale physical processes across a vast range of scales (here, from AGN accretion disks to cosmological scales). Our physical model predicts an intrinsic ionizing photon budget well above many of the estimates in the literature, meaning that helium reionization can comfortably be accomplished even with a relatively low escape fraction. We also make predictions for the AGN populations that are expected to be detected in future \emph{James Webb Space Telescope} surveys. △ Less

Submitted 27 September, 2021; originally announced September 2021.

Comments: 28 pages, 23 figures, accepted for publication in MNRAS

An empirical representation of a physical model for the ISM [CII], CO, and [CI] emission at redshift $1\leq z\leq9$

Authors: Shengqi Yang, Gergö Popping, Rachel S. Somerville, Anthony R. Pullen, Patrick C. Breysse, Abhishek S. Maniyar

Abstract: Sub-millimeter emission lines produced by the interstellar medium (ISM) are strong tracers of star formation and are some of the main targets of line intensity mapping (LIM) surveys. In this work we present an empirical multi-line emission model that simultaneously covers the mean, scatter, and correlations of [CII], CO J=1-0 to J=5-4, and [CI] lines in redshift range $1\leq z\leq9$. We assume the… ▽ More Sub-millimeter emission lines produced by the interstellar medium (ISM) are strong tracers of star formation and are some of the main targets of line intensity mapping (LIM) surveys. In this work we present an empirical multi-line emission model that simultaneously covers the mean, scatter, and correlations of [CII], CO J=1-0 to J=5-4, and [CI] lines in redshift range $1\leq z\leq9$. We assume the galaxy ISM line emission luminosity versus halo mass relations can be described by double power laws with redshift-dependent log normal scatter. The model parameters are then derived by fitting to the state of the art semi-analytic simulation results that have successfully reproduced multiple sub-millimeter line observations at $0\leq z\lesssim6$. We cross check the line emission statistics predicted by the semi-analytic simulation and our empirical model, finding that at $z\geq1$ our model reproduces the simulated line intensities with fractional error less than about 10%. The fractional difference is less than 25% for the power spectra. Grounded on physically-motivated and self-consistent galaxy simulations, this computationally efficient model will be helpful in forecasting ISM emission line statistics for upcoming LIM surveys. △ Less

Submitted 11 March, 2022; v1 submitted 17 August, 2021; originally announced August 2021.

Comments: 16 pages, 10 figures. Correct a typo in Eq 12. Add Figure 7

On estimating the cosmic molecular gas density from CO Line Intensity Mapping observations

Authors: Patrick C. Breysse, Shengqi Yang, Rachel S. Somerville, Anthony R. Pullen, Gergö Popping, Abhishek S. Maniyar

Abstract: The Millimeter-wave Intensity Mapping Experiment (mmIME) recently reported a detection of excess spatial fluctuations at a wavelength of 3 mm, which can be attributed to unresolved emission of several CO rotational transitions between $z\sim1-5$. We study the implications of this data for the high-redshift interstellar medium using a suite of state-of-the-art semianalytic simulations which have su… ▽ More The Millimeter-wave Intensity Mapping Experiment (mmIME) recently reported a detection of excess spatial fluctuations at a wavelength of 3 mm, which can be attributed to unresolved emission of several CO rotational transitions between $z\sim1-5$. We study the implications of this data for the high-redshift interstellar medium using a suite of state-of-the-art semianalytic simulations which have successfully reproduced many other sub-millimeter line observations across the relevant redshift range. We find that the semianalytic predictions are mildly in tension with the mmIME result, with a predicted CO power $\sim3.5σ$ below what was observed. We explore some simple modifications to the models which could resolve this tension. Increasing the molecular gas abundance at the relevant redshifts to $\sim10^8\ M_\odot\ \rm{Mpc}^{-3}$, a value well above that obtained from directly imaged sources, would resolve the discrepancy, as would assuming a CO-$H_2$ conversion factor $α_{\rm{CO}}$ of $\sim1.5\ M_{\odot}$ K$^{-1}$ $(\rm{km}/\rm{s})^{-1}$ pc$^{2}$, a value somewhat lower than is commonly assumed. We go on to demonstrate that these conclusions are quite sensitive to the detailed assumptions of our simulations, highlighting the need for more careful modeling efforts as more intensity mapping data become available. △ Less

Submitted 28 June, 2021; originally announced June 2021.

Comments: 15 pages, 8 figures, submitted to ApJ

The ALMA Spectroscopic Survey in the HUDF: A Search for [CII] Emitters at $6 \leq z \leq 8$

Authors: Bade D. Uzgil, Pascal A. Oesch, Fabian Walter, Manuel Aravena, Leindert Boogaard, Chris Carilli, Roberto Decarli, Tanio Díaz-Santos, Yoshi Fudamoto, Hanae Inami, Rychard Bouwens, Paulo C. Cortes, Pierre Cox, Emmanuele Daddi, Jorge González-López, Ivo Labbe, Gergö Popping, Dominik Riechers, Mauro Stefanon, Paul Van der Werf, Axel Weiss

Abstract: The ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS) Band 6 scan (212-272 GHz) covers potential [CII] emission in galaxies at $6\leq z \leq8$ throughout a 2.9 arcmin$^2$ area. By selecting on known Lyman-$α$ emitters (LAEs) and photometric dropout galaxies in the field, we perform targeted searches down to a 5$σ$ [CII] luminosity depth $L_{\mathrm{[CII]}}\sim2.0\times10^8$ L… ▽ More The ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS) Band 6 scan (212-272 GHz) covers potential [CII] emission in galaxies at $6\leq z \leq8$ throughout a 2.9 arcmin$^2$ area. By selecting on known Lyman-$α$ emitters (LAEs) and photometric dropout galaxies in the field, we perform targeted searches down to a 5$σ$ [CII] luminosity depth $L_{\mathrm{[CII]}}\sim2.0\times10^8$ L$_{\odot}$, corresponding roughly to star formation rates (SFRs) of $10$-$20$ M$_{\odot}$ yr$^{-1}$ when applying a locally calibrated conversion for star-forming galaxies, yielding zero detections. While the majority of galaxies in this sample are characterized by lower SFRs, the resulting upper limits on [CII] luminosity in these sources are consistent with the current literature sample of targeted ALMA observations of $z=6$-$7$ LAEs and Lyman-break galaxies (LBGs), as well as the locally calibrated relations between $L_{\mathrm{[CII]}}$ and SFR -- with the exception of a single [CII]-deficient, UV luminous LBG. We also perform a blind search for [CII]-bright galaxies that may have been missed by optical selections, resulting in an upper limit on the cumulative number density of [CII] sources with $L_{\mathrm{[CII]}}>2.0\times10^8$ L$_{\odot}$ ($5σ$) to be less than $1.8\times10^{-4}$ Mpc$^{-3}$ (90% confidence level). At this luminosity depth and volume coverage, we present an observed evolution of the [CII] luminosity function from $z=6$-$8$ to $z\sim0$ by comparing the ASPECS measurement to literature results at lower redshift. △ Less

Submitted 21 February, 2021; originally announced February 2021.

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

SIGAME v3: Gas Fragmentation in Post-processing of Cosmological Simulations for More Accurate Infrared Line Emission Modeling

Authors: Karen Pardos Olsen, Blakesley Burkhart, Mordecai-Mark Mac Low, Robin G. Treß, Thomas R. Greve, David Vizgan, Jay Motka, Josh Borrow, Gergö Popping, Romeel Davé, Rowan J. Smith, Desika Narayanan

Abstract: We present an update to the framework called SImulator of GAlaxy Millimeter/submillimeter Emission (SÍGAME). SÍGAME derives line emission in the far-infrared (FIR) for galaxies in particle-based cosmological hydrodynamics simulations by applying radiative transfer and physics recipes via a post-processing step after completion of the simulation. In this version, a new technique is developed to mod… ▽ More We present an update to the framework called SImulator of GAlaxy Millimeter/submillimeter Emission (SÍGAME). SÍGAME derives line emission in the far-infrared (FIR) for galaxies in particle-based cosmological hydrodynamics simulations by applying radiative transfer and physics recipes via a post-processing step after completion of the simulation. In this version, a new technique is developed to model higher gas densities by parametrizing the gas density probability distribution function (PDF) in higher resolution simulations for use as a look-up table, allowing for more adaptive PDFs than in previous work. SÍGAME v3 is tested on redshift z = 0 galaxies drawn from the SIMBA cosmological simulation for eight FIR emission lines tracing vastly different interstellar medium phases. Including dust radiative transfer with SKIRT and high resolution photo-ionization models with Cloudy, this new method is able to self-consistently reproduce observed relations between line luminosity and star formation rate in all cases, except for [NII]122, [NII]205 and [OI]63, the luminosities of which are overestimated by median factors of 1.6, 1.2 and 1.2 dex, respectively. We attribute the remaining disagreement with observations to the lack of precise attenuation of the interstellar light on subgrid scales (<200 pc). △ Less

Submitted 3 December, 2021; v1 submitted 4 February, 2021; originally announced February 2021.

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

The dust-continuum size of TNG50 galaxies at $z=1-5$: a comparison with the distribution of stellar light, stars, dust and H$_2$

Authors: Gergö Popping, Annalisa Pillepich, Gabriela Calistro Rivera, Sebastian Schulz, Lars Hernquist, Melanie Kaasinen, Federico Marinacci, Dylan Nelson, Mark Vogelsberger

Abstract: We present predictions for the extent of the dust-continuum emission of thousands of main-sequence galaxies drawn from the TNG50 simulation between $z=1-5$. To this aim, we couple the radiative transfer code SKIRT to the output of the TNG50 simulation and measure the dust-continuum half-light radius of the modeled galaxies, assuming a Milky Way dust type and a metallicity dependent dust-to-metal r… ▽ More We present predictions for the extent of the dust-continuum emission of thousands of main-sequence galaxies drawn from the TNG50 simulation between $z=1-5$. To this aim, we couple the radiative transfer code SKIRT to the output of the TNG50 simulation and measure the dust-continuum half-light radius of the modeled galaxies, assuming a Milky Way dust type and a metallicity dependent dust-to-metal ratio. The dust-continuum half-light radius at observed-frame 850 $μ$m is up to $\sim$75 per cent larger than the stellar half-mass radius, but significantly more compact than the observed-frame 1.6 $μ$m (roughly corresponding to H-band) half-light radius, particularly towards high redshifts: the compactness compared to the 1.6 $μ$m emission increases with redshift. This is driven by obscuration of stellar light from the galaxy centres, which increases the apparent extent of 1.6 $μ$m disk sizes relative to that at 850 $μ$m. The difference in relative extents increases with redshift because the observed-frame 1.6 $μ$m emission stems from ever shorter wavelength stellar emission. These results suggest that the compact dust-continuum emission observed in $z>1$ galaxies is not (necessarily) evidence of the buildup of a dense central stellar component. We also find that the dust-continuum half-light radius very closely follows the radius containing half the star formation in galaxies, indicating that single band dust-continuum emission is a good tracer of the location of (obscured) star formation. The dust-continuum emission is more compact than the H2 mass (for galaxies at $z\geq 2$) and the underlying dust mass. The dust emission strongly correlates with locations with the highest dust temperatures, which do not need to be the locations where most H$_2$ and/or dust is located. The presented results are a common feature of main-sequence galaxies. △ Less

Submitted 28 January, 2021; originally announced January 2021.

Comments: Submitted to MNRAS

The Evolution of the IR Luminosity Function and Dust-obscured Star Formation in the Last 13 Billion Years

Authors: J. A. Zavala, C. M. Casey, S. M. Manning, M. Aravena, M. Bethermin, K. I. Caputi, D. L. Clements, E. da Cunha, P. Drew, S. L. Finkelstein, S. Fujimoto, C. Hayward, J. Hodge, J. S. Kartaltepe, K. Knudsen, A. M. Koekemoer, A. S. Long, G. E. Magdis, A. W. S. Man, G. Popping, D. Sanders, N. Scoville, K. Sheth, J. Staguhn, S. Toft , et al. (3 additional authors not shown)

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

Submitted 12 January, 2021; originally announced January 2021.

Comments: Accepted for publication in The Astrophysical Journal

The prospects for observing [OIII] 52 micron emission from galaxies during the Epoch of Reionization

Authors: Shengqi Yang, Adam Lidz, Gergö Popping

Abstract: The [OIII] 88 $μ$m fine structure emission line has been detected into the Epoch of Reionization (EoR) from star-forming galaxies at redshifts $6 < z \lesssim 9$ with ALMA. These measurements provide valuable information regarding the properties of the interstellar medium (ISM) in the highest redshift galaxies discovered thus far. The [OIII] 88 $μ$m line observations leave, however, a degeneracy b… ▽ More The [OIII] 88 $μ$m fine structure emission line has been detected into the Epoch of Reionization (EoR) from star-forming galaxies at redshifts $6 < z \lesssim 9$ with ALMA. These measurements provide valuable information regarding the properties of the interstellar medium (ISM) in the highest redshift galaxies discovered thus far. The [OIII] 88 $μ$m line observations leave, however, a degeneracy between the gas density and metallicity in these systems. Here we quantify the prospects for breaking this degeneracy using future ALMA observations of the [\oiii] 52 $μ$m line. Among the current set of ten [OIII] 88 $μ$m emitters at $6 < z \lesssim 9$, we forecast 52 $μ$m detections (at 6-$σ$) in SXDF-NB1006-2, B14-6566, J0217-0208, and J1211-0118 within on-source observing times of 2-10 hours, provided their gas densities are larger than about $n_{\mathrm{H}} \gtrsim 10^2-10^3$ cm$^{-3}$. Other targets generally require much longer integration times for a 6-$σ$ detection. Either successful detections of the 52 $μ$m line, or reliable upper limits, will lead to significantly tighter constraints on ISM parameters. The forecasted improvements are as large as $\sim 3$ dex in gas density and $\sim 1$ dex in metallicity for some regions of parameter space. We suggest SXDF-NB1006-2 as a promising first target for 52 $μ$m line measurements. We discuss how such measurements will help in understanding the mass metallicity relationship during the EoR. △ Less

Submitted 29 March, 2021; v1 submitted 3 January, 2021; originally announced January 2021.

Comments: 5 pages, 2 figures. Added new section discussing the potential implications for the mass-metallicity relation

The Terahertz Intensity Mapper (TIM): a Next-Generation Experiment for Galaxy Evolution Studies

Authors: Joaquin Vieira, James Aguirre, C. Matt Bradford, Jeffrey Filippini, Christopher Groppi, Dan Marrone, Matthieu Bethermin, Tzu-Ching Chang, Mark Devlin, Oliver Dore, Jianyang Frank Fu, Steven Hailey Dunsheath, Gilbert Holder, Garrett Keating, Ryan Keenan, Ely Kovetz, Guilaine Lagache, Philip Mauskopf, Desika Narayanan, Gergo Popping, Erik Shirokoff, Rachel Somerville, Isaac Trumper, Bade Uzgil, Jonas Zmuidzinas

Abstract: Understanding the formation and evolution of galaxies over cosmic time is one of the foremost goals of astrophysics and cosmology today. The cosmic star formation rate has undergone a dramatic evolution over the course of the last 14 billion years, and dust obscured star forming galaxies (DSFGs) are a crucial component of this evolution. A variety of important, bright, and unextincted diagnostic l… ▽ More Understanding the formation and evolution of galaxies over cosmic time is one of the foremost goals of astrophysics and cosmology today. The cosmic star formation rate has undergone a dramatic evolution over the course of the last 14 billion years, and dust obscured star forming galaxies (DSFGs) are a crucial component of this evolution. A variety of important, bright, and unextincted diagnostic lines are present in the far-infrared (FIR) which can provide crucial insight into the physical conditions of galaxy evolution, including the instantaneous star formation rate, the effect of AGN feedback on star formation, the mass function of the stars, metallicities, and the spectrum of their ionizing radiation. FIR spectroscopy is technically difficult but scientifically crucial. Stratospheric balloons offer a platform which can outperform current instrument sensitivities and are the only way to provide large-area, wide bandwidth spatial/spectral mapping at FIR wavelengths. NASA recently selected TIM, the Terahertz Intensity Mapper, with the goal of demonstrating the key technical milestones necessary for FIR spectroscopy. The TIM instrument consists of an integral-field spectrometer from 240-420 microns with 3600 kinetic-inductance detectors (KIDs) coupled to a 2-meter low-emissivity carbon fiber telescope. In this paper, we will summarize plans for the TIM experiment's development, test and deployment for a planned flight from Antarctica. △ Less

Submitted 29 September, 2020; originally announced September 2020.

Comments: 8 pages, 5 figures, published in the 30th International Symposium on Space THz Technology (ISSTT2019), Gothenburg, Sweden, April 15-17, 2019

Report number: 31ST IEEE INTERNATIONAL SYMPOSIUM ON SPACE TERAHERTZ TECHNOLOGY

Multitracer cosmological line intensity mapping mock lightcone simulation

Authors: Shengqi Yang, Rachel S. Somerville, Anthony R. Pullen, Gergö Popping, Patrick C. Breysse, Abhishek S. Maniyar

Abstract: Sub-millimeter emission lines are important tracers of the cold gas and ionized environments of galaxies and are the targets for future line intensity mapping surveys. Physics-based simulations that predict multiple emission lines arising from different phases of the interstellar medium are crucial for constraining the global physical conditions of galaxies with upcoming LIM observations. In this… ▽ More Sub-millimeter emission lines are important tracers of the cold gas and ionized environments of galaxies and are the targets for future line intensity mapping surveys. Physics-based simulations that predict multiple emission lines arising from different phases of the interstellar medium are crucial for constraining the global physical conditions of galaxies with upcoming LIM observations. In this work we present a general framework for creating multitracer mock sub-millimeter line intensity maps based on physically grounded galaxy formation and sub-mm line emission models. We simulate a mock lightcone of 2 deg$^2$ over a redshift range $0\leq z\leq10$, comprising discrete galaxies and galaxy CII, CO, CI emission. We present simulated line intensity maps for two fiducial surveys with resolution and observational frequency windows representative of COMAP and EXCLAIM. We show that the star formation rate and line emission scaling relations predicted by our simulation significantly differ at low halo masses from widely used empirical relations, which are often calibrated to observations of luminous galaxies at lower redshifts. We show that these differences lead to significant changes in key summary statistics used in intensity mapping, such as the one point intensity probability density function and the power spectrum. It will be critical to use more realistic and complex models to forecast the ability of future line intensity mapping surveys to measure observables such as the cosmic star formation rate density. △ Less

Submitted 6 March, 2021; v1 submitted 24 September, 2020; originally announced September 2020.

Comments: 19 pages, 11 figures

The Evolution of the Baryons Associated with Galaxies Averaged over Cosmic Time and Space

Authors: Fabian Walter, Chris Carilli, Marcel Neeleman, Roberto Decarli, Gergo Popping, Rachel S. Somerville, Manuel Aravena, Frank Bertoldi, Leindert Boogaard, Pierre Cox, Elisabete da Cunha, Benjamin Magnelli, Danail Obreschkow, Dominik Riechers, Hans-Walter Rix, Ian Smail, Axel Weiss, Roberto J. Assef, Franz Bauer, Rychard Bouwens, Thierry Contini, Paulo C. Cortes, Emanuele Daddi, Tanio Diaz-Santo, Jorge Gonzalez-Lopez , et al. (8 additional authors not shown)

Abstract: We combine the recent determination of the evolution of the cosmic density of molecular gas (H_2) using deep, volumetric surveys, with previous estimates of the cosmic density of stellar mass, star formation rate and atomic gas (HI), to constrain the evolution of baryons associated with galaxies averaged over cosmic time and space. The cosmic HI and H_2 densities are roughly equal at z~1.5. The H_… ▽ More We combine the recent determination of the evolution of the cosmic density of molecular gas (H_2) using deep, volumetric surveys, with previous estimates of the cosmic density of stellar mass, star formation rate and atomic gas (HI), to constrain the evolution of baryons associated with galaxies averaged over cosmic time and space. The cosmic HI and H_2 densities are roughly equal at z~1.5. The H_2 density then decreases by a factor 6^{+3}_{-2} to today's value, whereas the HI density stays approximately constant. The stellar mass density is increasing continuously with time and surpasses that of the total gas density (HI and H_2) at redshift z~1.5. The growth in stellar mass cannot be accounted for by the decrease in cosmic H_2 density, necessitating significant accretion of additional gas onto galaxies. With the new H_2 constraints, we postulate and put observational constraints on a two step gas accretion process: (i) a net infall of ionized gas from the intergalactic/circumgalactic medium to refuel the extended HI reservoirs, and (ii) a net inflow of HI and subsequent conversion to H_2 in the galaxy centers. Both the infall and inflow rate densities have decreased by almost an order of magnitude since z~2. Assuming that the current trends continue, the cosmic molecular gas density will further decrease by about a factor of two over the next 5 Gyr, the stellar mass will increase by approximately 10%, and cosmic star formation activity will decline steadily toward zero, as the gas infall and accretion shut down. △ Less

Submitted 23 September, 2020; originally announced September 2020.

Comments: Accepted for publication in The Astrophysical Journal

The ALMA Spectroscopic Survey in the HUDF: Constraining the Molecular Content at $\log{(M_*/M_\odot)} \sim 9.5$ with CO stacking of MUSE detected $z\sim1.5$ Galaxies

Authors: Hanae Inami, Roberto Decarli, Fabian Walter, Axel Weiss, Chris Carilli, Manuel Aravena, Leindert Boogaard, Jorge González-López, Gergö Popping, Elisabete da Cunha, Roland Bacon, Franz Bauer, Thierry Contini, Paulo C. Cortes, Pierre Cox, Emanuele Daddi, Tanio Díaz-Santos, Melanie Kaasinen, Dominik A. Riechers, Jeff Wagg, Paul van der Werf, Lutz Wisotzki

Abstract: We report molecular gas mass estimates obtained from a stacking analysis of CO line emission in the ALMA Spectroscopic Survey (ASPECS) using the spectroscopic redshifts from the optical integral field spectroscopic survey by the Multi Unit Spectroscopic Explorer (MUSE) of the {\it Hubble} Ultra Deep Field (HUDF). Stacking was performed on subsets of the sample of galaxies classified by their stell… ▽ More We report molecular gas mass estimates obtained from a stacking analysis of CO line emission in the ALMA Spectroscopic Survey (ASPECS) using the spectroscopic redshifts from the optical integral field spectroscopic survey by the Multi Unit Spectroscopic Explorer (MUSE) of the {\it Hubble} Ultra Deep Field (HUDF). Stacking was performed on subsets of the sample of galaxies classified by their stellar mass and position relative to the main-sequence relation (on, above, below). Among all the CO emission lines, from \cotwoone to CO(6-5), with redshifts accessible via the ASPECS Band~3 and the MUSE data, \cotwoone provides the strongest constraints on the molecular gas content. We detect \cotwoone emission in galaxies down to stellar masses of $\log{(M_*/M_\odot)}=10.0$. Below this stellar mass, we present a new constraint on the molecular gas content of $z\sim1.5$ main-sequence galaxies by stacking based on the MUSE detections. We find that the molecular gas mass of main-sequence galaxies continuously decreases with stellar mass down to $\log{(M_*/M_\odot)}\approx9.0$. Assuming a metallicity-based CO--to--$\rm H_2$ conversion factor, the molecular gas-to-stellar mass ratio from $\log{(M_*/M_\odot)}\sim9.0$ to $\sim10.0$ does not seem to decrease as fast as for $\log{(M_*/M_\odot)}>10.0$, which is in line with simulations and studies at lower redshift. The inferred molecular gas density $ρ{\rm (H_2)}=(0.49\pm0.09)\times10^8\,{\rm M_\odot\,Mpc^{-3}}$ of MUSE-selected galaxies at $z\sim1.5$ is comparable with the one derived in the HUDF with a different CO selection. Using the MUSE data we recover most of the CO emission in our deep ALMA observations through stacking, demonstrating the synergy between volumetric surveys obtained at different wavebands. △ Less

Submitted 22 September, 2020; originally announced September 2020.

Comments: 26 pages, 12 figures, 6 tables accepted for publication in ApJ

The ALMA Spectroscopic Survey in the HUDF: Multi-band constraints on line luminosity functions and the cosmic density of molecular gas

Authors: Roberto Decarli, Manuel Aravena, Leindert Boogaard, Chris Carilli, Jorge González-López, Fabian Walter, Paulo C. Cortes, Pierre Cox, Elisabete da Cunha, Emanuele Daddi, Tanio Díaz-Santos, Jacqueline A. Hodge, Hanae Inami, Marcel Neeleman, Mladen Novak, Pascal Oesch, Gergö Popping, Dominik Riechers, Ian Smail, Bade Uzgil, Paul van der Werf, Jeff Wagg, Axel Weiss

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

Submitted 22 September, 2020; originally announced September 2020.

Comments: 23 pages, 11 figures, 5 tables. Accepted for publication in the Astrophysical Journal