Metal line emission around z<1 galaxies

Rajeshwari Dutta, Michele Fumagalli, Matteo Fossati, Marc Rafelski, Mitchell Revalski, Fabrizio Arrigoni Battaia, Valentina d'Odorico, Celine Peroux, Laura J Prichard, A.M Swinbank

► To cite this version:

Rajeshwari Dutta, Michele Fumagalli, Matteo Fossati, Marc Rafelski, Mitchell Revalski, et al. Metal line emission around z<1 galaxies. Astronomy & Astrophysics - A&A, 2024, 691, pp.A236. 10.1051/0004-6361/202450733. hal-04692567

HAL Id: hal-04692567

https://hal.science/hal-04692567v1

Submitted on 16 Jan 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Metal line emission around $z < 1$ galaxies

Rajeshwari Dutta<fnref n="1" />,*, Michele Fumagalli<fnref n="2" />,<fnref n="3" />, Matteo Fossati<fnref n="2" />,<fnref n="4" />, Marc Rafelski<fnref n="5" />,<fnref n="6" />, Mitchell Revalski<fnref n="5" />, Fabrizio Arrigoni Battaia<fnref n="7" />, Valentina D'Odorico<fnref n="3" />,<fnref n="8" />,<fnref n="9" />, Celine Péroux<fnref n="10" />,<fnref n="11" />, Laura J. Prichard<fnref n="5" />, and A. Mark Swinbank<fnref n="12" />

Received 15 May 2024 / Accepted 4 August 2024

## ABSTRACT

We characterize, for the first time, the average extended emission in multiple lines ([O II], [O III], and Hβ) around a statistical sample of 560 galaxies at $z \approx 0.25-0.85$. By stacking the Multi Unit Spectroscopic Explorer (MUSE) 3D data from two large surveys, the MUSE Analysis of Gas around Galaxies (MAGG) and the MUSE Ultra Deep Field (MUDF), we detect significant [O II] emission out to $\approx 40$ kpc, while [O III] and Hβ emission is detected out to $\approx 30$ kpc. Via comparisons with the nearby average stellar continuum emission, we find that the line emission at 20–30 kpc likely arises from the disk-halo interface. Combining our results with that of our previous study at $z \approx 1$, we find that the average [O II] surface brightness increases independently with redshift over $z \approx 0.4-1.3$ and with stellar mass over $M_* \approx 10^{6-12} M_\odot$, which is likely driven by the star formation rate as well as the physical conditions of the gas. By comparing the observed line fluxes with photoionization models, we find that the ionization parameter declines with distance, going from $\log q$ (cm s⁻<fnref n="1" />) $\approx 7.7$ at $\le 5$ kpc to $\approx 7.3$ at 20–30 kpc, which reflects a weaker radiation field in the outer regions of galaxies. The gas-phase metallicity shows no significant variation over 30 kpc, with a metallicity gradient of $\approx 0.003$ dex kpc⁻<fnref n="1" />, which indicates an efficient mixing of metals on these scales. Alternatively, there could be a significant contribution from shocks and diffuse ionized gas to the line emission in the outer regions.

**Key words.** galaxies: evolution – galaxies: halos – galaxies: ISM

## 1. Introduction

The majority of the baryons in the Universe are found in the form of gas in the intergalactic medium (IGM) and the circumgalactic medium (CGM) of galaxies. The gas acts as the fuel for the formation of stars in galaxies. The formation of stars leads to the creation of heavy elements or metals. Stellar winds and supernova explosions eject the metals out of stars into the interstellar medium (ISM) and even beyond into the CGM and the IGM. Galaxies thus form and evolve in complex ecosystems that are regulated by a cycle of gas flows (Somerville & Davé 2015; Tumlinson et al. 2017; Péroux & Howk 2020). This baryon cycle is believed to consist of: metal-poor gas accretion onto galaxies from filaments of the IGM; stellar feedback and active galactic nucleus-driven outflows of metal-enriched gas into the CGM; recycling of gas from the CGM back into the ISM through fountains; and transfer of gas between galaxies through interactions and mergers (Oppenheimer & Davé 2008; Dekel et al. 2009; Cresci et al. 2010; Kewley et al. 2010; Lilly et al. 2013; Fraternali 2017). All these processes leave their imprint on the distribution and

physical conditions of gas within and around galaxies, with the relative importance of each process evolving with cosmic time.

Strong rest-frame optical emission lines have long been used to trace the physical conditions of the ionized gas in the ISM. Several diagnostic line ratio diagrams have been developed to probe the chemical and ionization state of the ISM across redshifts (e.g., Kauffmann et al. 2003; Kewley et al. 2019; Maiolino & Mannucci 2019). Higher-redshift galaxies are found to have higher ionization parameters, harder ionizing spectra, and higher electron densities (e.g., Nakajima et al. 2013; Shirazi et al. 2014; Sanders et al. 2021). Furthermore, the metallicity of the gas in the ISM is found to decrease with increasing redshift at a fixed stellar mass (e.g., Erb et al. 2006; Maiolino et al. 2008; Sanders et al. 2021).

In the nearby Universe, the gas-phase metallicity in galaxies typically exhibits a negative gradient (i.e., the central regions are more metal-rich than the outskirts), which supports an inside-out picture of galaxy evolution (e.g., Zaritsky et al. 1994; Sánchez et al. 2014; Belfiore et al. 2017). Integral field unit (IFU) spectrographs, such as the Multi Unit Spectroscopic Explorer (MUSE; Bacon et al. 2010) and the K-band Multi

<footnote n="1">¹ IUCAA, Postbag 4 Ganeshkind, Pune 411007, India</footnote>

<footnote n="2">² Dipartimento di Fisica G. Occhialini, Università degli Studi di Milano Bicocca, Piazza della Scienza 3, 20126 Milano, Italy</footnote>

<footnote n="3">³ INAF – Osservatorio Astronomico di Trieste, Via G. B. Tiepolo 11, I-34143 Trieste, Italy</footnote>

<footnote n="4">⁴ INAF – Osservatorio Astronomico di Brera, Via Bianchi 46, 23087 Merate, (LC), Italy</footnote>

<footnote n="5">⁵ Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA</footnote>

<footnote n="6">⁶ Department of Physics and Astronomy, Johns Hopkins University, Baltimore, MD 21218, USA</footnote>

<footnote n="7">⁷ Max-Planck-Institut für Astrophysik, Karl-Schwarzschild-Str. 1, D-85748 Garching bei München, Germany</footnote>

<footnote n="8">⁸ Scuola Normale Superiore, P.zza dei Cavalieri, I-56126 Pisa, Italy</footnote>

<footnote n="9">⁹ IFPU-Institute for Fundamental Physics of the Universe, via Beirut 2, I-34151 Trieste, Italy</footnote>

<footnote n="10">¹⁰ European Southern Observatory, Karl-Schwarzschildstrasse 2, D-85748 Garching bei München, Germany</footnote>

<footnote n="11">¹¹ Aix Marseille Université, CNRS, LAM (Laboratoire d'Astrophysique de Marseille), UMR 7326, F-13388 Marseille, France</footnote>

<footnote n="12">¹² Centre for Extragalactic Astronomy, Department of Physics Durham University, South Road, Durham DH1 3LE, UK</footnote>

<footnote>* Corresponding author; rajeshwari.dutta@iucaa.in</footnote>
Object Spectrograph (K-MOS; Sharples et al. 2013) on the Very Large Telescope (VLT), and slitless grism spectrographs on the Hubble Space Telescope (HST) have enabled spatially resolved studies of the ISM up to $z \approx 3$. Such studies have revealed flatter metallicity gradients in $z \gtrsim 0.5$ galaxies with a large scatter and some galaxies even having positive gradients (e.g., Wuyts et al. 2016; Carton et al. 2018; Curti et al. 2020; Wang et al. 2020; Simons et al. 2021). Recently, James Webb Space Telescope (JWST) near-infrared (NIR) IFU observations found flat metallicity gradients in $z \approx 6-8$ galaxies (Venturi et al. 2024). The flattening of metallicity gradients could be due to several factors, including the inflow of metal-poor gas into the inner regions, outflows populating the outer regions with metals, recycling of metal-rich gas in the outer regions, and transfer of metals due to galaxy interactions and mergers (e.g., Cresci et al. 2010; Kewley et al. 2010; Rupke et al. 2010; Bresolin et al. 2012).

The above-cited studies probed the ionized gas within the stellar disks of galaxies (a few kiloparsecs, or a few effective radii). Beyond the stellar disks, galaxies are surrounded by halos of diffuse gas that extend out to a few hundred kiloparsecs (or a few virial radii). This is known as the CGM (Tumlinson et al. 2017), and it is relatively challenging to probe in emission due to the low gas density but is effectively probed in absorption against a bright background source such as a quasar. Along with the stellar and ISM properties of galaxies, the distribution and physical conditions of gas in the CGM can place crucial constraints on models of galaxy formation and evolution (e.g., Fumagalli et al. 2016; Liang et al. 2016; Turner et al. 2017; Oppenheimer et al. 2018; Weng et al. 2024). Based on several large observational campaigns carried out to characterize the CGM, we know that the CGM at $z \lesssim 4$ consists of metal-enriched, multiphase gas, which becomes more ionized with increasing distance from galaxies and which is influenced by both galaxy properties and environmental effects (e.g., Werk et al. 2014; Hamanowicz et al. 2020; Dutta et al. 2020, 2021; Galbiati et al. 2023; Qu et al. 2023).

Detecting emission in multiple lines from the diffuse gas in the extended stellar disk, disk-halo interface, and outer halo can provide us with additional and complementary insights into the distribution and physical conditions of the CGM. Sensitive optical IFU spectrographs such as MUSE and the Keck Cosmic Web Imager (KCWI; Morrissey et al. 2018) have facilitated detections of the low surface brightness (SB) halo gas at cosmological distances in emission. There have been a few detections of metal line emission extending up to tens of kiloparsecs around both individual galaxies and stacks of galaxies (e.g., Burchett et al. 2021; Zabl et al. 2021; Leclercq et al. 2022, 2024; Dutta et al. 2023; Guo et al. 2023).

Recently, by stacking the MUSE data of a sample of ≈600 galaxies, Dutta et al. (2023) detected, for the first time, the average Mg II and [O II] line emission up to ≈30–40 kpc from a general population of galaxy halos at $z \approx 1$. The shallower radial SB profile of Mg II up to ≈20 kpc compared to [O II] suggests that the resonant Mg II emission is affected by dust and radiative transfer effects, while the constant ratio of [O II] to Mg II over ≈20–40 kpc suggests a non-negligible in situ origin of the extended metal emission. At $z < 1$, the detection of extended emission in multiple lines (e.g., Hα, Hβ, [O II], and [O III]) around a few galaxies and quasars has enabled studies of the physical conditions and excitation mechanisms in the gas beyond the stellar disks (e.g., Johnson et al. 2018; Chen et al. 2019; Nielsen et al. 2023). Such studies help bridge our understanding of the gas within the stellar disks and the gas beyond it in the halo.

Here we extend the stacking analysis of Dutta et al. (2023) to $z < 1$ to investigate the average physical conditions in gas around a statistical sample of galaxies. The observations and the stacking procedure used in this work are explained in Sect. 2. The results of the stacked line emission and line ratios around galaxies are presented in Sect. 3. The results are summarized and discussed in Sect. 4. Throughout this work, we use a Planck 2015 cosmology with $H_0 = 67.7 \text{ km s}^{-1} \text{ Mpc}^{-1}$ and $\Omega_M = 0.307$ (Planck Collaboration XIII 2016).

## 2. Observations and analysis

This work is based upon MUSE observations from two large programs on the VLT: (1) MUSE Analysis of Gas around Galaxies (MAGG; Lofthouse et al. 2023), which consists of medium-deep, single pointing ($1 \times 1 \text{ arcmin}^2$) MUSE observations of 28 fields that are centered on quasars at $z \approx 3.2-4.5$; (2) MUSE Ultra Deep Field (MUDF; Fossati et al. 2019), which consists of very deep MUSE observations of a $1.5 \times 1.2 \text{ arcmin}^2$ region around a pair of quasars at $z \approx 3.2$. The MUSE exposure time in the MAGG survey is ≈4 h per field, except for two fields that have deeper exposure time of ≈10 h. The MUSE observations in the MUDF survey amount to a total of ≈143 h, with the exposure time decreasing from the center of the field to the outer regions. The average image quality of the MUSE data is ≈0.7 arcsec full-width at half-maximum (FWHM), which corresponds to $\le 5$ kpc at $z < 1$.

The observations and galaxy catalogs from these surveys are described in detail in Sect. 2 of Dutta et al. (2023). We selected all the galaxies in the MAGG and MUDF catalogs that lie within the redshift range $z = 0.25-0.85$. This selection ensured that we got simultaneous coverage of the [O II]λλ3727,3729, Hβ λ4863, [O III]λ4960 and λ5008 emission lines in the MUSE spectra. The total sample consists of 560 galaxies, 480 from the MAGG survey, and 80 from the MUDF survey. The total MUSE exposure time of the galaxies used for stacking in this work is ≈6264 h.

To obtain the physical properties of the galaxies such as stellar masses and star formation rates (SFRs), the MUSE spectra and photometry were jointly fit with stellar population synthesis models using the Monte Carlo Spectro-Photometric Fitter (MC-SPF; Fossati et al. 2018). For the MUDF galaxy sample, HST optical and NIR photometry in five bands (Revalski et al. 2023) and HAWK-I K-band photometry were used in addition to the MUSE spectra and photometry. We refer to Sect. 2 of Dutta et al. (2020) for details of the spectral energy distribution (SED) fitting process. Briefly, MC-SPF adopts the Bruzual & Charlot (2003) models at solar metallicity, the Chabrier (2003) initial mass function, nebular emission lines from the models of Byler et al. (2018), and the dust attenuation law of Calzetti et al. (2000).

The distributions of redshifts and stellar masses of the galaxy sample used for stacking in this work are shown in Fig. 1. The median redshift of the sample is ≈0.6. The typical uncertainty in the redshift estimates from MUSE spectra is ≈60 km s⁻¹. The galaxy sample spans the stellar mass range of ≈10⁶⁻¹¹ M☉, with a median stellar mass of ≈10⁹ M☉. The typical uncertainty in the estimates of stellar mass from SED fitting is ≈0.1–0.2 dex.

We followed the same procedure to stack the MUSE data of the galaxies as described in Sect. 3 of Dutta et al. (2023). In brief, for each galaxy we first extracted a sub-cube of area $200 \times 200 \text{ kpc}^2$ centered on the galaxy spanning 100 Å in rest wavelength around the emission lines of interest, after subtracting the continuum emission and masking out all the other continuum sources. We performed both mean and median stacking of all the sub-cubes in the rest-frame wavelengths of the galaxies.
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Sellar mass versus redshift of the sample of galaxies used for stacking. The galaxies from the MAGG survey are marked as light purple circles, and the galaxies from the MUDF survey are marked as dark purple squares. The size of the symbols is proportional to the MUSE exposure time of the galaxies. The top and right panels show the histograms of redshifts and stellar masses, respectively, in light purple for the MAGG galaxies, and in dark purple for the MUDF galaxies. The dashed lines mark the median values of the MAGG galaxies, and the dotted lines mark the median values of the MUDF galaxies.</figcaption></figure>

The results from both the stacks are found to be consistent within $1\sigma$. Here we present the results based on median stacking, which is more robust against outliers. To subtract the continuum, we fit a low-order spline to the continuum emission around the lines for each spaxel in the cube. We repeated the above procedure 100 times with repetition to obtain the sample variance from the 16th and 84th percentiles of the bootstrapped sample.

We corrected for the Milky Way extinction following Schlafly & Finkbeiner (2011) and the Fitzpatrick (1999) extinction curve. We applied a correction for dust to the fluxes of each galaxy before stacking using the dust extinction derived by MC-SPF and the extinction curve of Calzetti et al. (2000). We checked that there is no significant difference in the results if instead we correct for dust after stacking using the median dust extinction of the sample. We note that there likely are radial gradients in the dust attenuation within a galaxy (e.g., Wuyts et al. 2012; Nelson et al. 2016; Tacchella et al. 2018; Wilman et al. 2020), which we do not take into account in this work because spatially resolved dust attenuation distribution in high-z galaxies is still not well understood, particularly at the larger spatial scales probed in this work. For the median stellar mass of this sample, the dust attenuation is found to be weak and flat on average in the ISM of $z \approx 1$ galaxies (Nelson et al. 2016). To correct the $H\beta$ flux for the underlying stellar absorption, we increased the flux by 3% based on the typical average correction found in the literature (e.g., Cullen et al. 2021; Sanders et al. 2021; Curti et al. 2023; Stephenson et al. 2024).

To generate pseudo-narrowband (NB) images of the line emission, we summed the stacked cube from $-250$ to $250$ km s$^{-1}$ in rest-frame velocity around the $H\beta$ and [O III] lines, and from $-250$ to $500$ km s$^{-1}$ around the $\lambda3727$ line of the [O II] doublet to cover both the lines of the doublet. This velocity range is found to encompass all the emission in the stacked spectra (see Fig. 2). For reference, the median virial velocity expected for

the host halos of the galaxies in this sample is $\approx 100$ km s$^{-1}$. To make a comparison with the line NB images, we generated NB images of the stellar continuum emission by averaging the NB images over two velocity windows around $\pm 1500$ km s$^{-1}$ of each line with the same velocity width as for the line NB. The line and continuum NB images were used to estimate the average SB radial profiles. We did not take the inclination and orientation of the galaxies into account during stacking because we did not have HST imaging available for the full sample. Therefore, the stacked emission at any projected separation represents the azimuthally averaged emission of the galaxy sample at that separation.

# 3. Results

## 3.1. Line emission around galaxies

First we characterized the average emission around $z \approx 0.6$ galaxies from the [O II], $H\beta$, and [O III] lines. We extracted the 1D spectra from the 3D median stacked cubes by summing up the flux over different annular regions. Figure 2 shows the spectra extracted within 0–5 kpc, 5–10 kpc, 10–20 kpc, 20–30 kpc, and 30–40 kpc around the [O II] $\lambda\lambda3727, 3729$ doublet, $H\beta$, and [O III] $\lambda5008$ lines. We note that the [O III] $\lambda4960$ line, not shown here, is found to have 1/3 the flux of the [O III] $\lambda5008$ line as per theoretical expectations. The [O II] line emission is detected out to $\approx 40$ kpc at $3\sigma$, while the $H\beta$ and [O III] emission lines are detected out to $\approx 30$ kpc at $3\sigma$. It can be seen that, as expected, the flux of the emission lines decreases going outward from the center. We fit Gaussian profiles to the four emission lines with the constraint that the velocity centroid and velocity width of all the lines are equal. We fit a double Gaussian profile to the [O II] doublet, allowing the ratio of the line fluxes, $\lambda3729/\lambda3727$, to vary between the low density limit of 0.35 and the high density limit of 1.5. We find that the line fluxes decrease by a factor of $\approx 10–20$ going from the central 5 kpc region to the outer 20–30 kpc annular region. In addition to the decrease in flux, the emission lines become broader in the outer regions, with the velocity dispersion (not corrected for line spread function) going from $\approx 75$ km s$^{-1}$ within 5 kpc to $\approx 96$ km s$^{-1}$ at 20–30 kpc. The [O II] doublet ratio increases from $\approx 1.3$ within 5 kpc to the low density limit of 1.5 at 30–40 kpc.

The top row of Fig. 3 shows the NB images of the line emission in the median stacked cube. The NB images of the continuum emission around the lines are shown in the second row for comparison. The combined NB image of the [O III] $\lambda4960$ and $\lambda5008$ lines is shown. The [O II] line emission is the most spatially extended, followed by the [O III] and $H\beta$ emission. To quantify the average radial profile and extent of the emission, the azimuthally averaged SB in annuli of radius 10 kpc, as obtained from the NB images, are shown in Fig. 4. The [O II] SB profile is more radially extended than that of the continuum, while the $H\beta$ and [O III] SB profiles are similar to that of the continuum. The continuum emission near the [O II] line extends out to $\approx 20$ kpc at a $3\sigma$ SB level of $\approx 10^{-20}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$. This suggests that we are probing the stellar disk or ISM up to $\approx 20$ kpc. Beyond that at 20–30 kpc, the [O II] emission is significantly enhanced compared to the continuum, while $H\beta$ and [O III] are marginally enhanced. The line emission in this region is most likely originating from the disk-halo interface, where the ISM is transitioning into the CGM. In the case of [O II], we detect significant emission out to 30–40 kpc, which likely originates from the CGM.
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Spectra of different emission lines (labeled at the top of each column) extracted from the median stacked cubes of the sample. From top to bottom, the spectra represent the total flux in annuli in radius ranges of 0–5 kpc, 5–10 kpc, 10–20 kpc, 20–30 kpc, and 30–40 kpc. The 16th and 84th percentiles of the spectra from bootstrapping analysis are shown as vertical bars. The Gaussian fits to the lines are shown as dotted red lines. The velocity windows used to construct the pseudo-NB images and SB profiles are marked with vertical dashed lines.</figcaption></figure>

## 3.2. Evolution of [O II] emission

Based on a similar stacking analysis of higher-redshift galaxies ($z \approx 0.7–1.5$) in the MAGG and MUDF surveys, Dutta et al. (2023) found that the [O II] line emission becomes brighter and more spatially extended with increasing stellar mass and redshift, and suggested that this was likely due to higher SFRs at higher stellar masses and redshifts. To further investigate the evolution of [O II] emission with redshift, we combined our sample with that presented in Dutta et al. (2023), and conducted two control experiments. First, we formed four subsamples in the redshift ranges, $0.25 \le z < 0.5$, $0.5 \le z < 0.7$, $0.7 \le z < 1.0$, and $1.0 \le z < 1.5$, which are matched in stellar mass within a factor of two, such that the maximum difference between the cumulative distributions of the stellar mass in each of the redshift bins is $\le 0.1$ and the $p$-value is $\ge 0.8$ based on two-sided Kolmogorov-Smirnov test. In this way, we were able to form a stellar mass-matched sample of 70 galaxies in each of the four redshift bins. Next, similar to above, we formed subsamples in the four redshift bins that are matched within a factor of two in both stellar mass and SFR. This led to a matched sample of 34 galaxies in each redshift bin.

In the left panel of Fig. 5 we plot the average [O II] SB, corrected for cosmological dimming, as a function of redshift, for both the matched samples. The average values are estimated in two regions: within a circular aperture of radius 20 kpc that that probes the ISM and within an annular aperture between radii

20 and 30 kpc that likely probes the disk-halo interface. For the inner region, the average SB is found to increase by a factor of $\approx 20$ from $z \approx 0.4$ to $z \approx 1.3$ in the case of the stellar mass-matched samples, while in the case of the stellar mass- and SFR-matched samples, the dependence of SB on redshift is weaker, with the average SB increasing by a factor of $\approx 4$ from $z \approx 0.4$ to $z \approx 1.3$. For the outer region, the trend with redshift again becomes weaker for the stellar mass- and SFR-matched samples, although there is no significant difference between the two matched samples. This indicates that the higher SFR of galaxies at higher redshifts is likely a major factor behind the enhanced [O II] emission in the ISM at high redshifts. However, there are likely additional factors, such as gas density, that contribute to the enhanced extended [O II] emission at high redshifts. Indeed, studies have found that the density of the warm ionized gas is higher at $z \approx 2–3$ compared to local galaxies that are matched in stellar mass and star formation activities (e.g., Shirazi et al. 2014; Davies et al. 2021).

Furthermore, we performed similar control experiments as above to investigate the dependence of [O II] emission on stellar mass using the combined sample of this work and Dutta et al. (2023). Firstly, we formed four subsamples in the stellar mass bins, $M_* = 10^{6–8}$, $10^{8–9}$, $10^{9–10}$, and $10^{10–12} M_{\odot}$, which are matched in redshift within 0.3, such that the $p$-value from Kolmogorov-Smirnov test is $\ge 0.8$. This led to a redshift matched sample of 51 galaxies in each of the stellar mass bins. Secondly, we formed four subsamples in the above stellar mass bins that
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Pseudo-NB images of the median stacked sample smoothed using a Gaussian kernel of standard deviation 0.3 arcsec. The top row shows, from left to right, the NB images of the [O II] ($\lambda$3727 + $\lambda$3729), H$\beta$, and [O III] ($\lambda$4960 + $\lambda$5008) line emission. The bottom row shows, from left to right, the NB images of the continuum emission near the [O II], H$\beta$, and [O III] lines. The red contour is plotted at three times the rms noise in the NB images.</figcaption></figure>

are matched within a factor of two in SFR in addition to being matched in redshift as above, leading to samples of 25 galaxies in each mass bin. We then performed the [O II] line emission stacking for these subsamples. The average [O II] SB values in the inner 20 kpc and outer 20–30 kpc regions are shown in the right panel of Fig. 5 for the matched samples. In the inner region, the [O II] SB increases by a factor of ≈8 from $M_*= 10^{6-8} M_\odot$ to $10^{10-12} M_\odot$ for the redshift matched samples, and by a factor of ≈2 when matched additionally in SFR. In the outer region, the increasing trend with stellar mass is consistent within the large uncertainties for both the matched samples. Therefore, the dependence of the [O II] emission on stellar mass is likely driven by a combination of dependence on star formation activity and other physical conditions in the gas. We note that in the above analysis, as mentioned in Sect. 2, we corrected the fluxes for dust before stacking using the stellar extinction obtained from SED fitting, assuming that the nebular and stellar extinction follow each other.

### 3.3. Line ratios around galaxies

Having detected emission from multiple lines in the stacked cube, we next investigated the average physical conditions in the ionized gas around $z \approx 0.6$ galaxies using different line ratios. In particular, the O32 line ratio, defined as $[\text{O III}]\lambda5008/[\text{O II}]\lambda3727, 3729$, is an indicator of the degree of ionization of the gas, while the R23 line ratio, defined as

([O III]$\lambda$4960, 5008+[O II]$\lambda$3727, 3939)/H$\beta$, is sensitive to the gas-phase metallicity. The O32 versus R23 diagram has been used to study the ionization and metallicity in the ISM using integrated spectra of galaxies. Star-forming galaxies generally follow a trend in the O32–R23 diagram, going from low R23 and O32 values at high metallicity to high R23 and O32 values at low metallicity, although R23 begins to decrease at very low metallicity (≤0.1 Z$_\odot$; Kewley et al. 2019). We note that according to the mass-metallicity relation, the median stellar mass of the sample is approximately at the turn-over metallicity of R23.

Thanks to stacking of the MUSE 3D data, here we were able to spatially resolve the above line ratios and study how they vary on average around a statistical sample of star-forming galaxies. We note that physical conditions such as metallicity can vary with stellar mass and redshift (Tremonti et al. 2004; Maiolino & Mannucci 2019). Therefore, we created similar stacks in bins of stellar mass and redshift below and above the median values, and also by weighting with the stellar mass, to confirm that the trends with physical separation we find below remain similar over the stellar mass and redshift range probed here.

To compute the line ratios, we integrated the flux in the spectra extracted from the stacked cube over the velocity window marked in Fig. 2. For the [O II] doublet, we considered the total flux of both the lines. The error in the integrated flux was obtained by propagating the error in the spectra. If the integrated flux was detected at less than 2σ, we considered it an upper limit.
<figure><img src="image_5.png" /></figure>

**Fig. 4.** SB of the median stacked sample as a function of projected separation from the galaxy center. The panels from left to right show the SB profiles of the [O II] (λ3727 + λ3729), Hβ, and [O III] (λ4960 + λ5008) emission. The SB is estimated as the azimuthally averaged value in annuli of radius 10 kpc from the pseudo-NB images shown in Fig. 3. In each panel, the solid blue line shows the SB profile of the line emission, the dotted black line shows the SB profile of the continuum emission near the line (averaged over two velocity windows with the same width as that of the line), the dashed light blue line shows the 3σ upper limit on the line SB, and the dashed gray line shows the 3σ upper limit on the continuum SB. The 16th and 84th percentiles of the SB profiles from bootstrapping analysis are shown as vertical bars. Note that the average image quality of the MUSE data is ≤5 kpc FWHM at z ≤ 1.

<figure><img src="image_6.png" /></figure>

**Fig. 5.** Evolution of [O II] SB with redshift and stellar mass. Left: Average SB of [O II] emission, corrected for cosmological dimming, in four different redshift bins. The purple circles are for samples that are matched in stellar mass, whereas the orange squares are for samples that are matched in both stellar mass and SFR. Right: Average SB of [O II] emission in four different stellar mass bins. The purple circles are for samples that are matched in redshift, whereas the orange squares are for samples that are matched in both redshift and SFR. In both panels, the filled symbols represent the average values within a circular aperture of radius 20 kpc, and the open symbols the average values within an annular region between radii 20 and 30 kpc.

The left panel of Fig. 6 shows the O32 versus R23 line ratios estimated from the stacked spectra in different annuli (Fig. 2). The O32 ratio is found to decrease going from the central 0–5 kpc bin (log O32 ≈ −0.1) to the outer 20–30 kpc bin (log O32 ≈ −0.5). The R23 ratio is found to increase from the central 0–5 kpc bin (log R23 ≈ 0.8) to the 10–20 kpc bin (log R23 ≈ 0.9), and then decrease in the 20–30 kpc bin (log R23 ≈ 0.75).

For comparison, we also plot the O32 versus R23 ratios from integrated spectra of lower- and higher-redshift galaxies. The median value of the ratios in the sample of z ≈ 0.1 galaxies from the Sloan Digital Sky Survey (SDSS; York et al. 2000) are obtained using the Max-Planck-Institute for Astrophysics – Johns Hopkins University (MPA–JHU) catalog<fnref n="1" />. The line ratios at z ≈ 2.3 and z ≈ 3.3 are obtained from the composite spectra of star-forming galaxies in different stellar mass bins

(≈2 × 10⁹−4 × 10¹⁰ M☉) in the MOSFIRE Deep Evolution Field (MOSDEF) survey (Sanders et al. 2021), with the line ratios decreasing with increasing stellar mass. The stacked line ratios in the central region, within 20 kpc, are similar to the average line ratios found in the ISM of z ≈ 2–3 galaxies with stellar mass, log M* ≈ 10.2–10.6 M☉. Furthermore, the stacked line ratios in the central region are consistent with the trend of R23 increasing at a fixed O32, and O32 decreasing at a fixed R23, from z ≈ 0 to z ≈ 3. On the other hand, the stacked line ratios at 20–30 kpc are more similar to those found in the ISM of local SDSS galaxies. This indicates a transition in the physical conditions or the ionization mechanism of the warm ionized gas from the inner to the outer regions (i.e., from the ISM to the disk-halo interface and the CGM).

To investigate this further, we compared our results with photoionization models. Strong line ratios – calibrated either theoretically via comparison with photoionization models or empir-

<footnote n="1">¹ https://www.mpa.garching.mpg.de/SDSS/DR7/</footnote>
<figure><img src="image_7.png" /><figcaption>**Fig. 6.** Stacked line ratios around galaxies. Left: O32 versus R23 diagram. The line ratios estimated from the stacked spectra (Fig. 2) are shown as circles that are color-coded by the average projected separation from the galaxy center. The annular bins (in kpc) corresponding to the ratios are also labeled next to each circle. Arrows denote $2\sigma$ limits on the line ratios. The median value of the line ratios of $z \approx 0.1$ galaxies from SDSS is shown as the dashed black line. Line ratios from composite spectra of galaxies in the MOSDEF survey (Sanders et al. 2021) are shown as downward orange triangles for the $z \approx 2.3$ sample and as purple upward triangles for the $z \approx 3.3$ sample. Center: Ionization parameter versus metallicity (symbols are same as in the left panel) estimated from the emission line fluxes in the stacked spectra and photoionization models using the code IZI. The solar metallicity is demarcated by the dotted line. Right: O23 versus R23 diagram of the stacked sample (symbols are same as in the left panel) compared with results from the “shock only” (solid lines) and “shock+precursor” models (dotted lines) from Allen et al. (2008). Shock models are shown for several magnetic parameters and shock velocities, as denoted in the plot.</figcaption></figure>

ically via comparison with estimates from direct method in sources when possible – have been used extensively to infer the ionization parameter and chemical abundance of the photoionized gas in the ISM (Maiolino & Mannucci 2019; Kewley et al. 2019). However, the O32 ratio, typically used as an estimator of the ionization parameter, is influenced by the gas-phase metallicity and gas pressure. Similarly, the R23 ratio, which is commonly used to determine chemical abundances, is sensitive to the ionization parameter and gas pressure, and is moreover double-valued, having both a low and a high abundance branch. On the other hand, Bayesian methods that use theoretical photoionization models to fit emission line fluxes can be used to simultaneously derive the ionization parameter ($q$) and gas-phase metallicity ($Z$), although they likely still suffer from some of the same limitations as the individual line ratio calibrations.

Here, we used the code IZI (Blanc et al. 2015; Mingozzi et al. 2020) to simultaneously infer $q$ and $Z$ in a Bayesian method by comparing the emission line fluxes in the stacked spectra with photoionization models. We adopted the models of Levesque et al. (2010) that were calculated using the high mass-loss tracks, a constant star formation history, a Salpeter IMF with a $100 M_{\odot}$ upper cutoff, an electron density of $100 \text{ cm}^{-3}$, and an age of 6 Myr. We verified that the radial trends of $q$ and $Z$ remain similar using the Levesque et al. (2010) models with different parameters or different photoionization models (e.g., Kewley et al. 2001; Dopita et al. 2013). We note that the photoionization models assume that the density of the emitting gas is constant across the scales of interest here, which may not be valid. The central panel of Fig. 6 shows the $q$ and $Z$ values estimated by IZI, which also takes into account upper limits on the line fluxes. The ionization parameter is found to decrease from $\log q$ (cm s$^{-1}$) $\approx$7.7 to $\approx$7.3 on going outward from the central 0–5 kpc region to the outer 20–30 kpc region, most likely due to a weaker radiation field in the outskirts of galaxies. While there appears to be a transition in the line ratios from the inner 20 kpc to the outer 20–30 kpc region in the O32–R23 diagram, there is no significant variation in the gas-phase metallicity between the inner and outer regions based on photoionization models, with $12 + \log(O/H)^2$ going from $\approx$8.3 at $\le$5 kpc to $\approx$8.4 at 20–30 kpc. Using the empirical calibrations of Curti et al. (2017), we get a

similar trend of metallicity, with $12 + \log(O/H) \approx 8.4$ at $\le$5 kpc and $\approx$8.5 at 20–30 kpc.

The above analysis assumes that the warm ionized gas is photoionized. While this assumption is likely valid in the ISM, it may not necessarily hold in the extended disk or CGM, where shocks and turbulence could also ionize the gas. To investigate alternate ionization mechanisms, we compared the line ratios from the stacked spectra with shock models. The right panel of Fig. 6 compares the O32 and R23 line ratios from the stacked spectra with shock models assuming a Large Magellanic Cloud abundance and a pre-shock density of $1 \text{ cm}^{-3}$ for several magnetic parameters and shock velocities from Allen et al. (2008). The “shock only” models include contribution from radiative shocks only, while the “shock+precursor” models include contribution from the H II region ahead of the shock front as well. The line ratios in the inner 5 kpc region are consistent with the “shock+precursor” models with magnetic parameter, $Bn^{-1/2} = 2–4 \mu \text{ G cm}^{3/2}$, and velocity $200 \text{ km s}^{-1}$. This suggests that photoionization is the dominant ionization mechanism in the central region. On the other hand, at 20–30 kpc, the line ratios are consistent with the “shock only” models with magnetic parameter, $Bn^{-1/2} = 2–4 \mu \text{ G cm}^{3/2}$, and velocity $200 \text{ km s}^{-1}$. This indicates that outer regions could be affected by radiative shocks arising due to gas flows and interactions, although the observed velocity dispersions of the emission lines are less than the required shock velocities.

# 4. Discussion and summary

We have presented the average line emission around a statistical sample of 560 galaxies at $z \approx 0.25–0.85$ based on stacking of the MUSE 3D data of two large surveys, MAGG and MUDF. The redshift range was selected to have simultaneous coverage of the [O II] $\lambda\lambda3727, 3729$, H$\beta$ $\lambda$4863, [O III] $\lambda$4960, and $\lambda$5008 emission lines. The average [O II] line emission is detected out to $\approx$40 kpc, and the average [O III] and H$\beta$ emission is detected out to $\approx$30 kpc. For comparison, the average stellar continuum emission is detected out to $\approx$20 kpc. Therefore, the extended emission likely probes the disk-halo interface or the transition region between the ISM and the CGM.

Combining our sample with that of Dutta et al. (2023), who conducted a similar stacking analysis at higher redshifts, we find that the [O II] emission increases independently with both

<footnote n="2">² Solar metallicity, $12 + \log(O/H) = 8.69$ (Asplund et al. 2009).</footnote>
redshift over $z \approx 0.4–1.3$ and with stellar mass over $M_* \approx 10^{6–12} M_\odot$. Based on a control analysis, we find that the enhancement in the average [O II] SB in the inner 20 kpc region tracing the ISM is driven to a large extent by the higher SFRs at higher redshifts and stellar masses. On the other hand, additional factors, such as the evolution of physical conditions in the warm ionized gas, likely contribute to the observed enhancement in the [O II] SB in the outer 20–30 kpc region with redshift and stellar mass.

To investigate the average physical conditions in the extended disks or the disk-halo interface of galaxies, we analyzed the emission line ratios in the stacked spectra. Both the O32 and R23 line ratios are found to decrease from the center out to 30 kpc. The line ratios at <20 kpc occupy a similar region in the O32–R23 diagram as those obtained from integrated spectra probing the ISM of $z \approx 2–3$ galaxies. On the other hand, the line ratios at >20 kpc are shifted in the O32–R23 diagram, which indicates a transition in the physical conditions and/or ionization mechanisms of the warm ionized gas from the ISM to the CGM.

We simultaneously estimated the ionization parameter and gas-phase metallicity under the assumption that the gas is photoionized. For this purpose, we used the Bayesian code izi to compare the emission line fluxes in the stacked spectra with the photoionization models of Levesque et al. (2010). The ionization parameter exhibits a declining trend from the central region out to 30 kpc. This is consistent with a picture in which the number of H II regions decreases with distance from galaxy center, resulting in a weaker radiation field in the outskirts. On the other hand, the metallicity does not show any significant trend, remaining more or less similar up to 30 kpc. As discussed in Sect. 1, several mechanisms acting simultaneously are likely influencing the distribution of metals within and around galaxies. For example, the mixing of metals due to the balance between the inflow of metal-poor gas, outflow of metal-rich gas, and gas transfer due to mergers could lead to a flat metallicity gradient. Several studies have found that the metallicity gradients in the ISM of $z \ge 0.5$ galaxies are generally flat with a large scatter (e.g., Wuyts et al. 2016; Carton et al. 2018; Curti et al. 2020; Wang et al. 2020; Simons et al. 2021; Venturi et al. 2024 see however Møller & Christensen 2020). Studies of the CGM metallicity have also indicated that the CGM is complex, with gas outflows, accretion, and recycling all at play (e.g., Lehner et al. 2019; Pointon et al. 2019; Guo et al. 2020). Extending such studies into the disk-halo interface, we find that the average metallicity gradient is flat ($\approx 0.003 \text{ dex kpc}^{-1}$) out to 30 kpc, indicating an efficient mixing of metals at these larger scales.

We point out that other ionization mechanisms in addition to photoionization by H II regions might be altering the observed line ratios and metallicity measurements in the outer regions. For example, the gas could be ionized by radiative shocks arising due to stellar outflows, supernovae, gas accretion from the CGM, or interactions (e.g., Fossati et al. 2016; Epinat et al. 2018; Chen et al. 2019; Franchetto et al. 2021). By comparing the line ratios from the stacked spectra with those from the shock models of Allen et al. (2008), we find that the line ratios at 20–30 kpc can be explained by shock velocities of $\approx 200 \text{ km s}^{-1}$. However, it may be difficult to reconcile the observed narrow line widths ($\le 100 \text{ km s}^{-1}$) with the required shock velocities, although slow shock models with a velocity dispersion of $\approx 100 \text{ km s}^{-1}$ have been found to produce line ratios consistent with active galactic nucleus-like excitation due to shocks (Rich et al. 2011). In addition, there could be a contribution from the diffuse ionized gas to the line emission in

the outer regions, leading to line ratios significantly different from those produced by H II regions (e.g., Zhang et al. 2017; Vale Asari et al. 2019; Mannucci et al. 2021). The diffuse ionized gas could be ionized by photons leaking out of H II regions, hot evolved stars, and shocks. In the future, additional multi-wavelength observations such as NIR observations of H$\alpha$, [N II], and [S II] emission lines would facilitate distinguishing between different ionization mechanisms.

We note that we stacked the emission of galaxies at random disk orientations, washing out any azimuthally dependent signatures of metallicity gradients if present, such as those due to biconical outflows of metal-rich gas from the center of galaxies or the accretion of metal-poor gas along the plane of the galaxy disks (e.g., Bouché et al. 2013; Péroux et al. 2020). In Dutta et al. (2023), the average [O II] emission at $z \approx 1$ was found to be $\approx 30\%$ more extended along the major axis than along the minor axis, while the average Mg II emission at $z \approx 1$ was found to be enhanced up to 10 kpc along the minor axis in the stacking analysis of Guo et al. (2023). In the current sample, the number of galaxies with orientation measurements from high-spatial-resolution HST imaging is not sufficiently large to probe extended emission from multiple lines in stacking.

Nevertheless, the average metallicity gradient in a statistical sample can still place useful constraints on galaxy formation models. For example, comparing different feedback prescriptions in cosmological hydrodynamical simulations, Gibson et al. (2013) found that an enhanced feedback model gives rise to flat metallicity gradients in galaxies. Recently, using EAGLE simulations, Tissera et al. (2022) found that the median metallicity gradient is close to zero at all redshifts, with mergers and gas accretion both influencing the gradients, while Hemler et al. (2021) found predominantly negative metallicity gradients at all redshifts using TNG-50 simulations. This shows that the gradient can be a useful tool for discriminating between different feedback models. The flat metallicity gradient ($\approx 0.003 \text{ dex kpc}^{-1}$) out to 30 kpc found in this work supports a model in which feedback and interactions regulate the gradients. Future detections of extended emission in multiple lines from the gas around a larger sample of individual galaxies, as well as an extension of the stacking analysis to higher redshifts using IFU data from the JWST and the Extremely Large Telescope, are expected to constrain galaxy formation models even further by characterizing the radial gradients of the physical, chemical, and ionization conditions of gas from the ISM to the CGM.

**Acknowledgements.** We thank the anonymous referee for their useful comments. RD thanks Palle Møller for helpful comments. This work has been supported by Fondazione Cariplo, grant No 2018–2329. This work is based on observations collected at the European Organisation for Astronomical Research in the Southern Hemisphere under ESO programme IDs 197.A-0384 and 1100.A-0528. Based on observations with the NASA/ESA Hubble Space Telescope obtained, from the Data Archive at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Incorporated, under NASA contract NAS 5-26555. Support for Program numbers 15637, and 15968 were provided through grants from the STScI under NASA contract NAS 5-26555. This work used the DiRAC Data Centric system at Durham University, operated by the Institute for Computational Cosmology on behalf of the STFC DiRAC HPC Facility (www.dirac.ac.uk). This equipment was funded by BIS National E-infrastructure capital grant ST/K00042X/1, STFC capital grants ST/H008519/1 and ST/K00087X/1, STFC DiRAC Operations grant ST/K003267/1 and Durham University. DiRAC is part of the National E-Infrastructure.

## References

Allen, M. G., Groves, B. A., Dopita, M. A., Sutherland, R. S., & Kewley, L. J. 2008, ApJS, 178, 20
Aslund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481

Bacon, R., Accardo, M., Adjali, L., et al. 2010, Proc. SPIE, 7735, 773508

Belfiore, F., Maiolino, R., Tremonti, C., et al. 2017, MNRAS, 469, 151

Blanc, G. A., Kewley, L., Vogt, F. P. A., & Dopita, M. A. 2015, ApJ, 798, 99

Bouché, N., Murphy, M. T., Kacprzak, G. G., et al. 2013, Science, 341, 50

Bresolin, F., Kennicutt, R. C., & Ryan-Weber, E. 2012, ApJ, 750, 122

Bruzual, G., & Charlot, S. 2003, MNRAS, 344, 1000

Burchett, J. N., Rubin, K. H. R., Prochaska, J. X., et al. 2021, ApJ, 909, 151

Byler, N., Dalcanton, J. J., Conroy, C., et al. 2018, ApJ, 863, 14

Calzetti, D., Armus, L., Bohlin, R. C., et al. 2000, ApJ, 533, 682

Carton, D., Brinchmann, J., Contini, T., et al. 2018, MNRAS, 478, 4293

Chabrier, G. 2003, PASP, 115, 763

Chen, H.-W., Boettcher, E., Johnson, S. D., et al. 2019, ApJ, 878, L33

Cresci, G., Mannucci, F., Maiolino, R., et al. 2010, Nature, 467, 811

Cullen, F., Shapley, A. E., McLure, R. J., et al. 2021, MNRAS, 505, 903

Curti, M., Cresci, G., Mannucci, F., et al. 2017, MNRAS, 465, 1384

Curti, M., Maiolino, R., Cirasuolo, M., et al. 2020, MNRAS, 492, 821

Curti, M., D'Eugenio, F., Carniani, S., et al. 2023, MNRAS, 518, 425

Davies, R. L., Förster Schreiber, N. M., Genzel, R., et al. 2021, ApJ, 909, 78

Dekel, A., Birnboim, Y., Engel, G., et al. 2009, Nature, 457, 451

Dopita, M. A., Sutherland, R. S., Nicholls, D. C., Kewley, L. J., & Vogt, F. P. A. 2013, ApJS, 208, 10

Dutta, R., Fumagalli, M., Fossati, M., et al. 2020, MNRAS, 499, 5022

Dutta, R., Fumagalli, M., Fossati, M., et al. 2021, MNRAS, 508, 4573

Dutta, R., Fossati, M., Fumagalli, M., et al. 2023, MNRAS, 522, 535

Epinat, B., Contini, T., Finley, H., et al. 2018, A&A, 609, A40

Erb, D. K., Shapley, A. E., Pettini, M., et al. 2006, ApJ, 644, 813

Fitzpatrick, E. L. 1999, PASP, 111, 63

Fossati, M., Fumagalli, M., Boselli, A., et al. 2016, MNRAS, 455, 2028

Fossati, M., Mendel, J. T., Boselli, A., et al. 2018, A&A, 614, A57

Fossati, M., Fumagalli, M., Lofthouse, E. K., et al. 2019, MNRAS, 490, 1451

Franchetto, A., Tonnesen, S., Poggianti, B. M., et al. 2021, ApJ, 922, L6

Fraternali, F. 2017, Astrophys. Space Sci. Lib., 430, 323

Fumagalli, M., O'Meara, J. M., & Prochaska, J. X. 2016, MNRAS, 455, 4100

Galbiati, M., Fumagalli, M., Fossati, M., et al. 2023, MNRAS, 524, 3474

Gibson, B. K., Pilkington, K., Brook, C. B., Stinson, G. S., & Bailin, J. 2013, A&A, 554, A47

Guo, Y., Maiolino, R., Jiang, L., et al. 2020, ApJ, 898, 26

Guo, Y., Bacon, R., Bouché, N. F., et al. 2023, Nature, 624, 53

Hamanowicz, A., Péroux, C., Zwaan, M. A., et al. 2020, MNRAS, 492, 2347

Hemler, Z. S., Torrey, P., Qi, J., et al. 2021, MNRAS, 506, 3024

Johnson, S. D., Chen, H.-W., Straka, L. A., et al. 2018, ApJ, 869, L1

Kauffmann, G., Heckman, T. M., Tremonti, C., et al. 2003, MNRAS, 346, 1055

Kewley, L. J., Dopita, M. A., Sutherland, R. S., Heisler, C. A., & Trevena, J. 2001, ApJ, 556, 121

Kewley, L. J., Rupke, D., Zahid, H. J., Geller, M. J., & Barton, E. J. 2010, ApJ, 721, L48

Kewley, L. J., Nicholls, D. C., & Sutherland, R. S. 2019, ARA&A, 57, 511

Leclercq, F., Verhamme, A., Epinat, B., et al. 2022, A&A, 663, A11

Leclercq, F., Chisholm, J., King, W., et al. 2024, A&A, 687, A73

Lehner, N., Wotta, C. B., Howk, J. C., et al. 2019, ApJ, 887, 5

Levesque, E. M., Kewley, L. J., & Larson, K. L. 2010, AJ, 139, 712

Liang, C. J., Kravtsov, A. V., & Agertz, O. 2016, MNRAS, 458, 1164

Lilly, S. J., Carollo, C. M., Pipino, A., Renzini, A., & Peng, Y. 2013, ApJ, 772, 119

Lofthouse, E. K., Fumagalli, M., Fossati, M., et al. 2023, MNRAS, 518, 305

Maiolino, R., & Mannucci, F. 2019, A&ARv, 27, 3

Maiolino, R., Nagao, T., Grazian, A., et al. 2008, A&A, 488, 463

Mannucci, F., Belfiore, F., Curti, M., et al. 2021, MNRAS, 508, 1582

Mingozzi, M., Belfiore, F., Cresci, G., et al. 2020, A&A, 636, A42

Møller, P., & Christensen, L. 2020, MNRAS, 492, 4805

Morrissey, P., Matuszewski, M., Martin, D. C., et al. 2018, ApJ, 864, 93

Nakajima, K., Ouchi, M., Shimasaku, K., et al. 2013, ApJ, 769, 3

Nelson, E. J., van Dokkum, P. G., Momcheva, I. G., et al. 2016, ApJ, 817, L9

Nielsen, N. M., Fisher, D. B., Kacprzak, G. G., et al. 2023, ArXiv e-prints [arXiv:2311.00856]

Oppenheimer, B. D., & Davé, R. 2008, MNRAS, 387, 577

Oppenheimer, B. D., Schaye, J., Crain, R. A., Werk, J. K., & Richings, A. J. 2018, MNRAS, 481, 835

Péroux, C., & Howk, J. C. 2020, ARA&A, 58, 363

Péroux, C., Nelson, D., van de Voort, F., et al. 2020, MNRAS, 499, 2462

Planck Collaboration XIII. 2016, A&A, 594, A13

Pointon, S. K., Kacprzak, G. G., Nielsen, N. M., et al. 2019, ApJ, 883, 78

Qu, Z., Chen, H.-W., Rudie, G. C., et al. 2023, MNRAS, 524, 512

Revalski, M., Rafelski, M., Fumagalli, M., et al. 2023, ApJS, 265, 40

Rich, J. A., Kewley, L. J., & Dopita, M. A. 2011, ApJ, 734, 87

Rupke, D. S. N., Kewley, L. J., & Barnes, J. E. 2010, ApJ, 710, L156

Sánchez, S. F., Rosales-Ortega, F. F., Iglesias-Páramo, J., et al. 2014, A&A, 563, A49

Sanders, R. L., Shapley, A. E., Jones, T., et al. 2021, ApJ, 914, 19

Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103

Sharples, R., Bender, R., Agudo Berbel, A., et al. 2013, The Messenger, 151, 21

Shirazi, M., Brinchmann, J., & Rahmati, A. 2014, ApJ, 787, 120

Simons, R. C., Papovich, C., Momcheva, I., et al. 2021, ApJ, 923, 203

Somerville, R. S., & Davé, R. 2015, ARA&A, 53, 51

Stephenson, H. M. O., Stott, J. P., Cullen, F., et al. 2024, MNRAS, 527, 7891

Tacchella, S., Carollo, C. M., Förster Schreiber, N. M., et al. 2018, ApJ, 859, 56

Tissera, P. B., Rosas-Guevara, Y., Sillero, E., et al. 2022, MNRAS, 511, 1667

Tremonti, C. A., Heckman, T. M., Kauffmann, G., et al. 2004, ApJ, 613, 898

Tumlinson, J., Peeples, M. S., & Werk, J. K. 2017, ARA&A, 55, 389

Turner, M. L., Schaye, J., Crain, R. A., et al. 2017, MNRAS, 471, 690

Vale Asari, N., Couto, G. S., Cid Fernandes, R., et al. 2019, MNRAS, 489, 4721

Venturi, G., Carniani, S., Parlanti, E., et al. 2024, A&A, in press, https://doi.org/10.1051/0004-6361/202449855

Wang, X., Jones, T. A., Treu, T., et al. 2020, ApJ, 900, 183

Weng, S., Péroux, C., Ramesh, R., et al. 2024, MNRAS, 527, 3494

Werk, J. K., Prochaska, J. X., Tumlinson, J., et al. 2014, ApJ, 792, 8

Wilman, D. J., Fossati, M., Mendel, J. T., et al. 2020, ApJ, 892, 1

Wuys, S., Förster Schreiber, N. M., Genzel, R., et al. 2012, ApJ, 753, 114

Wuys, E., Wisnioski, E., Fossati, M., et al. 2016, ApJ, 827, 74

York, D. G., Adelman, J., & Anderson, J. E. 2000, AJ, 120, 1579

Zabl, J., Bouché, N. F., Wisotzki, L., et al. 2021, MNRAS, 507, 4294

Zaritsky, D., Kennicutt, R. C., Jr., Huchra, J. P., et al. 1994, ApJ, 420, 87

Zhang, K., Yan, R., Bundy, K., et al. 2017, MNRAS, 466, 3217