Showing 1–50 of 251 results for author: Ellison, S

The MAGPI Survey: Evolution of radial trends in star formation activity across cosmic time

Authors: Marcie Mun, Emily Wisnioski, Andrew J. Battisti, J. Trevor Mendel, Sara L. Ellison, Edward N. Taylor, Claudia D. P. Lagos, Katherine E. Harborne, Caroline Foster, Scott M. Croom, Sabine Bellstedt, Stefania Barsanti, Anshu Gupta, Lucas M. Valenzuela, Qian-Hui Chen, Kathryn Grasha, Tamal Mukherjee, Hye-Jin Park, Piyush Sharda, Sarah M. Sweet, Rhea-Silvia Remus, Tayyaba Zafar

Abstract: Using adaptive optics with the Multi-Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT), the Middle Ages Galaxy Properties with Integral Field Spectroscopy (MAGPI) survey allows us to study the spatially resolved Universe at a crucial time of ~4 Gyr ago ($z$ ~ 0.3) when simulations predict the greatest diversity in evolutionary pathways for galaxies. We investigate the radial tre… ▽ More Using adaptive optics with the Multi-Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT), the Middle Ages Galaxy Properties with Integral Field Spectroscopy (MAGPI) survey allows us to study the spatially resolved Universe at a crucial time of ~4 Gyr ago ($z$ ~ 0.3) when simulations predict the greatest diversity in evolutionary pathways for galaxies. We investigate the radial trends in the star formation (SF) activity and luminosity-weighted stellar ages as a function of offset from the star-forming main sequence (SFMS) for a total of 294 galaxies. Using both H$α$ emission and the 4000 Angstrom break (i.e., D4000) as star formation rate (SFR) tracers, we find overall flat radial profiles for galaxies lying on and above the SFMS, suggestive of physical processes that enhance/regulate SF throughout the entire galaxy disc. However, for galaxies lying below the SFMS, we find positive gradients in SF suggestive of inside-out quenching. Placing our results in context with results from other redshift regimes suggests an evolution in radial trends at $z$ ~ 0.3 for SF galaxies above the SFMS, from uniformly enhanced SF at $z$ ~ 1 and $z$ ~ 0.3 to centrally enhanced SF at $z$ ~ 0 (when averaged over a wide range of mass). We also capture higher local SFRs for galaxies below the SFMS compared to that of $z$ ~ 0, which can be explained by a larger population of quenched satellites in the local Universe and/or different treatments of limitations set by the D4000-sSFR relation. △ Less

Submitted 24 April, 2024; originally announced April 2024.

Comments: 19 pages, 15 figures, 4 tables, accepted for publication in MNRAS

The ALMaQUEST Survey XV: The Dependence of the Molecular-to-Atomic Gas Ratios on Resolved Optical Diagnostics

Authors: Niankun Yu, Zheng Zheng, Chao-Wei Tsai, Pei Zuo, Sara L. Ellison, David V. Stark, Di Li, Jingwen Wu, Karen L. Masters, Ting Xiao, Yinghui Zheng, Zongnan Li, Kai Zhang, Hongying Chen, Shu Liu, Sihan Jiao, Fanyi Meng

Abstract: The atomic-to-molecular gas conversion is a critical step in the baryon cycle of galaxies, which sets the initial conditions for subsequent star formation and influences the multi-phase interstellar medium. We compiled a sample of 94 nearby galaxies with observations of multi-phase gas contents by utilizing public H I, CO, and optical IFU data from the MaNGA survey together with new FAST H I obser… ▽ More The atomic-to-molecular gas conversion is a critical step in the baryon cycle of galaxies, which sets the initial conditions for subsequent star formation and influences the multi-phase interstellar medium. We compiled a sample of 94 nearby galaxies with observations of multi-phase gas contents by utilizing public H I, CO, and optical IFU data from the MaNGA survey together with new FAST H I observations. In agreement with previous results, our sample shows that the global molecular-to-atomic gas ratio ($R_{\rm mol} \equiv$ log $M_{\rm H_2}/M_{\rm H\ I}$) is correlated with the global stellar mass surface density $μ_*$ with a Kendall's $τ$ coefficient of 0.25 and $p < 10^{-3}$, less tightly but still correlated with stellar mass and NUV$-$ r color, and not related to the specific star formation rate (sSFR). The cold gas distribution and kinematics inferred from the H I and CO global profile asymmetry and shape do not significantly rely on $R_{\rm mol}$. Thanks to the availability of kpc-scale observations of MaNGA, we decompose galaxies into H II, composite, and AGN-dominated regions by using the BPT diagrams. With increasing $R_{\rm mol}$, the fraction of H II regions within 1.5 effective radius decreases slightly; the density distribution in the spatially resolved BPT diagram also changes significantly, suggesting changes in metallicity and ionization states. Galaxies with high $R_{\rm mol}$ tend to have high oxygen abundance, both at one effective radius with a Kendall's $τ$ coefficient of 0.37 ($p < 10^{-3}$) and their central regions. Among all parameters investigated here, the oxygen abundance at one effective radius has the strongest relation with global $R_{\rm mol}$, but the dependence of gas conversion on gas distribution and galaxy ionization states is weak. △ Less

Submitted 28 March, 2024; originally announced March 2024.

Comments: Accepted by SCPMA. We sincerely appreciate the constructive suggestions and comments from Lihwai Lin

Does the Fundamental Metallicity Relation Evolve with Redshift? I: The Correlation Between Offsets from the Mass-Metallicity Relation and Star Formation Rate

Authors: Alex M. Garcia, Paul Torrey, Sara Ellison, Kathryn Grasha, Lars Hernquist, Henry R. M. Zovaro, Qian-Hui Chen, Z. S. Hemler, Lisa J. Kewley, Erica J. Nelson, Ruby J. Wright

Abstract: The scatter about the mass-metallicity relation (MZR) has a correlation with the star formation rate (SFR) of galaxies. The lack of evidence of evolution in correlated scatter at $z\lesssim2.5$ leads many to refer to the relationship between mass, metallicity, and SFR as the Fundamental Metallicity Relation (FMR). Yet, recent high-redshift (z>3) JWST observations have challenged the fundamental (i… ▽ More The scatter about the mass-metallicity relation (MZR) has a correlation with the star formation rate (SFR) of galaxies. The lack of evidence of evolution in correlated scatter at $z\lesssim2.5$ leads many to refer to the relationship between mass, metallicity, and SFR as the Fundamental Metallicity Relation (FMR). Yet, recent high-redshift (z>3) JWST observations have challenged the fundamental (i.e., redshift-invariant) nature of the FMR. In this work, we show that the cosmological simulations Illustris, IllustrisTNG, and EAGLE all predict MZRs that exhibit scatter with a secondary dependence on SFR up to $z=8$. We introduce the concept of a "strong" FMR, where the strength of correlated scatter does not evolve with time, and a "weak" FMR, where there is some time evolution. We find that each simulation analysed has a weak FMR -- there is non-negligible evolution in the strength of the correlation with SFR. Furthermore, we show that the scatter is reduced an additional ~10-40% at $z\gtrsim3$ when using a weak FMR, compared to assuming a strong FMR. These results highlight the importance of avoiding coarse redshift binning when assessing the FMR. △ Less

Submitted 10 May, 2024; v1 submitted 13 March, 2024; originally announced March 2024.

Comments: 10 pages, 5 figures (Figure 5 is summary). Accepted to MNRAS

Interacting galaxies in the IllustrisTNG simulations -- VI: Reconstructed orbits, close encounters and mergers

Authors: David R. Patton, Lawrence Faria, Maan H. Hani, Paul Torrey, Sara L. Ellison, Shivani D. Thakur, Raven I. Westlake

Abstract: Cosmological simulations have been used to study interacting galaxies as a function of galaxy pair separation, enabling comparisons with observational studies of galaxy pairs. The study of interacting galaxies as a function of time (i.e. merger stage) has mostly been limited to high resolution merger simulations, due to the poor time sampling available in cosmological simulations. Building on an e… ▽ More Cosmological simulations have been used to study interacting galaxies as a function of galaxy pair separation, enabling comparisons with observational studies of galaxy pairs. The study of interacting galaxies as a function of time (i.e. merger stage) has mostly been limited to high resolution merger simulations, due to the poor time sampling available in cosmological simulations. Building on an earlier study of galaxy pairs in the IllustrisTNG cosmological simulations, we reconstruct the orbits of galaxy pairs involving massive galaxies ($M_* > 10^{10}M_{\odot}$) at redshifts of $0 \leq z < 1$, using a novel kinematic interpolation scheme to model the orbits in between the IllustrisTNG snapshots (which are separated by 162 Myr on average). We assess the accuracy of these interpolations using a pre-existing suite of merger simulations, and find that kinematic interpolations provide a remarkable improvement in accuracy compared with interpolations that use only radial separations or 3D positions. We find that nearly 90 per cent of the closest pairs ($r < 25$ kpc) have had a pericentre encounter within the past Gyr. Many of these close pairs are found on rapidly shrinking orbits, and roughly 85 per cent of these pairs will merge within 1 Gyr. However, approximately 3 per cent of these close pairs appear to be flyby systems that will never merge. These reconstructed orbits will be used in future studies to investigate how and when galaxy properties change during close encounters and mergers between galaxies in IllustrisTNG. △ Less

Submitted 27 February, 2024; originally announced February 2024.

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

The ALMaQUEST Survey XIII: Understanding radial trends in star formation quenching via the relative roles of gas availability and star formation efficiency

Authors: Hsi-An Pan, Lihwai Lin, Sara L. Ellison, Mallory D. Thorp, Sebastian F. Sanchez, Asa F. L. Bluck, Francesco Belfiore, Joanna M. Piotrowska, Jillian M. Scudder, William M. Baker

Abstract: Star formation quenching is one of the key processes that shape the evolution of galaxies. In this study, we investigate the changes in molecular gas and star formation properties as galaxies transit from the star-forming main sequence to the passive regime. Our analysis reveals that as galaxies move away from the main sequence towards the green valley the radial profile of specific star formation… ▽ More Star formation quenching is one of the key processes that shape the evolution of galaxies. In this study, we investigate the changes in molecular gas and star formation properties as galaxies transit from the star-forming main sequence to the passive regime. Our analysis reveals that as galaxies move away from the main sequence towards the green valley the radial profile of specific star formation rate surface density ($Σ_\mathrm{sSFR}$) is suppressed compared with main sequence galaxies out to a galactocentric radius of 1.5 $R_{e}$ ($\sim$ 7 kpc for our sample). By combining radial profiles of gas fraction ($f_\mathrm{gas}$) and star formation efficiency (SFE), we can discern the underlying mechanism that determines $Σ_\mathrm{sSFR}$ at different galactocentric radii. Analysis of relative contributions of $f_\mathrm{gas}$ and SFE to $Σ_\mathrm{sSFR}$ uncovers a diverse range of quenching modes. Star formation in approximately half of our quenching galaxies is primarily driven by a single mode (i.e. either $f_\mathrm{gas}$ or SFE), or a combination of both. A collective analysis of all galaxies reveals that the reduction in star formation within the central regions ($R$ $<$ 0.5 $R_{e}$) is primarily attributable to a decrease in SFE. Conversely, in the disk regions ($R$ $>$ 0.5 $R_{e}$), both $f_\mathrm{gas}$ and SFE contribute to the suppression of star formation. Our findings suggest that multiple quenching mechanisms may be at play in our sample galaxies, and even within a single galaxy. We also compare our observational outcomes with those from galaxy simulations and discuss the implications of our data. △ Less

Submitted 11 February, 2024; originally announced February 2024.

Comments: 27 pages, 8 figures, 1 table. Accepted for publication in The Astrophysical Journal

Interacting galaxies in the IllustrisTNG simulations -- VII: The connection between the most luminous active galactic nuclei and galaxy interactions

Authors: Shoshannah Byrne-Mamahit, David R. Patton, Sara L. Ellison, Robert Bickley, Leonardo Ferreira, Maan Hani, Salvatore Quai, Scott Wilkinson

Abstract: We investigate the connection between the most luminous active galactic nuclei (AGN), galaxy pairs, and post-mergers in the IllustrisTNG simulation. We select galaxy pairs and post-mergers with a mass ratio between 1:10 $< μ<$ 1:1 and a redshift between $0<z<1$. We compare the incidence of luminous AGN in pairs with matched non-pair controls, finding that AGN with luminosity… ▽ More We investigate the connection between the most luminous active galactic nuclei (AGN), galaxy pairs, and post-mergers in the IllustrisTNG simulation. We select galaxy pairs and post-mergers with a mass ratio between 1:10 $< μ<$ 1:1 and a redshift between $0<z<1$. We compare the incidence of luminous AGN in pairs with matched non-pair controls, finding that AGN with luminosity $L_{\mathrm{bol}}>10^{44}$ erg/s occur in $\sim $26\% of paired galaxies with a companion within 20 kpc, compared with $\sim $14\% in controls (a fractional excess of $\sim$2). The enhancement of AGN in galaxy pairs is luminosity dependent, with the highest excess (a factor of $\sim6\pm2$ at the closest separations) for AGN with $L_{\mathrm{bol}}>10^{45}$ erg/s. Additionally, pairs exhibit a modest yet statistically significant excess of luminous AGN up to $\sim$150 kpc in separation. For pairs which merge between $0<z<1$, AGN fractions are elevated between 1.5 Gyr before and after coalescence, with the highest excess closest to coalescence. Our results indicate that pre-coalescence interactions drive excesses of luminous AGN, but that luminous AGN in galaxy pairs are not ubiquitous. Finally, we investigate what fraction of AGN can be associated with an interaction (either having a companion within 100 kpc or a merger within the last 500 Myr). For AGN with $L_{\mathrm{bol}}>10^{45}$ erg/s, $\sim$55\% are interacting, compared with a 30\% interaction fraction in AGN with $10^{44}<L_{\mathrm{bol}}<10^{44.5}$ erg/s. Our results support a picture in which interactions play a dominant role in (but are not the sole cause of) triggering the most luminous AGN. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: 16 pages, 9 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

X-ray AGNs with SRG/eROSITA: Multi-wavelength observations reveal merger triggering and post-coalescence circumnuclear blowout

Authors: Robert W. Bickley, Sara L. Ellison, Mara Salvato, Samir Salim, David R. Patton, Andrea Merloni, Shoshannah Byrne-Mamahit, Leonardo Ferreira, Scott Wilkinson

Abstract: Major mergers between galaxies are predicted to fuel their central supermassive black holes (SMBHs), particularly after coalescence. However, determining the prevalence of active galactic nuclei (AGNs) in mergers remains a challenge, because AGN diagnostics are sensitive to details of the central structure (e.g., nuclear gas clouds, geometry and orientation of a dusty torus) that are partly decoup… ▽ More Major mergers between galaxies are predicted to fuel their central supermassive black holes (SMBHs), particularly after coalescence. However, determining the prevalence of active galactic nuclei (AGNs) in mergers remains a challenge, because AGN diagnostics are sensitive to details of the central structure (e.g., nuclear gas clouds, geometry and orientation of a dusty torus) that are partly decoupled from SMBH accretion. X-rays, expected to be ubiquitous among accreting systems, are detectable through non-Compton-thick screens of obscuring material, and thus offer the potential for a more complete assessment of AGNs in mergers. But, extant statistical X-ray studies of AGNs in mergers have been limited by either sparse, heterogeneous, or shallow on-sky coverage. We use new X-ray observations from the first SRG/eROSITA all-sky data release to characterize the incidence, luminosity, and observability of AGNs in mergers. Combining machine learning and visual classification, we identify 923 post-mergers in Dark Energy Camera Legacy Survey (DECaLS) imaging and select 4,565 interacting galaxy pairs (with separations <120 kpc and mass ratios within 1:10) from the Sloan Digital Sky Survey. We find that galaxies with X-ray AGNs are 2.0+/-0.24 times as likely to be identified as post-mergers compared to non-AGN controls, and that post-mergers are 1.8+/-0.1 times as likely to host an X-ray AGN as non-interacting controls. A multi-wavelength census of X-ray, optical, and mid-IR-selected AGNs suggests a picture wherein the underlying AGN fraction increases during pair-phase interactions, that galaxy pairs within ~20 kpc become heavily obscured, and that the obscuration often clears post-coalescence. △ Less

Submitted 30 January, 2024; originally announced January 2024.

Comments: 20 pages, 14 figures. Submitted to MNRAS, post-merger catalog will be available through the journal upon acceptance

The limitations (and potential) of non-parametric morphology statistics for post-merger identification

Authors: Scott Wilkinson, Sara L. Ellison, Connor Bottrell, Robert W. Bickley, Shoshannah Byrne-Mamahit, Leonardo Ferreira, David R. Patton

Abstract: Non-parametric morphology statistics have been used for decades to classify galaxies into morphological types and identify mergers in an automated way. In this work, we assess how reliably we can identify galaxy post-mergers with non-parametric morphology statistics. Low-redshift (z<0.2), recent (t_post-merger < 200 Myr), and isolated (r > 100 kpc) post-merger galaxies are drawn from the Illustris… ▽ More Non-parametric morphology statistics have been used for decades to classify galaxies into morphological types and identify mergers in an automated way. In this work, we assess how reliably we can identify galaxy post-mergers with non-parametric morphology statistics. Low-redshift (z<0.2), recent (t_post-merger < 200 Myr), and isolated (r > 100 kpc) post-merger galaxies are drawn from the IllustrisTNG100-1 cosmological simulation. Synthetic r-band images of the mergers are generated with SKIRT9 and degraded to various image qualities, adding observational effects such as sky noise and atmospheric blurring. We find that even in perfect quality imaging, the individual non-parametric morphology statistics fail to recover more than 55% of the post-mergers, and that this number decreases precipitously with worsening image qualities. The realistic distributions of galaxy properties in IllustrisTNG allow us to show that merger samples assembled using individual morphology statistics are biased towards low mass, high gas fraction, and high mass ratio. However, combining all of the morphology statistics together using either a linear discriminant analysis or random forest algorithm increases the completeness and purity of the identified merger samples and mitigates bias with various galaxy properties. For example, we show that in imaging similar to that of the 10-year depth of the Legacy Survey of Space and Time (LSST), a random forest can identify 89% of mergers with a false positive rate of 17%. Finally, we conduct a detailed study of the effect of viewing angle on merger observability and find that there may be an upper limit to merger recovery due to the orientation of merger features with respect to the observer. △ Less

Submitted 24 January, 2024; originally announced January 2024.

Comments: 32 pages, 21 figures Accepted for publication by MNRAS

Interplay of Stellar and Gas-Phase Metallicities: Unveiling Insights for Stellar Feedback Modeling with Illustris, IllustrisTNG, and EAGLE

Authors: Alex M. Garcia, Paul Torrey, Kathryn Grasha, Lars Hernquist, Sara Ellison, Henry R. M. Zovaro, Z. S. Hemler, Erica J. Nelson, Lisa J. Kewley

Abstract: The metal content of galaxies provides a window into their formation in the full context of the cosmic baryon cycle. In this study, we examine the relationship between stellar mass and stellar metallicity (${\rm MZ}_*{\rm R}$) in the hydrodynamic simulations Illustris, TNG, and EAGLE to understand the global properties of stellar metallicities within the feedback paradigm employed by these simulat… ▽ More The metal content of galaxies provides a window into their formation in the full context of the cosmic baryon cycle. In this study, we examine the relationship between stellar mass and stellar metallicity (${\rm MZ}_*{\rm R}$) in the hydrodynamic simulations Illustris, TNG, and EAGLE to understand the global properties of stellar metallicities within the feedback paradigm employed by these simulations. Interestingly, we observe significant variations in the overall normalization and redshift evolution of the ${\rm MZ}_*{\rm R}$ across the three simulations. However, all simulations consistently demonstrate a tertiary dependence on the specific star formation rate (sSFR) of galaxies. This finding parallels the relationship seen in both simulations and observations between stellar mass, gas-phase metallicity, and some proxy of galaxy gas content (e.g., SFR, gas fraction, atomic gas mass). Since we find this correlation exists in all three simulations, each employing a sub-grid treatment of the dense, star-forming interstellar medium (ISM) to simulate smooth stellar feedback, we interpret this result as a fairly general feature of simulations of this kind. Furthermore, with a toy analytic model, we propose that the tertiary correlation in the stellar component is sensitive to the extent of the ``burstiness'' of feedback within galaxies. △ Less

Submitted 11 March, 2024; v1 submitted 22 January, 2024; originally announced January 2024.

Comments: 15 pages, 9 figures, Accepted to MNRAS

Stellar Populations With Optical Spectra: Deep Learning vs. Popular Spectrum Fitting Codes

Authors: Joanna Woo, Dan Walters, Finn Archinuk, S. M. Faber, Sara L. Ellison, Hossen Teimoorinia, Kartheik Iyer

Abstract: We compare the performance of several popular spectrum fitting codes (Firefly, starlight, pyPipe3D and pPXF), and a deep-learning convolutional neural network (StarNet), in recovering known stellar population properties (mean stellar age, stellar metallicity, stellar mass-to-light ratio M*/L_r and the internal E(B-V)) of simulated galaxy spectra in optical wavelengths. Our mock spectra are constru… ▽ More We compare the performance of several popular spectrum fitting codes (Firefly, starlight, pyPipe3D and pPXF), and a deep-learning convolutional neural network (StarNet), in recovering known stellar population properties (mean stellar age, stellar metallicity, stellar mass-to-light ratio M*/L_r and the internal E(B-V)) of simulated galaxy spectra in optical wavelengths. Our mock spectra are constructed from star-formation histories from the IllustrisTNG100-1 simulation. These spectra mimic the Sloan Digital Sky Survey (SDSS) through a novel method of including the noise, sky residuals and emission lines taken directly from SDSS. We find that StarNet vastly outperforms all conventional codes in both speed and recovery of stellar population properties (error scatter < 0.08 dex, average biases < 0.02 dex for all tested quantities), but it requires an appropriate training set. Of the non-machine-learning codes, pPXF was a factor of 3-4 times faster than the other codes, and was the best in recovering stellar population properties (error scatter of < 0.11 dex, average biases < 0.08 dex). However, the errors and biases are strongly dependent on both true and predicted values of stellar age and metallicity, and signal-to-noise ratio. The biases of all codes can approach 0.15 dex in stellar ages, metallicities and log M*/L_r , but remain < 0.05 for E(B-V). Using unrealistic Gaussian noise in the construction of mock spectra will underestimate the errors in the metallicities by a factor of two or more, and mocks without emission lines will underestimate the errors in stellar age and M*/L_r by a factor of two. △ Less

Submitted 25 April, 2024; v1 submitted 22 January, 2024; originally announced January 2024.

Comments: Accepted to MNRAS

The ALMaQUEST Survey XII: Dense Molecular Gas as traced by HCN and HCO$^{+}$ in Green Valley Galaxies

Authors: Lihwai Lin, Hsi-An Pan, Sara L. Ellison, Nanase Harada, Maria J. Jimenez-Donaire, K. Decker French, William M. Baker, Bau-Ching Hsieh, Yusei Koyama, Carlos Lopez-Coba, Tomonari Michiyama, Kate Rowlands, Sebastian F. Sanchez, Mallory Thorp

Abstract: We present ALMA observations of two dense gas tracers, HCN(1-0) and HCO$^{+}$(1-0), for three galaxies in the green valley and two galaxies on the star-forming main sequence with comparable molecular gas fractions as traced by the CO(1-0) emissions, selected from the ALMaQUEST survey. We investigate whether the deficit of molecular gas star formation efficiency (SFE$_{\rm mol}$) that leads to the… ▽ More We present ALMA observations of two dense gas tracers, HCN(1-0) and HCO$^{+}$(1-0), for three galaxies in the green valley and two galaxies on the star-forming main sequence with comparable molecular gas fractions as traced by the CO(1-0) emissions, selected from the ALMaQUEST survey. We investigate whether the deficit of molecular gas star formation efficiency (SFE$_{\rm mol}$) that leads to the low specific star formation rate in these green valley galaxies is due to a lack of dense gas (characterized by the dense gas fraction $f_{\rm dense}$) or the low star formation efficiency of dense gas (SFE$_{\rm dense}$). We find that SFE$_{\rm mol}$ as traced by the CO emissions, when considering both star-forming and retired spaxels together, is tightly correlated with SFE$_{\rm dense}$ and depends only weakly on $f_{\rm dense}$. The specific star formation rate (sSFR) on kpc scales is primarily driven by SFE$_{\rm mol}$ and SFE$_{\rm dense}$, followed by the dependence on $f_{\rm mol}$, and is least correlated with $f_{\rm dense}$ or the dense-to-stellar mass ratio ($R_{\rm dense}$). When compared with other works in the literature, we find that our green valley sample shows lower global SFE$_{\rm mol}$ as well as lower SFE$_{\rm dense}$ while exhibiting similar dense gas fractions when compared to star-forming and starburst galaxies. We conclude that the star formation of the 3 green valley galaxies with a normal abundance of molecular gas is suppressed mainly due to the reduced SFE$_{\rm dense}$ rather than the lack of dense gas. △ Less

Submitted 11 January, 2024; originally announced January 2024.

Comments: 20 pages, 13 figures, ApJ accepted

The ALMaQUEST Survey XIV: do radial molecular gas flows affect the star-forming ability of barred galaxies?

Authors: Lucy M. Hogarth, Amélie Saintonge, Tim A. Davis, Sara L. Ellison, Lihwai Lin, Carlos López-Cobá, Hsi-An Pan, Mallory D. Thorp

Abstract: We investigate whether barred galaxies are statistically more likely to harbour radial molecular gas flows and what effect those flows have on their global properties. Using 46 galaxies from the ALMA-MaNGA QUEnching and STar formation (ALMaQUEST) survey, we identify galaxies hosting optical bars using a combination of the morphological classifications in Galaxy Zoo 2 and HyperLEDA. In order to det… ▽ More We investigate whether barred galaxies are statistically more likely to harbour radial molecular gas flows and what effect those flows have on their global properties. Using 46 galaxies from the ALMA-MaNGA QUEnching and STar formation (ALMaQUEST) survey, we identify galaxies hosting optical bars using a combination of the morphological classifications in Galaxy Zoo 2 and HyperLEDA. In order to detect radial molecular gas flows, we employ full 3D kinematic modelling of the ALMaQUEST CO(1-0) datacubes. By combining our bar classifications with our radial bar-driven flow detections, we find that galaxies classed as barred are statistically more likely to host large-scale radial gas motions compared to their un-barred and edge-on galaxy counterparts. Moreover, the majority of barred galaxies require multi-component surface brightness profiles in their best-fit models, indicative of the presence of resonance systems. We find that galaxies classed as barred with radial bar-driven flows ("barred + radial flow" subset) have significantly suppressed global star-formation efficiencies compared to barred galaxies without radial bar-driven flows and galaxies in the other morphological sub-samples. Our "barred + radial flow" subset galaxies also possess consistently centrally concentrated molecular gas distributions, with no indication of depleted gas mass fractions, suggesting that gas exhaustion is not the cause of their suppressed star formation. Furthermore, these objects have higher median gas mass surface densities in their central 1 kpc, implying that a central gas enhancements do not fuel central starbursts in these objects. We propose that dynamical effects, such as shear caused by large-scale inflows of gas, act to gravitationally stabilise the inner gas reservoirs. △ Less

Submitted 4 February, 2024; v1 submitted 22 December, 2023; originally announced December 2023.

Comments: 20 pages, 20 figures

The ALMaQUEST Survey XI: A strong but non-linear relationship between star formation and dynamical equilibrium pressure

Authors: Sara L. Ellison, Hsi-An Pan, Asa F. L. Bluck, Mark R. Krumholz, Lihwai Lin, Leslie Hunt, Edvige Corbelli, Mallory D. Thorp, Jorge Barrera-Ballesteros, Sebastian F. Sanchez, Jillian M. Scudder, Salvatore Quai

Abstract: We present the extended ALMA MaNGA QUEnching and STar formation survey, a combination of the original 46 ALMaQUEST galaxies plus new ALMA observations for a further 20 interacting galaxies. Three well-studied scaling relations are fit to the 19,999 star-forming spaxels in the extended sample, namely the resolved Schmidt-Kennicutt (rSK) relation, the resolved star forming main sequence (rSFMS) and… ▽ More We present the extended ALMA MaNGA QUEnching and STar formation survey, a combination of the original 46 ALMaQUEST galaxies plus new ALMA observations for a further 20 interacting galaxies. Three well-studied scaling relations are fit to the 19,999 star-forming spaxels in the extended sample, namely the resolved Schmidt-Kennicutt (rSK) relation, the resolved star forming main sequence (rSFMS) and the resolved molecular gas main sequence (rMGMS). We additionally investigate the relationship between the dynamical equilibrium pressure (PDE) and star formation rate surface density (Sigma_SFR), which we refer to as the resolved PDE (rPDE) relation. Contrary to previous studies that have focussed on normal star-forming galaxies and found an approximately linear rPDE relation, the presence of more vigourously star-forming galaxies in the extended ALMaQUEST sample reveals a marked turnover in the relation at high pressures. Although the scatter around the linear fit to the rPDE relation is similar to the other three relations, a random forest analysis, which can extract non-linear dependences, finds that PDE is unambiguously more important than either Sigma_H2 or Sigma_star for predicting Sigma_SFR. We compare the observed rPDE relation to the prediction of the pressure-regulated feedback-modulated (PRFM) model of star formation, finding that galaxies residing on the global SFMS do indeed closely follow the rPDE relation predicted by the PRFM theory. However, galaxies above and below the global SFMS show significant deviations from the model. Galaxies with high SFR are instead consistent with models that include other contributions to turbulence in addition to the local star formation feedback. △ Less

Submitted 5 December, 2023; originally announced December 2023.

Comments: Accepted for publication in MNRAS

$Σ_{\mathrm{SFR}}$-M* Diagram: A Valuable Galaxy Evolution Diagnostic to Complement (s)SFR-M* Diagrams

Authors: Samir Salim, Sandro Tacchella, Chandler Osborne, S. M. Faber, Janice C. Lee, Sara L. Ellison

Abstract: The specific star formation rate (sSFR) is commonly used to describe the level of galaxy star formation (SF) and to select quenched galaxies. However, being a relative measure of the young-to-old population, an ambiguity in its interpretation may arise because a small sSFR can be either because of a substantial previous mass build up, or because SF is low. We show, using large samples spanning 0 <… ▽ More The specific star formation rate (sSFR) is commonly used to describe the level of galaxy star formation (SF) and to select quenched galaxies. However, being a relative measure of the young-to-old population, an ambiguity in its interpretation may arise because a small sSFR can be either because of a substantial previous mass build up, or because SF is low. We show, using large samples spanning 0 < z < 2, that the normalization of SFR by the physical extent over which SF is taking place (i.e., SFR surface density, $Σ_{\mathrm{SFR}}$) overcomes this ambiguity. $Σ_{\mathrm{SFR}}$ has a strong physical basis, being tied to the molecular gas density and the effectiveness of stellar feedback, so we propose $Σ_{\mathrm{SFR}}$-M* as an important galaxy evolution diagram to complement (s)SFR-M* diagrams. Using the $Σ_{\mathrm{SFR}}$-M* diagram we confirm the Schiminovich et al. (2007) result that the level of SF along the main sequence today is only weakly mass dependent - high-mass galaxies, despite their redder colors, are as active as blue, low-mass ones. At higher redshift, the slope of the "$Σ_{\mathrm{SFR}}$ main sequence" steepens, signaling the epoch of bulge build-up in massive galaxies. We also find that $Σ_{\mathrm{SFR}}$ based on the optical isophotal radius more cleanly selects both the starbursting and the spheroid-dominated (early-type) galaxies than sSFR. One implication of our analysis is that the assessment of the inside-out vs. outside-in quenching scenarios should consider both sSFR and $Σ_{\mathrm{SFR}}$ radial profiles, because ample SF may be present in bulges with low sSFR (red color). △ Less

Submitted 17 October, 2023; originally announced October 2023.

Comments: 16 pages. Accepted to ApJ. Comments on content or relevant missing references welcome

VERTICO and IllustrisTNG: The spatially resolved effects of environment on galactic gas

Authors: Adam R. H. Stevens, Toby Brown, Benedikt Diemer, Annalisa Pillepich, Lars Hernquist, Dylan Nelson, Yannick M. Bahé, Alessandro Boselli, Timothy A. Davis, Pascal J. Elahi, Sara L. Ellison, María J. Jiménez-Donaire, Ian D. Roberts, Kristine Spekkens, Vicente Villanueva, Adam B. Watts, Christine D. Wilson, Nikki Zabel

Abstract: It has been shown in previous publications that the TNG100 simulation quantitatively reproduces the observed reduction in each of the total atomic and total molecular hydrogen gas for galaxies within massive halos, i.e.~dense environments. In this Letter, we study how well TNG50 reproduces the resolved effects of a Virgo-like cluster environment on the gas surface densities of satellite galaxies w… ▽ More It has been shown in previous publications that the TNG100 simulation quantitatively reproduces the observed reduction in each of the total atomic and total molecular hydrogen gas for galaxies within massive halos, i.e.~dense environments. In this Letter, we study how well TNG50 reproduces the resolved effects of a Virgo-like cluster environment on the gas surface densities of satellite galaxies with $m_* > \! 10^9\,{\rm M}_\odot$ and ${\rm SFR} \! > 0.05\,{\rm M}_\odot\,{\rm yr}^{-1}$. We select galaxies in the simulation that are analogous to those in the HERACLES and VERTICO surveys, and mock-observe them to the common specifications of the data. Although TNG50 does not quantitatively match the observed gas surface densities in the centers of galaxies, the simulation does qualitatively reproduce the trends of gas truncation and central density suppression seen in VERTICO in both HI and H$_2$. This result promises that modern cosmological hydrodynamic simulations can be used to reliably model the post-infall histories of cluster satellite galaxies. △ Less

Submitted 11 October, 2023; originally announced October 2023.

Comments: 9 pages, 4 figures, accepted in ApJL

VERTICO VII: Environmental quenching caused by suppression of molecular gas content and star formation efficiency in Virgo Cluster galaxies

Authors: Toby Brown, Ian D. Roberts, Mallory Thorp, Sara L. Ellison, Nikki Zabel, Christine D. Wilson, Yannick M. Bahé, Dhruv Bisaria, Alberto D. Bolatto, Alessandro Boselli, Aeree Chung, Luca Cortese, Barbara Catinella, Timothy A. Davis, María J. Jiménez-Donaire, Claudia D. P. Lagos, Bumhyun Lee, Laura C. Parker, Rory Smith, Kristine Spekkens, Adam R. H. Stevens, Vicente Villanueva, Adam B. Watts

Abstract: We study how environment regulates the star formation cycle of 33 Virgo Cluster satellite galaxies on 720 parsec scales. We present the first resolved star-forming main sequence for cluster galaxies, dividing the sample based on their global HI properties and comparing to a control sample of field galaxies. HI-poor cluster galaxies have reduced star formation rate (SFR) surface densities with resp… ▽ More We study how environment regulates the star formation cycle of 33 Virgo Cluster satellite galaxies on 720 parsec scales. We present the first resolved star-forming main sequence for cluster galaxies, dividing the sample based on their global HI properties and comparing to a control sample of field galaxies. HI-poor cluster galaxies have reduced star formation rate (SFR) surface densities with respect to both HI-normal cluster and field galaxies (0.5 dex), suggesting that mechanisms regulating the global HI content are responsible for quenching local star formation. We demonstrate that the observed quenching in HI-poor galaxies is caused by environmental processes such as ram pressure stripping (RPS) simultaneously reducing molecular gas surface density and star formation efficiency (SFE), compared to regions in HI-normal systems (by 0.38 and 0.22 dex, respectively). We observe systematically elevated SFRs that are driven by increased molecular gas surface densities at fixed stellar mass surface density in the outskirts of early-stage RPS galaxies, while SFE remains unchanged with respect to the field sample. We quantify how RPS and starvation affect the star formation cycle of inner and outer galaxy discs as they are processed by the cluster. We show both are effective quenching mechanisms with the key difference being that RPS acts upon the galaxy outskirts while starvation regulates the star formation cycle throughout disc, including within the truncation radius. For both processes, the quenching is caused by a simultaneous reduction in molecular gas surface densities and SFE at fixed stellar mass surface density. △ Less

Submitted 21 August, 2023; originally announced August 2023.

Comments: 17 pages, 1 table, 5 figures, accepted for publication in ApJ

VERTICO VI: Cold-gas asymmetries in Virgo cluster galaxies

Authors: Ian D. Roberts, Toby Brown, Nikki Zabel, Christine D. Wilson, Aeree Chung, Laura C. Parker, Dhruv Bisaria, Alessandro Boselli, Barbara Catinella, Ryan Chown, Luca Cortese, Timothy A. Davis, Sara Ellison, Maria Jesus Jimenez-Donaire, Bumhyun Lee, Rory Smith, Kristine Spekkens, Adam R. H. Stevens, Mallory Thorp, Vincente Villanueva, Adam B. Watts, Charlotte Welker, Hyein Yoon

Abstract: We analyze cold-gas distributions in Virgo cluster galaxies using resolved CO(2-1) (tracing molecular hydrogen, H2) and HI observations from the Virgo Environment Traced In CO (VERTICO) and the VLA Imaging of Virgo in Atomic Gas (VIVA) surveys. From a theoretical perspective, it is expected that environmental processes in clusters will have a stronger influence on diffuse atomic gas compared to th… ▽ More We analyze cold-gas distributions in Virgo cluster galaxies using resolved CO(2-1) (tracing molecular hydrogen, H2) and HI observations from the Virgo Environment Traced In CO (VERTICO) and the VLA Imaging of Virgo in Atomic Gas (VIVA) surveys. From a theoretical perspective, it is expected that environmental processes in clusters will have a stronger influence on diffuse atomic gas compared to the relatively dense molecular gas component, and that these environmental perturbations can compress the cold interstellar medium in cluster galaxies leading to elevated star formation. In this work we observationally test these predictions for star-forming satellite galaxies within the Virgo cluster. We divide our Virgo galaxy sample into HI-normal, HI-tailed, and HI-truncated classes and show, unsurprisingly, that the HI-tailed galaxies have the largest quantitative HI asymmetries. We also compare to a control sample of non-cluster galaxies and find that Virgo galaxies, on average, have HI asymmetries that are 40 +/- 10 per cent larger than the control. There is less separation between control, HI-normal, HI-tailed, and HI-truncated galaxies in terms of H2 asymmetries, and on average, Virgo galaxies have H2 asymmetries that are only marginally (20 +/- 10 per cent) larger than the control sample. We find a weak correlation between HI and H2 asymmetries over our entire sample, but a stronger correlation for those specific galaxies being strongly impacted by environmental perturbations. Finally, we divide the discs of the HI-tailed Virgo galaxies into a leading half and trailing half according to the observed tail direction. We find evidence for excess molecular gas mass on the leading halves of the disc. This excess molecular gas on the leading half is accompanied by an excess in star formation rate such that the depletion time is, on average, unchanged. △ Less

Submitted 24 May, 2023; originally announced May 2023.

Comments: 15 pages, 8 figures, 1 table, accepted for publication in A&A

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

Post-Starburst Properties of Post-Merger Galaxies

Authors: Wenhao Li, Preethi Nair, Kate Rowlands, Karen Masters, David Stark, Niv Drory, Sara Ellison, Jimmy Irwin, Shobita Satyapal, Amy Jones, William Keel, Kavya Mukundan, Zachary Tu

Abstract: Post-starburst galaxies (PSBs) are transition galaxies showing evidence of recent rapid star formation quenching. To understand the role of galaxy mergers in triggering quenching, we investigate the incidence of PSBs and resolved PSB properties in post-merger galaxies using both SDSS single-fiber spectra and MaNGA resolved IFU spectra. We find post-mergers have a PSB excess of 10 - 20 times that r… ▽ More Post-starburst galaxies (PSBs) are transition galaxies showing evidence of recent rapid star formation quenching. To understand the role of galaxy mergers in triggering quenching, we investigate the incidence of PSBs and resolved PSB properties in post-merger galaxies using both SDSS single-fiber spectra and MaNGA resolved IFU spectra. We find post-mergers have a PSB excess of 10 - 20 times that relative to their control galaxies using single-fiber PSB diagnostics. A similar excess of ~ 19 times is also found in the fraction of central (C)PSBs and ring-like (R)PSBs in post-mergers using the resolved PSB diagnostic. However, 60% of the CPSBs + RPSBs in both post-mergers and control galaxies are missed by the single-fiber data. By visually inspecting the resolved PSB distribution, we find that the fraction of outside-in quenching is 7 times higher than inside-out quenching in PSBs in post-mergers while PSBs in control galaxies do not show large differences in these quenching directions. In addition, we find a marginal deficit of HI gas in PSBs relative to non-PSBs in post-mergers using the MaNGA-HI data. The excesses of PSBs in post-mergers suggest that mergers play an important role in triggering quenching. Resolved IFU spectra are important to recover the PSBs missed by single-fiber spectra. The excess of outside-in quenching relative to inside-out quenching in post-mergers suggests that AGN are not the dominant quenching mechanism in these galaxies, but that processes from the disk (gas inflows/consumption and stellar feedback) play a more important role. △ Less

Submitted 12 May, 2023; originally announced May 2023.

Comments: Accepted in MNRAS on May 12 2023, 19 pages, 15 figures, 4 tables

Evidence of First Stars-enriched Gas in High-redshift Absorbers

Authors: A. Saccardi, S. Salvadori, V. D'Odorico, G. Cupani, M. Fumagalli, T. A. M. Berg, G. D. Becker, S. Ellison, S. Lopez

Abstract: The first stars were born from chemically pristine gas. They were likely massive, and thus they rapidly exploded as supernovae, enriching the surrounding gas with the first heavy elements. In the Local Group, the chemical signatures of the first stellar population were identified among low-mass, long-lived, very metal-poor ([Fe/H]<-2) stars, characterized by high abundances of carbon over iron ([C… ▽ More The first stars were born from chemically pristine gas. They were likely massive, and thus they rapidly exploded as supernovae, enriching the surrounding gas with the first heavy elements. In the Local Group, the chemical signatures of the first stellar population were identified among low-mass, long-lived, very metal-poor ([Fe/H]<-2) stars, characterized by high abundances of carbon over iron ([C/Fe]>+0.7): the so-called carbon-enhanced metal-poor stars. Conversely, a similar carbon excess caused by first-star pollution was not found in dense neutral gas traced by absorption systems at different cosmic time. Here we present the detection of 14 very metal-poor, optically thick absorbers at redshift z~3-4. Among these, 3 are carbon-enhanced and reveal an overabundance with respect to Fe of all the analyzed chemical elements (O, Mg, Al, and Si). Their relative abundances show a distribution with respect to [Fe/H] that is in very good agreement with those observed in nearby very metal-poor stars. All the tests we performed support the idea that these C-rich absorbers preserve the chemical yields of the first stars. Our new findings suggest that the first-star signatures can survive in optically thick but relatively diffuse absorbers, which are not sufficiently dense to sustain star formation and hence are not dominated by the chemical products of normal stars. △ Less

Submitted 3 May, 2023; originally announced May 2023.

Interacting galaxies in the IllustrisTNG simulations -- V: Comparing the influence of star-forming vs. passive companions

Authors: Westley Brown, David R. Patton, Sara L. Ellison, Lawrence Faria

Abstract: We study interacting galaxy pairs in the TNG100-1 and TNG300-1 cosmological simulations using previously generated closest companion samples. We study the specific star formation rates (sSFR) of massive ($10^{10} M_{\odot} < M_* < 10^{12} M_{\odot}$) galaxies at $z \leq 0.2$ as a function of separation from the closest companion galaxy. We split our sample based on whether the companion galaxy is… ▽ More We study interacting galaxy pairs in the TNG100-1 and TNG300-1 cosmological simulations using previously generated closest companion samples. We study the specific star formation rates (sSFR) of massive ($10^{10} M_{\odot} < M_* < 10^{12} M_{\odot}$) galaxies at $z \leq 0.2$ as a function of separation from the closest companion galaxy. We split our sample based on whether the companion galaxy is star-forming or passive. We find that galaxies with close star-forming companions have sSFRs that are enhanced (on average) by a factor of $2.9 \pm 0.3$ in TNG100-1 and $2.27 \pm 0.06$ in TNG300-1 compared to controls, with enhancements present out to separations of $\sim 300$ kpc. Galaxies with passive companions in TNG300-1 exhibit mild sSFR suppression ($\sim12$ percent) at 100-300 kpc and small sSFR enhancements at separations below 50 kpc. sSFR suppression is strongest in pairs where the galaxy's stellar mass is more than 2 times that of its passive companion. By generating a stellar mass-matched ("twinned") sample in TNG300-1, we show that differences in sSFR trends between companion types are not a result of intrinsic stellar mass differences in star-forming vs. passive galaxies. We compare with an analogous sample of galaxy pairs from SDSS, finding consistent results between observations and simulations. Overall, we find that star-forming galaxies show enhanced sSFRs regardless of companion type, but that galaxies with close passive companions are more likely to be passive themselves. △ Less

Submitted 27 April, 2023; originally announced April 2023.

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

Galaxy interactions are the dominant trigger for local type 2 quasars

Authors: J. C. S. Pierce, C. Tadhunter, C. Ramos Almeida, P. Bessiere, J. V. Heaton, S. L. Ellison, G. Speranza, Y. Gordon, C. O'Dea, L. Grimmett, L. Makrygianni

Abstract: The triggering mechanism for the most luminous, quasar-like active galactic nuclei (AGN) remains a source of debate, with some studies favouring triggering via galaxy mergers, but others finding little evidence to support this mechanism. Here, we present deep Isaac Newton Telescope/Wide Field Camera imaging observations of a complete sample of 48 optically-selected type 2 quasars $-$ the QSOFEED s… ▽ More The triggering mechanism for the most luminous, quasar-like active galactic nuclei (AGN) remains a source of debate, with some studies favouring triggering via galaxy mergers, but others finding little evidence to support this mechanism. Here, we present deep Isaac Newton Telescope/Wide Field Camera imaging observations of a complete sample of 48 optically-selected type 2 quasars $-$ the QSOFEED sample (L$_{\rm [OIII]}>$10$^{8.5}$$L_{\odot}$; $z < 0.14$). Based on visual inspection by eight classifiers, we find clear evidence that galaxy interactions are the dominant triggering mechanism for quasar activity in the local universe, with 65$^{+6}_{-7}$ per cent of the type 2 quasar hosts showing morphological features consistent with galaxy mergers or encounters, compared with only 22$^{+5}_{-4}$ per cent of a stellar-mass- and redshift-matched comparison sample of non-AGN galaxies $-$ a 5$σ$ difference. The type 2 quasar hosts are a factor 3.0$^{+0.5}_{-0.8}$ more likely to be morphologically disturbed than their matched non-AGN counterparts, similar to our previous results for powerful 3CR radio AGN of comparable [OIII] emission-line luminosity and redshift. In contrast to the idea that quasars are triggered at the peaks of galaxy mergers as the two nuclei coalesce, and only become visible post-coalescence, the majority of morphologically-disturbed type 2 quasar sources in our sample are observed in the pre-coalescence phase (61$^{+8}_{-9}$ per cent). We argue that much of the apparent ambiguity that surrounds observational results in this field is a result of differences in the surface brightness depths of the observations, combined with the effects of cosmological surface brightness dimming. △ Less

Submitted 27 March, 2023; originally announced March 2023.

Comments: Accepted for publication in MNRAS. 17 pages, 7 figures

VERTICO V: The environmentally driven evolution of the inner cold gas discs of Virgo cluster galaxies

Authors: Adam B. Watts, Luca Cortese, Barbara Catinella, Toby Brown, Christine D. Wilson, Nikki Zabel, Ian D. Roberts, Timothy A. Davis, Mallory Thorp, Aeree Chung, Adam R. H. Stevens, Sara L. Ellison, Kristine Spekkens, Laura C. Parker, Yannick M. Bahé, Vicente Villanueva, María Jiménez-Donaire, Dhruv Bisaria, Alessandro Boselli, Alberto D. Bolatto, Bumhyun Lee

Abstract: The quenching of cluster satellite galaxies is inextricably linked to the suppression of their cold interstellar medium (ISM) by environmental mechanisms. While the removal of neutral atomic hydrogen (HI) at large radii is well studied, how the environment impacts the remaining gas in the centres of galaxies, which are dominated by molecular gas, is less clear. Using new observations from the Virg… ▽ More The quenching of cluster satellite galaxies is inextricably linked to the suppression of their cold interstellar medium (ISM) by environmental mechanisms. While the removal of neutral atomic hydrogen (HI) at large radii is well studied, how the environment impacts the remaining gas in the centres of galaxies, which are dominated by molecular gas, is less clear. Using new observations from the Virgo Environment traced in CO survey (VERTICO) and archival HI data, we study the HI and molecular gas within the optical discs of Virgo cluster galaxies on 1.2-kpc scales with spatially resolved scaling relations between stellar ($Σ_{\star}$), HI ($Σ_\mathrm{HI}$), and molecular gas ($Σ_\mathrm{mol}$) surface densities. Adopting HI deficiency as a measure of environmental impact, we find evidence that, in addition to removing the HI at large radii, the cluster processes also lower the average $Σ_\mathrm{HI}$ of the remaining gas even in the central 1.2 kpc. The impact on molecular gas is comparatively weaker than on the HI, and we show that the lower $Σ_\mathrm{mol}$ gas is removed first. In the most HI-deficient galaxies, however, we find evidence that environmental processes reduce the typical $Σ_\mathrm{mol}$ of the remaining gas by nearly a factor of 3. We find no evidence for environment-driven elevation of $Σ_\mathrm{HI}$ or $Σ_\mathrm{mol}$ in HI-deficient galaxies. Using the ratio of $Σ_\mathrm{mol}$-to-$Σ_\mathrm{HI}$ in individual regions, we show that changes in the ISM physical conditions, estimated using the total gas surface density and midplane hydrostatic pressure, cannot explain the observed reduction in molecular gas content. Instead, we suggest that direct stripping of the molecular gas is required to explain our results. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Comments: Accepted for publication in PASA. Main text is 19 pages including 12 figures and 3 tables, plus 1 appendix. A 2.5 min, high-level summary can be found at https://youtu.be/7djMmVEpDVc

A Multiwavelength Study of Active Galactic Nuclei in Post-Merger Remnants

Authors: Wenhao Li, Preethi Nair, Jimmy Irwin, Sara Ellison, Shobita Satyapal, Niv Drory, Amy Jones, William Keel, Karen Masters, David Stark, Russell Ryan, Kavya Mukundan

Abstract: We investigate the role of galaxy mergers in triggering AGN in the nearby Universe. Our analysis is based on a sample of 79 post-merger remnant galaxies with deep X-ray observations from Chandra/XMM-Newton capable of detecting a low-luminosity AGN of > 10^40.5 erg s^-1. This sample is derived from a visually classified, volume-limited sample of 807 post-mergers identified in the Sloan Digital Sky… ▽ More We investigate the role of galaxy mergers in triggering AGN in the nearby Universe. Our analysis is based on a sample of 79 post-merger remnant galaxies with deep X-ray observations from Chandra/XMM-Newton capable of detecting a low-luminosity AGN of > 10^40.5 erg s^-1. This sample is derived from a visually classified, volume-limited sample of 807 post-mergers identified in the Sloan Digital Sky Survey Data Release 14 with log M*/M_sun > 10.5 and 0.02 < z < 0.06. We find the X-ray AGN fraction in this sample is 55.7% +\- 5.6% compared to 23.6% +\- 2.8% for a mass and redshift matched non-interacting control sample. The multi-wavelength AGN fraction (identified as an AGN in one of X-ray, IR, radio or optical diagnostics) for post-mergers is 76.6% +\- 4.8% compared to 39.1% +\- 3.2% for controls. Thus post-mergers exhibit a high overall AGN fraction with an excess between 2 - 4 depending on the AGN diagnostics used. In addition, we find most optical, IR, and radio AGN are also identified as X-ray AGN while a large fraction of X-ray AGN are not identified in any other diagnostic. This highlights the importance of deep X-ray imaging to identify AGN. We find the X-ray AGN fraction of post-mergers is independent of stellar mass above log M*/M_sun > 10.5 unlike the trend seen in control galaxies. Overall, our results show that post-merger galaxies are a good tracer of the merger-AGN connection and strongly support the theoretical expectations that mergers trigger AGN. △ Less

Submitted 22 February, 2023; v1 submitted 15 January, 2023; originally announced January 2023.

Comments: Accepted for publication in ApJ, 27 pages, 16 figures, 3 tables

Journal ref: ApJ 944 168 (2023)

AGN in post-mergers from the Ultraviolet Near Infrared Optical Northern Survey

Authors: Robert W. Bickley, Sara L. Ellison, David R. Patton, Scott Wilkinson

Abstract: The kinematic disturbances associated with major galaxy mergers are known to produce gas inflows, which in turn may trigger accretion onto the supermassive black holes (SMBH) of the participant galaxies. While this effect has been studied in galaxy pairs, the frequency of active galactic nuclei (AGN) in fully coalesced post-merger systems is poorly constrained due to the limited size or impurity o… ▽ More The kinematic disturbances associated with major galaxy mergers are known to produce gas inflows, which in turn may trigger accretion onto the supermassive black holes (SMBH) of the participant galaxies. While this effect has been studied in galaxy pairs, the frequency of active galactic nuclei (AGN) in fully coalesced post-merger systems is poorly constrained due to the limited size or impurity of extant post-merger samples. Previously, we combined convolutional neural network (CNN) predictions with visual classifications to identify a highly pure sample of 699 post-mergers in deep r-band imaging. In the work presented here, we quantify the frequency of AGN in this sample using three metrics: optical emission lines, mid-infrared (mid- IR) colour, and radio detection of low-excitation radio galaxies (LERGs). We also compare the frequency of AGN in post-mergers to that in a sample of spectroscopically identified galaxy pairs. We find that AGN identified by narrow-line optical emission and mid-IR colour have an increased incidence rate in post-mergers, with excesses of ~4 over mass- and redshift-matched controls. The optical and mid-IR AGN excesses in post-mergers exceed the values found for galaxy pairs, indicating that AGN activity in mergers peaks after coalescence. Conversely, we recover no significant excess of LERGs in post-mergers or pairs. Finally, we find that the [OIII] luminosity (a proxy for SMBH accretion rate) in post-mergers that host an optical AGN is ~0.3 dex higher on average than in non-interacting galaxies with an optical AGN, suggesting that mergers generate higher accretion rates than secular triggering mechanisms. △ Less

Submitted 9 January, 2023; originally announced January 2023.

Comments: 13 pages, 5 figures

The interconnection between galaxy mergers, AGN activity and rapid quenching of star formation in simulated post-merger galaxies

Authors: Salvatore Quai, Shoshannah Byrne-Mamahit, Sara L. Ellison, David R. Patton, Maan H. Hani

Abstract: We investigate the role of galaxy mergers on supermassive black hole (SMBH) accretion and star formation quenching in three state-of-the-art cosmological simulations with contrasting physics models: EAGLE, Illustris and IllustrisTNG. We find that recently coalesced 'post-mergers' in all three simulations have elevated SMBH accretion rates by factors of ~2-5. However, rapid (within 500 Myr of coale… ▽ More We investigate the role of galaxy mergers on supermassive black hole (SMBH) accretion and star formation quenching in three state-of-the-art cosmological simulations with contrasting physics models: EAGLE, Illustris and IllustrisTNG. We find that recently coalesced 'post-mergers' in all three simulations have elevated SMBH accretion rates by factors of ~2-5. However, rapid (within 500 Myr of coalescence) quenching of star formation is rare, with incidence rates of 0.4% in Illustris, 4.5% in EAGLE and 10% in IllustrisTNG. The rarity of quenching in post-mergers results from substantial gas reservoirs that remain intact after the merger. The post-mergers that do successfully quench tend to be those that had both low pre-merger gas fractions as well as those that experience the largest gas losses. Although rare, the recently quenched fraction of post-mergers is still elevated compared to a control sample of non-mergers by factors of two in IllustrisTNG and 11 in EAGLE. Conversely, quenching is rarer in Illustris post-mergers than in their control. Recent observational results by Ellison et al. have found rapid quenching to be at least 30 times more common in post-mergers, a significantly higher excess than found in any of the simulations. Our results, therefore, indicate that whilst merger-induced SMBH accretion is a widespread prediction of the simulations, its link to quenching depends sensitively on the physics models, and that none of the subgrid models of the simulations studied here can fully capture the connection between mergers and rapid quenching seen in observations. △ Less

Submitted 20 December, 2022; originally announced December 2022.

Comments: This article has been accepted for publication in Monthly Notices of the Royal Astronomical Society

PHANGS-JWST First Results: A Global and Moderately Resolved View of Mid-Infrared and CO Line Emission from Galaxies at the Start of the JWST Era

Authors: Adam K. Leroy, Alberto D. Bolatto, Karin Sandstrom, Erik Rosolowsky, Ashley. T. Barnes, F. Bigiel, Médéric Boquien, Jakob S. den Brok, Yixian Cao, Jérémy Chastenet, Mélanie Chevance, I-Da Chiang, Ryan Chown, Dario Colombo, Sara L. Ellison, Eric Emsellem, Kathryn Grasha, Jonathan D. Henshaw, Annie Hughes, Ralf S. Klessen, Eric W. Koch, Jaeyeon Kim, Kathryn Kreckel, J. M. Diederik Kruijssen, Kirsten L. Larson , et al. (19 additional authors not shown)

Abstract: We explore the relationship between mid-infrared (mid-IR) and CO rotational line emission from massive star-forming galaxies, which is one of the tightest scalings in the local universe. We assemble a large set of unresolved and moderately ($\sim 1$ kpc) spatially resolved measurements of CO (1-0) and CO (2-1) intensity, $I_{\rm CO}$, and mid-IR intensity, $I_{\rm MIR}$, at 8, 12, 22, and 24$μ$m.… ▽ More We explore the relationship between mid-infrared (mid-IR) and CO rotational line emission from massive star-forming galaxies, which is one of the tightest scalings in the local universe. We assemble a large set of unresolved and moderately ($\sim 1$ kpc) spatially resolved measurements of CO (1-0) and CO (2-1) intensity, $I_{\rm CO}$, and mid-IR intensity, $I_{\rm MIR}$, at 8, 12, 22, and 24$μ$m. The $I_{\rm CO}$ vs. $I_{\rm MIR}$ relationship is reasonably described by a power law with slopes $0.7{-}1.2$ and normalization $I_{\rm CO} \sim 1$ K km s$^{-1}$ at $I_{\rm MIR} \sim 1$ MJy sr$^{-1}$. Both the slopes and intercepts vary systematically with choice of line and band. The comparison between the relations measured for CO~(1-0) and CO (2-1) allow us to infer that $R_{21} \propto I_{\rm MIR}^{0.2}$, in good agreement with other work. The $8μ$m and $12μ$m bands, with strong PAH features, show steeper CO vs. mid-IR slopes than the $22μ$m and $24μ$m, consistent with PAH emission arising not just from CO-bright gas but also from atomic or CO-dark gas. The CO-to-mid-IR ratio correlates with global galaxy stellar mass ($M_\star$) and anti-correlates with SFR/$M_\star$. At $\sim 1$ kpc resolution, the first four PHANGS-JWST targets show CO to mid-IR relationships that are quantitatively similar to our larger literature sample, including showing the steep CO-to-mid-IR slopes for the JWST PAH-tracing bands, although we caution that these initial data have a small sample size and span a limited range of intensities. △ Less

Submitted 27 December, 2022; v1 submitted 19 December, 2022; originally announced December 2022.

Comments: 29 pages, 13 figures, key quantitative results in Table 3, Accepted as part of a PHANGS-JWST Focus Issue to appear in ApJ

Interacting galaxies in the IllustrisTNG simulations -- IV: Enhanced Supermassive Black Hole Accretion Rates in Post-Merger Galaxies

Authors: Shoshannah Byrne-Mamahit, Maan Hani, Sara Ellison, Salvatore Quai, David Patton

Abstract: We present an analysis of the instantaneous supermassive black hole (SMBH) accretion rates in a collection of 1563 post-merger galaxies drawn from the IllustrisTNG simulation. Our sample consists of galaxies that have experienced a merger in the last simulation snapshot (within ~160 Myrs of coalescence) in the redshift range 0<z<1, with merger stellar mass ratios >1:10 and post-merger stellar mass… ▽ More We present an analysis of the instantaneous supermassive black hole (SMBH) accretion rates in a collection of 1563 post-merger galaxies drawn from the IllustrisTNG simulation. Our sample consists of galaxies that have experienced a merger in the last simulation snapshot (within ~160 Myrs of coalescence) in the redshift range 0<z<1, with merger stellar mass ratios >1:10 and post-merger stellar masses > $10^{10} M_{\odot}$. We find that, on average, the accretion rates of the post-mergers are ~1.7 times higher than in a control sample and that post-mergers are 3-4 times more likely to experience a luminous active galactic nuclei (AGN) phase than isolated galaxies. SMBH accretion rate enhancements persist for ~2 Gyrs after coalescence, significantly exceeding the ~500 Myr lifetime of star formation rate enhancements. We find that the presence of simultaneous enhancements in both the star formation and SMBH accretion rates depends on both the mass ratio of the merger and on the gas mass of the post-merger galaxy. Despite these accretion rate enhancements, only ~35% of post-mergers experience a luminous AGN ($L_{bol}>10^{44}$ erg/s) within 500 Myrs after coalescence, and fewer than 10\% achieve a luminosity in excess of $L_{bol}>10^{45}$ erg/s. Moreover, only ~10\% of the highest luminosity ($L_{bol}>10^{45}$ erg/s) AGN in the IllustrisTNG galaxy sample are recent mergers. Our results are therefore consistent with a picture in which mergers can (but don't always) trigger AGN activity, but where the majority of galaxies hosting high luminosity AGN are not recent mergers. △ Less

Submitted 14 December, 2022; originally announced December 2022.

Comments: This article has been accepted for publication in Monthly Notices of the Royal Astronomical Society

VERTICO III: The Kennicutt-Schmidt relation in Virgo cluster galaxies

Authors: M. J. Jiménez-Donaire, T. Brown, C. D. Wilson, I. D. Roberts, N. Zabel, S. L. Ellison, M. Thorp, V. Villanueva, R. Chown, D. Bisaria, A. D. Bolatto, A. Boselli, B. Catinella, A. Chung, L. Cortese, T. A. Davis, C. D. P. Lagos, B. Lee, L. C. Parker, K. Spekkens, A. R. H. Stevens, J. Sun

Abstract: In this VERTICO science paper we aim to study how the star formation process depends on galactic environment and gravitational interactions in the context of galaxy evolution. We explore the scaling relation between the star formation rate (SFR) surface density and the molecular gas surface density, also known as the Kennicutt-Schmidt (KS) relation, in a subsample of Virgo cluster spiral galaxies.… ▽ More In this VERTICO science paper we aim to study how the star formation process depends on galactic environment and gravitational interactions in the context of galaxy evolution. We explore the scaling relation between the star formation rate (SFR) surface density and the molecular gas surface density, also known as the Kennicutt-Schmidt (KS) relation, in a subsample of Virgo cluster spiral galaxies. We use new ACA and TP observations from the VERTICO-ALMA Large Program at 720pc resolution to resolve the molecular gas content, as traced by the 12CO(2-1) transition, across the disks of 37 spiral galaxies in the Virgo cluster. In combination with archival observations, we estimate the parameters of the KS relation for the entire ensemble of galaxies, and within individual galaxies. We find the KS slope for the entire population to be N=0.97+/-0.07, with a characteristic molecular gas depletion time of 1.86Gyr for our full sample, in agreement with previous work in isolated star-forming galaxies. In individual galaxies, we find KS slopes ranging between 0.69 and 1.40, and typical star formation efficiencies (SFE) that can vary from galaxy to galaxy by a factor of ~4. These galaxy-to-galaxy variations account for ~0.20dex in scatter in the ensemble KS relation, which is characterized by a 0.42dex scatter. We find that the HI-deficient galaxies in the Virgo cluster show a steeper resolved KS relation and lower molecular gas efficiencies than HI-normal cluster galaxies. While the molecular gas content in Virgo cluster galaxies appears to behave similarly to that in isolated galaxies, our VERTICO sample shows that cluster environments play a key role in regulating star formation. The environmental mechanisms affecting the HI galaxy content also have a direct impact in the SFE of molecular gas in cluster galaxies, leading to longer depletion times in HI-deficient members. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 671, A3 (2023)

A New Physical Picture for AGNs Lacking Optical Emission Lines

Authors: Christopher J. Agostino, Samir Salim, Sara L. Ellison, Robert W. Bickley, S. M. Faber

Abstract: In this work, we use ~500 low-redshift (z ~ 0.1) X-ray AGNs observed by XMM-Newton and SDSS to investigate the prevalence and nature of AGNs that apparently lack optical emission lines (``optically dull AGNs''). Although 1/4 of spectra appear absorption-line dominated in visual assessment, line extraction with robust continuum subtraction from the MPA/JHU catalog reveals usable [OIII] measurements… ▽ More In this work, we use ~500 low-redshift (z ~ 0.1) X-ray AGNs observed by XMM-Newton and SDSS to investigate the prevalence and nature of AGNs that apparently lack optical emission lines (``optically dull AGNs''). Although 1/4 of spectra appear absorption-line dominated in visual assessment, line extraction with robust continuum subtraction from the MPA/JHU catalog reveals usable [OIII] measurements in 98% of the sample, allowing us to study [OIII]-underluminous AGNs together with more typical AGNs in the context of the L$_{\mathrm{[OIII]}}$--L$_{X}$ relation. We find that ``optically dull AGNs'' do not constitute a distinct population of AGNs. Instead, they are the [OIII]-underluminous tail of a single, unimodal L$_{\mathrm{[OIII]}}$--L$_{X}$ relation that has substantial scatter (0.6 dex). We find the degree to which an AGN is underluminous in [OIII] correlates with the specific SFR or D$_{4000}$ index of the host, which are both linked to the molecular gas fraction. Thus the emerging physical picture for the large scatter seems to involve the gas content of the narrow-line region. We find no significant role for previously proposed scenarios for the presence of optically dull AGNs, such as host dilution or dust obscuration. Despite occasionally weak lines in SDSS spectra, >80% of X-ray AGNs are identified as such with the BPT diagram. >90% are classified as AGNs based only on [NII]/H$α$, providing more complete AGN samples when [OIII] or H$β$ are weak. X-ray AGNs with LINER spectra obey essentially the same \lxo\ relation as Seyfert 2s, suggesting their line emission is produced by AGN activity. △ Less

Submitted 14 November, 2022; originally announced November 2022.

Comments: 21 pages, 14 figures, submitted to ApJ. Comments welcome

LoTSS Jellyfish Galaxies IV: Enhanced Star Formation on the Leading Half of Cluster Galaxies and Gas Compression in IC3949

Authors: Ian D. Roberts, Maojin Lang, Daria Trotsenko, Ashley Bemis, Sara L. Ellison, Lihwai Lin, Hsi-An Pan, Alessandro Ignesti, Sarah Leslie, Reinout J. van Weeren

Abstract: With MaNGA integral field spectroscopy, we present a resolved analysis of star formation for 29 jellyfish galaxies in nearby clusters, identified from radio continuum imaging taken by the Low Frequency Array. Simulations predict enhanced star formation on the "leading half" of galaxies undergoing ram pressure stripping, and in this work we report observational evidence for this elevated star forma… ▽ More With MaNGA integral field spectroscopy, we present a resolved analysis of star formation for 29 jellyfish galaxies in nearby clusters, identified from radio continuum imaging taken by the Low Frequency Array. Simulations predict enhanced star formation on the "leading half" of galaxies undergoing ram pressure stripping, and in this work we report observational evidence for this elevated star formation. The dividing line (through the galaxy center) that maximizes this star formation enhancement is systematically tied to the observed direction of the ram pressure stripped tail, suggesting a physical connection between ram pressure and this star formation enhancement. We also present a case study on the distribution of molecular gas in one jellyfish galaxy from our sample, IC3949, using ALMA CO J=1-0, HCN J=1-0, and HCO$^+$ J=1-0 observations from the ALMaQUEST survey. The $\mathrm{H_2}$ depletion time (as traced by CO) in IC3949 ranges from $\sim\!1\,\mathrm{Gyr}$ in the outskirts of the molecular gas disk to $\sim\!11\,\mathrm{Gyr}$ near the galaxy center. IC3949 shows a clear region of enhanced star formation on the leading half of the galaxy where the average depletion time is $\sim\!2.7\,\mathrm{Gyr}$, in line with the median value for the galaxy on the whole. Dense gas tracers, HCN and HCO$^+$, are only detected at the galaxy center and on the leading half of IC3949. Our results favour a scenario in which ram pressure compresses the interstellar medium, promoting the formation of molecular gas that in turn fuels a localized increase of star formation. △ Less

Submitted 28 October, 2022; originally announced October 2022.

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

Resolved Molecular Gas Observations of MaNGA Post-starbursts Reveal a Tumultuous Past

Authors: Justin Atsushi Otter, Kate Rowlands, Katherine Alatalo, Ho-Hin Leung, Vivienne Wild, Yuanze Luo, Andreea O. Petric, Elizaveta Sazonova, David V. Stark, Timothy Heckman, Timothy A. Davis, Sara Ellison, K. Decker French, William Baker, Asa F. L. Bluck, Lauranne Lanz, Lihwai Lin, Charles Liu, Carlos López Cobá, Karen L. Masters, Preethi Nair, Hsi-an Pan, Rogemar A. Riffel, Jillian M. Scudder, Adam Smercina , et al. (2 additional authors not shown)

Abstract: Post-starburst galaxies (PSBs) have recently and rapidly quenched their star-formation, thus they are an important way to understand how galaxies transition from star-forming late-types to quiescent early-types. The recent discovery of large cold gas reservoirs in PSBs calls into question the theory that galaxies must lose their gas to become quiescent. Optical Integral Field Spectroscopy (IFS) su… ▽ More Post-starburst galaxies (PSBs) have recently and rapidly quenched their star-formation, thus they are an important way to understand how galaxies transition from star-forming late-types to quiescent early-types. The recent discovery of large cold gas reservoirs in PSBs calls into question the theory that galaxies must lose their gas to become quiescent. Optical Integral Field Spectroscopy (IFS) surveys have revealed two classes of PSBs: central PSBs with central quenching regions and ring PSBs with quenching in their outskirts. We analyze a sample of 13 nearby (z < 0.1) PSBs with spatially resolved optical IFS data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey and matched resolution Atacama Large (sub-)Millimeter Array (ALMA) observations of $^{12}$CO(1-0). Disturbed stellar kinematics in 7/13 of our PSBs and centrally concentrated molecular gas is consistent with a recent merger for most of our sample. In galaxies without merger evidence, alternate processes may funnel gas inwards and suppress star-formation, which may include outflows, stellar bars, and minor mergers or interactions. The star-formation efficiencies of the post-starburst regions in nearly half our galaxies are suppressed while the gas fractions are consistent with star-forming galaxies. AGN feedback may drive this stabilization, and we observe AGN-consistent emission in the centers of 5/13 galaxies. Finally, our central and ring PSBs have similar properties except the ionized and molecular gas in central PSBs is more disturbed. Overall, the molecular gas in our PSBs tends to be compact and highly disturbed, resulting in concentrated gas reservoirs unable to form stars efficiently. △ Less

Submitted 21 October, 2022; originally announced October 2022.

Comments: Accepted for publication in The Astrophysical Journal

VERTICO IV: Environmental Effects on the Gas Distribution and Star Formation Efficiency of Virgo Cluster Spirals

Authors: Vicente Villanueva, Alberto D. Bolatto, Stuart Vogel, Tobias Brown, Christine D. Wilson, Nikki Zabel, Sara Ellison, Adam R. H. Stevens, Maria Jesus Jimenez Donaire, Kristine Spekkens, Mallory Thorp, Timothy A. Davis, Laura C. Parker, Ian D. Roberts, Dhruv Bisaria, Alessandro Boselli, Barbara Catinella, Aeree Chung, Luca Cortese, Bumhyun Lee, Adam Watts

Abstract: We measure the molecular-to-atomic gas ratio, $R_{\rm mol}$, and the star formation rate (SFR) per unit molecular gas mass, SFE$_{\rm mol}$, in 38 nearby galaxies selected from the Virgo Environment Traced in CO (VERTICO) survey. We determine their scale-lengths for the molecular and stellar components and find a roughly 3:5 ratio between them compared to $\sim$1:1 in field galaxies, indicating th… ▽ More We measure the molecular-to-atomic gas ratio, $R_{\rm mol}$, and the star formation rate (SFR) per unit molecular gas mass, SFE$_{\rm mol}$, in 38 nearby galaxies selected from the Virgo Environment Traced in CO (VERTICO) survey. We determine their scale-lengths for the molecular and stellar components and find a roughly 3:5 ratio between them compared to $\sim$1:1 in field galaxies, indicating that the CO emission is more centrally concentrated than the stars. We compute $R_{\rm mol}$ as a function of different physical quantities. While the spatially-resolved $R_{\rm mol}$ on average decreases with increasing radius, we find that the mean molecular-to-atomic gas ratio within the stellar effective radius $R_{\rm e}$, $R_{\rm mol}(r<R_{\rm e})$, shows a systematic increase with the level of H$_{\rm I}$, truncation and/or asymmetry (H$_{\rm I}$ perturbation). Analysis of the molecular- and the atomic-to-stellar mass ratios within $R_{\rm e}$, $R^{\rm mol}_{\star}(r<R_{\rm e})$ and $R^{\rm atom}_{\star}(r<R_{\rm e})$, shows that VERTICO galaxies have increasingly lower $R^{\rm atom}_{\star}(r<R_{\rm e})$ for larger levels of H$_{\rm I}$perturbation (compared to field galaxies matched in stellar mass), but no significant change in $R^{\rm mol}_{\star}(r<R_{\rm e})$. We also measure a clear systematic decrease of the SFE$_{\rm mol}$ within $R_{\rm e}$, SFE$_{\rm mol}(r<R_{\rm e})$, with increasingly perturbed H$_{\rm I}$. Therefore, compared to galaxies from the field, VERTICO galaxies are more compact in CO emission in relation to their stellar distribution, but increasingly perturbed atomic gas increases their $R_{\rm mol}$ and decreases the efficiency with which their molecular gas forms stars. (abridged) △ Less

Submitted 1 November, 2022; v1 submitted 11 October, 2022; originally announced October 2022.

The metallicity's fundamental dependence on both local and global galactic quantities

Authors: William M. Baker, Roberto Maiolino, Francesco Belfiore, Mirko Curti, Asa F. L. Bluck, Lihwai Lin, Sara L. Ellison, Mallory Thorp, Hsi-An Pan

Abstract: We study the scaling relations between gas-phase metallicity, stellar mass surface density ($Σ_*$), star formation rate surface density ($Σ_{SFR}$), and molecular gas surface density ($Σ_{H_2}$) in local star-forming galaxies on scales of a kpc. We employ optical integral field spectroscopy from the MaNGA survey, and ALMA data for a subset of MaNGA galaxies. We use Partial Correlation Coefficients… ▽ More We study the scaling relations between gas-phase metallicity, stellar mass surface density ($Σ_*$), star formation rate surface density ($Σ_{SFR}$), and molecular gas surface density ($Σ_{H_2}$) in local star-forming galaxies on scales of a kpc. We employ optical integral field spectroscopy from the MaNGA survey, and ALMA data for a subset of MaNGA galaxies. We use Partial Correlation Coefficients and Random Forest regression to determine the relative importance of local and global galactic properties in setting the gas-phase metallicity. We find that the local metallicity depends primarily on $Σ_*$ (the resolved mass-metallicity relation, rMZR), and has a secondary anti-correlation with $Σ_{SFR}$ (i.e. a spatially-resolved version of the 'Fundamental Metallicity Relation', rFMR). We find that $Σ_{H_2}$ is less important than $Σ_{SFR}$ in determining the local metallicity. This result indicates that gas accretion, resulting in local metallicity dilution and local boosting of star formation, is unlikely to be the primary origin of the rFMR. The local metallicity depends also on the global properties of galaxies. We find a strong dependence on the total stellar mass ($M_*$) and a weaker (inverse) dependence on the total SFR. The global metallicity scaling relations, therefore, do not simply stem out of their resolved counterparts; global properties and processes, such as the global gravitational potential well, galaxy-scale winds and global redistribution/mixing of metals, likely contribute to the local metallicity, in addition to local production and retention. △ Less

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

Comments: 22 pages, 24 figures, Accepted for publication in MNRAS

Galaxy mergers can rapidly shut down star formation

Authors: Sara L. Ellison, Scott Wilkinson, Joanna Woo, Ho-Hin Leung, Vivienne Wild, Robert W. Bickley, David R. Patton, Salvatore Quai, Stephen Gwyn

Abstract: Galaxy mergers trigger both star formation and accretion onto the central supermassive black hole. As a result of subsequent energetic feedback processes, it has long been proposed that star formation may be promptly extinguished in galaxy merger remnants. However, this prediction of widespread, rapid quenching in late stage mergers has been recently called into question with modern simulations an… ▽ More Galaxy mergers trigger both star formation and accretion onto the central supermassive black hole. As a result of subsequent energetic feedback processes, it has long been proposed that star formation may be promptly extinguished in galaxy merger remnants. However, this prediction of widespread, rapid quenching in late stage mergers has been recently called into question with modern simulations and has never been tested observationally. Here we perform the first empirical assessment of the long-predicted end phase in the merger sequence. Based on a sample of ~500 post-mergers identified from the Ultraviolet Near Infrared Optical Northern Survey (UNIONS), we show that the frequency of post-merger galaxies that have rapidly shutdown their star formation following a previous starburst is 30-60 times higher than expected from a control sample of non-merging galaxies. No such excess is found in a sample of close galaxy pairs, demonstrating that mergers can indeed lead to a rapid halt to star formation, but that this process only manifests after coalescence. △ Less

Submitted 15 September, 2022; originally announced September 2022.

Comments: Accepted for publication in MNRAS Letters

The ALMaQUEST Survey X: What powers merger induced star formation?

Authors: Mallory D. Thorp, Sara L. Ellison, Hsi-An Pan, Lihwai Lin, David R. Patton, Asa F. L. Bluck, Dan Walters, Jillian M. Scudder

Abstract: Galaxy mergers are known to trigger both extended and central star formation. However, what remains to be understood is whether this triggered star formation is facilitated by enhanced star formation efficiencies, or an abundance of molecular gas fuel. This work presents spatially resolved measurements of CO emission collected with the Atacama Large Millimetre Array (ALMA) for 20 merging galaxies… ▽ More Galaxy mergers are known to trigger both extended and central star formation. However, what remains to be understood is whether this triggered star formation is facilitated by enhanced star formation efficiencies, or an abundance of molecular gas fuel. This work presents spatially resolved measurements of CO emission collected with the Atacama Large Millimetre Array (ALMA) for 20 merging galaxies (either pairs or post-mergers) selected from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. Eleven additional merging galaxies are selected from the ALMA MaNGA QUEnching and STar formation (ALMaQUEST) survey, resulting in a set of 31 mergers at various stages of interaction and covering a broad range of star formation rates (SFR). We investigate galaxy-to-galaxy variations in the resolved Kennicutt-Schmidt relation (rKS: $Σ_{H_2}$ vs. $Σ_{SFR}$), the resolved molecular gas main sequence (rMGMS: $Σ_{\star}$ vs. $Σ_{H_2}$), and the resolved star-forming main sequence (rSFMS: $Σ_{\star}$ vs. $Σ_{SFR}$). We quantify offsets from these resolved relations to determine if star formation rate, molecular gas fraction, and/or star formation efficiency (SFE) is enhanced in different regions of an individual galaxy. By comparing offsets in all three parameters we can discern whether gas fraction or SFE powers an enhanced $Σ_{SFR}$. We find that merger-induced star formation can be driven by a variety of mechanisms, both within a galaxy and between different mergers, regardless of interaction stage. △ Less

Submitted 12 August, 2022; originally announced August 2022.

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

The merger fraction of post-starburst galaxies in UNIONS

Authors: Scott Wilkinson, Sara L. Ellison, Connor Bottrell, Robert W. Bickley, Stephen Gwyn, Jean-Charles Cuillandre, Vivienne Wild

Abstract: Post-starburst (PSB) galaxies are defined as having experienced a recent burst of star formation, followed by a prompt truncation in further activity. Identifying the mechanism(s) causing a galaxy to experience a post-starburst phase therefore provides integral insight into the causes of rapid quenching. Galaxy mergers have long been proposed as a possible post-starburst trigger. Effectively testi… ▽ More Post-starburst (PSB) galaxies are defined as having experienced a recent burst of star formation, followed by a prompt truncation in further activity. Identifying the mechanism(s) causing a galaxy to experience a post-starburst phase therefore provides integral insight into the causes of rapid quenching. Galaxy mergers have long been proposed as a possible post-starburst trigger. Effectively testing this hypothesis requires a large spectroscopic galaxy survey to identify the rare PSBs as well as high quality imaging and robust morphology metrics to identify mergers. We bring together these critical elements by selecting PSBs from the overlap of the Sloan Digital Sky Survey and the Canada-France Imaging Survey and applying a suite of classification methods: non-parametric morphology metrics such as asymmetry and Gini-M20, a convolutional neural network trained to identify post-merger galaxies, and visual classification. This work is therefore the largest and most comprehensive assessment of the merger fraction of PSBs to date. We find that the merger fraction of PSBs ranges from 19% to 42% depending on the merger identification method and details of the PSB sample selection. These merger fractions represent an excess of 3-46x relative to non-PSB control samples. Our results demonstrate that mergers play a significant role in generating PSBs, but that other mechanisms are also required. However, applying our merger identification metrics to known post-mergers in the IllustrisTNG simulation shows that ~70% of recent post-mergers (<200 Myr) would not be detected. Thus, we cannot exclude the possibility that nearly all post-starburst galaxies have undergone a merger in their recent past. △ Less

Submitted 8 July, 2022; originally announced July 2022.

Comments: 21 pages, 16 figures

The ALMaQUEST survey. VIII. What causes the velocity discrepancy between CO and H$α$ rotation curves in galaxies?

Authors: Yung-Chau Su, Lihwai Lin, Hsi-An Pan, Carlos López Cobá, Bau-Ching Hsieh, Sebastián F. Sánchez, Mallory D. Thorp, Martin Bureau, Sara L. Ellison

Abstract: We compare the CO(1-0) and H$α$ kinematics in 34 nearby galaxies, selected from the ALMaQUEST and EDGE-CALIFA surveys. We use 3-D Barolo, a 3D tilted ring model, to derive the CO and H$α$ rotation curves. Before comparing rotation curves in the 34 nearby galaxies, we found systematics between the MaNGA and the CALIFA data using eight MaNGA-CALIFA overlapping galaxies. We assume the rotation curves… ▽ More We compare the CO(1-0) and H$α$ kinematics in 34 nearby galaxies, selected from the ALMaQUEST and EDGE-CALIFA surveys. We use 3-D Barolo, a 3D tilted ring model, to derive the CO and H$α$ rotation curves. Before comparing rotation curves in the 34 nearby galaxies, we found systematics between the MaNGA and the CALIFA data using eight MaNGA-CALIFA overlapping galaxies. We assume the rotation curves based on the MaNGA data are accurate and made the corresponding correction to the CALIFA data. Our result shows that $\sim$56% (19/34) of our galaxies present slower H$α$ rotation curves compared to the CO rotation curves, with a median value of 6.5 km/s. The remaining galaxies (15/34) show consistent CO-H$α$ rotation velocity within uncertainties. As a result, the Ha rotation may underestimate the total dynamical mass by ~6% for a circular velocity of 200 km/s (the median value in our sample). Furthermore, the velocity difference between the CO and H$α$ rotation velocity is found to correlate with the velocity dispersion difference between CO and H$α$, suggesting that the gas pressure plays a role in the velocity discrepancy. After incorporating the effect of pressure support due to the turbulent gas motion to our sample, the median value of the velocity differences decreases to 1.9 km/s, which in turn reduces the underestimation of dynamical mass to $\sim$2%. Finally, we also investigate the role that the extra-planar diffuse ionized gas (eDIG) plays in the CO-H$α$ velocity discrepancy. △ Less

Submitted 14 June, 2022; v1 submitted 13 June, 2022; originally announced June 2022.

Star formation characteristics of CNN-identified post-mergers in the Ultraviolet Near Infrared Optical Northern Survey (UNIONS)

Authors: Robert W. Bickley, Sara L. Ellison, David R. Patton, Connor Bottrell, Stephen Gwyn, Michael J. Hudson

Abstract: The importance of the post-merger epoch in galaxy evolution has been well-documented, but post-mergers are notoriously difficult to identify. While the features induced by mergers can sometimes be distinctive, they are frequently missed by visual inspection. In addition, visual classification efforts are highly inefficient because of the inherent rarity of post-mergers (~1% in the low-redshift Uni… ▽ More The importance of the post-merger epoch in galaxy evolution has been well-documented, but post-mergers are notoriously difficult to identify. While the features induced by mergers can sometimes be distinctive, they are frequently missed by visual inspection. In addition, visual classification efforts are highly inefficient because of the inherent rarity of post-mergers (~1% in the low-redshift Universe), and non-parametric statistical merger selection methods do not account for the diversity of post-mergers or the environments in which they appear. To address these issues, we deploy a convolutional neural network (CNN) which has been trained and evaluated on realistic mock observations of simulated galaxies from the IllustrisTNG simulations, to galaxy images from the Canada France Imaging Survey (CFIS), which is part of the Ultraviolet Near Infrared Optical Northern Survey (UNIONS). We present the characteristics of the galaxies with the highest CNN-predicted post-merger certainties, as well as a visually confirmed subset of 699 post-mergers. We find that post-mergers with high CNN merger probabilities (p(x)>0.8) have an average star formation rate that is 0.1 dex higher than a mass- and redshift-matched control sample. The SFR enhancement is even greater in the visually confirmed post-merger sample, a factor of two higher than the control sample. △ Less

Submitted 27 May, 2022; originally announced May 2022.

Comments: 14 pages, 9 figures, visually confirmed post-merger catalog available via MNRAS

Galaxy Pairs in the Sloan Digital Sky Survey XV: Properties of Ionised Outflows

Authors: William Matzko, Shobita Satyapal, Sara L. Ellison, Remington O. Sexton, Nathan J. Secrest, Gabriela Canalizo, Laura Blecha, David R. Patton, Jillian M. Scudder

Abstract: Powerful outflows are thought to play a critical role in galaxy evolution and black hole growth. We present the first large-scale systematic study of ionised outflows in paired galaxies and post-mergers compared to a robust control sample of isolated galaxies. We isolate the impact of the merger environment to determine if outflow properties depend on merger stage. Our sample contains $\sim$4,000… ▽ More Powerful outflows are thought to play a critical role in galaxy evolution and black hole growth. We present the first large-scale systematic study of ionised outflows in paired galaxies and post-mergers compared to a robust control sample of isolated galaxies. We isolate the impact of the merger environment to determine if outflow properties depend on merger stage. Our sample contains $\sim$4,000 paired galaxies and $\sim$250 post-mergers in the local universe ($0.02 \leq z \leq 0.2$) from the SDSS DR 7 matched in stellar mass, redshift, local density of galaxies, and [OIII] $λ$5007 luminosity to a control sample of isolated galaxies. By fitting the [OIII] $λ$5007 line, we find ionised outflows in $\sim$15 per cent of our entire sample. Outflows are much rarer in star-forming galaxies compared to AGN, and outflow incidence and velocity increase with [OIII] $λ$5007 luminosity. Outflow incidence is significantly elevated in the optical+mid-infrared selected AGN compared to purely optical AGN; over 60 per cent show outflows at the highest luminosities ($L_{\mathrm{[OIII] \lambda5007}}$ $\gtrsim$ 10$^{42}$ erg s$^{-1}$), suggesting mid-infrared AGN selection favours galaxies with powerful outflows, at least for higher [OIII] $λ$5007 luminosities. However, we find no statistically significant difference in outflow incidence, velocity, and luminosity in mergers compared to isolated galaxies, and there is no dependence on merger stage. Therefore, while interactions are predicted to drive gas inflows and subsequently trigger nuclear star formation and accretion activity, when the power source of the outflow is controlled for, the merging environment has no further impact on the large-scale ionised outflows as traced by [OIII] $\lambda5007$. △ Less

Submitted 26 May, 2022; originally announced May 2022.

Comments: 18 pages, 15 figures, to be published in Monthly Notices of the Royal Astronomical Society

VERTICO II: effects of HI-identified environmental mechanisms on molecular gas

Authors: Nikki Zabel, Toby Brown, Christine D. Wilson, Timothy A. Davis, Luca Cortese, Laura C. Parker, Alessandro Boselli, Barbara Catinella, Ryan Chown, Aeree Chung, Tirna Deb, Sara L. Ellison, María J. Jiménez-Donaire, Bumhyun Lee, Ian D. Roberts, Kristine Spekkens, Adam R. H. Stevens, Mallory Thorp, Stephanie Tonnesen, Vicente Villanueva

Abstract: In this VERTICO early science paper we explore in detail how environmental mechanisms, identified in HI, affect the resolved properties of molecular gas reservoirs in cluster galaxies. The molecular gas is probed using ALMA ACA (+TP) observations of 12CO(2-1) in 51 spiral galaxies in the Virgo cluster (of which 49 are detected), all of which are included in the VIVA HI survey. The sample spans a s… ▽ More In this VERTICO early science paper we explore in detail how environmental mechanisms, identified in HI, affect the resolved properties of molecular gas reservoirs in cluster galaxies. The molecular gas is probed using ALMA ACA (+TP) observations of 12CO(2-1) in 51 spiral galaxies in the Virgo cluster (of which 49 are detected), all of which are included in the VIVA HI survey. The sample spans a stellar mass range of 9 < log M*/Msol < 11. We study molecular gas radial profiles, isodensity radii, and surface densities as a function of galaxy HI deficiency and morphology. There is a weak correlation between global HI and H2 deficiencies, and resolved properties of molecular gas correlate with HI deficiency: galaxies that have large HI deficiencies have relatively steep and truncated molecular gas radial profiles, which is due to the removal of low-surface density molecular gas on the outskirts. Therefore, while the environmental mechanisms observed in HI also affect molecular gas reservoirs, there is only a moderate reduction of the total amount of molecular gas. △ Less

Submitted 12 January, 2023; v1 submitted 11 May, 2022; originally announced May 2022.

Comments: Published in ApJ. 22 pages, 6 figures, 1 table, 1 appendix. Erratum accepted for publication in ApJ

The evolution of the Si IV content in the Universe from the epoch of reionization to cosmic noon

Authors: V. D'Odorico, K. Finlator, S. Cristiani, G. Cupani, S. Perrotta, F. Calura, M. Cènturion, G. Becker, T. A. M. Berg, S. Lopez, S. Ellison, E. Pomante

Abstract: We investigate the abundance and distribution of metals in the high-redshift intergalactic medium and circum-galactic medium through the analysis of a sample of almost 600 SiIV absorption lines detected in high and intermediate resolution spectra of 147 quasars. The evolution of the number density of SiIV lines, the column density distribution function and the cosmic mass density are studied in th… ▽ More We investigate the abundance and distribution of metals in the high-redshift intergalactic medium and circum-galactic medium through the analysis of a sample of almost 600 SiIV absorption lines detected in high and intermediate resolution spectra of 147 quasars. The evolution of the number density of SiIV lines, the column density distribution function and the cosmic mass density are studied in the redshift interval 1.7 <= z <= 6.2 and for log N(SiIV) >= 12.5. All quantities show a rapid increase between z~6 and z< 5 and then an almost constant behaviour to z~2 in very good agreement with what is already observed for CIV absorption lines. The present results are challenging for numerical simulations: when simulations reproduce our SiIV results, they tend to underpredict the properties of CIV, and when the properties of CIV are reproduced, the number of strong SiIV lines (log N(SiIV) > 14) is overpredicted. △ Less

Submitted 24 February, 2022; originally announced February 2022.

Comments: 14 pages, 8 figures, accepted for publication in MNRAS

Quenching Timescales in the IllustrisTNG Simulation

Authors: Dan Walters, Joanna Woo, Sara L. Ellison

Abstract: The timescales for galaxy quenching offer clues to its underlying physical drivers. We investigate central galaxy quenching timescales in the IllustrisTNG 100-1 simulation, their evolution over time, and the pre-quenching properties of galaxies that predict their quenching timescales. Defining quenching duration $τ_q$ as the time between crossing sSFR thresholds, we find that $\sim$40% of galaxies… ▽ More The timescales for galaxy quenching offer clues to its underlying physical drivers. We investigate central galaxy quenching timescales in the IllustrisTNG 100-1 simulation, their evolution over time, and the pre-quenching properties of galaxies that predict their quenching timescales. Defining quenching duration $τ_q$ as the time between crossing sSFR thresholds, we find that $\sim$40% of galaxies quench rapidly with $τ_q<$1 Gyr, but a substantial tail of galaxies can take up to 10 Gyr to quench. Furthermore, 29% of galaxies that left the star forming main sequence (SFMS) more than 2 Gyr ago never fully quench by $z=0$. While the median $τ_q$ is fairly constant with epoch, the rate of galaxies leaving the SFMS increases steadily over cosmic time, with the rate of slow quenchers being dominant around $z\sim2$ to 0.7. Compared to fast quenchers ($τ_q<$1 Gyr), slow-quenching galaxies ($τ_q>$1 Gyr) were more massive, had more massive black holes, had larger stellar radii and accreted gas with higher specific angular momentum (AM) prior to quenching. These properties evolve little by $z=0$, except for the accreting gas AM for fast quenchers, which reaches the same high AM as the gas in slow quenchers. By $z=0$, slow quenchers also have residual star formation in extended gas rings. Using the expected relationship between stellar age gradient and $τ_q$ for inside-out quenching we find agreement with MaNGA IFU observations. Our results suggest the accreting gas AM and potential well depth determine the quenching timescale. △ Less

Submitted 31 January, 2022; originally announced February 2022.

Comments: 18 pages, 18 figures, accepted for publication in MNRAS

The quenching of galaxies, bulges, and disks since cosmic noon: A machine learning approach for identifying causality in astronomical data

Authors: Asa F. L. Bluck, Roberto Maiolino, Simcha Brownson, Christopher J. Conselice, Sara L. Ellison, Joanna M. Piotrowska, Mallory D. Thorp

Abstract: We present an analysis of the quenching of star formation in galaxies, bulges, and disks throughout the bulk of cosmic history, from $z=2-0$. We utilise observations from the SDSS and MaNGA at low redshifts. We complement these data with observations from CANDELS at high redshifts. Additionally, we compare the observations to detailed predictions from the LGalaxies semi-analytic model. To analyse… ▽ More We present an analysis of the quenching of star formation in galaxies, bulges, and disks throughout the bulk of cosmic history, from $z=2-0$. We utilise observations from the SDSS and MaNGA at low redshifts. We complement these data with observations from CANDELS at high redshifts. Additionally, we compare the observations to detailed predictions from the LGalaxies semi-analytic model. To analyse the data, we developed a machine learning approach utilising a Random Forest classifier. We first demonstrate that this technique is extremely effective at extracting causal insight from highly complex and inter-correlated model data, before applying it to various observational surveys. Our primary observational results are as follows: At all redshifts studied in this work, we find bulge mass to be the most predictive parameter of quenching, out of the photometric parameter set (incorporating bulge mass, disk mass, total stellar mass, and $B/T$ structure). Moreover, we also find bulge mass to be the most predictive parameter of quenching in both bulge and disk structures, treated separately. Hence, intrinsic galaxy quenching must be due to a stable mechanism operating over cosmic time, and the same quenching mechanism must be effective in both bulge and disk regions. Despite the success of bulge mass in predicting quenching, we find that central velocity dispersion is even more predictive (when available in spectroscopic data sets). In comparison to the LGalaxies model, we find that all of these observational results may be consistently explained through quenching via preventative `radio-mode' active galactic nucleus (AGN) feedback. Furthermore, many alternative quenching mechanisms (including virial shocks, supernova feedback, and morphological stabilisation) are found to be inconsistent with our observational results and those from the literature. △ Less

Submitted 19 January, 2022; originally announced January 2022.

Comments: Accepted to A&A; 40 pages; 20 figures

Journal ref: A&A 659, A160 (2022)

ALMaQUEST -- VII: Star Formation Scaling Relations of Green Valley Galaxies

Authors: Lihwai Lin, Sara L. Ellison, Hsi-An Pan, Mallory Thorp, Po-Chieh Yu, Francesco Belfiore, Bau-Ching Hsieh, Roberto Maiolino, S. Ramya, Sebastian Sanchez, Yung-Chau Su

Abstract: We utilize the ALMA-MaNGA QUEnch and STar formation (ALMaQUEST) survey to investigate the kpc-scale scaling relations, presented as the resolved star forming main sequence (rSFMS: $Σ_{\rm SFR}$ vs. $Σ_{*}$), the resolved Schmidt-Kennicutt relation (rSK: $Σ_{\rm SFR}$ vs. $Σ_{\rm H_{2}}$), and the resolved molecular gas main sequence (rMGMS: $Σ_{\rm H_{2}}$ vs. $Σ_{*}$), for 11478 star-forming and… ▽ More We utilize the ALMA-MaNGA QUEnch and STar formation (ALMaQUEST) survey to investigate the kpc-scale scaling relations, presented as the resolved star forming main sequence (rSFMS: $Σ_{\rm SFR}$ vs. $Σ_{*}$), the resolved Schmidt-Kennicutt relation (rSK: $Σ_{\rm SFR}$ vs. $Σ_{\rm H_{2}}$), and the resolved molecular gas main sequence (rMGMS: $Σ_{\rm H_{2}}$ vs. $Σ_{*}$), for 11478 star-forming and 1414 retired spaxels (oversampled by a factor of $\sim20$) located in 22 green valley (GV) and 12 main sequence (MS) galaxies. For a given galaxy type (MS or GV), the retired spaxels are found to be offset from the sequences formed by the star-forming spaxels on the rSFMS, rSK, and rMGMS planes, toward lower absolute values of sSFR, SFE, and $f_{\rm H_{2}}$ by $\sim$ 1.1, 0.6, and 0.5 dex. The scaling relations for GV galaxies are found to be distinct from that of the MS galaxies, even if the analyses are restricted to the star-forming spaxels only. It is found that for star-forming spaxels, sSFR, SFE, and $f_{\rm H_{2}}$ in GV galaxies are reduced by $\sim$0.36, 0.14, and 0.21 dex, respectively, compared to those in MS galaxies. Therefore, the suppressed sSFR/SFE/$f_{\rm H_{2}}$ in GV galaxies are associated with not only an increased proportion of retired regions in GV galaxies but also a depletion of these quantities in star-forming regions. Finally, the reduction of SFE and $f_{\rm H_{2}}$ in GV galaxies relative to MS galaxies is seen in both bulge and disk regions (albeit with larger uncertainties), suggesting that statistically, quenching in the GV population may persist from the inner to the outer regions. △ Less

Submitted 14 January, 2022; originally announced January 2022.

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

The ALMaQUEST Survey IX: The nature of the resolved star forming main sequence

Authors: William M. Baker, Roberto Maiolino, Asa F. L. Bluck, Lihwai Lin, Sara L. Ellison, Francesco Belfiore, Hsi-An Pan, Mallory Thorp

Abstract: We investigate the nature of the scaling relations between the surface density of star formation rate ($Σ_{\rm SFR}$), stellar mass ($Σ_*$), and molecular gas mass ($Σ_{\rm H_2}$), aiming at distinguishing between the relations that are primary, i.e. more fundamental, and those which are instead an indirect by-product of the other relations. We use the ALMaQUEST survey and analyse the data by usin… ▽ More We investigate the nature of the scaling relations between the surface density of star formation rate ($Σ_{\rm SFR}$), stellar mass ($Σ_*$), and molecular gas mass ($Σ_{\rm H_2}$), aiming at distinguishing between the relations that are primary, i.e. more fundamental, and those which are instead an indirect by-product of the other relations. We use the ALMaQUEST survey and analyse the data by using both partial correlations and Random Forest regression techniques. We unambiguously find that the strongest intrinsic correlation is between $Σ_{\rm SFR}$ and $Σ_{\rm H_2}$ (i.e. the resolved Schmidt-Kennicutt relation), followed by the correlation between $Σ_{\rm H_2}$ and $Σ_*$ (resolved Molecular Gas Main Sequence, rMGMS). Once these two correlations are taken into account, we find that there is no evidence for any intrinsic correlation between $Σ_{\rm SFR}$ and $Σ_*$, implying that SFR is entirely driven by the amount of molecular gas, while its dependence on stellar mass (i.e. the resolved Star Forming Main Sequence, rSFMS) simply emerges as a consequence of the relationship between molecular gas and stellar mass. △ Less

Submitted 10 January, 2022; originally announced January 2022.

Comments: Accepted for publication in MNRAS, 7 pages, 4 figures

The combined and respective roles of imaging and stellar kinematics in identifying galaxy merger remnants

Authors: Connor Bottrell, Maan Hani, Hossen Teimoorinia, David R. Patton, Sara L. Ellison

Abstract: One of the central challenges to establishing the role of mergers in galaxy evolution is the selection of pure and complete merger samples in observations. In particular, while large and reasonably pure interacting galaxy pair samples can be obtained with relative ease via spectroscopic criteria, automated selection of post-coalescence merger remnants is restricted to the physical characteristics… ▽ More One of the central challenges to establishing the role of mergers in galaxy evolution is the selection of pure and complete merger samples in observations. In particular, while large and reasonably pure interacting galaxy pair samples can be obtained with relative ease via spectroscopic criteria, automated selection of post-coalescence merger remnants is restricted to the physical characteristics of remnants alone. Furthermore, such selection has predominantly focused on imaging data -- whereas kinematic data may offer a complimentary basis for identifying merger remnants. Therefore, we examine the theoretical utility of both the morphological and kinematic features of merger remnants in distinguishing galaxy merger remnants from other galaxies. Deep classification models are calibrated and evaluated using idealized synthetic images and line-of-sight stellar velocity maps of a heterogeneous population of galaxies and merger remnants from the TNG100 cosmological hydrodynamical simulation. We show that even idealized stellar kinematic data has limited utility compared to imaging and under-performs by $2.1\%\pm0.5\%$ in completeness and $4.7\%\pm0.4\%$ in purity for our fiducial model architecture. Combining imaging and stellar kinematics offers a small boost in completeness (by $1.8\%\pm0.4\%$, compared to $92.7\%\pm0.2\%$ from imaging alone) but no change in purity ($0.1\%\pm0.3\%$ improvement compared to $92.7\%\pm0.2\%$, evaluated with equal numbers of merger remnant and non-remnant control galaxies). Classification accuracy of all models is particularly sensitive to physical companions at separations $\lesssim40$ kpc and to time-since-coalescence. Taken together, our results show that the stellar kinematic data has little to offer in compliment to imaging for merger remnant identification in a heterogeneous galaxy population. △ Less

Submitted 10 January, 2022; originally announced January 2022.

Comments: MNRAS accepted; 20 pages, 12 figures and 2 tables; data available upon request

Do AGN triggering mechanisms vary with radio power? II. The importance of mergers as a function of radio power and optical luminosity

Authors: J. C. S. Pierce, C. N. Tadhunter, Y. Gordon, C. Ramos Almeida, S. L. Ellison, C. O'Dea, L. Grimmett, L. Makrygianni, P. S. Bessiere, P. Doña Girón

Abstract: Investigation of the triggering mechanisms of radio AGN is important for improving our general understanding of galaxy evolution. In the first paper in this series, detailed morphological analysis of high-excitation radio galaxies (HERGs) with intermediate radio powers suggested that the importance of triggering via galaxy mergers and interactions increases strongly with AGN radio power and weakly… ▽ More Investigation of the triggering mechanisms of radio AGN is important for improving our general understanding of galaxy evolution. In the first paper in this series, detailed morphological analysis of high-excitation radio galaxies (HERGs) with intermediate radio powers suggested that the importance of triggering via galaxy mergers and interactions increases strongly with AGN radio power and weakly with optical emission-line luminosity. Here, we use an online classification interface to expand our morphological analysis to a much larger sample of 155 active galaxies (3CR radio galaxies, radio-intermediate HERGs and Type 2 quasars) that covers a broad range in both 1.4 GHz radio power and [OIII]$λ$5007 emission-line luminosity. All active galaxy samples are found to exhibit excesses in their rates of morphological disturbance relative to 378 stellar-mass- and redshift-matched non-active control galaxies classified randomly and blindly alongside them. These excesses are highest for the 3CR HERGs (4.7$σ$) and Type 2 quasar hosts (3.7$σ$), supporting the idea that galaxy mergers provide the dominant triggering mechanism for these subgroups. When the full active galaxy sample is considered, there is clear evidence to suggest that the enhancement in the rate of disturbance relative to the controls increases strongly with [OIII]$λ$5007 emission-line luminosity but not with 1.4 GHz radio power. Evidence that the dominant AGN host types change from early-type galaxies at high radio powers to late-type galaxies at low radio powers is also found, suggesting that triggering by secular, disk-based processes holds more importance for lower-power radio AGN. △ Less

Submitted 4 November, 2021; originally announced November 2021.

Comments: Accepted for publication in MNRAS. 21 pages, 11 figures. Supplementary information will be available online with the published journal article (also available from the corresponding author by request)

VERTICO: The Virgo Environment Traced In CO Survey

Authors: Toby Brown, Christine D. Wilson, Nikki Zabel, Timothy A. Davis, Alessandro Boselli, Aeree Chung, Sara L. Ellison, Claudia D. P. Lagos, Adam R. H. Stevens, Luca Cortese, Yannick M. Bahé, Dhruv Bisaria, Alberto D. Bolatto, Claire R. Cashmore, Barbara Catinella, Ryan Chown, Benedikt Diemer, Pascal J. Elahi, Maan H. Hani, María J. Jiménez-Donaire, Bumhyun Lee, Katya Leidig, Angus Mok, Karen Pardos Olsen, Laura C. Parker , et al. (11 additional authors not shown)

Abstract: We present the Virgo Environment Traced in CO (VERTICO) survey, a new effort to map $^{12}$CO($2-1$), $^{13}$CO($2-1$), and C$^{18}$O($2-1$) in 51 Virgo Cluster galaxies with the Atacama Compact Array, part of the Atacama Large Millimeter/submillimeter Array (ALMA). The primary motivation of VERTICO is to understand the physical mechanisms that perturb molecular gas disks, and therefore star forma… ▽ More We present the Virgo Environment Traced in CO (VERTICO) survey, a new effort to map $^{12}$CO($2-1$), $^{13}$CO($2-1$), and C$^{18}$O($2-1$) in 51 Virgo Cluster galaxies with the Atacama Compact Array, part of the Atacama Large Millimeter/submillimeter Array (ALMA). The primary motivation of VERTICO is to understand the physical mechanisms that perturb molecular gas disks, and therefore star formation and galaxy evolution, in dense environments. This first paper contains an overview of VERTICO's design and sample selection, $^{12}$CO($2-1$) observations, and data reduction procedures. We characterize global $^{12}$CO($2-1$) fluxes and molecular gas masses for the 49 detected VERTICO galaxies, provide upper limits for the two non-detections, and produce resolved $^{12}$CO($2-1$) data products (median resolution $= 8^{\prime\prime} \approx 640~{\rm pc}$). Azimuthally averaged $^{12}$CO($2-1$) radial intensity profiles are presented along with derived molecular gas radii. We demonstrate the scientific power of VERTICO by comparing the molecular gas size--mass scaling relation for our galaxies with a control sample of field galaxies, highlighting the strong effect that radius definition has on this correlation. We discuss the drivers of the form and scatter in the size--mass relation and highlight areas for future work. VERTICO is an ideal resource for studying the fate of molecular gas in cluster galaxies and the physics of environment-driven processes that perturb the star formation cycle. Upon public release, the survey will provide a homogeneous legacy dataset for studying galaxy evolution in our closest cluster. △ Less

Submitted 1 November, 2021; originally announced November 2021.

Comments: 68 pages, 13 Figures, 2 Figure Sets, Accepted for publication in ApJS, Online FITS versions of Tables 1, 2, and 3 are available with the journal publication

HI absorption at z~0.7 against the lobe of the powerful radio galaxy PKS 0409-75

Authors: Elizabeth K. Mahony, James R. Allison, Elaine M. Sadler, Sara L. Ellison, Sui Ann Mao, Raffaella Morganti, Vanessa A. Moss, Amit Seta, Clive N. Tadhunter, Simon Weng, Matthew T. Whiting, Hyein Yoon, Martin Bell, John D. Bunton, Lisa Harvey-Smith, Amy Kimball, Bärbel S. Koribalski, Max A. Voronkov

Abstract: We present results from a search for the HI 21-cm line in absorption towards 16 bright radio sources with the 6-antenna commissioning array of the Australian Square Kilometre Array Pathfinder (ASKAP). Our targets were selected from the 2-Jy sample, a flux-limited survey of the southern radio sky with extensive multi-wavelength follow-up. Two sources were detected in HI absorption including a new d… ▽ More We present results from a search for the HI 21-cm line in absorption towards 16 bright radio sources with the 6-antenna commissioning array of the Australian Square Kilometre Array Pathfinder (ASKAP). Our targets were selected from the 2-Jy sample, a flux-limited survey of the southern radio sky with extensive multi-wavelength follow-up. Two sources were detected in HI absorption including a new detection towards the bright FRII radio galaxy PKS 0409-75 at a redshift of z=0.674. The HI absorption line is blueshifted by ~3300 km/s compared to the optical redshift of the host galaxy of PKS 0409-75 at z=0.693. Deep optical imaging and spectroscopic follow-up with the GMOS instrument on the Gemini-South telescope reveal that the HI absorption is associated with a galaxy in front of the southern radio lobe with a stellar mass of $3.2 - 6.8 \times 10^{11}M_\odot$, a star-formation rate of $\sim 1.24 M_\odot$ yr$^{-1}$ and an estimated HI column density of $2.16\times10^{21}$ cm$^{-2}$, assuming a spin temperature of $T_{\rm spin}=500$ K and source covering factor of $C_{\rm f}=0.3$. Using polarisation measurements of PKS 0409-75 from the literature we estimate the magnetic field of the absorbing galaxy to be ~14.5$μ$G, consistent with field strengths observed in nearby spiral galaxies, but larger than expected for an elliptical galaxy. Results from this pilot study can inform future surveys as new wide-field telescopes allow us to search for 21-cm HI absorption towards all bright radio sources as opposed to smaller targeted samples. △ Less

Submitted 19 October, 2021; originally announced October 2021.

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

Gas flows in galaxy mergers: supersonic turbulence in bridges, accretion from the circumgalactic medium, and metallicity dilution

Authors: Martin Sparre, Joseph Whittingham, Mitali Damle, Maan H. Hani, Philipp Richter, Sara L. Ellison, Christoph Pfrommer, Mark Vogelsberger

Abstract: In major galaxy mergers, the orbits of stars are violently perturbed, and gas is torqued to the centre, diluting the gas metallicity and igniting a starburst. In this paper, we study the gas dynamics in and around merging galaxies using a series of cosmological magneto-hydrodynamical (MHD) zoom-in simulations. We find that the gas bridge connecting the merging galaxies pre-coalescence is dominated… ▽ More In major galaxy mergers, the orbits of stars are violently perturbed, and gas is torqued to the centre, diluting the gas metallicity and igniting a starburst. In this paper, we study the gas dynamics in and around merging galaxies using a series of cosmological magneto-hydrodynamical (MHD) zoom-in simulations. We find that the gas bridge connecting the merging galaxies pre-coalescence is dominated by turbulent pressure, with turbulent Mach numbers peaking at values of 1.6-3.3. This implies that bridges are dominated by supersonic turbulence, and are thus ideal candidates for studying the impact of extreme environments on star formation. We also find that gas accreted from the circumgalactic medium (CGM) during the merger significantly contributes (27-51 per cent) to the star formation rate (SFR) at the time of coalescence and drives the subsequent reignition of star formation in the merger remnant. Indeed, 19-53 per cent of the SFR at z=0 originates from gas belonging to the CGM prior the merger. Finally, we investigate the origin of the metallicity-diluted gas at the centre of merging galaxies. We show that this gas is rapidly accreted onto the galactic centre with a time-scale much shorter than that of normal star-forming galaxies. This explains why coalescing galaxies are not well-captured by the fundamental metallicity relation. △ Less

Submitted 7 October, 2021; originally announced October 2021.

Comments: 15 pages, 13 figures, submitted to MNRAS