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

The effect of galaxy interactions on star formation rates in the COLIBRE simulations

Authors: Quinten W. E. van Zegveld, Evgenii Chaikin, Joop Schaye, David R. Patton, Sara L. Ellison, Alejandro Benítez-Llambay, Filip Huško, Robert J. McGibbon, Sylvia Ploeckinger, Alexander J. Richings, Matthieu Schaller

Abstract: Observations and theory indicate that galaxy interactions enhance star formation rates (SFRs). However, the degree of enhancement and its dependence on the properties of the interacting galaxies vary across different studies. In this work, we use the COLIBRE simulations of galaxy formation to investigate the effect of interactions on the SFRs of star-forming galaxies at redshift $z\approx0$. The C… ▽ More Observations and theory indicate that galaxy interactions enhance star formation rates (SFRs). However, the degree of enhancement and its dependence on the properties of the interacting galaxies vary across different studies. In this work, we use the COLIBRE simulations of galaxy formation to investigate the effect of interactions on the SFRs of star-forming galaxies at redshift $z\approx0$. The COLIBRE simulations capture the multiphase nature of the interstellar medium and have volumes up to $200^3$ and $400^3$ cMpc$^3$ at m6 (gas and dark-matter particle mass $\sim10^6~\mathrm{M_\odot}$) and m7 ($\sim10^7~\mathrm{M_\odot}$) resolutions, respectively. After constructing samples of interacting galaxies (with mass ratios $>0.1$) and isolated controls, matched in stellar mass, large- and small-scale environment, and redshift, we show that the average specific SFR (sSFR) of interacting galaxies is enhanced by up to a factor of $\approx2$ for separations of $\approx10$ kpc. The enhancement decreases with pair separation but remains significant out to $\approx200$ kpc. The enhancement increases with increasing numerical resolution, is more pronounced in the central regions of galaxies, and decreases with increasing stellar mass at fixed separation. Mergers with higher mass ratios induce stronger sSFR enhancement. We compare our results with observational data from the SDSS, finding good agreement in the dependence of the mean sSFR enhancement on separation, but underpredicting its normalisation by a factor of $\approx2$. Finally, we show that the pre-merger sSFR enhancement of resolved interactions accounts for $\approx2$ per cent of the $z\approx0$ cosmic SFR density. △ Less

Submitted 12 August, 2026; originally announced August 2026.

Comments: 21 pages, 20 figures; submitted to MNRAS

MAUVE-MUSE: When Metallicity Follows or Fights Star Formation -- A Mass-Dependent Inversion in Virgo Galaxies

Authors: Rongjun Huang, Luca Cortese, Barbara Catinella, Luke J. M. Davies, Toby Brown, Andrei Ristea, Alessandro Boselli, Andrew J. Battisti, Vicente Villanueva, Kristine Spekkens, Sara L. Ellison, Daniel A. Dale, Sabine Thater, Amirnezam Amiri

Abstract: Although globally-integrated studies often find that, at fixed stellar mass, high star formation rate (SFR) galaxies are relatively metal-poor while lower-SFR systems are more metal-rich, the corresponding coupling between gas-phase metallicity ($Z_{\rm gas}$) and star formation on sub-galactic scales remains poorly constrained. In this study, we analyse 14 Virgo spirals from the MAUVE-MUSE survey… ▽ More Although globally-integrated studies often find that, at fixed stellar mass, high star formation rate (SFR) galaxies are relatively metal-poor while lower-SFR systems are more metal-rich, the corresponding coupling between gas-phase metallicity ($Z_{\rm gas}$) and star formation on sub-galactic scales remains poorly constrained. In this study, we analyse 14 Virgo spirals from the MAUVE-MUSE survey to revisit the resolved mass-metallicity relation (rMZR) and its secondary dependence on SFR surface density ($Σ_\mathrm{SFR}$) at $\sim 100$\,pc scales. We construct co-spatial maps of stellar mass surface density ($Σ_*$), $Σ_\mathrm{SFR}$, and oxygen abundance. MAUVE-MUSE galaxies follow a standard rMZR, but when binned by $Σ_*$, we find a mass-dependent inversion in the $Z_\mathrm{gas}$-$Σ_\mathrm{SFR}$ relation using O3N2 calibrations: the standard anti-correlation is confined to low-$Σ_*$ bins, while high-$Σ_*$ regions show a positive correlation, inverting at $\log_{10}(Σ_*/M_\odot\,\mathrm{kpc}^{-2})\simeq 7.5$-8.0. Correlated and anti-correlated \ion{H}{ii} regions coexist within the same discs; the mass dependence emerges only when grouping spaxels by $Σ_*$. We develop a spatially resolved gas-regulator model showing this $Z_\mathrm{gas}$-$Σ_\mathrm{SFR}$ (anti-)correlation arises from competition between star-formation-driven and gas-supply-driven variability. This framework naturally extrapolates to the integrated scenario, providing a unified explanation for resolved and global relations. However, the presence and strength of the $Z_\mathrm{gas}$-$Σ_\mathrm{SFR}$ (anti-)correlation depend strongly on the metallicity indicator used, highlighting the challenge of disentangling physical secondary trends within metallicity scaling relations. △ Less

Submitted 29 May, 2026; originally announced May 2026.

Comments: 20 pages, 7 figues, accepted by MNRAS

A Comprehensive Analysis of WISE Mid-Infrared Colors for Obscured AGN Selection

Authors: Anika Goel, Samir Salim, Sara L. Ellison, Shobita Satyapal, Sheyda Salehirad, Robert W. Bickley, Christopher J. Agostino

Abstract: In this paper, we investigate the robustness of WISE mid-IR color selection (W1-W2) for identifying obscured (Type 2) active galactic nuclei (AGNs) at low redshift (z<0.3), using a sample of ~360,000 SDSS galaxies classified via emission lines into Seyfert 2 (Sy2), LINER, and star-forming (BPT-SF) galaxies. We find that the K-correction is essential to remove non-AGN contamination, and once applie… ▽ More In this paper, we investigate the robustness of WISE mid-IR color selection (W1-W2) for identifying obscured (Type 2) active galactic nuclei (AGNs) at low redshift (z<0.3), using a sample of ~360,000 SDSS galaxies classified via emission lines into Seyfert 2 (Sy2), LINER, and star-forming (BPT-SF) galaxies. We find that the K-correction is essential to remove non-AGN contamination, and once applied the simple W1-W2>0.5 selection emerges as optimal in terms of purity and completeness of AGN selection. However, we confirm that even this lenient cut selects only ~13% of Sy2 galaxies and that achieving W1-W2>0.5 requires AGN contributing >75% of the total infrared luminosity, which is uncommon. Although mid-IR-selected Sy2s tend to be luminous, the high [OIII] luminosity does not guarantee red W1-W2 (nor does any other tested global or NLR-scale parameter), suggesting the critical role of obscuration on smaller scales. <1% of BPT-SF systems (but making ~20% of all mid-IR selected galaxies) exhibit W1-W2>0.5 colors. Such colors cannot be reproduced by models of star-heated dust alone. Red BPT-SFs tend to have higher W4 luminosities than expected from SF, indicating true AGNs. Intriguingly, mid-IR AGNs in massive bulges ($M_{\mathrm{bulge}} \gtrsim 10^{10} M_{\odot}$) predominantly (84%) manifest themselves as BPT-AGNs, whereas those in low-mass bulges ($\lesssim 10^{10} M_{\odot}$) mostly (60%) manifest as BPT-SF. This BPT-AGN vs.\ BPT-SF dichotomy does not extend to total stellar mass. We conclude that although the mid-IR AGN selection is incomplete, its strength lies in identifying optically inconspicuous AGNs with low-mass bulges, regardless of the total mass. △ Less

Submitted 31 March, 2026; originally announced April 2026.

Comments: Submitted to ApJ. Comments welcome

EMBERS I: Low redshift post-starburst galaxies are frequently depleted in molecular gas relative to star forming progenitors

Authors: Ben F. Rasmussen, María Jesús Jiménez-Donaire, Sara L. Ellison, Vivienne Wild, Kate Rowlands, Qifeng Huang, Jing Wang, Dong Yang, Scott Wilkinson, Blake Ledger, Toby Brown, Ho-Hin Leung, Shoshannah Byrne-Mamahit

Abstract: The cold gas content of post-starburst galaxies (PSBs) provides important insight into the mechanisms that drive rapid quenching, but a multiphase assessment of both the atomic and molecular gas in PSBs does not yet exist. We introduce the Ensemble of Multiphase Baryons Evolving in Rapidly-quenching Systems, or EMBERS, a homogeneously selected, nearly mass- and redshift-complete survey of the glob… ▽ More The cold gas content of post-starburst galaxies (PSBs) provides important insight into the mechanisms that drive rapid quenching, but a multiphase assessment of both the atomic and molecular gas in PSBs does not yet exist. We introduce the Ensemble of Multiphase Baryons Evolving in Rapidly-quenching Systems, or EMBERS, a homogeneously selected, nearly mass- and redshift-complete survey of the global atomic (HI) and molecular gas (H2) in PSBs, observed with the Five Hundred-metre Aperture Spherical Telescope (FAST) and the Institut de radioastronomie millimetrique (IRAM) 30m telescope. We present new CO(1-0) observations for 52 PSBs with the IRAM 30m, which, combined with 9 archival observations, gives a total H2 sample of 61, of which 58/61 have ancillary HI measurements. We detect CO(1-0) in 34/61 galaxies, corresponding to molecular gas fractions (fH2 = MH2/M*) ranging from two to 250 per cent. By comparing with a stellar-mass matched star-forming (SF) control sample from xCOLD GASS, we find that PSBs on average are 0.3-0.6 dex depleted in H2. However, considering both HI and H2, individual PSBs host diverse gas reservoirs ranging from gas-rich in both phases, elevated in one phase, or gas-poor, the latter of which is common at lower stellar mass. The existence of gas-normal and gas-depleted PSBs in both phases suggests that some PSBs may rejuvenate their star formation, but the rapid shutdown of star formation in others is likely terminal. Despite this diversity, the majority of EMBERS PSBs are gas-poor compared to SF controls, with the typical PSB hosting gas reservoirs intermediate to those found in star-forming and quenched galaxies. △ Less

Submitted 27 February, 2026; originally announced March 2026.

Comments: 27 pages, 9 figures, Accepted to The Monthly Notices of the Royal Astronomical Society

Journal ref: Mon Not R Astron Soc (2026)

Merger Driven or Internal Evolution? A New Morphological Study of Tidal Disruption Event Host Galaxies

Authors: Janet N. Y. Chang, Connor Bottrell, Lixin Dai, Rudrani Kar Chowdhury, Meng Gu, Renbin Yan, Leonardo Ferreira, Sara L. Ellison, Scott Wilkinson, Thomas de Boer

Abstract: The host galaxies of tidal disruption events (TDEs) show enhanced central stellar concentration and are preferentially found in poststarburst and green valley populations. This connection has led to the proposal that TDE host galaxies likely have gone through recent mergers. We conduct a new morphological study of 14 TDE host galaxies, using the r-band images from the Sloan Digital Sky Survey (SDS… ▽ More The host galaxies of tidal disruption events (TDEs) show enhanced central stellar concentration and are preferentially found in poststarburst and green valley populations. This connection has led to the proposal that TDE host galaxies likely have gone through recent mergers. We conduct a new morphological study of 14 TDE host galaxies, using the r-band images from the Sloan Digital Sky Survey (SDSS), Dark Energy Camera Legacy Survey, and Ultraviolet Near-Infrared Optical Northern Survey, with the images from the latter two surveys having much higher depth and resolution than SDSS. We examine galaxy structures using conventional methods and also apply diagnostics of merger activity from a suite of machine learning models. Consistent with previous studies, our results show that TDE host galaxies are ~16% more centrally concentrated when compared to non-TDE-host controls. However, surprisingly, TDE hosts lack any indication of significant recent merger activity from both morphological analysis and the machine learning merger classifier. Instead, our results reveal that TDE host galaxies in the green valley are approximately 1.5-3 times more likely to have bar-like or ringlike structures compared to their controls. Based on these results, we propose that bar-driven secular evolution, instead of mergers, likely dominates the recent evolution of the TDE hosts found in the green valley, which can simultaneously explain their distinctive nuclear properties and enhanced TDE rates. △ Less

Submitted 8 June, 2026; v1 submitted 6 February, 2026; originally announced February 2026.

Comments: 15 pages, 10 figures

The ALMaQUEST Survey XVII: Unveiling Multiple Quenching Pathways in Green Valley Galaxies via Molecular Gas and Quenching Timescale Analyses

Authors: Lihwai Lin, Po-Feng Wu, Mallory D. Thorp, Asa F. L. Bluck, Hsi-An Pan, Sara L. Ellison, Kate Rowlands, Justin Atsushi Otter, Sebastián F. Sánchez

Abstract: Statistically, green valley (GV) galaxies exhibit lower molecular gas fractions ($f_{gas}$) and reduced star formation efficiency (SFE) compared to star-forming galaxies. However, it remains unclear whether quenching is primarily driven by one factor or results from a combination of mechanisms in individual GV galaxies. In this study, we address this question by examining the spatial distributions… ▽ More Statistically, green valley (GV) galaxies exhibit lower molecular gas fractions ($f_{gas}$) and reduced star formation efficiency (SFE) compared to star-forming galaxies. However, it remains unclear whether quenching is primarily driven by one factor or results from a combination of mechanisms in individual GV galaxies. In this study, we address this question by examining the spatial distributions of star formation and molecular gas in 28 GVs selected from the ALMaQUEST survey and additional literature samples. For each galaxy, we identify regions with suppressed specific star formation rate (sSFR) and measure $Δf_{gas}$ and $Δ$SFE-offsets from the resolved scaling relations of the star-forming main sequence galaxies. By comparing the fraction of regions with negative $Δf_{gas}$ and $Δ$SFE, we classify 35.7$\pm$13.2\% (57.1$\pm$17.9\%) of GV galaxies as $f_{gas}$-driven, 39.3$\pm$14.0\% (39.3$\pm$14.0\%) as SFE-driven, and 25.0$\pm$10.6\% (3.6$\pm$3.6\%) as mixed mode when adopting a fixed (variable) CO-to-$\rm H_{2}$ conversion factor ($α_{CO}$). These results indicate that GVs undergo quenching through multiple pathways. As sSFR decreases from the main sequence to the green valley, we observe a transition toward predominantly SFE-driven quenching, possibly linked to internal processes such as morphological quenching or AGN activity. We further estimate the quenching timescale ($τ_{decay}$), defined as the time from the peak SFR to 1/e (approximately 37\%) of its value, using integrated MaNGA spectra. SFE-driven quenching is typically associated with short $τ_{decay}$ , while $f_{gas}$-driven quenching shows a broader range. Overall, 75\% of GVs exhibit $τ_{decay}$ shorter than 1 Gyr, suggesting that quenching in most GVs proceeds rapidly, challenging purely slow-quenching scenarios like starvation. △ Less

Submitted 14 January, 2026; originally announced January 2026.

Comments: 20 pages, 9 figures + 1 figure set, accepted for publication in ApJ

Journal ref: ApJ, 999, 263 (2026)

Galaxy Mergers in UNIONS -- II: Predicting Timescales in the Post-Merger Regime

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

Abstract: Galaxy mergers are critical events that influence galaxy evolution by driving processes such as enhanced star formation, quenching, and active galactic nucleus (AGN) activity. However, constraining the timescales over which these processes occur in the post-merger phase has remained a significant challenge. This study extends the MUlti-Model Merger Identifier (\textsc{Mummi}) framework to predict… ▽ More Galaxy mergers are critical events that influence galaxy evolution by driving processes such as enhanced star formation, quenching, and active galactic nucleus (AGN) activity. However, constraining the timescales over which these processes occur in the post-merger phase has remained a significant challenge. This study extends the MUlti-Model Merger Identifier (\textsc{Mummi}) framework to predict post-merger timescales ($T_{PM}$) for galaxies, leveraging machine learning models trained on realism-enhanced mock observations derived from the IllustrisTNG simulations. By classifying post-merger galaxies into four temporal bins spanning 0 to 1.76 Gyr after coalescence, \textsc{Mummi} achieves time classification accuracies exceeding 70 per cent. We apply this framework to the Ultraviolet Near Infrared Optical Northern Survey (UNIONS), yielding a catalog of 8,716 post-merger galaxies with $T_{PM}$ predictions and stellar masses $\log(M_*/M_\odot) \geq 10$ at redshifts 0.03 < z < 0.3. These results provide a robust methodology to connect galaxy interaction timescales with physical processes, enabling detailed studies of galaxy evolution in the post-merger regime. △ Less

Submitted 12 January, 2026; originally announced January 2026.

Comments: Accepted for publication in MNRAS; 12 pages, 10 figures, data available at https://github.com/astroferreira/MUMMI_UNIONS

Galaxy evolution in the post-merger regime. IV -- The long-term effect of mergers on galactic stellar mass growth and distribution

Authors: Sara L. Ellison, Leonardo Ferreira

Abstract: Galaxy mergers are known to trigger bursts of central star formation, which should therefore lead to stellar mass growth in their inner regions. However, observational measurements of this `burst mass fraction' are scant. Here, we assemble a large sample of ~14,000 post-coalescence galaxies that have recently completed their merger-induced star formation, and compare various measurements of centra… ▽ More Galaxy mergers are known to trigger bursts of central star formation, which should therefore lead to stellar mass growth in their inner regions. However, observational measurements of this `burst mass fraction' are scant. Here, we assemble a large sample of ~14,000 post-coalescence galaxies that have recently completed their merger-induced star formation, and compare various measurements of central stellar mass with a matched control sample. Specifically, we quantify (at fixed redshift, star formation rate and total stellar mass) the stellar mass enhancement within a fixed angular aperture (Delta M_{star,fibre}) and in the galactic bulge (Delta M_{star,bulge}), finding burst mass fractions of 10 -- 20%. 61 galaxies in our sample are at z<0.05 and have integral field unit data from the Mapping Galaxies at Apache Point (MaNGA) survey, allowing further kpc-scale assessment of excess stellar mass and radial gradients. When assessed within apertures defined in units of kpc we again find a ~ 15 -- 20% excess of stellar mass in the central regions of the post-mergers compared with matched controls. However, within apertures defined in units of effective radius this stellar mass enhancement increases to 40%, suggesting that the relative structure/size of the galaxy is important for regulating the location of the merger induced star formation. Moreover, we find that these stellar mass enhancements are spatially extended, out to ~7 kpc or around 1 R/R_e, although the small sample size of the MaNGA overlap limits our radial sampling. Our work represents the first direct measurement of merger-induced stellar mass growth that is independent of stellar population modelling, or fitting light profiles, demonstrating significant and extended mass build-up in late stage post-mergers. △ Less

Submitted 23 April, 2026; v1 submitted 26 November, 2025; originally announced November 2025.

Comments: Published the Open Journal of Astrophysics

SDSS-ALMA Legacy Value Archival Gas Exploration (SALVAGE) -- I: global star formation is governed by central (not global) molecular gas

Authors: Scott Wilkinson, Toby Brown, Chiara Circosta, Sara L. Ellison, Blake Ledger, Samuel D. Fielder

Abstract: Star-forming galaxies form tight relations between their stellar mass, star-formation rate, and molecular gas reservoir on global and resolved scales. On the path to quiescence, the exchange between gas and stars must inevitably be broken. Understanding the mechanisms governing star formation and quenching therefore requires observations of both the stellar and molecular gas components. To this en… ▽ More Star-forming galaxies form tight relations between their stellar mass, star-formation rate, and molecular gas reservoir on global and resolved scales. On the path to quiescence, the exchange between gas and stars must inevitably be broken. Understanding the mechanisms governing star formation and quenching therefore requires observations of both the stellar and molecular gas components. To this end, we have assembled a sample of 277 galaxies ($0.02 \lesssim z \lesssim 0.25$) with semi-resolved optical and millimetre $^{12}$CO(1-0) data, wherein the properties of the inner $\thicksim$2 kpc can be distinguished from the outer regions. This effort was made possible by the Sloan Digital Sky Survey (SDSS) catalogues and the maturing archive of the Atacama Large (sub-)Millimetre Array (ALMA). We call this dataset the SDSS-ALMA Legacy Value Archival Gas Exploration (SALVAGE). In this work, we leverage SALVAGE to provide a semi-resolved perspective on global scaling relations and why some galaxies deviate from them. In agreement with previous work, we find that the offset of a galaxy from the global star-forming main sequence (SFMS) is driven by its inner star formation rate. With the relative inner and outer distributions of molecular gas fraction and star formation efficiency, we investigate whether the central star formation driving global changes is due to fuel availability or efficiency. We find that the position of a galaxy within the SFMS is largely due to the inner star-formation efficiency, while departure from the SFMS is driven by availability of central gas. The central few kpc are thus the most consequential region for galaxy evolution at low redshift. △ Less

Submitted 24 November, 2025; v1 submitted 10 November, 2025; originally announced November 2025.

Comments: 21 pages, 11 figures. Accepted by MNRAS. Data available at https://www.canfar.net/storage/vault/list/AstroDataCitationDOI/CISTI.CANFAR/25.0077/data

Interacting galaxies in the IllustrisTNG simulations - IX: Mini mergers trigger AGN in cosmological simulations

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

Abstract: Galaxy mergers are transformative events that can cause gaseous inflows capable of triggering active galactic nuclei (AGN). Previous studies of AGN in simulations have mainly focused on major interactions (i.e. between approximately equal mass galaxies), which produce the strongest inflows and, therefore, would be the most likely to trigger AGN activity. However, minor interactions are far more co… ▽ More Galaxy mergers are transformative events that can cause gaseous inflows capable of triggering active galactic nuclei (AGN). Previous studies of AGN in simulations have mainly focused on major interactions (i.e. between approximately equal mass galaxies), which produce the strongest inflows and, therefore, would be the most likely to trigger AGN activity. However, minor interactions are far more common and may still enhance accretion onto supermassive black holes. We present an analysis of post-merger galaxies from the IllustrisTNG simulation with stellar mass ratios of $μ>$1:100. We select post-mergers from the TNG50-1 simulation, from redshifts $0\leq z< 1$, with stellar masses greater than $10^{10}M_{\odot}$. We find an AGN excess in post-mergers with a stellar mass ratio as low as 1:40. The AGN excess is mass ratio and luminosity dependent, with 1.2-1.6 times more AGN found in post-mergers of 1:40$\leq μ< $1:4 than in matched non-merger controls, and as many as 6 times more AGN found in major $μ\geq$1:4 post-mergers. The AGN excess is long lived, between 500 Myr to 1 Gyr post-coalescence, across all of the mass ratio regimes. We demonstrate that the most luminous AGN in the simulation overwhelmingly occur in either post-mergers or pairs (with $μ\geq $1:40). Finally, we demonstrate that mini mergers are likely to be overlooked in observational studies due to the weakness of features usually associated with recent merger activity, such as tidal streams and shells, making it challenging to completely account for merger-induced AGN activity even in deep galaxy surveys. △ Less

Submitted 14 October, 2025; originally announced October 2025.

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

Decoding the Single-peaked HI Spectra of Low Redshift Post-starburst Galaxies

Authors: Qifeng Huang, Jing Wang, Sara L. Ellison, Zezhong Liang, Xuchen Lin, Dong Yang

Abstract: Recent observations with the Five-hundred-meter Aperture Spherical Telescope (FAST) have revealed abundant reservoirs of neutral hydrogen (HI) in low redshift post-starburst galaxies (PSBs), raising the question of why star formation ceases rapidly in these systems. In this study, we present a detailed analysis of the shape of the integrated HI spectra of 67 PSBs. We find that PSBs exhibit signifi… ▽ More Recent observations with the Five-hundred-meter Aperture Spherical Telescope (FAST) have revealed abundant reservoirs of neutral hydrogen (HI) in low redshift post-starburst galaxies (PSBs), raising the question of why star formation ceases rapidly in these systems. In this study, we present a detailed analysis of the shape of the integrated HI spectra of 67 PSBs. We find that PSBs exhibit significantly higher HI spectral concentration values ($K$) compared to a matched sample from xGASS, and are more comparable to those of starburst galaxies. By extending our analysis to spatially resolved HI data from THINGS and ATLAS$^{\rm 3D}$, we show that both centrally concentrated HI distributions and dynamically unsettled HI can effectively increase $K$, while non-axisymmetric structures only contribute to the scatter of the $K$ distribution. Distinguishing between central concentration and dynamically unsettled gas as the origin of high $K$ can be achieved by measuring the spectral asymmetry ($A_{\rm F}$), making the $K$-$A_{\rm F}$ plane a powerful diagnostic tool for identifying galaxies with unsettled HI using integrated spectra alone. Based on their location in the $K$-$A_{\rm F}$ plane, we find that most PSBs are not dominated by unsettled HI, but rather exhibit elevated central gas concentration. Both modes of gas redistribution in PSBs may eventually contribute to their quenching. △ Less

Submitted 26 August, 2025; originally announced August 2025.

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

Unmixed metals: Variations in the enrichment of z~4 sub-damped Lyman alpha systems

Authors: Trystyn A. M. Berg, Louise A. Welsh, Ryan J. Cooke, Lise Christensen, Valentina D'Odorico, Sara L. Ellison, Sebastian Lopez

Abstract: The chemical abundance patterns of near-pristine objects provide important constraints on the properties of the first generations of stars in the Universe. We present the chemical abundances of five very metal-poor ([M/H]<-2.5) sub damped Lyman alpha systems (subDLAs) covering the redshift range $3.6<z<4.3$, identified with the XQ-100 survey. We find that the subDLAs in our sample show consistent… ▽ More The chemical abundance patterns of near-pristine objects provide important constraints on the properties of the first generations of stars in the Universe. We present the chemical abundances of five very metal-poor ([M/H]<-2.5) sub damped Lyman alpha systems (subDLAs) covering the redshift range $3.6<z<4.3$, identified with the XQ-100 survey. We find that the subDLAs in our sample show consistent chemical abundance patterns (in particular [C/O], [Al/O], and [Fe/O]) with those of very metal-poor DLAs. Based on Voigt profile fitting, the chemical abundance ratios [C/O], [Al/O], and [Si/O] of individual velocity components in at least three of the subDLAs shows some intrinsic scatter. In order to verify these chemical inhomogeneities in absorption components, we present a novel method for computing ionization corrections (ICs) on a component-by-component basis and show that ICs alone cannot explain the variations in [C/O], [Al/O], and [Si/O] between components of the same absorber at ~2 sigma significance. Comparing the observed abundance ratios to the simulated core-collapse supernovae yields of early stellar populations, we find that all individual components of the subDLAs appears to be enriched by progenitor masses of <30 M_sol. The observed inhomogeneities between components can be reproduced by differences in the progenitor mass or supernova explosion energy. As such, the observed chemical inhomogeneities between components can be explained by poorly mixed gas from different nucleosynthetic events. △ Less

Submitted 15 July, 2025; originally announced July 2025.

Comments: 21 pages, 10 figures, accepted for publication in MNRAS

UNIONS: The Ultraviolet Near-Infrared Optical Northern Survey

Authors: Stephen Gwyn, Alan W. McConnachie, Jean-Charles Cuillandre, Ken C. Chambers, Eugene A. Magnier, Michael J. Hudson, Masamune Oguri, Hisanori Furusawa, Hendrik Hildebrandt, Raymond Carlberg, Sara L. Ellison, Junko Furusawa, Raphaël Gavazzi, Rodrigo Ibata, Yannick Mellier, Ken Osato, H. Aussel, Lucie Baumont, Manuel Bayer, Olivier Boulade, Patrick Côté, David Chemaly, Cail Daley, Pierre-Alain Duc, A. Ellien , et al. (64 additional authors not shown)

Abstract: The Ultraviolet Near-Infrared Optical Northern Survey (UNIONS) is a "collaboration of collaborations" that is using the Canada-France-Hawai'i Telescope, the Pan-STARRS telescopes, and the Subaru Observatory to obtain $ugriz$ images of a core survey region of 6250 deg$^2$ of the northern sky. The $10σ$ point source depth of the data, as measured within a 2-arcsecond diameter aperture, are… ▽ More The Ultraviolet Near-Infrared Optical Northern Survey (UNIONS) is a "collaboration of collaborations" that is using the Canada-France-Hawai'i Telescope, the Pan-STARRS telescopes, and the Subaru Observatory to obtain $ugriz$ images of a core survey region of 6250 deg$^2$ of the northern sky. The $10σ$ point source depth of the data, as measured within a 2-arcsecond diameter aperture, are $[u,g,r,i,z] = [23.7, 24.5, 24.2, 23.8, 23.3]$\ in AB magnitudes. UNIONS is addressing some of the most fundamental questions in astronomy, including the properties of dark matter, the growth of structure in the Universe from the very smallest galaxies to large-scale structure, and the assembly of the Milky Way. It is set to become the major ground-based legacy survey for the northern hemisphere for the next decade and provides an essential northern complement to the static-sky science of the Vera C. Rubin Observatory's Legacy Survey of Space and Time. UNIONS supports the core science mission of the {\it Euclid} space mission by providing the data necessary in the northern hemisphere for the calibration of the wavelength dependence of the {\it Euclid} point-spread function and derivation of photometric redshifts in the North Galactic Cap. This region contains the highest quality sky for {\it Euclid}, with low backgrounds from the zodiacal light, stellar density, extinction, and emission from Galactic cirrus. Here, we describe the UNIONS survey components, science goals, data products, and the current status of the overall program. △ Less

Submitted 17 March, 2025; originally announced March 2025.

Comments: submitted to AJ

Metallicity Gradients in Modern Cosmological Simulations I: Tension Between Smooth Stellar Feedback Models and Observations

Authors: Alex M. Garcia, Paul Torrey, Aniket Bhagwat, Ruby J. Wright, Qian-hui Chen, Kathryn Grasha, Sophia Ridolfo, Z. S. Hemler, Arnab Sarkar, Priyanka Chakraborty, Erica J. Nelson, Ryan L. Sanders, Tiago Costa, Mark Vogelsberger, Lisa J. Kewley, Sara L. Ellison, Lars Hernquist

Abstract: The metallicity of galaxies, and its variation with galactocentric radius, provides key insights into the formation histories of galaxies and the physical processes driving their evolution. In this work, we analyze the radial metallicity gradients of star forming galaxies in the EAGLE, Illustris, IllustrisTNG, and SIMBA cosmological simulations across a broad mass (… ▽ More The metallicity of galaxies, and its variation with galactocentric radius, provides key insights into the formation histories of galaxies and the physical processes driving their evolution. In this work, we analyze the radial metallicity gradients of star forming galaxies in the EAGLE, Illustris, IllustrisTNG, and SIMBA cosmological simulations across a broad mass ($10^{8.0}M_\odot\leq M_\star \lesssim10^{12.0}M_\odot$) and redshift ($0\leq z\leq8$) range. We find that all simulations predict strong negative (i.e., radially decreasing) metallicity gradients at early cosmic times, likely due to their similar treatments of relatively smooth stellar feedback not providing sufficient mixing to quickly flatten gradients. The strongest redshift evolution occurs in galaxies with stellar masses of $10^{10.0}-10^{11.0}M_\odot$, while galaxies with stellar masses $< 10^{10}M_\odot$ and $>10^{11}M_\odot$ exhibit weaker redshift evolution. Our results of negative gradients at high-redshift contrast with the many positive and flat gradients in the $1<z<4$ observational literature. At $z>6$, the negative gradients observed with JWST and ALMA are flatter than those in simulations, albeit with closer agreement than at lower redshift. Overall, we suggest that these smooth stellar feedback galaxy simulations may not sufficiently mix their metal content radially, and that either stronger stellar feedback or additional subgrid turbulent metal diffusion models may be required to better reproduce observed metallicity gradients. △ Less

Submitted 29 July, 2025; v1 submitted 5 March, 2025; originally announced March 2025.

Comments: 16 pages, 4 figures, + appendices. Accepted to ApJ

Low redshift post-starburst galaxies host abundant HI reservoirs

Authors: Sara L. Ellison, Qifeng Huang, Dong Yang, Jing Wang, Vivienne Wild, Ben F. Rasmussen, Maria-Jesus Jimenez-Donaire, Kate Rowlands, Scott Wilkinson, Toby Brown, Ho-Hin Leung

Abstract: Studying the gas content of post-starburst (PSB) galaxies can provide valuable clues regarding the process of fast quenching. Although previous works have studied the molecular gas content of PSBs, only a handful of HI measurements exist. Here, we present new Five hundred metre Aperture Spherical Telescope (FAST) 21cm observations of 44 PSBs, leading to 43 detections or sensitive upper limits of H… ▽ More Studying the gas content of post-starburst (PSB) galaxies can provide valuable clues regarding the process of fast quenching. Although previous works have studied the molecular gas content of PSBs, only a handful of HI measurements exist. Here, we present new Five hundred metre Aperture Spherical Telescope (FAST) 21cm observations of 44 PSBs, leading to 43 detections or sensitive upper limits of HI, which we combine with 25 archival measurements, for a total sample of 68 PSB MHI measurements. HI is detected in 57/68 galaxies, with HI masses ranging from MHI ~10^8.5 up to 10^10 Msun and gas fractions (fHI = MHI/M*) from a few percent up to almost 30 percent. Post-starbursts therefore retain ample atomic gas reservoirs, despite no longer forming stars. By comparing with a stellar mass-matched sample of star-forming galaxies in xGASS, we find that PSBs have, on average, gas fractions lower by ~0.2-0.4 dex, consistent with a mild reduction compared with their progenitor population. However, PSBs show a diversity of HI properties; about half have HI gas masses within the expected scatter of the star-forming population with the remaining 50 per cent up to a factor of 10 more gas-poor. Compared with galaxies in the green valley, about two thirds of PSBs have gas fractions within the expected range, with the remaining third up to a factor of 10 more gas-rich. Our results demonstrate that quenching in PSBs is not the result of wholesale removal of the atomic gas reservoir and that the population has atomic gas fractions that span the range from star-forming to green valley galaxies. We find no correlation between HI gas mass and time since burst; even galaxies a Gyr past their burst can remain HI-normal. The significant gas reservoirs remaining in many PSBs leaves open the possibility for future rekindling of star formation. △ Less

Submitted 6 July, 2025; v1 submitted 4 March, 2025; originally announced March 2025.

Comments: Published in the Open Journal of Astrophysics

Unions with UNIONS: Using galaxy-galaxy lensing to probe galaxy mergers

Authors: Isaac Cheng, Jack Elvin-Poole, Michael J. Hudson, Ruxin Barré, Sara L. Ellison, Robert W. Bickley, Thomas J. L. de Boer, Sébastien Fabbro, Leonardo Ferreira, Sacha Guerrini, Hendrik Hildebrandt, Martin Kilbinger, Alan W. McConnachie, Ludovic van Waerbeke, Anna Wittje

Abstract: We use galaxy-galaxy lensing to investigate how the dark matter (DM) haloes and stellar content of galaxies with $0.012 \leq z \leq 0.32$ and $10 \leq \log_{10}(M_\star/\mathrm{M}_\odot) \leq 12$ change as a result of the merger process. To this end, we construct two samples of galaxies obtained from the Ultraviolet Near Infrared Optical Northern Survey (UNIONS), comprising 1 623 post-mergers and… ▽ More We use galaxy-galaxy lensing to investigate how the dark matter (DM) haloes and stellar content of galaxies with $0.012 \leq z \leq 0.32$ and $10 \leq \log_{10}(M_\star/\mathrm{M}_\odot) \leq 12$ change as a result of the merger process. To this end, we construct two samples of galaxies obtained from the Ultraviolet Near Infrared Optical Northern Survey (UNIONS), comprising 1 623 post-mergers and $\sim$30 000 non-merging controls, that live in low-density environments to use as our lenses. These samples are weighted to share the same distributions of stellar mass, redshift, and geometric mean distance to a galaxy's three nearest neighbours to ensure differences in the lensing signal are due to the merger process itself. We do not detect a statistically significant difference in the excess surface density profile of post-mergers and non-merging controls with current data. Fitting haloes composed of a point-like stellar mass component and an extended DM structure described by a Navarro-Frenk-White profile to the lensing measurements yields, for both samples, halo masses of $M_\text{halo} \sim 4\times10^{12}\,\mathrm{M}_\odot$ and a moderately negative correlation between $M_\text{halo}$ and concentration $c$. This allows us to rule out, at the 95% confidence level, merger-induced starbursts in which more than 60% of the stellar mass is formed in the burst. The application of our methods to upcoming surveys that are able to provide samples $\sim$10$\times$ larger than our current catalogue are expected to detect the weak-lensing signatures of mergers and further constrain their properties. △ Less

Submitted 3 September, 2025; v1 submitted 1 February, 2025; originally announced February 2025.

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

Star Formation Rates, Metallicities, and Stellar Masses on kpc-scales in TNG50

Authors: Jia Qi, Alex M. Garcia, Davis Robinson, Paul Torrey, Jorge Moreno, Kara N. Green, Aaron S. Evans, Z. S. Hemler, Lars Hernquist, Sara L. Ellison

Abstract: Integral field units (IFU) have extended our knowledge of galactic properties to kpc (or, sometimes, even smaller) patches of galaxies. These scales are where the physics driving galaxy evolution (feedback, chemical enrichment, etc.) take place. Quantifying the spatially-resolved properties of galaxies, both observationally and theoretically, is therefore critical to our understanding of galaxy ev… ▽ More Integral field units (IFU) have extended our knowledge of galactic properties to kpc (or, sometimes, even smaller) patches of galaxies. These scales are where the physics driving galaxy evolution (feedback, chemical enrichment, etc.) take place. Quantifying the spatially-resolved properties of galaxies, both observationally and theoretically, is therefore critical to our understanding of galaxy evolution. To this end, we investigate spatially-resolved scaling relations within galaxies of $M_\star>10^{9.0}$ at $z=0$ in IllustrisTNG. We examine both the resolved star-forming main sequence (rSFMS) and the resolved mass-metallicity relation (rMZR) using $1~{\rm kpc}\times1~{\rm kpc}$ maps. We find that the rSFMS in IllustrisTNG is well-described by a power-law, but is significantly shallower than the observed rSFMS. However, the disagreement between the rSFMS of IllustrisTNG and observations is likely driven by an overestimation of AGN feedback in IllustrisTNG for the higher mass hosts. Conversely, the rMZR for IllustrisTNG has very good agreement with observations. Furthermore, we argue that the rSFMS is an indirect result of the Schmidt-Kennicutt (SK) law and local gas relation, which are both independent of host galaxy properties. Finally, we expand upon a localized leaky-box model to study the evolution of idealized spaxels and find that it provides a good description of these resolved relations. The degree of agreement, however, between idealized spaxels and simulated spaxels depends on the `net' outflow rate for the spaxel, and the IllustrisTNG scaling relations indicate a preference for a low net outflow rate. △ Less

Submitted 5 September, 2025; v1 submitted 30 January, 2025; originally announced January 2025.

Comments: Submitted to ApJ; Comments Welcome

The CIViL* Survey: The Discovery of a C IV Dichotomy in the CGM of L* Galaxies

Authors: Samantha L. Garza, Jessica K. Werk, Trystyn A. M. Berg, Yakov Faerman, Benjamin D. Oppenheimer, Rongmon Bordoloi, Sara L. Ellison

Abstract: This paper investigates C IV absorption in the circumgalactic medium (CGM) of L* galaxies and its relationship with galaxy star formation rates. We present new observations from the C IV in L* survey (CIViL*; PID$\#$17076) using the Hubble Space Telescope/Cosmic Origins Spectrograph. By combining these measurements with archival C IV data (46 observations total), we estimate detection fractions fo… ▽ More This paper investigates C IV absorption in the circumgalactic medium (CGM) of L* galaxies and its relationship with galaxy star formation rates. We present new observations from the C IV in L* survey (CIViL*; PID$\#$17076) using the Hubble Space Telescope/Cosmic Origins Spectrograph. By combining these measurements with archival C IV data (46 observations total), we estimate detection fractions for star-forming (sSFR $>$ 10$^{-11}$ yr$^{-1}$) and passive galaxies (sSFR $\leq$ 10$^{-11}$ yr$^{-1}$) to be 72$_{-18}^{+14}$\% [21/29] and 23$_{-15}^{+27}$\% [3/13], respectively. This indicates a significant dichotomy in C IV presence between L* star-forming and passive galaxies, with over 99% confidence. This finding aligns with Tumlinson et al. (2011), which noted a similar dichotomy in O VI absorption. Our results imply a substantial carbon reservoir in the CGM of L* galaxies, suggesting a minimum carbon mass of $\gtrsim$ 3.03 $\times$ 10$^{6}$ M$_{\odot}$ out to 120 kpc. Together, these findings highlight a strong connection between star formation in galaxies and the state of their CGM, providing insight into the mechanisms governing galaxy evolution. △ Less

Submitted 16 December, 2024; originally announced December 2024.

Comments: 14 pages, 4 figures, ApJ Letters, accepted

Galaxy evolution in the post-merger regime. III -- The triggering of active galactic nuclei peaks immediately after coalescence

Authors: Sara L. Ellison, Leonardo Ferreira, Robert Bickley, Tess Grindlay, Samir Salim, Shoshannah Byrne-Mamahit, Shobita Satyapal, David R. Patton, Jillian M. Scudder

Abstract: Galaxy mergers have been shown to trigger AGN in the nearby universe, but the timescale over which this process happens remains unconstrained. The Multi-Model Merger Identifier (MUMMI) machine vision pipeline has been demonstrated to provide reliable predictions of time post-merger (T_PM) for galaxies selected from the Ultraviolet Near Infrared and Optical Northern Survey (UNIONS) up to T_PM=1.76… ▽ More Galaxy mergers have been shown to trigger AGN in the nearby universe, but the timescale over which this process happens remains unconstrained. The Multi-Model Merger Identifier (MUMMI) machine vision pipeline has been demonstrated to provide reliable predictions of time post-merger (T_PM) for galaxies selected from the Ultraviolet Near Infrared and Optical Northern Survey (UNIONS) up to T_PM=1.76 Gyr after coalescence. By combining the post-mergers identified in UNIONS with pre-coalescence galaxy pairs, we can study the triggering of AGN throughout the merger sequence. AGN are identified using a range of complementary metrics: mid-IR colours, narrow emission lines and broad emission lines, which can be combined to provide insight into the demographics of dust and luminosity of the AGN population. Our main results are: 1) Regardless of the metric used, we find that the peak AGN excess (compared with a matched control sample) occurs immediately after coalescence, at 0 < T_PM < 0.16 Gyr. 2) The excess of AGN is observed until long after coalescence; both the mid-IR selected AGN and broad line AGN are more common than in the control sample even in the longest time bin of our sample (0.96 < T_PM < 1.76 Gyr). 3) The AGN excess is larger for more luminous and bolometrically dominant AGN, and we find that AGN in post-mergers are generally more luminous than secularly triggered events. 4) A deficit of broad line AGN in the pre-merger phase, that evolves into an excess in post-mergers is consistent with evolution of the covering fraction of nuclear obscuring material. Before coalescence, tidally triggered inflows increase the covering fraction of nuclear dust; in the post-merger regime feedback from the AGN clears (at least some of) this material. 5) The statistical peak in the triggering of starbursts occurs contemporaneously with AGN, within 0.16 Gyr of coalescence. △ Less

Submitted 2 February, 2025; v1 submitted 3 December, 2024; originally announced December 2024.

Comments: Accepted for publication in the Open Journal of Astrophysics

The MAGPI Survey: radial trends in star formation across different cosmological simulations in comparison with observations at $z \sim$ 0.3

Authors: Marcie Mun, Emily Wisnioski, Katherine E. Harborne, Claudia D. P. Lagos, Lucas M. Valenzuela, Rhea-Silvia Remus, J. Trevor Mendel, Andrew J. Battisti, Sara L. Ellison, Caroline Foster, Matias Bravo, Sarah Brough, Scott M. Croom, Tianmu Gao, Kathryn Grasha, Anshu Gupta, Yifan Mai, Anilkumar Mailvaganam, Eric G. M. Muller, Gauri Sharma, Sarah M. Sweet, Edward N. Taylor, Tayyaba Zafar

Abstract: We investigate the internal and external mechanisms that regulate and quench star formation (SF) in galaxies at $z \sim 0.3$ using MAGPI observations and the EAGLE, Magneticum, and IllustrisTNG cosmological simulations. Using SimSpin to generate mock observations of simulated galaxies, we match detection/resolution limits in star formation rates and stellar mass, along with MAGPI observational det… ▽ More We investigate the internal and external mechanisms that regulate and quench star formation (SF) in galaxies at $z \sim 0.3$ using MAGPI observations and the EAGLE, Magneticum, and IllustrisTNG cosmological simulations. Using SimSpin to generate mock observations of simulated galaxies, we match detection/resolution limits in star formation rates and stellar mass, along with MAGPI observational details including the average point spread function and pixel scale. While we find a good agreement in the slope of the global star-forming main sequence (SFMS) between MAGPI observations and all three simulations, the slope of the resolved SFMS does not agree within 1 $-$ 2$σ$. Furthermore, in radial SF trends, good agreement between observations and simulations exists only for galaxies far below the SFMS, where we capture evidence for inside-out quenching. The simulations overall agree with each other between $\sim1.5-4 \ R_{\rm e}$ but show varying central suppression within $R \sim 1.5 \ R_{\rm e}$ for galaxies on and below the SFMS, attributable to different AGN feedback prescriptions. All three simulations show similar dependencies of SF radial trends with environment. Central galaxies are subject to both internal and external mechanisms, showing increased SF suppression in the centre with increasing halo mass, indicating AGN feedback. Satellite galaxies display increasing suppression in the outskirts as halo mass increases, indicative of environmental processes. These results demonstrate the power of spatially resolved studies of galaxies; while global properties align, radial profiles reveal discrepancies between observations and simulations and their underlying physics. △ Less

Submitted 26 November, 2024; originally announced November 2024.

Comments: 20 pages, 10 figures, submitted to MNRAS

Galaxy evolution in the post-merger regime. II -- Post-merger quenching peaks within 500 Myr of coalescence

Authors: Sara L. Ellison, Leonardo Ferreira, Vivienne Wild, Scott Wilkinson, Kate Rowlands, David R. Patton

Abstract: (Abridged) Mechanisms for quenching star formation in galaxies remain hotly debated, with galaxy mergers an oft-proposed pathway. In Ellison et al. (2022) we tested this scenario by quantifying the fraction of recently and rapidly quenched post-starbursts (PSBs) in a sample of post-merger galaxies identified in the Ultraviolet Near Infrared Optical Northern Survey (UNIONS). With our recent develop… ▽ More (Abridged) Mechanisms for quenching star formation in galaxies remain hotly debated, with galaxy mergers an oft-proposed pathway. In Ellison et al. (2022) we tested this scenario by quantifying the fraction of recently and rapidly quenched post-starbursts (PSBs) in a sample of post-merger galaxies identified in the Ultraviolet Near Infrared Optical Northern Survey (UNIONS). With our recent development of the Multi-Model Merger Identifier (MUMMI) neural network ensemble (Ferreira et al. 2024a,b), we are now additionally able to predict the time since coalescence (T_PM) for the UNIONS post-merger galaxies up to T_PM = 1.8 Gyr, allowing us to further dissect the merger sequence and measure more precisely when quenching occurs. Based on a sample of 5927 z<0.3 post-mergers identified in UNIONS, we find that the post-coalescence population evolves from one dominated by star-forming (and starbursting) galaxies at 0 < T_PM < 0.16 Gyr, through to a population that is dominated by quenched galaxies by T_PM ~ 1.5 Gyr. We find a PSB excess throughout the post-merger regime, but with a clear peak at 0.16 < T_PM < 0.48 Gyr. In this post-merger time range PSBs are more common than in control galaxies by factors of 30-100, an excess that drops sharply at longer times since merger. We also quantify the fraction of PSBs that are mergers and find that the majority (75%) of classically selected E+A are identified as mergers, with a lower merger fraction (60%) amongst PCA selected PSBs. Our results demonstrate that 1) galaxy-galaxy interactions can lead to rapid post-merger quenching within 0.5 Gyr of coalescence, 2) the majority of (but not all) PSBs at low z are linked to mergers and 3) quenching pathways are diverse, with different PSB selection techniques likely identifying galaxies quenched by different physical processes with an additional dependence on stellar mass. △ Less

Submitted 20 December, 2024; v1 submitted 8 October, 2024; originally announced October 2024.

Comments: Accepted for publication in the Open Journal of Astrophysics. 15 pages, 11 figures. For a temporal view of star formation in the post-merger regime see companion paper Ferreira et al. 2024c

Galaxy evolution in the post-merger regime I -- Most merger-induced in-situ stellar mass growth happens post-coalescence

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

Abstract: Galaxy mergers can enhance star formation rates throughout the merger sequence, with this effect peaking around the time of coalescence. However, owing to a lack of information about their time of coalescence, post-mergers could only previously be studied as a single, time-averaged population. We use timescale predictions of post-coalescence galaxies in the UNIONS survey, based on the Multi-Model… ▽ More Galaxy mergers can enhance star formation rates throughout the merger sequence, with this effect peaking around the time of coalescence. However, owing to a lack of information about their time of coalescence, post-mergers could only previously be studied as a single, time-averaged population. We use timescale predictions of post-coalescence galaxies in the UNIONS survey, based on the Multi-Model Merger Identifier deep learning framework (\textsc{Mummi}) that predicts the time elapsed since the last merging event. For the first time, we capture a complete timeline of star formation enhancements due to galaxy mergers by combining these post-merger predictions with data from pre-coalescence galaxy pairs in SDSS. Using a sample of $564$ galaxies with $M_* \geq 10^{10} M_\odot$ at $0.005 < z < 0.3$ we demonstrate that: 1) galaxy mergers enhance star formation by, on average, up to a factor of two; 2) this enhancement peaks within 500 Myr of coalescence; 3) enhancements continue for up to 1~Gyr after coalescence; and 4) merger-induced star formation significantly contributes to galaxy mass assembly, with galaxies increasing their final stellar masses by, $10\%$ to $20\%$ per merging event, producing on average $\log(M_*/M_\odot) = {9.56_{-0.19}^{+0.13}}$ more mass than non-interacting star-forming galaxies solely due to the excess star formation. △ Less

Submitted 8 October, 2024; originally announced October 2024.

Comments: Letter Submitted to MNRAS, 6 pages, 4 figures. For a temporal view of quenching in the post-merger regime see companion paper Ellison et al. 2024

The effect of image quality on galaxy merger identification with deep learning

Authors: Robert W. Bickley, Scott Wilkinson, Leonardo Ferreira, Sara L. Ellison, Connor Bottrell, Debarpita Jyoti

Abstract: Studies have shown that the morphologies of galaxies are substantially transformed following coalescence after a merger, but post-mergers are notoriously difficult to identify, especially in imaging that is shallow or low-resolution. We train convolutional neural networks (CNNs) to identify simulated post-merger galaxies in a range of image qualities, modelled after five real surveys: the Sloan Di… ▽ More Studies have shown that the morphologies of galaxies are substantially transformed following coalescence after a merger, but post-mergers are notoriously difficult to identify, especially in imaging that is shallow or low-resolution. We train convolutional neural networks (CNNs) to identify simulated post-merger galaxies in a range of image qualities, modelled after five real surveys: the Sloan Digital Sky Survey (SDSS), the Dark Energy Camera Legacy Survey (DECaLS), the Canada-France Imaging Survey (CFIS), the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP), and the Legacy Survey of Space and Time (LSST). Holding constant all variables other than imaging quality, we present the performance of the CNNs on reserved test set data for each image quality. The success of CNNs on a given dataset is found to be sensitive to both imaging depth and resolution. We find that post-merger recovery generally increases with depth, but that limiting 5 sigma point-source depths in excess of ~25 mag, similar to what is achieved in CFIS, are only marginally beneficial. Finally, we present the results of a cross-survey inference experiment, and find that CNNs trained on a given image quality can sometimes be applied to different imaging data to good effect. The work presented here therefore represents a useful reference for the application of CNNs for merger searches in both current and future imaging surveys. △ Less

Submitted 25 September, 2024; originally announced September 2024.

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

Are all starbursts equal? Star-formation-rate profiles in merger versus secular starbursts

Authors: Mallory D. Thorp, Sara L. Ellison, Ana Galicia

Abstract: Galaxy interactions can trigger drastic changes in the resolved star-forming properties of their constituents, but it remains unclear as to whether those changes are discernible from secular starburst triggers. In this Letter we investigate whether or not post-merger galaxies create unique star-forming trends on a kiloparsec scale. We present radial trends in star-formation-rate (SFR) surface dens… ▽ More Galaxy interactions can trigger drastic changes in the resolved star-forming properties of their constituents, but it remains unclear as to whether those changes are discernible from secular starburst triggers. In this Letter we investigate whether or not post-merger galaxies create unique star-forming trends on a kiloparsec scale. We present radial trends in star-formation-rate (SFR) surface density ($Σ_{SFR}$) for 150 post-merger galaxies with moderate to extremely heightened global SFRs using observations from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. We juxtapose these profiles with those of noninteracting galaxies (excluding both galaxy pairs and post-merger galaxies) with similarly enhanced global SFRs. Post-merger galaxies have a much stronger central starburst than isolated galaxies with similar global star-formation enhancements. Indeed, isolated starburst galaxies (SBs) lack a marked central enhancement and instead show a fairly uniform enhancement in $Σ_{SFR}$ with radius. Moreover, the difference in central star formation between post-merger galaxies and noninteracting galaxies is more radially extended and pronounced when the global enhancement in star formation is larger. We conclude that post-merger galaxies create a unique signature in their resolved star-forming properties that is distinct from secular processes that can trigger similar global SFR enhancements. △ Less

Submitted 23 September, 2024; originally announced September 2024.

Comments: Accepted for publication in A&A Letters. 5 pages, 2 figures (+ Appendix 3 pages, 4 figures)

The infrared luminosity of retired and post-starburst galaxies: A cautionary tale for star formation rate measurements

Authors: Vivienne Wild, Natalia Vale Asari, Kate Rowlands, Sara L. Ellison, Ho-Hin Leung, Christy Tremonti

Abstract: In galaxies with significant ongoing star formation there is an impressively tight correlation between total infrared luminosity (L$_{TIR}$) and H$α$ luminosity (L$_{Hα}$), when H$α$ is properly corrected for stellar absorption and dust attenuation. This long-standing result gives confidence that both measurements provide accurate estimates of a galaxy's star formation rate (SFR), despite their di… ▽ More In galaxies with significant ongoing star formation there is an impressively tight correlation between total infrared luminosity (L$_{TIR}$) and H$α$ luminosity (L$_{Hα}$), when H$α$ is properly corrected for stellar absorption and dust attenuation. This long-standing result gives confidence that both measurements provide accurate estimates of a galaxy's star formation rate (SFR), despite their differing origins. To test the extent to which this holds in galaxies with lower specific SFR (sSFR=SFR/M$_{gal}$), we combine optical spectroscopy from the SDSS with multi-wavelength (FUV to FIR) photometric observations from GAMA. We find that L$_{TIR}$/L$_{Hα}$ increases steadily with decreasing H$α$ equivalent width (W$_{Hα}$, a proxy for sSFR), indicating that both luminosities cannot provide a valid measurement of SFR in galaxies below the canonical star-forming sequence. For both `retired galaxies' and `post-starburst galaxies', L$_{TIR}$/L$_{Hα}$ can be up to a factor of 30 larger than for star-forming galaxies. The smooth change in L$_{TIR}$/L$_{Hα}$, irrespective of star formation history, ionisation or heating source, dust temperature or other properties, suggests that the value of L$_{TIR}$/L$_{Hα}$ is determined by the balance between star-forming regions and ambient interstellar medium contributing to both L$_{TIR}$ and H$α$. It is not a result of the differing timescales of star formation that these luminosities probe. While L$_{Hα}$ can only be used to estimate the SFR for galaxies with $W_{Hα}$ > 3A (sSFR $\gtrsim 10^{-11.5}$/yr), we argue that the mid- and far-infrared can only be used to estimate the SFR of galaxies on the star-forming sequence with W$_{Hα}$ >10A (sSFR $\gtrsim 10^{-10.5}$/yr). We find no evidence for dust obscured star-formation in local post-starburst galaxies. △ Less

Submitted 8 January, 2025; v1 submitted 13 September, 2024; originally announced September 2024.

Comments: Accepted in the Open Journal of Astrophysics. 15 pages, 5 figures

Journal ref: Open Journal of Astrophysics, Volume 8, 9th January 2025

The molecular gas content throughout the low-z merger sequence

Authors: Mark T. Sargent, S. L. Ellison, J. T. Mendel, A. Saintonge, D. Cs. Molnár, T. Schwandt, J. M. Scudder, G. Violino

Abstract: Exploiting IRAM 30m CO spectroscopy, we find that SDSS post-merger galaxies display gas fractions and depletion times enhanced by 25-50%, a mildly higher CO excitation, and standard molecular-to-atomic gas ratios, compared to non-interacting galaxies with similar redshift, stellar mass ($M_{\star}$) and star-formation rate (SFR). To place these results in context, we compile further samples of int… ▽ More Exploiting IRAM 30m CO spectroscopy, we find that SDSS post-merger galaxies display gas fractions and depletion times enhanced by 25-50%, a mildly higher CO excitation, and standard molecular-to-atomic gas ratios, compared to non-interacting galaxies with similar redshift, stellar mass ($M_{\star}$) and star-formation rate (SFR). To place these results in context, we compile further samples of interacting or starbursting galaxies, from pre-coalescence kinematic pairs to post-starbursts, carefully homogenising gas mass, $M_{\star}$ and SFR measurements in the process. We explore systematics by duplicating our analysis for different SFR and $M_{\star}$ estimators, finding good qualitative agreement in general. Gas fractions and depletion times are enhanced in interacting pairs, albeit by less than for post-mergers. Among all samples studied, gas fraction and depletion time enhancements appear largest in young (a few 100 Myr) post-starbursts. While there is only partial overlap between post-mergers and post-starbursts, this suggests that molecular gas reservoirs are boosted throughout most stages of galaxy interactions, plausibly due to torque-driven inflows of halo gas and gas compression. The gas fraction and depletion time offsets of mergers and post-starbursts anti-correlate with distance from the galaxy main sequence $Δ({\rm MS})$, evidencing the role of SFE in driving the high SFRs of the strongest starbursts. Post-starbursts display the steepest dependency of gas fraction and SFE-offsets on $Δ({\rm MS})$, with an evolving normalisation that reflects gas reservoir depletion over time. Our multi-sample analysis paints a coherent picture of the starburst-merger connection throughout the low-z merger sequence. It reconciles contradictory literature findings by highlighting that gas fraction enhancements and SFE variations both play their part in merger-driven star formation. △ Less

Submitted 31 March, 2025; v1 submitted 10 September, 2024; originally announced September 2024.

Comments: 48 pages, 21 figures. Accepted for publication in the Open Journal of Astrophysics

The First Large Absorption Survey in HI (FLASH): II. Pilot Survey data release and first results

Authors: Hyein Yoon, Elaine M. Sadler, Elizabeth K. Mahony, J. N. H. S. Aditya, James R. Allison, Marcin Glowacki, Emily F. Kerrison, Vanessa A. Moss, Renzhi Su, Simon Weng, Matthew Whiting, O. Ivy Wong, Joseph R. Callingham, Stephen J. Curran, Jeremy Darling, Alastair C. Edge, Sara L. Ellison, Kimberly L. Emig, Lilian Garratt-Smithson, Gordon German, Kathryn Grasha, Baerbel S. Koribalski, Raffaella Morganti, Tom Oosterloo, Céline Péroux , et al. (19 additional authors not shown)

Abstract: The First Large Absorption Survey in HI (FLASH) is a large-area radio survey for neutral hydrogen in the redshift range 0.4<z<1.0, using the 21cm HI absorption line as a probe of cold neutral gas. FLASH uses the ASKAP radio telescope and is the first large 21cm absorption survey to be carried out without any optical preselection of targets. We use an automated Bayesian line-finding tool to search… ▽ More The First Large Absorption Survey in HI (FLASH) is a large-area radio survey for neutral hydrogen in the redshift range 0.4<z<1.0, using the 21cm HI absorption line as a probe of cold neutral gas. FLASH uses the ASKAP radio telescope and is the first large 21cm absorption survey to be carried out without any optical preselection of targets. We use an automated Bayesian line-finding tool to search through large datasets and assign a statistical significance to potential line detections. Two Pilot Surveys, covering around 3000 deg$^2$ of sky, were carried out in 2019-22 to test and verify the strategy for the full FLASH survey. The processed data from these Pilot Surveys (spectral-line cubes, continuum images, and catalogues) are available online. Here, we describe the FLASH spectral-line and continuum data and discuss the quality of the HI spectra and the completeness of our automated line search. Finally, we present a set of 30 new HI absorption lines that were robustly detected in the Pilot Surveys. These lines span a wide range in HI optical depth, including three lines with a peak optical depth $τ>1$, and appear to be a mixture of intervening and associated systems. Interestingly, around two-thirds of the lines found in this untargeted sample are detected against sources with a peaked-spectrum radio continuum, which are only a minor (5-20%) fraction of the overall radio-source population. The overall detection rate for HI absorption lines in the Pilot Surveys (0.3 to 0.5 lines per ASKAP field) is a factor of two below the expected value. One possible reason for this is the presence of a range of spectral-line artefacts in the Pilot Survey data that have now been mitigated and are not expected to recur in the full FLASH survey. A future paper will discuss the host galaxies of the HI absorption systems identified here. △ Less

Submitted 22 May, 2025; v1 submitted 13 August, 2024; originally announced August 2024.

Comments: Accepted for publication in PASA. 44 pages, 24 figures, 7 tables

Galaxy Mergers in UNIONS -- I: A Simulation-driven Hybrid Deep Learning Ensemble for Pure Galaxy Merger Classification

Authors: Leonardo Ferreira, Robert W. Bickley, Sara L. Ellison, David R. Patton, Shoshannah Byrne-Mamahit, Scott Wilkinson, Connor Bottrell, Sébastien Fabbro, Stephen D. J. Gwyn, Alan McConnachie

Abstract: Merging and interactions can radically transform galaxies. However, identifying these events based solely on structure is challenging as the status of observed mergers is not easily accessible. Fortunately, cosmological simulations are now able to produce more realistic galaxy morphologies, allowing us to directly trace galaxy transformation throughout the merger sequence. To advance the potential… ▽ More Merging and interactions can radically transform galaxies. However, identifying these events based solely on structure is challenging as the status of observed mergers is not easily accessible. Fortunately, cosmological simulations are now able to produce more realistic galaxy morphologies, allowing us to directly trace galaxy transformation throughout the merger sequence. To advance the potential of observational analysis closer to what is possible in simulations, we introduce a supervised deep learning Convolutional Neural Network (CNN) and Vision Transformer (ViT) hybrid framework, Mummi (MUlti Model Merger Identifier). Mummi is trained on realism-added synthetic data from IllustrisTNG100-1, and is comprised of a multi-step ensemble of models to identify mergers and non-mergers, and to subsequently classify the mergers as interacting pairs or post-mergers. To train this ensemble of models, we generate a large imaging dataset of 6.4 million images targeting UNIONS with RealSimCFIS. We show that Mummi offers a significant improvement over many previous machine learning classifiers, achieving 95% pure classifications even at Gyr long timescales when using a jury-based decision making process, mitigating class imbalance issues that arise when identifying real galaxy mergers from $z=0$ to $0.3$. Additionally, we can divide the identified mergers into pairs and post-mergers at 96% success rate. We drastically decrease the false positive rate in galaxy merger samples by 75%. By applying Mummi to the UNIONS DR5-SDSS DR7 overlap, we report a catalog of 13,448 high confidence galaxy merger candidates. Finally, we demonstrate that Mummi produces powerful representations solely using supervised learning, which can be used to bridge galaxy morphologies in simulations and observations. △ Less

Submitted 25 July, 2024; originally announced July 2024.

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

Correction to Conversions between gas-phase metallicities in MaNGA

Authors: Jillian M. Scudder, Sara L. Ellison, Loubna El Meddah El Idrissi, Henry Poetrodjojo

Abstract: We present a brief correction to Scudder et al. (2021) due to an error in the O3N2 based metallicity calibrations presented in that work. Conclusions are unchanged, but metallicity values shift by $\sim$0.03 dex, with polynomials for affected conversions shifted by the same amount. We present updated materials here. We present a brief correction to Scudder et al. (2021) due to an error in the O3N2 based metallicity calibrations presented in that work. Conclusions are unchanged, but metallicity values shift by $\sim$0.03 dex, with polynomials for affected conversions shifted by the same amount. We present updated materials here. △ Less

Submitted 22 July, 2024; originally announced July 2024.

Comments: 4 pages, 3 figures. Supplementary materials available online. Accepted for publication in MNRAS

Does the Fundamental Metallicity Relation Evolve with Redshift? II: The Evolution in Normalisation of the Mass-Metallicity Relation

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

Abstract: The metal content of galaxies is a direct probe of the baryon cycle. A hallmark example is the relationship between a galaxy's stellar mass, star formation rate (SFR), and gas-phase metallicity: the Fundamental Metallicity Relation (FMR). While low-redshift ($z\lesssim4$) observational studies suggest that the FMR is redshift-invariant, recent high-$z$ JWST data indicate deviations from the FMR es… ▽ More The metal content of galaxies is a direct probe of the baryon cycle. A hallmark example is the relationship between a galaxy's stellar mass, star formation rate (SFR), and gas-phase metallicity: the Fundamental Metallicity Relation (FMR). While low-redshift ($z\lesssim4$) observational studies suggest that the FMR is redshift-invariant, recent high-$z$ JWST data indicate deviations from the FMR established at low-$z$. In this study, we utilize the FMR to predict the evolution of the normalisation of the mass-metallicity relation (MZR) using the cosmological simulations Illustris, IllustrisTNG, EAGLE, and SIMBA. Our findings demonstrate that a $z = 0$ calibrated FMR struggles to predict the evolution in the MZR of each simulation. To quantify the divergence of the predictions, we introduce the concepts of a ``static'' FMR, where the role of the SFR in setting the normalization of the MZR does not change with redshift, and a ``dynamic'' FMR, where the role of SFR evolves over time. We find static FMRs in SIMBA and dynamic FMRs in Illustris, IllustrisTNG and EAGLE. We suggest that the differences between these models likely points to the subtle differences in the implementation of the baryon cycle. Moreover, we echo recent JWST results at $z > 4$ by finding significant offsets from the FMR in IllustrisTNG and EAGLE, suggesting that the observed FMR may have a similar dynamic trend as these simulations. Overall, our findings imply that the current FMR framework neglects important time variations of these simulations' baryon cycles. △ Less

Submitted 14 November, 2024; v1 submitted 8 July, 2024; originally announced July 2024.

Comments: Accepted to MNRAS. 17 pages, 8 Figures. See also paper one of series (arXiv:2403.08856)

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

Report number: 2024SCPMA..6799811Y

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 11 August, 2024; v1 submitted 30 January, 2024; originally announced January 2024.

Comments: 20 pages, 14 figures. Accepted for publication in MNRAS, post-merger catalog will be available via the journal article

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

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

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

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