Search for: All records

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

Abstract The metallicity of galaxies, and its variation with galactocentric radius, provides key insights into the formation histories of galaxies and the physical processes driving their evolution. In this work, we analyze the radial metallicity gradients of star-forming galaxies in the EAGLE, Illustris, IllustrisTNG, and SIMBA cosmological simulations across broad mass (10^8.0M_⊙≤M_⋆ ≲ 10^12.0M_⊙) and redshift (0 ≤z≤ 8) ranges. We find that all simulations predict strong negative (i.e., radially decreasing) metallicity gradients at early cosmic times, likely due to their similar treatments of relatively smooth stellar feedback not providing sufficient mixing to quickly flatten gradients. The strongest redshift evolution occurs in galaxies with stellar masses of 10^10.0–10^11.0M_⊙, while galaxies with stellar mass < 10¹⁰M_⊙and >10¹¹M_⊙exhibit weaker redshift evolution. Our result of negative gradients at high redshift contrast with the many positive and flat gradients in the 1 < z < 4 observational literature. Atz > 6, the negative gradients observed with JWST and the Atacama Large Millimeter/submillimeter Array are flatter than those in simulations, albeit with closer agreement than at lower redshift. Overall, we suggest that these smooth stellar feedback galaxy simulations may not sufficiently mix their metal content radially, and that either stronger stellar feedback or additional subgrid turbulent metal diffusion models may be required to better reproduce observed metallicity gradients. 

A brief correction to a publication published in 2020, which created a 0.03 dex offset in published calibrations. 

ABSTRACT The scatter about the mass-metallicity relation (MZR) has a correlation with the star formation rate (SFR) of galaxies. The lack of evidence of evolution in correlated scatter at z ≲ 2.5 leads many to refer to the relationship between mass, metallicity, and SFR as the Fundamental Metallicity Relation (FMR). Yet, recent high-redshift (z > 3) JWST observations have challenged the fundamental (i.e. redshift-invariant) nature of the FMR. In this work, we show that the cosmological simulations Illustris, IllustrisTNG, and Evolution and Assembly of GaLaxies and their Environment (EAGLE) all predict MZRs that exhibit scatter with a secondary dependence on SFR up to z = 8. We introduce the concept of a ‘strong’ FMR, where the strength of correlated scatter does not evolve with time, and a ‘weak’ FMR, where there is some time evolution. We find that each simulation analysed has a statistically significant weak FMR – there is non-negligible evolution in the strength of the correlation with SFR. Furthermore, we show that the scatter is reduced an additional ∼10–40 per cent at z ≳ 3 when using a weak FMR, compared to assuming a strong FMR. These results highlight the importance of avoiding coarse redshift binning when assessing the FMR. 

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

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

ABSTRACT Powerful outflows are thought to play a critical role in galaxy evolution and black hole growth. We present the first large-scale systematic study of ionized 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 ∼4000 paired galaxies and ∼250 post-mergers in the local universe (0.02 ≤ z ≤ 0.2) from the Sloan Digital Sky Survey Data Release 7 (SDSS DR 7) matched in stellar mass, redshift, local density of galaxies, and [O iii] λ5007 luminosity to a control sample of isolated galaxies. By fitting the [O iii] λ5007 line, we find ionized outflows in ∼15 per cent of our entire sample. Outflows are much rarer in star-forming galaxies compared to active galactic nuclei (AGNs), and outflow incidence and velocity increase with [O iii] λ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]~\lambda 5007}}\, \gtrsim$$ 1042 erg s−1), suggesting mid-infrared AGN selection favours galaxies with powerful outflows, at least for higher [O iii] λ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 ionized outflows as traced by [O iii] λ5007. 

Abstract Galaxy mergers are predicted to trigger accretion onto the central supermassive black holes, with the highest rates occurring during final coalescence. Previously, we have shown elevated rates of both optical and mid-IR selected active galactic nuclei (AGN) in post-mergers, but to date the prevalence of X-ray AGN has not been examined in the same systematic way. We present XMM-Newton data of 43 post-merger galaxies selected from the Sloan Digital Sky Survey along with 430 non-interacting control galaxies matched in stellar mass, redshift, and environment in order to test for an excess of hard X-ray (2–10 keV) emission in post-mergers attributable to triggered AGN. We find 2 X-ray detections in the post-mergers ($$4.7^{+9.3}_{-3.8}\%$$) and 9 in the controls ($$2.1^{+1.5}_{-1.0}\%$$), an excess of $$2.22^{+4.44}_{-2.22}$$, where the confidence intervals are 90%. While we therefore do not find statistically significant evidence for an X-ray AGN excess in post-mergers (p = 0.26), we find a factor of ∼17 excess of mid-IR AGN in our sample, consistent with past work and inconsistent with the observed X-ray excess (p = 2.7 × 10−4). Dominant, luminous AGN are therefore more frequent in post-mergers, and the lack of a comparable excess of 2–10 keV X-ray AGN suggests that AGN in post-mergers are more likely to be heavily obscured. Our results are consistent with the post-merger stage being characterised by enhanced AGN fueling, heavy AGN obscuration, and more intrinsically luminous AGN, in line with theoretical predictions. 

ABSTRACT We present the measured gas-phase metal column densities in 155 sub-damped Ly α 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 ≤ zabs ≤ 4.3. Using all available column densities of the ionic species investigated (mainly C iv, Si ii, Mg ii, Si iv, Al ii, Fe ii, C ii, and O i; 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 H i-weighted gas-phase metallicity evolution of subDLAs is consistent with damped Ly α systems (DLAs). When ionization corrections are included, subDLAs are systematically more metal poor than DLAs (between ≈0.5σ and ≈3σ significance) by up to ≈1.0 dex over the redshift range 3 ≤ zabs ≤ 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. 

ABSTRACT We investigate the spatial structure and evolution of star formation and the interstellar medium (ISM) in interacting galaxies. We use an extensive suite of parsec-scale galaxy-merger simulations (stellar mass ratio = 2.5:1), which employs the ‘Feedback In Realistic Environments-2’ model (fire-2). This framework resolves star formation, feedback processes, and the multiphase structure of the ISM. We focus on the galaxy-pair stages of interaction. We find that close encounters substantially augment cool (H i) and cold-dense (H2) gas budgets, elevating the formation of new stars as a result. This enhancement is centrally concentrated for the secondary galaxy, and more radially extended for the primary. This behaviour is weakly dependent on orbital geometry. We also find that galaxies with elevated global star formation rate (SFR) experience intense nuclear SFR enhancement, driven by high levels of either star formation efficiency (SFE) or available cold-dense gas fuel. Galaxies with suppressed global SFR also contain a nuclear cold-dense gas reservoir, but low SFE levels diminish SFR in the central region. Concretely, in the majority of cases, SFR enhancement in the central kiloparsec is fuel-driven (55 per cent for the secondary, 71 per cent for the primary) – while central SFR suppression is efficiency-driven (91 per cent for the secondary, 97 per cent for the primary). Our numerical predictions underscore the need of substantially larger, and/or merger-dedicated, spatially resolved galaxy surveys – capable of examining vast and diverse samples of interacting systems – coupled with multiwavelength campaigns aimed to capture their internal ISM structure.