Search

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 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

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

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.

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

Quantitative morphologies, such as asymmetry and concentration, have long been used as an effective way to assess the distribution of galaxy starlight in large samples. Application of such quantitative indicators to other data products could provide a tool capable of capturing the 2-dimensional distribution of a range of galactic properties, such as stellar mass or star-formation rate maps. In this work, we utilize galaxies from the Illustris and IllustrisTNG simulations to assess the applicability of concentration and asymmetry indicators to the stellar mass distribution in galaxies. Specifically, we test whether the intrinsic values of concentration and asymmetry (measured directly from the simulation stellar mass particle maps) are recovered after the application of measurement uncertainty and a point spread function (PSF). We find that random noise has a non-negligible systematic effect on asymmetry that scales inversely with signal-to-noise, particularly at signal-to-noise less than 100. We evaluate different methods to correct for the noise contribution to asymmetry at very low signal-to-noise, where previous studies have been unable to explore due to systematics. We present algebraic corrections for noise and resolution to recover the intrinsic morphology parameters. Using Illustris as a comparison dataset, we evaluate the robustness of these fits in the presence of a different physics model, and confirm these correction methods can be applied to other datasets. Lastly, we provide estimations for the uncertainty on different correction methods at varying signal-to-noise and resolution regimes., Comment: 19 pages, 18 figures, accepted for publication in MNRAS

In order to investigate the role of gas in the demise of star formation on kpc-scales, we compare the resolved molecular gas main sequence (rMGMS: Sigma_* vs Sigma_H2) of star-forming regions to the sequence of `retired' regions that have ceased to form new stars. Using data from the ALMaQUEST survey, we find that retired spaxels form a rMGMS that is distinct from that of star-forming spaxels, offset to lower Sigma_H2 at fixed Sigma_* by a factor of ~5. We study the rMGMS of star-forming and retired spaxels on a galaxy-by-galaxy basis for eight individual ALMaQUEST galaxies. Six of these galaxies have their retired spaxels concentrated within the central few kpc. Molecular gas is detected in 40-100% of retired spaxels in the eight galaxies in our sample. Both the star-forming and retired rMGMS show a diversity in normalization from galaxy-to-galaxy. However, in any given galaxy, the rMGMS for retired regions is found to be distinct from the star-forming sequence and gas fractions of retired spaxels are up to an order of magnitude lower than the star-forming spaxels. We conclude that quenching is associated with a depletion (but not absence) of molecular gas via a mechanism that typically begins in the centre of the galaxy., Comment: MNRAS, in press

Using a sample of ~15,000 kpc-scale star-forming spaxels in 28 galaxies drawn from the ALMA-MaNGA QUEnching and STar formation (ALMaQUEST) survey, we investigate the galaxy-to-galaxy variation of the `resolved' Schmidt-Kennicutt relation (rSK; Sigma_H2 - Sigma_SFR), the `resolved' star forming main sequence (rSFMS; Sigma_* - Sigma_SFR) and the `resolved' molecular gas main sequence (rMGMS; Sigma_* - Sigma_H2). The rSK relation, rSFMS and rMGMS all show significant galaxy-to-galaxy variation in both shape and normalization, indicating that none of these relations is universal between galaxies. The rSFMS shows the largest galaxy-to-galaxy variation and the rMGMS the least. By defining an `offset' from the average relations, we compute a Delta_rSK, Delta_rSFMS, Delta_rMGMS for each galaxy, to investigate correlations with global properties. We find the following correlations with at least 2 sigma significance: the rSK is lower (i.e. lower star formation efficiency) in galaxies with higher M_*, larger Sersic index and lower specific SFR (sSFR); the rSFMS is lower (i.e. lower sSFR) in galaxies with higher M_* and larger Sersic index; the rMGMS is lower (i.e. lower gas fraction) in galaxies with lower sSFR. In the ensemble of all 15,000 data points, the rSK relation and rMGMS show equally tight scatters and strong correlation coefficients, compared with a larger scatter and weaker correlation in the rSFMS. Moreover, whilst there is no correlation between Delta_rSK and Delta_rMGMS in the sample, the offset of a galaxy's rSFMS does correlate with both of the other two offsets. Our results therefore indicate that the rSK and rMGMS are independent relations, whereas the rSFMS is a result of their combination., Comment: MNRAS, in press

The ALMaQUEST (ALMA-MaNGA QUEnching and STar formation) survey is a program with spatially-resolved $^{12}$CO(1-0) measurements obtained with the Atacama Large Millimeter Array (ALMA) for 46 galaxies selected from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) DR15 optical integral-field spectroscopic survey. The aim of the ALMaQUEST survey is to investigate the dependence of star formation activity on the cold molecular gas content at kpc scales in nearby galaxies. The sample consists of galaxies spanning a wide range in specific star formation rate (sSFR), including starburst (SB), main-sequence (MS), and green valley (GV) galaxies. In this paper, we present the sample selection and characteristics of the ALMA observations, and showcase some of the key results enabled by the combination of spatially-matched stellar populations and gas measurements. Considering the global (aperture-matched) stellar mass, molecular gas mass, and star formation rate of the sample, we find that the sSFR depends on both the star formation efficiency (SFE) and the molecular gas fraction ($f_{\rm H_{2}}$), although the correlation with the latter is slightly weaker. Furthermore, the dependence of sSFR on the molecular gas content (SFE or $f_{\rm H_{2}}$) is stronger than that on either the atomic gas fraction or the molecular-to-atomic gas fraction, albeit with the small HI sample size. On kpc scales, the variations in both SFE and $f_{\rm H_{2}}$ within individual galaxies can be as large as 1-2 dex thereby demonstrating that the availability of spatially-resolved observations is essential to understand the details of both star formation and quenching processes., Comment: 20 pages. ApJS accepted. Complete figure sets of Figures 7 & 8 are available in the online journal

We investigate how star formation quenching proceeds within central and satellite galaxies using spatially resolved spectroscopy from the SDSS-IV MaNGA DR15. We adopt a complete sample of star formation rate surface densities ($\Sigma_{\rm SFR}$), derived in Bluck et al. (2020), to compute the distance at which each spaxel resides from the resolved star forming main sequence ($\Sigma_{\rm SFR} - \Sigma_*$ relation): $\Delta \Sigma_{\rm SFR}$. We study galaxy radial profiles in $\Delta \Sigma_{\rm SFR}$, and luminosity weighted stellar age (${\rm Age_L}$), split by a variety of intrinsic and environmental parameters. Via several statistical analyses, we establish that the quenching of central galaxies is governed by intrinsic parameters, with central velocity dispersion ($\sigma_c$) being the most important single parameter. High mass satellites quench in a very similar manner to centrals. Conversely, low mass satellite quenching is governed primarily by environmental parameters, with local galaxy over-density ($\delta_5$) being the most important single parameter. Utilising the empirical $M_{BH}$ - $\sigma_c$ relation, we estimate that quenching via AGN feedback must occur at $M_{BH} \geq 10^{6.5-7.5} M_{\odot}$, and is marked by steeply rising $\Delta \Sigma_{\rm SFR}$ radial profiles in the green valley, indicating `inside-out' quenching. On the other hand, environmental quenching occurs at over-densities of 10 - 30 times the average galaxy density at z$\sim$0.1, and is marked by steeply declining $\Delta \Sigma_{\rm SFR}$ profiles, indicating `outside-in' quenching. Finally, through an analysis of stellar metallicities, we conclude that both intrinsic and environmental quenching must incorporate significant starvation of gas supply., Comment: Accepted for publication in MNRAS. 40 pages; 24 figures

Starburst galaxies have elevated star formation rates (SFRs) for their stellar mass. In Ellison et al. (2018) we used integral field unit (IFU) maps of star formation rate surface density (Sigma_SFR) and stellar mass surface density (Sigma_*) to show that starburst galaxies in the local universe are driven by SFRs that are preferentially boosted in their central regions. Here, we present molecular gas maps obtained with the Atacama Large Millimeter Array (ALMA) observatory for 12 central starburst galaxies at z~0 drawn from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. The ALMA and MaNGA data are well matched in spatial resolution, such that the ALMA maps of molecular gas surface density (Sigma_H2) can be directly compared with MaNGA maps at kpc-scale resolution. The combination of Sigma_H2, Sigma_* and Sigma_SFR at the same resolution allow us to investigate whether central starbursts are driven primarily by enhancements in star formation efficiency (SFE) or by increased gas fractions. By computing offsets from the resolved Kennicutt-Schmidt relation (Sigma_H2 vs. Sigma_SFR) and the molecular gas main sequence (Sigma_* vs. Sigma_H2), we conclude that the primary driver of the central starburst is an elevated SFE. We also show that the enhancement in Sigma_SFR is accompanied by a dilution in O/H, consistent with a triggering that is induced by metal poor gas inflow. These observational signatures are found in both undisturbed (9/12 galaxies in our sample) and recently merged galaxies, indicating that both interactions and secular mechanisms contribute to central starbursts., Comment: Accepted for publication in MNRAS

Using a sample of 11,478 spaxels in 34 galaxies with molecular gas, star formation and stellar maps taken from the ALMA-MaNGA QUEnching and STar formation (ALMaQUEST) survey, we investigate the parameters that correlate with variations in star formation rates on kpc scales. We use a combination of correlation statistics and an artificial neural network to quantify the parameters that drive both the absolute star formation rate surface density (Sigma_SFR), as well as its scatter around the resolved star forming main sequence (Delta Sigma_SFR). We find that Sigma_SFR is primarily regulated by molecular gas surface density (Sigma_H2) with a secondary dependence on stellar mass surface density (Sigma_*), as expected from an `extended Kennicutt-Schmidt relation'. However, Delta Sigma_SFR is driven primarily by changes in star formation efficiency (SFE), with variations in gas fraction playing a secondary role. Taken together, our results demonstrate that whilst the absolute rate of star formation is primarily set by the amount of molecular gas, the variation of star formation rate above and below the resolved star forming main sequence (on kpc scales) is primarily due to changes in SFE., Comment: Accepted for publication in MNRAS

We present an analysis of star formation and quenching in the SDSS-IV MaNGA-DR15, utilising over 5 million spaxels from $\sim$3500 local galaxies. We estimate star formation rate surface densities ($\Sigma_{\rm SFR}$) via dust corrected $H\alpha$ flux where possible, and via an empirical relationship between specific star formation rate (sSFR) and the strength of the 4000 Angstrom break (D4000) in all other cases. We train a multi-layered artificial neural network (ANN) and a random forest (RF) to classify spaxels into `star forming' and `quenched' categories given various individual (and groups of) parameters. We find that global parameters (pertaining to the galaxy as a whole) perform collectively the best at predicting when spaxels will be quenched, and are substantially superior to local/ spatially resolved and environmental parameters. Central velocity dispersion is the best single parameter for predicting quenching in central galaxies. We interpret this observational fact as a probable consequence of the total integrated energy from AGN feedback being traced by the mass of the black hole, which is well known to correlate strongly with central velocity dispersion. Additionally, we train both an ANN and RF to estimate $\Sigma_{\rm SFR}$ values directly via regression in star forming regions. Local/ spatially resolved parameters are collectively the most predictive at estimating $\Sigma_{\rm SFR}$ in these analyses, with stellar mass surface density at the spaxel location ($\Sigma_*$) being by far the best single parameter. Thus, quenching is fundamentally a global process but star formation is governed locally by processes within each spaxel., Comment: Accepted for publication in MNRAS. 45 pages; 20 figures

The origin of the star forming main sequence ( i.e., the relation between star formation rate and stellar mass, globally or on kpc-scales; hereafter SFMS) remains a hotly debated topic in galaxy evolution. Using the ALMA-MaNGA QUEnching and STar formation (ALMaQUEST) survey, we show that for star forming spaxels in the main sequence galaxies, the three local quantities, star-formation rate surface density (\sigsfr), stellar mass surface density (\sigsm), and the \h2~mass surface density (\sigh2), are strongly correlated with one another and form a 3D linear (in log) relation with dispersion. In addition to the two well known scaling relations, the resolved SFMS (\sigsfr~ vs. \sigsm) and the Schmidt-Kennicutt relation (\sigsfr~ vs. \sigh2; SK relation), there is a third scaling relation between \sigh2~ and \sigsm, which we refer to as the `molecular gas main sequence' (MGMS). The latter indicates that either the local gas mass traces the gravitational potential set by the local stellar mass or both quantities follow the underlying total mass distributions. The scatter of the resolved SFMS ($\sigma \sim 0.25$ dex) is the largest compared to those of the SK and MGMS relations ($\sigma \sim$ 0.2 dex). A Pearson correlation test also indicates that the SK and MGMS relations are more strongly correlated than the resolved SFMS. Our result suggests a scenario in which the resolved SFMS is the least physically fundamental and is the consequence of the combination of the SK and the MGMS relations., Comment: 4 figures, 1 table, accepted by ApJL

Large galaxy surveys have demonstrated that galaxy-galaxy mergers can dramatically change the morphologies, star formation rates (SFRs) and metallicities of their constituents. However, most statistical studies have been limited to the measurement of global quantities, through large fibres or integrated colours. In this work, we present the first statistically significant study of spatially resolved star formation and metallicity profiles using integral field spectroscopy, using a sample of $\sim$20,000 spaxels in 36 visually selected post-merger galaxies from the Mapping Nearby Galaxies with Apache Point Observatory (MaNGA) survey. By measuring offsets from SFR and metallicity scaling relations on a spaxel-by-spaxel basis, we are able to quantify where in the galaxy these properties are most affected by the interaction. We find that the SFR enhancements are generally centrally peaked, by a factor of 2.5 on average, in agreement with predictions from simulations. However, there is considerable variation in the SFR behaviour in the galactic outskirts, where both enhancement and suppression are seen. The median SFR remains enhanced by 0.1 dex out to at least 1.9 $R_{e}$. The metallicity is also affected out to these large radii, typically showing a suppression of $\sim$-0.04dex., Comment: 6 pages, 4 figures, accepted for publication in MNRAS Letters