Search

In galaxies with significant ongoing star formation there is an impressively tight correlation between total infrared luminosity (L$_{TIR}$) and H$\alpha$ luminosity (L$_{H\alpha}$), when H$\alpha$ 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/Mgal, where Mgal is the stellar mass), we combine optical spectroscopy from the Sloan Digital Sky Survey (SDSS) with multi-wavelength (FUV to FIR) photometric observations from the Galaxy And Mass Assembly survey (GAMA). We find that L$_{TIR}$/L$_{H\alpha}$increases steadily with decreasing H$\alpha$ equivalent width (W$_{H\alpha}$, 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\alpha}$ can be up to a factor of 30 larger than for star-forming galaxies. The smooth change in L$_{TIR}$/L$_{H\alpha}$, irrespective of star formation history, ionisation or heating source, dust temperature or other properties, suggests that the value of L$_{TIR}$/L$_{H\alpha}$ is given by the balance between star-forming regions and ambient interstellar medium contributing to both L$_{TIR}$ and L$_{H\alpha}$. While L$_{H\alpha}$ can only be used to estimate the SFR for galaxies with W$_{H\alpha}$ > 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, and in particular only for galaxies with W$_{H\alpha}$ >10 A (sSFR $\gtrsim 10^{-10.5}$/yr). We find no evidence for dust obscured star-formation in post-starburst galaxies., Comment: Submitted to the Open Journal of Astrophysics. Comments welcome. 14 pages, 4 figures

Exploiting IRAM 30 m 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. Molecular gas fractions and depletion times are enhanced in interacting pairs, albeit 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 their distance from the galaxy main sequence $\Delta({\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 $\Delta({\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 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., Comment: 46 pages, 21 figures. Submitted to the Open Journal of Astrophysics

The First Large Absorption Survey in HI (FLASH) is a large-area radio survey for neutral hydrogen in the redshift range 0.41$, and appear to be a mixture of intervening and associated systems. 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. There are several possible reasons for this, but one likely factor 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., Comment: 46 pages, 25 figures, 10 tables. Submitted to PASA

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., Comment: 23 pages, 16 figures. Accepted for publication in MNRAS

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., Comment: 4 pages, 3 figures. Supplementary materials available online. Accepted for publication in MNRAS

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 JWST data indicate deviations from this model. 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 Illustris and SIMBA and dynamic FMRs in 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 be dynamic as well. Overall, our findings imply that the current FMR framework neglects important variations in the baryon cycle through cosmic time., Comment: Submitted to MNRAS. 17 pages, 8 Figures. See also paper one of series (arXiv:2403.08856). Comments welcome

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$\alpha$ 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., Comment: 19 pages, 15 figures, 4 tables, accepted for publication in MNRAS

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 $\mu_*$ with a Kendall's $\tau$ 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 $\tau$ 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., Comment: Accepted by SCPMA. We sincerely appreciate the constructive suggestions and comments from Lihwai Lin

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., Comment: 15 pages, 16 figures. Accepted for publication in MNRAS

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 ($\Sigma_\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 $\Sigma_\mathrm{sSFR}$ at different galactocentric radii. Analysis of relative contributions of $f_\mathrm{gas}$ and SFE to $\Sigma_\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., Comment: 27 pages, 8 figures, 1 table. Accepted for publication in The Astrophysical Journal

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 $< \mu <$ 1:1 and a redshift between $010^{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 $010^{45}$ erg/s, $\sim$55\% are interacting, compared with a 30\% interaction fraction in AGN with $10^{44}

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., Comment: 20 pages, 14 figures. Accepted for publication in MNRAS, post-merger catalog will be available via the journal article

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., Comment: 32 pages, 21 figures Accepted for publication by MNRAS

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., Comment: Accepted to MNRAS

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., Comment: 20 pages, 13 figures, ApJ accepted

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., Comment: 20 pages, 20 figures

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., Comment: Accepted for publication in MNRAS

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, $\Sigma_{\mathrm{SFR}}$) overcomes this ambiguity. $\Sigma_{\mathrm{SFR}}$ has a strong physical basis, being tied to the molecular gas density and the effectiveness of stellar feedback, so we propose $\Sigma_{\mathrm{SFR}}$-M* as an important galaxy evolution diagram to complement (s)SFR-M* diagrams. Using the $\Sigma_{\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 "$\Sigma_{\mathrm{SFR}}$ main sequence" steepens, signaling the epoch of bulge build-up in massive galaxies. We also find that $\Sigma_{\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 $\Sigma_{\mathrm{SFR}}$ radial profiles, because ample SF may be present in bulges with low sSFR (red color)., Comment: 16 pages. Accepted to ApJ. Comments on content or relevant missing references welcome

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., Comment: 9 pages, 4 figures, accepted in ApJL

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., Comment: 17 pages, 1 table, 5 figures, accepted for publication in ApJ

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., Comment: 16 pages, 12 figures. Accepted for publication in MNRAS

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$\sigma$ 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., Comment: Accepted for publication in MNRAS. 17 pages, 7 figures

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 (${\Sigma}_{\star}$), HI (${\Sigma}_\mathrm{HI}$), and molecular gas (${\Sigma}_\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 ${\Sigma}_\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 ${\Sigma}_\mathrm{mol}$ gas is removed first. In the most HI-deficient galaxies, however, we find evidence that environmental processes reduce the typical ${\Sigma}_\mathrm{mol}$ of the remaining gas by nearly a factor of 3. We find no evidence for environment-driven elevation of ${\Sigma}_\mathrm{HI}$ or ${\Sigma}_\mathrm{mol}$ in HI-deficient galaxies. Using the ratio of ${\Sigma}_\mathrm{mol}$-to-${\Sigma}_\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., Comment: 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

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., Comment: 13 pages, 5 figures

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., Comment: This article has been accepted for publication in Monthly Notices of the Royal Astronomical Society

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$\mu$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\mu$m and $12\mu$m bands, with strong PAH features, show steeper CO vs. mid-IR slopes than the $22\mu$m and $24\mu$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., Comment: 29 pages, 13 figures, key quantitative results in Table 3, Accepted as part of a PHANGS-JWST Focus Issue to appear in ApJ

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

We describe the scientific goals and survey design of the First Large Absorption Survey in HI (FLASH), a wide field survey for 21-cm line absorption in neutral atomic hydrogen (HI) at intermediate cosmological redshifts. FLASH will be carried out with the Australian Square Kilometre Array Pathfinder (ASKAP) radio telescope and is planned to cover the sky south of $\delta \approx +40$deg at frequencies between 711.5 and 999.5MHz. At redshifts between $z = 0.4$ and $1.0$ (look back times of 4 - 8Gyr), the HI content of the Universe has been poorly explored due to the difficulty of carrying out radio surveys for faint 21-cm line emission and, at ultra-violet wavelengths, space-borne searches for Damped Lyman-$\alpha$ absorption in quasar spectra. The ASKAP wide field of view and large spectral bandwidth, in combination with a radio-quiet site, will enable a search for absorption lines in the radio spectra of bright continuum sources over 80% of the sky. This survey is expected to detect at least several hundred intervening 21-cm absorbers, and will produce an HI-absorption-selected catalogue of galaxies rich in cool, star-forming gas, some of which may be concealed from optical surveys. Likewise, at least several hundred associated 21-cm absorbers are expected to be detected within the host galaxies of radio sources at $0.4 < z < 1.0$, providing valuable kinematical information for models of gas accretion and jet-driven feedback in radio-loud active galactic nuclei. FLASH will also detect OH 18-cm absorbers in diffuse molecular gas, megamaser OH emission, radio recombination lines, and stacked HI emission., Comment: 31 pages, 10 figures, 6 tables, accepted for publication in PASA

We present polynomial conversions between each of 11 different strong line gas-phase metallicity calibrations, each based on $\sim$ 1.1 million star-forming spaxels in the public Sloan Digital Sky Survey (SDSS) Data Release 15 (DR15) Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. For this sample, which is $\sim$ 20 times larger than previous works, we present 5th order polynomial fits for each of 110 possible calibration conversions, for both Small Magellanic Cloud (SMC)-type and Milky Way (MW)-type dust corrections. The typical $2\sigma$ scatter around our polynomial fits is 0.1 dex; we present the range over which the metallicities are valid. Conversions between metallicities which rely on the same set of line ratios, or a heavily shared set of emission lines, have reduced scatter in their conversions relative to those conversions with little overlap in required emission lines. Calibration conversions with less consistent sets of emission lines also have increased galaxy-to-galaxy variability, and this variability can account for up to 35% of the total scatter. We also compare our conversions to previous work with the single fibre SDSS DR7 spectra along with higher spatial resolution data from the TYPHOON Integral Field Spectroscopy survey, resulting in comparison samples with spatial resolutions from several kpc down to $\sim$100 pc. Our metallicity conversions, obtained with the large sample of MaNGA, are robust against the influence of diffuse ionized gas, redshift, effective radius and spatial blurring, and are therefore consistent across both integrated spectra and the high resolution integral field spectroscopy data., Comment: MNRAS accepted. 20 pages, 14 figures, 4 tables. Supplementary material available online