Researcher profile

Sara L. Ellison

Sara L. Ellison contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
22works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2026arXiv

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

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

preprint2022arXiv

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

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

preprint2022arXiv

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

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

preprint2022arXiv

Quenching Timescales in the IllustrisTNG Simulation

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

preprint2022arXiv

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

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.

preprint2022arXiv

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

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

preprint2022arXiv

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

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

preprint2022arXiv

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

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

preprint2022arXiv

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

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.

preprint2022arXiv

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

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&#39; 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.

preprint2021arXiv

A re-assessment of strong line metallicity conversions in the machine learning era

Strong line metallicity calibrations are widely used to determine the gas phase metallicities of individual HII regions and entire galaxies. Over a decade ago, based on the Sloan Digital Sky Survey Data Release 4 (SDSS DR4), Kewley \& Ellison published the coefficients of third-order polynomials that can be used to convert between different strong line metallicity calibrations for global galaxy spectra. Here, we update the work of Kewley \& Ellison in three ways. First, by using a newer data release (DR7), we approximately double the number of galaxies used in polynomial fits, providing statistically improved polynomial coefficients. Second, we include in the calibration suite five additional metallicity diagnostics that have been proposed in the last decade and were not included by Kewley \& Ellison. Finally, we develop a new machine learning approach for converting between metallicity calibrations. The random forest algorithm is non-parametric and therefore more flexible than polynomial conversions, due to its ability to capture non-linear behaviour in the data. The random forest method yields the same accuracy as the (updated) polynomial conversions, but has the significant advantage that a single model can be applied over a wide range of metallicities, without the need to distinguish upper and lower branches in $R_{23}$ calibrations. The trained random forest is made publicly available for use in the community.

preprint2021arXiv

Sub-damped Lyman alpha systems in the XQ-100 survey II -- Chemical evolution at 2.4<z<4.3

We present the measured gas-phase metal column densities in 155 sub-damped Lyman alpha systems (subDLAs) with the aim to investigate the contribution of subDLAs to the chemical evolution of the Universe. The sample was identified within the absorber-blind XQ-100 quasar spectroscopic survey over the redshift range 2.4<=z<=4.3. Using all available column densities of the ionic species investigated (mainly CIV, SiII, MgII, SiIV, AlII, FeII, CII, and OI; in order of decreasing detection frequency), we estimate the ionization-corrected gas-phase metallicity of each system using Markov Chain Monte Carlo techniques to explore a large grid of Cloudy ionization models. Without accounting for ionization and dust depletion effects, we find that the HI-weighted gas-phase metallicity evolution of subDLAs are consistent with damped Lyman alpha systems (DLAs). When ionization corrections are included, subDLAs are systematically more metal-poor than DLAs (between ~0.5 sigma and ~3 sigma significance) by up to ~1.0 dex over the redshift range 3<=z<=4.3. The correlation of gas-phase [Si/Fe] with metallicity in subDLAs appears to be consistent with that of DLAs, suggesting that the two classes of absorbers have a similar relative dust depletion pattern. As previously seen for Lyman limit systems, the gas-phase [C/O] in subDLAs remains constantly solar for all metallicities indicating that both subDLAs and Lyman limit systems could trace carbon-rich ejecta, potentially in circumgalactic environments.

preprint2020arXiv

A successful search for intervening 21 cm HI absorption in galaxies at 0.4 < z <1.0 with the Australian Square Kilometre Array Pathfinder (ASKAP)

We have used the Australian Square Kilometre Array Pathfinder (ASKAP) radio telescope to search for intervening 21 cm neutral hydrogen (HI) absorption along the line of sight to 53 bright radio continuum sources. Our observations are sensitive to HI column densities typical of Damped Lyman Alpha absorbers (DLAs) in cool gas with an HI spin temperature below about 300-500 K. The six-dish Boolardy Engineering Test Array (BETA) and twelve-antenna Early Science array (ASKAP-12) covered a frequency range corresponding to redshift $0.4<z<1.0$ and $0.37<z<0.77$ respectively for the HI line. Fifty of the 53 radio sources observed have reliable optical redshifts, giving a total redshift path $Δz$ = 21.37. This was a spectroscopically-untargeted survey, with no prior assumptions about the location of the lines in redshift space. Four intervening HI lines were detected, two of them new. In each case, the estimated HI column density lies above the DLA limit for HI spin temperatures above 50-80 K, and we estimate a DLA number density at redshift $z\sim0.6$ of $n(z)=0.19\substack{+0.15 \\ -0.09}$. This value lies somewhat above the general trend of $n(z)$ with redshift seen in optical DLA studies. Although the current sample is small, it represents an important proof of concept for the much larger 21cm First Large Absorption Survey in HI (FLASH) project to be carried out with the full 36-antenna ASKAP telescope, probing a total redshift path $Δz\sim\,50,000$.

preprint2020arXiv

Galaxy mergers in EAGLE do not induce a significant amount of black hole growth yet do increase the rate of luminous AGN

We investigate the connection between galaxy--galaxy mergers and enhanced black hole (BH) growth using the cosmological hydrodynamical EAGLE simulation. We do this via three methods of analysis, investigating: the merger fraction of AGN, the AGN fraction of merging systems and the AGN fraction of galaxies with close companions. In each case, we find an increased abundance of AGN within merging systems relative to control samples of inactive or isolated galaxies (by up to a factor of $\approx 3$ depending on the analysis method used), confirming that mergers are enhancing BH accretion rates for at least a subset of the galaxy population. The greatest excess of AGN triggered via a merger are found in lower mass ($M_* \sim 10^{10}$~\Msol) gas rich ($f_{\mathrm{gas}} > 0.2$) central galaxies with lower mass BHs ($M_{\mathrm{BH}} \sim 10^{7}$~\Msol) at lower redshifts ($z<1$). We find no enhancement of AGN triggered via mergers in more massive galaxies ($M_* \gtrsim 10^{11}$~\Msol). The enhancement of AGN is not uniform throughout the phases of a merger, and instead peaks within the early \emph{remnants} of merging systems (typically lagging $\approx 300$~Myr post-coalescence of the two galaxies at $z=0.5$). We argue that neither major ($M_{\mathrm{*,1}} / M_{\mathrm{*,2}} \geq \frac{1}{4}$) nor minor mergers ($\frac{1}{10} < M_{\mathrm{*,1}} / M_{\mathrm{*,2}} < \frac{1}{4}$) are statistically relevant for enhancing BH masses globally. Whilst at all redshifts the galaxies experiencing a merger have accretion rates that are on average 2--3 times that of isolated galaxies, the majority of mass that is accreted onto BHs occurs outside the periods of a merger. We compute that on average no more than 15\% of a BHs final day mass comes from the enhanced accretion rates triggered via a merger

preprint2020arXiv

Galaxy pairs in the Sloan Digital Sky Survey -- XIV. Galaxy mergers do not lie on the Fundamental Metallicity Relation

In recent observational studies, star-forming galaxies have been shown to follow a relation often dubbed the fundamental metallicity relation (FMR). This relation links the stellar mass of a galaxy with its star formation rate (SFR) and its gas-phase metallicity. Specifically, the FMR predicts that galaxies, at a given stellar mass, exhibit lower metallicities for higher SFRs. This trend is qualitatively consistent with observations of galaxy pairs, which have been robustly shown to experience increasing gas-phase metallicity dilution and enhanced star formation activity with decreasing projected separation. In this work, we show that, despite the qualitative consistency with FMR expectations, the observed O/H dilution in galaxy pairs of the Sloan Digital Sky Survey is stronger than what is predicted by the FMR. We conclude that the evolutionary phase of galaxies interacting with companions is not encoded in the FMR, and thus, mergers constitute a clearly defined population of outliers. We find that galaxies in pairs are consistent with the FMR only when their separation is larger than 110 kpc. Finally, we also quantify the local environment of the pairs using the number of galaxy neighbours within $2\, \mathrm{Mpc}$, $N_2$, and the projected separation to the second closest galaxy, $r_2$. We find that pairs are more sensitive to a second companion than to the local galaxy density, displaying less elevated SFRs with smaller values of $r_2$.

preprint2020arXiv

Interacting galaxies in the IllustrisTNG simulations -- I: Triggered star formation in a cosmological context

We use the IllustrisTNG cosmological hydrodynamical simulations to investigate how the specific star formation rates (sSFRs) of massive galaxies $(M_* > 10^{10} M_\odot)$ depend on the distance to their closest companions. We estimate sSFR enhancements by comparing with control samples that are matched in redshift, stellar mass, local density and isolation, and we restrict our analysis to pairs with stellar mass ratios of 0.1 to 10. At small separations (~15 kpc), the mean sSFR is enhanced by a factor of 2.0 $\pm$ 0.1 in the flagship (110.7 Mpc)$^3$ simulation (TNG100-1). Statistically significant enhancements extend out to 3D separations of 280 kpc in the (302.6 Mpc)$^3$ simulation (TNG300-1). We find similar trends in the EAGLE and Illustris simulations, although their sSFR enhancements are lower than those in TNG100-1 by about a factor of two. Enhancements in IllustrisTNG galaxies are seen throughout the redshift range explored ($0 \leq z < 1$), with the strength of the enhancements decreasing with increasing redshift for galaxies with close companions. In order to more closely compare with observational results, we separately consider 2D projected distances between galaxies in IllustrisTNG. We detect significant sSFR enhancements out to projected separations of 260 kpc in TNG300-1, with projection effects diluting the size of the enhancements by about 20 per cent below 50 kpc. We find similar sSFR enhancements in TNG100-1 and Sloan Digital Sky Survey galaxies, with enhancements extending out to projected separations of about 150 kpc for star-forming galaxies at $z < 0.2$. Finally, by summing over all separations, we estimate that the presence of closest companions boosts the average sSFR of massive galaxies in TNG100-1 by 14.5 per cent.

preprint2020arXiv

Interacting galaxies in the IllustrisTNG simulations -- II: Star formation in the post-merger stage

Galaxy mergers are a major evolutionary transformation whose effects are borne out by a plethora of observations and numerical simulations. However, most previous simulations have used idealised, isolated, binary mergers and there has not been significant progress on studying statistical samples of galaxy mergers in large cosmological simulations. We present a sample of 27,691 post-merger (PM) galaxies ($0\le z \le 1$) identified from IllustrisTNG: a cosmological, large box, magneto-hydrodynamical simulation suite. The PM sample spans a wide range of merger and galaxy properties ($M_\star$, $μ$, $f_\mathrm{gas}$). We demonstrate that star forming (SF) PMs exhibit enhanced star formation rates (SFRs) on average by a factor of $\sim 2$, while the passive PMs show no statistical enhancement. We find that the SFR enhancements: (1) show no dependence on redshift, (2) anti-correlate with the PM&#39;s stellar mass, and (3) correlate with the gas fraction of the PM&#39;s progenitors. However, SF PMs show stronger enhancements which may indicate other processes being at play (e.g., gas phase, feedback efficiency). Although the SFR enhancement correlates mildly with the merger mass ratio, the more abundant minor mergers ($0.1 \le μ< 0.3 $) still contribute $\sim 50\%$ of the total SFR enhancement. By tracing the PM sample forward in time, we find that galaxy mergers can drive significant SFR enhancements which decay over $\sim 0.5$ Gyr independent of the merger mass ratio, although the decay timescale is dependent on the simulation resolution. The strongest merger-driven starburst galaxies evolve to be passive/quenched on faster timescales than their controls.

preprint2020arXiv

The ALMaQUEST Survey: II. What drives central starbursts at z~0?

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.

preprint2020arXiv

Variations in the slope of the resolved star-forming main sequence: a tool for constraining the mass of star-forming regions

The correlation between galaxies&#39; integrated stellar masses and star formation rates (the `star formation main sequence&#39;; SFMS) is a well-established scaling relation. Recently, surveys have found a relationship between the star formation rate and stellar mass surface densities on kpc and sub-kpc scales (the `resolved SFMS&#39;; rSFMS). In this work, we demonstrate that the rSFMS emerges naturally in FIRE-2 zoom-in simulations of Milky Way-mass galaxies. We make SFR and stellar mass maps of the simulated galaxies at a variety of spatial resolutions and star formation averaging time-scales and fit the rSFMS using multiple methods from the literature. While the absolute value of the SFMS slope depends on the fitting method, the slope is steeper for longer star formation time-scales and lower spatial resolutions regardless of the fitting method employed. We present a toy model that quantitatively captures the dependence of the simulated galaxies&#39; rSFMS slope on spatial resolution and use it to illustrate how this dependence can be used to constrain the characteristic mass of star-forming clumps.

preprint2019arXiv

Are galactic star formation and quenching governed by local, global or environmental phenomena?

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 ($Σ_{\rm SFR}$) via dust corrected $Hα$ 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&#39; and `quenched&#39; 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 $Σ_{\rm SFR}$ values directly via regression in star forming regions. Local/ spatially resolved parameters are collectively the most predictive at estimating $Σ_{\rm SFR}$ in these analyses, with stellar mass surface density at the spaxel location ($Σ_*$) 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.

preprint2019arXiv

Probing Structure in Cold Gas at $z \lesssim 1$ with Gravitationally Lensed Quasar Sight Lines

Absorption spectroscopy of gravitationally lensed quasars (GLQs) enables study of spatial variations in the interstellar and/or circumgalactic medium of foreground galaxies. We report observations of 4 GLQs, each with two images separated by 0.8-3.0&#34;, that show strong absorbers at redshifts 0.4$<$$z_{abs}$$<$1.3 in their spectra, including some at the lens redshift with impact parameters 1.5-6.9 kpc. We measure H I Lyman lines along two sight lines each in five absorbers (10 sight lines in total) using HST STIS, and metal lines using Magellan Echellette or Sloan Digital Sky Survey. Our data have doubled the lens galaxy sample with measurements of H I column densities ($N_{\rm H I}$) and metal abundances along multiple sight lines. Our data, combined with the literature, show no strong correlation between absolute values of differences in $N_{\rm H I}$, $N_{\rm Fe II}$, or [Fe/H] and the sight line separations at the absorber redshifts for separations of 0-8 kpc. The estimated abundance gradients show a tentative anti-correlation with abundances at galaxy centers. Some lens galaxies show inverted gradients, possibly suggesting central dilution by mergers or infall of metal-poor gas. [Fe/H] measurements and masses estimated from GLQ astrometry suggest the lens galaxies lie below the total mass-metallicity relation for early-type galaxies as well as measurements for quasar-galaxy pairs and gravitationally lensed galaxies at comparable redshifts. This difference may arise in part from the dust depletion of Fe. Higher resolution measurements of H and metals (especially undepleted elements) for more GLQ absorbers and accurate lens redshifts are needed to confirm these trends.

preprint2019arXiv

The ALMaQUEST Survey: III. Scatter in the resolved star forming main sequence is primarily due to variations in star formation efficiency

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&#39;. 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.