Open Access
Issue
A&A
Volume 686, June 2024
Article Number A81
Number of page(s) 15
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/202348399
Published online 30 May 2024

© The Authors 2024

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1. Introduction

The evolution of galaxies is associated with global features, such as morphology and kinematics, as well as properties that drive the evolution of the stellar populations that they have (e.g., chemical composition, star formation rates (SFRs), initial mass, and luminosity functions). Physical phenomena such as inflow and outflow processes are also known to be key in driving the chemical evolution of galaxies.

Seminal works in the 70s showed that supernovae (SNe) explosions and stellar winds play a significant role in the chemical evolution. Supernovae explosions bring enough energy to ionize and accelerate particles in the interstellar medium, heating the environment and expelling the gas component through galactic winds. This process is more efficient in low-mass galaxies, where reaching the escape velocity is easier because of their shallow potential wells (Larson 1974; Lequeux et al. 1979). In addition, inflows of pristine gas can dilute metals, increasing star formation and leading to the synthesis of new heavy elements (Finlator & Davé 2008). The balance between these processes naturally produces the mass-metallicity relation (MZR), according to which high-mass galaxies present high gas-phase metallicities.

In the last two decades, the MZR has shown itself to be a powerful tool to constrain the chemical evolution of galaxies. Tremonti et al. (2004), using a large sample of star-forming (SF) galaxies, constructed the MZR of the local Universe, showing statistically robust evidence of the existence of this relation and confirming the suggestion introduced by Larson (1974) that low abundances and a tight (±0.1 dex) MZR are driven by efficient gas loss. On the other hand, the MZR tends to flatten for massive galaxies with M ≥ 3 × 1010M (Kauffmann et al. 2004) because of galactic downsizing, where the closed-box-like behavior begins (Cowie et al. 1996; Neistein et al. 2006; Sánchez et al. 2013; Spitoni et al. 2020). The MZR also evolves in redshift: both the gas-phase metal content and stellar masses of galaxies increase across cosmic time, where the higher-redshift MZR (z < 3.5) presents a higher downshift with respect to the local Universe (Savaglio et al. 2005; Erb et al. 2006; Maiolino et al. 2008; Moustakas et al. 2011; Zahid et al. 2014; Gao et al. 2018). Evidence of smooth metallicity gradients and the weak dependence on the SFRs in SF galaxies at z ∼ 3 show that processes involving material flows are stronger at higher redshifts, naturally resulting in a shallower slope across cosmic time (Troncoso et al. 2014). On the other hand, the formation processes of galaxies in the hierarchical growth scenario and the key role of SNe feedback, which acts as the main regulator of the star formation and morphology of the galaxies, proves that mergers depending on the gas fractions of their galaxies drive changes in the evolution of the MZR (Tissera et al. 2005; de Rossi et al. 2007; Kobayashi et al. 2007; Scannapieco et al. 2008).

The internal phenomena of SF galaxies that drive the MZR have been extensively studied (see Maiolino & Mannucci 2019 for a detailed review). However, a similar investigation of the physical processes intervening in the chemical evolution of galaxies in dense environments such as groups and galaxy clusters is missing. Therefore, constructing the MZRs in a galaxy cluster allows the exploration of the chemical evolution of its members in a broad range of masses at a given redshift, under a variety of dynamical and environmental effects.

Clusters of galaxies present deep potential wells originating in the interaction of smaller-scale structures or the evolution of overdensities in the early Universe (Press & Schechter 1974; Springel et al. 2005). In the former scenario, the “preprocessing” can affect galaxies’ chemical evolution as well as their morphology and kinematics (Fujita 2004). In a recent work, Kim et al. (2023) showed that in galaxy clusters with ongoing dynamical activity, early infallers will stop the star formation activity due to strong strangulation processes, while the recent infallers will suffer from ram-pressure stripping (RPS) when galaxies interact with the intracluster medium (ICM). Moreover, the scenario will be even more chaotic when infall makes physical mechanisms such as accretion, tidal interactions, and mergers more likely (Barton et al. 2000; Alonso et al. 2007; Ellison et al. 2008, 2019).

Up to date, the general picture of the role of the physical phenomena driving the evolution of galaxies in clusters is still unclear. However, the virial mass of galaxy clusters plays a crucial role, underlying differences in the chemical content of cluster members with respect to field galaxies (Petropoulou et al. 2012). This allows the exploration of environmental mechanisms that affect the chemical evolution of galaxies and to investigate the origin of the quenching of galaxies in more detail. For example, in high-mass-density environments, galaxies found at the cluster center or group center will experience enhancements in the gas-phase abundances up to ∼0.05 dex (Ellison et al. 2009). This could be produced (i) by a metal-poor inflow suppression toward the galactic center, (ii) preprocessed gas from the ICM, (iii) tidal interactions producing starbursts, (iv) strangulation, and (v) RPS (Gunn & Gott 1972) acting on galaxies with nonzero metallicity gradients (Ellison et al. 2009; Petropoulou et al. 2011, 2012; Gupta et al. 2016). On the other hand, although galaxy clusters with M ∼ 1014M also contain SF galaxies suffering environmental effects, these do not seem to be strong enough to produce significant changes in the chemical content of the cluster’s members. However, a weak metallicity dependence on clustercentric radii is seen (Petropoulou et al. 2012; Gupta et al. 2016; Vaughan et al. 2020; Lara-López et al. 2022), showing the relevance of galaxy cluster masses in the exploration of the phenomena that shape the chemical evolution of galaxies.

In this work, we study the role of galaxy clusters in the MZR, by using SF galaxies at redshift z ∼ 0.3 belonging to AC114 (ABELL S1077). AC114 (α = 22h58m52s.3, δ = −34° 46′55″; J2000) has been widely studied near the cluster core, rcore = 0.27 h−1 Mpc (Couch & Sharples 1987). The cluster core shows strong lensing arcs produced by the central dominant (cD) galaxy of AC114 (Smail et al. 1995; Natarajan et al. 1998; Campusano et al. 2001). AC114 hosts a higher fraction of blue, late-type galaxies than lower-redshift galaxy clusters (Couch et al. 1998; Sereno et al. 2010), which is a typical symptom of the Butcher-Oemler effect (Butcher & Oemler 1978). Roughly 60% of galaxies are found outside the core region, and its red sequence is separated from the blue galaxies (Proust et al. 2015). A soft extended X-ray tail originates close to the cluster core, extending ∼1.5′ toward the southeast (∼400 kpc), suggesting that the substructure could be part of an ongoing merger event (De Filippis et al. 2004).

This is the third paper in a series. In Paper I, Proust et al. (2015) revisited the dynamical state of AC114. They find that the large velocity dispersion of 1893 82 + 73 km s 1 $ 1893_{-82}^{+73}\, \mathrm{km \,s^{-1}} $ and the shape of its infalling pattern is indicative of AC114 being in a radial relaxation phase, showing a very elongated radial filament structure spanning 12 000 km s−1. From the dynamical activity, they find that AC114 is a massive galaxy cluster with a virial mass of Mvir = (4.3  ±  0.7)×1015h−1M. In Paper II, Saviane et al. (2023), by studying stellar populations of passive galaxies and emission line galaxies (ELGs), find that the stellar MZR is steeper than the local Universe analogue. Moreover, the stellar evolution of galaxies in the center is faster than that found in those on the outskirts of AC114, where the latter are still forming stars, a signal that strangulation processes are stronger in the central region of AC114.

In this work, we study the AC114-MZR through a chemical analysis of the gas-phase component of its SF galaxies. In Sect. 2, we present the spectroscopic and photometric observations used, together with the data reduction. In Sect. 3, we describe the data analysis, including the redshift, mass, and metallicity measurements and estimates. In Sect. 4, we present the MZRs constructed for the foreground, cluster members, and background galaxies, together with their properties. In Sects. 5 and 6 we discuss the results that we obtain and presents our conclusions.

In this paper, we adopt the following cosmological parameters: H0 = 69.6 km s−1 Mpc−1, ΩM = 0.286, and ΩΛ = 0.714 in a flat Universe according to Wright (2006).

2. Observations and data reduction

2.1. VIMOS spectroscopic data

In this work, we use the same spectroscopic data as Proust et al. (2015) and Saviane et al. (2023), hereafter P15 and S23. The data was obtained with the VIsible MultiObject Spectrograph (VIMOS, Le Fèvre et al. 2003) mounted on the ESO Very Large Telescope (VLT), under the program 083.A-0566 (PI: I. Saviane), between August 16 and September 25, 2009. VIMOS was a visible (360–1000 nm) wide field imager and multi-object spectrograph mounted on the Nasmyth focus B of UT3 Melipal. The instrument is made of four identical arms with a field of view (FOV) of 7′×8′ with a 0.205″ pixel size and a gap between each quadrant of ∼2′. Seven exposures each for the HR-red and MR grisms (630–870 nm, R = 2500 and 500–1000 nm, R = 580, respectively, for a 1″ slit) were obtained, for a total observing time of ∼14 h.

The galaxy selection was made from the pre-imaging frames of the cluster. Nonstellar objects in the region were selected by eye in order to punch the maximum number of slits into each of the four quadrants. Such a visual inspection allowed us to discriminate between extended objects and stars.

The VIMOS Pipeline in the EsoReflex environment1 (Freudling et al. 2013) was used to perform bias subtraction, flat normalization and cosmic ray removal in each scientific frame. Then, sky subtraction was done, taking the median of the sky spectrum located between the galaxy spectra in the two-dimensional frames. The spectral flux was calibrated by using two standard stars from Hamuy et al. (1992, 1994): the F-type LTT1788 (V = 13.16, B − V = 0.47) and DA-type LTT7987 (V = 12.23, B − V = 0.05).

Once the galaxy spectra were reduced, we proceeded with emission line corrections to get proper flux measurements. We subtracted stellar continuum spectra by fitting Chebyshev polynomials with the Specutils package (Astropy-Specutils Development Team 2019). Emission line fluxes from collisional de-excitation (O+,O++,N+,S++) were then measured with a Gaussian fitting. To correct for dust and reddening effects, Hα and Hβ fluxes were measured by Gaussian fitting as well, and the observed hydrogen line ratio (FHα/FHβ) was compared with the theoretical expectations. The H line ratio, IHα/IHβ= 2.85, was adopted from case B (Te ≃ 10 000 K) following Hummer & Storey (1987). The reddening constant was calculated as CHβ = [log(IHα/IHβ) − log(FHα/FHβ)]/[f(λ) − f(Hβ)], where f(λ) = ⟨A(λ)/A(V)⟩ is the extinction law at a given wavelength of the modified version of Cardelli et al. (1989) in the optical domain, obtained from O’Donnell (1994). We also assumed that RV = 3.1. Finally, emission line measurements were corrected by applying Iλ/IHβ = (Fλ/FHβ)×10CHβ[f(λ)−fHβ].

Recombination H lines can be affected by underlying stellar absorption, producing changes in the metallicity estimations due to the use of Balmer lines in the most common calibrators. However, the spectra of our sample are contaminated by fringing patterns in λ > 7500 Å, so it is difficult to make an appropriate correction for the stellar absorption (at z ∼ 0.3, the Hα emission line falls at ∼ 8500 Å).

To evaluate a proper correction of this effect, we employed the Hopkins et al. (2003) correction, which was tested in galaxies of the SDSS survey. The authors show that an appropriate stellar absorption correction is about EW = 1.3 Å, which in this work is translated in a flux correction of Fλ ∼ 0.001 and ∼0.002 × 10−16 erg s−1 cm−2 on average for Hα and Hβ, respectively. Therefore, these corrections were applied to the emission line measurements of the Balmer lines.

Our sample consists of 184 spectra of galaxies, where 25 (13%) are not reported in earlier works with the same data, because in this paper the observations were reduced with the latest version (4.1.8) of the VIMOS pipeline. The sample is listed in Table A.2, with galaxies labeled according to the quadrant and slit of extraction. Additional galaxies are listed in Table A.2. with the syntax #quadrant-A#. Our sample is shown on top of an image of AC114 in Fig. 1.

thumbnail Fig. 1.

Observed galaxies (colored circles) on top of an image taken from the DECaLS DR10 viewer for the AC114 cluster. Quadrants Q1, Q2, Q3, and Q4 are counterclockwise, with Q4 in the bottom right clump.

Applying the same criteria as in S23, a manual inspection of the reduced spectra was conducted to classify the sample into passive galaxies and ELGs, based on the presence or absence of the [OII]λλ3727, [OIII]λλ4959, 5007, and/or Hγ, Hβ emission lines. Figure 2 shows the difference between a passive galaxy and an ELG in the sample. We identified 103 (56%) passive and 81 (44%) ELGs. The top panel shows the spectrum of the passive galaxy Q4-31 (#quadrant #slit). Only absorption features can be identified in the spectrum, such as CaII H & K, the CH G band, the Mg triplet, and the Na D doublet, shown with red lines from left to right, respectively. The bottom panel of Fig. 2 shows the spectrum of the ELG Q4-14. In this case, only emission lines are present, such as [OII]λλ3727, Hγ,[OIII]λ4363, Hβ, [OIII]λλ4959, 5007, Hα, [NII]λ6584, and [SII]λλ6716, 6731 shown with blue lines from left to right, respectively.

thumbnail Fig. 2.

Reduced spectra of a passive galaxy (top panel) and an ELG (bottom panel). Red lines indicate absorption features, whereas blue lines indicate emission features.

Most of the galaxies are in Q4 (62 galaxies, bottom right in Fig. 1), which covers the central region of AC114. The other 35, 41, and 46 galaxies are grouped in the Q1, Q2, and Q3 quadrants, respectively (counterclockwise). The clumps are produced by the VIMOS CDD-chips gap of ∼2′.

2.2. Optical and near-infrared photometric data

To construct the MZR, photometric data is needed to compute the stellar masses of the galaxies. Optical photometry was found for the VIMOS FoV of AC114 in the DECaLS DR10 Survey (see Dey et al. 2019 for an overview of the DESI Legacy Imaging Surveys). The DECaLS program made use of the Dark Energy Camera (DECam) on the Blanco 4m telescope, located at the Cerro Tololo Inter-American Observatory. DR10 covers the South Galactic Cap region at Dec ≤ 34° in the g, r, z, i bands. 121 (66%) and 139 (76%) galaxies in our sample have available photometry in the g and i bands, respectively. However, no photometry is available in the r and z bands.

The top panel of Fig. 3 shows the distribution of apparent magnitudes for both g and i band passes with blue and orange lines, respectively. Galaxies observed in the g band cover a magnitude range from 18.83 to 25.10, with a clear peak at 22.46, while the respective i band magnitudes cover a range of 16.67 to 23.77, with three prominent clumps centered on 22.24, 20.76, and 19.32.

thumbnail Fig. 3.

Top panel: distribution of photometric data collected from DECaLS DR10 survey. The blue line indicates the distribution of g-band magnitudes and the orange line indicates the distribution of i-band magnitudes. Bottom panel: distribution of photometric data collected by the VIKING DR4 survey. The blue line indicates the distribution of H-band magnitudes and the orange line indicates the distribution of Ks-band magnitudes.

Bell & de Jong (2001), studying stellar mass-to-light ratios in optical and near-infrared (NIR) band passes, found that NIR luminosities are less sensitive to the past star formation history, smoothing out variations with respect to the color. Therefore, they concluded that NIR luminosities lead to more accurate stellar masses. With this motivation, we also searched for available NIR photometric data, which we found in the VIKING DR4 ESO survey. The VISTA Kilo-degree Infrared Galaxy Survey (VIKING)2 is a wide-area (covering a final area of 1350 square degrees), intermediate-depth (5σ detection limit J ∼ 21 on Vega system) NIR imaging survey in the five broadband filters, Z, Y, J, H, and Ks. The sky coverage is at high galactic latitudes and includes two main stripes of ∼70 × 70° each: one in the South Galactic cap near Dec ≃ −30° and one near Dec ∼ 0°. A cone search of 25′ was conducted, centered on the cluster, and 121 (66%) of the galaxies in our data set were found. Near-infrared photometry in VIKING DR4 was only available in the H and Ks bands.

The bottom panel of Fig. 3 shows the distribution of apparent magnitudes for our galaxy sample for both H and Ks bands, shown with blue and orange lines, respectively. Galaxies observed in the H and Ks bands cover a range of magnitudes from 14.90 to 20.60 and from 14.10 to 19.79, respectively.

As is discussed in P15 and S23, our galaxy sample reaches luminosities as faint as R ∼ 21, but it starts to suffer incompleteness at R ∼ 20. We searched in the galaxy sample for objects with values R ∼ 20 and found that the galaxy sample is incomplete at g ∼ 21.4 and i ∼ 19.8 mag in the optical bands, whereas the same is seen at H ∼ 17.7 and Ks ∼ 17.3 mag in the NIR band passes. The photometry is used to compute the stellar masses via mass-to-light ratios in the galaxy sample. The incompleteness starts to be significant in galaxies with ∼1010M.

The cluster members in the galaxy sample are well distributed in terms of mass, spanning values from ∼109M to ∼1011M, as is expected for a construction of the MZR. Therefore, the incompleteness does not affect our results.

3. Data analysis

3.1. Redshift estimations

The galaxies’ redshift, z, was measured by cross-correlation (Tonry & Davis 1979) between the VIMOS spectra and synthetic ones at z = 0 in the wavelength range from 3700 Å to 8800 Å. In this work, we made use of the Bagpipes models3 (Carnall et al. 2019) to generate a standard passive and a SF synthetic spectrum, shown in Fig. 4 with blue and orange, respectively. We considered a galaxy whose spectrum presents a evolved stellar continuum with relevant absorption lines typically found in old stellar populations to represent a passive galaxy. We achieved this, generating a synthetic spectrum of a galaxy with stellar mass M* = 1010M and 10 Gyr since the star formation began by using the exponential-τ model. On the other hand, a SF galaxy has a small stellar continuum with prominent emission lines: we calculated the spectrum of a galaxy with a stellar mass of M* = 108M, 1 Gyr after the star formation began, by using the exponential-τ model, and an ionization parameter, log(U) = −2.5.

thumbnail Fig. 4.

Synthetic spectra reproduced with the Bagpipes models. The blue spectrum represents a passive galaxy. The orange spectrum represents an SF galaxy.

Once the synthetic spectra were ready to use, we cross-correlated them with the VIMOS spectra. We used the IRAF4 task “fxcor” to estimate the redshift of galaxies. The galaxies of the sample have ∼0.1 < z < 0.7, with uncertainties of 0.0014 on average. The redshift distribution of the sample is shown in Fig. 5. Emission line galaxy and passive galaxy distributions are shown with blue and red lines, respectively. The distribution reveals that passive galaxies are found at the cluster redshift, meaning 0.28 < z < 0.34 according to P15. In this region, 58 (73%) of a total of 80 galaxies are passive, while 22 (27%) are ELGs. Outside the cluster redshift, it is clear that there is an additional structure that is dominated by ELGs centered on z ∼ 0.41. According to P15, this background structure has a probability of ∼60% of being a group.

thumbnail Fig. 5.

Redshift distribution of the galaxy sample. Upper panel: distribution of passive galaxies (red line) and ELGs (blue line). Lower panel: redshift distribution of the complete sample.

3.2. Stellar mass estimates

3.2.1. Synthetic photometry

The masses of the galaxies were obtained by using stellar mass-to-light ratios. The data acquired from DECaLS DR10 give us photometry for 66% and 76% of the galaxies in the sample in the g and i bands, respectively. For galaxies without photometry, we obtained synthetic photometry by integrating the spectral flux of the galaxies. Optical stellar mass estimates are done via the theoretical mass-to-light ratio in the i band. Therefore, synthetic photometry was estimated for this filter only because it has available photometry for more galaxies than the g band does.

We made use of the PYPHOT5 (Fouesneau 2022) Python package. As a first step, we calculated the synthetic photometry of the LTT7987 standard star. We obtain magnitudes underestimated by ΔMi = Mstandard, i − Mpyphot, i = 0.16. Therefore, this amount was applied as a correction when the artificial magnitudes were computed. In addition, the 1″ slits used in the VIMOS observations do not entirely cover the projected area for most of the galaxies of the sample, which is ∼3″. Therefore, we expect that synthetic photometry will present overestimated apparent magnitudes (underestimated luminosities). In fact, this is seen as an offset in the top panel of Fig. 6, which presents a comparison between the artificial magnitudes and those collected from DECaLS DR10. Synthetic magnitudes are ∼1.49 larger in the i band, with a scatter of σ = 0.46. Hence, overestimations were corrected by subtracting this offset from the artificial magnitudes. The resulting correction is shown in the bottom panel of Fig. 6. We adopted the scatter, σi, as a statistical error of the artificial magnitudes. This is similar to the optical magnitudes from DECaLS, which present a mean uncertainty of 0.47 mag.

thumbnail Fig. 6.

Comparison between DECaLS DR10 (on the x axis) and those from artificial magnitudes (on the y axis). Colors are the same as in Fig. 5. The gray line is an identity line. Top panel: comparison between synthetic and DECaLS magnitudes in the i band. Bottom panel: comparison between synthetic and DECaLS magnitudes in the i band with the applied correction.

3.2.2. Optical-based stellar masses

To compute the absolute magnitudes, we followed the same procedure reported in S23, which converts apparent magnitudes to absolute ones, with i − Mi = 5log DL − 5. Here, DL is the luminosity distance, which is a function of the redshift. A relation was computed for a few select values following Wright (2006). Interpolating a second-degree polynomial it was found that the luminosity distances expressed in Mpc obey DL = 2245.7 z2 + 4682.4 z − 20.838, valid only in z ∈ [0.05, 0.7] (see their Fig. 11, S23).

Absolute magnitudes in the i band were converted to luminosities, with L = 10−0.4(Mi − Mi, ⊙) L. Following S23, for ELGs and passive galaxies we adopted ages of 2 Gyr and 10 Gyr, respectively. By using the theoretical age-mass-to-light relation from E-MILES models (Vazdekis et al. 2012, 2016; Ricciardelli et al. 2012), the corresponding ages gives: M/Li = 0.75  ±  0.20 (ELGs) and 2.48  ±  0.52 (passive). These models adopt evolutionary tracks from Girardi et al. (2000), assuming a sub-solar metallicity of [M/H]= − 0.4 and a Chabrier initial mass function (IMF). Finally, mass-to-light ratios together with the luminosities were converted into the stellar masses of the galaxies. Uncertainties were propagated according to the expression δ f = i ( f / x i ) 2 δ x i 2 $ \delta f=\sqrt{\sum_{i}(\partial f/\partial x_{i})^{2}\delta x_{i}^{2}} $, where f is the propagated function, xi, and δxi are the parameters with their respective uncertainties. Optical-based mass estimations shows mean uncertainties of ∼2%. Table A.3 presents absolute magnitudes in the optical band used in this work, together with the respective mass estimations. The distribution of masses is shown in the upper panel of Fig. 7. Masses of ELGs are distributed between ∼108 and ∼1011M, whereas those of passive ones lie from ∼109 to ∼1012M. From the upper panel of Fig. 7 it is easy to see that, in general, ELGs are less massive than passive galaxies: masses are more concentrated at ∼109.5 and ∼1010.5M for the ELGs and passive components, respectively.

thumbnail Fig. 7.

Central panel: comparison between i and Ks-based masses. Blue and red dots represent ELGs and passive galaxies, respectively. The gray line corresponds to the 1:1 line ratio. Upper panel: distribution of optical luminosity based stellar masses. The black line shows the mass distribution of the complete sample, whereas the blue and red areas indicate the distribution of ELGs and passive galaxies, respectively. Right panel: distribution of near-infrared luminosity based stellar masses.

3.2.3. Near-infrared-based stellar masses

Given that we were able to obtain NIR luminosities for 66% of the galaxy sample, optical-based and NIR-based stellar masses can be compared. Bell & de Jong (2001), by using spectro-photometric models of galaxy evolution, found that there are substantial variations in the stellar mass-to-light ratios within and among galaxies as a function of the color, showing variations of a factor of three and seven in optical band passes and two in the NIR ones. These findings drive to a natural conclusion that NIR-based mass estimations are the best way to obtain the stellar masses of galaxies, because recent star formation can generate optical light that dominates the SED, leading to underestimates of the stellar mass, whereas recent star formation does not affect much NIR light.

E-MILES models were used one more time to compute mass-to-light ratios as a function of the age, resulting in M/LKs = 0.26  ±  0.09 (ELGs) and M/LKs = 0.91  ±  0.13 (passive). Ks absolute magnitudes, together with NIR-based stellar mass estimations, are presented in Table A.3. The mass distribution is also shown in the right panel of Fig. 7. Emission line galaxies span values from ∼109 to ∼1011M, and passives ones from ∼1010 and ∼1012M. Despite the lack of NIR photometry for 34% of the sample, the feature that passive galaxies are more massive in general than ELGs is clearly seen.

The central panel of Fig. 7 shows a comparison between optical-based and NIR-based stellar masses. The comparison shows that the mass estimations agree, following a 1:1 relation, for galaxies with M* < 1010.5, where the sample is mostly dominated by ELGs. On the other hand, at masses M* > 1010.5, galaxies tend to fall above the 1:1 relation by 0.23 orders of magnitude on average. This region is mostly dominated by passive galaxies. However, this offset is within the uncertainty of the mass estimates, which is σ = 0.29 mag. Therefore, it does not change our conclusions about the MZR construction.

3.3. Metallicity estimations

3.3.1. Classification of star-forming galaxies

In Sect. 2.1 we define ELGs by the presence of emission lines in their spectra. However, the ionization sources that trigger emission lines can be produced by active galactic nuclei (AGNs) and/or SF activity. To discriminate the actual nature of our galaxies, we made use of the emission line ratios, [NII]λ6584/Hα and [OIII]λ5007/Hβ, in the Baldwin-Philips-Terlevich (BPT) diagram (Baldwin et al. 1981). In fact, the BPT diagram is a useful tool with which to identify the dominant ionization source of emission line galaxies. Figure 8 shows the ELGs of the sample in the diagram. The red dots are cluster members and the black diamonds correspond to foreground and background galaxies. The dashed black line is the Kewley limit (KL). Galaxies found below the KL have dominant ionization from HII regions, where ongoing star formation is present. However, ELGs can have ionization from both HII regions and AGN activity, i.e. they are composite galaxies. Therefore, we used the Kauffmann et al. (2003) demarcation limit to select pure SF galaxies and also study the possible effects of composite galaxies in the MZR. Composite galaxies are shown in Fig. 8 with open blue diamonds. Two composite galaxies are present in AC114, two belong to the foreground, and three belong to the background galaxies.

thumbnail Fig. 8.

Baldwin-Philips-Terlevich diagram of the sample. The dashed black line is the KL at z = 0 (Kewley et al. 2001), which excludes pure AGNs. The solid black line is the Kauffmann et al. (2003) limit, which excludes composite galaxies (ionization from star formation and AGNs). The red points are galaxies at the cluster redshift. Black diamonds correspond to foreground and background galaxies. According to both demarcation limits, all ELGs are SF galaxies. Data points with blue open diamonds are galaxies classified as composite galaxies, with two of them being composite galaxies at the cluster redshift and five of them found outside the cluster redshift.

The fact that there is not a significant number of composite galaxies contaminating the galaxy sample allows us to make a direct and homogeneous comparison of the MZR for AC114, together with the foreground and background structures, without any other kind of contamination regarding the ionization sources.

3.3.2. Gas-phase metallicities

The first half of the construction of the MZR was done with optical and NIR luminosities to estimate the stellar mass of the galaxies. We now continued to estimate the metallicities of the SF galaxies in the sample.

The best technique to estimate gas-phase metallicities is through measuring intensities of emission line ratios sensitive to the electron temperature, Te, together with the computation of line emissivities. This technique is the so-called “direct method” (Peimbert 1967; Aller 1984). However, Te-sensitive emission line ratios involve lines such as the auroral [OIII]λ4363, [NII]λ5755, or [SII]λλ4068, 4076. These lines are ∼3 − 4 orders of magnitude fainter than Balmer emission lines (Maiolino & Mannucci 2019). Therefore, auroral lines can be well measured in metal-poor or high-temperature environments. Otherwise, auroral lines get lost in the noise of the galactic spectra.

Some alternative methods of estimating metallicities have been proposed; the so-called “strong-line methods”, calibrated on Te-based oxygen abundances or photoionization models (e.g., Pilyugin 2001; Kobulnicky & Kewley 2004; Edmunds & Pagel 1984; Ferland et al. 1998; Pérez-Montero 2015). In this work, we chose the well-known strong-line O3N2 index, defined as O 3 N 2 = log ( [ OIII ] λ 5007 / H β [ NII ] λ 6584 / H α ) $ \mathrm{O3N2}= \log{\left( \frac{\mathrm{[OIII]}\lambda5007/\mathrm{H}\beta}{\mathrm{[NII]}\lambda6584/\mathrm{H}\alpha} \right)} $, which correlates linearly with oxygen abundances in this way: 12 + log(O/H) = 8.73 − 0.32 × O3N2 (Pettini & Pagel 2004). There exist at least ten different calibrators in the literature regarding gas-phase metallicity estimations (see Kewley & Ellison 2008 for a detailed analysis of the most common calibrators and comparisons between them). Despite the significant differences between some calibrators (up to 0.7 dex), the authors also provide consistent conversions between different calibrators within uncertainties of 0.2 dex. Those will also be used in Sect. 4 to make comparisons with the MZRs reported in the literature.

We chose to estimate metallicities with the O3N2 index for several reasons: (i) this index is calibrated with Te-based oxygen abundances, which is the most reliable method despite its intrinsic biases such as temperature fluctuations, underestimations in the high metallicity range, and the way of estimating the electron temperature (Kobulnicky & Skillman 1996; Pilyugin & Grebel 2016; Cameron et al. 2023). (ii) The relation between O3N2 and oxygen abundances is not bi-valued, at odds with, for instance, the R23 index; (iii) the emission lines used in the O3N2 calibrator are closer together (Hβ − [OIII]λ5007 in ∼ 150 Å and Hα − [NII]λ6584 at ∼ 20 Å); hence, there are no significant noise changes in the observed spectra, the internal absorption does not affect the measurements, and the flux calibration is also less of a concern. (iv) Our constructed MZRs with field-foreground and field-background galaxies show oxygen abundances estimates consistent with those expected at their respective redshifts, despite the intrisic bias of the O3N2 index with the ionization.

Oxygen abundances of the sample of galaxies are shown in Table A.4. The estimated oxygen abundances have values from 8.10  ± 0.06 < 12 + log(O/H)O3N2 < 8.77  ± 0.19 for 54 (67%) SF galaxies. The other SF galaxies do not have all these four emission lines to compute their respective metallicities.

4. Results

The reported statistical parameters relative to the analysis presented in this section are listed in Table A.1.

4.1. The MZR of AC114, foreground and background galaxies

We constructed the MZRs with 54 SF galaxies for three subsamples labeled: foreground (0.1 < z < 0.28) with 19 SF galaxies, AC114 (0.28 ≤ z ≤ 0.34) with 20 SF members, and background (0.34 < z < 0.7) with 15 SF galaxies. The redshift boundaries used were adopted from P15. In Fig. 9, the MZRs are compared with the local Universe MZR from Tremonti et al. (2004). To do the comparison, we used the Kewley & Ellison (2008) conversion to get the local Universe MZR in O3N2-based oxygen abundances. The left panel of Fig. 9 shows the MZRs of the foreground and background galaxies. Foreground and background galaxies’ oxygen abundances correlate with the stellar mass, showing the existence of foreground and background MZRs (hereafter fMZR and bMZR). The non-parametric Spearman correlation test was performed for galaxies with < 1σ scatter, showing correlation coefficients of r = 0.97 (p-value 2.6 × 10−8) and r = 0.90 (p-value 1.5 × 10−4). The 1σ scatter is σ = 0.118 and σ = 0.120, showing that both the fMZR and the bMZR are tightly correlated, respectively. Additionally, they are downshifted by 0.09 and 0.20 dex in metallicities with respect to the local MZR, a clear sign of evolution with redshift of the global MZR: at a given mass, higher-redshift galaxies will present lower metal content.

thumbnail Fig. 9.

Mass-metallicity relations constructed with the SF galaxies of the sample compared with the local Universe MZR from Tremonti et al. (2004). The left panel shows the foreground MZR (fMZR) with blue dots and background MZR (bMZR) with red dots. Lines correspond to linear fits. The shaded area corresponds to the 1σ scatter for the fMZR and bMZR. Right panel: MZR with the SF members of AC114 (cMZR). The black line and shaded area are the linear fit and the respective 1σ scatter. The dashed blue and red lines are the fMZR and bMZR linear fits. The open circles (left) and open circles (right) show galaxies whose stellar mass was estimated using synthetic photometry.

The shallower slope of the fMZR (0.236  ±  0.003) with respect to the bMZR (0.312  ±  0.009) is produced mainly by two physical phenomena: (i) galaxies grow hierarchically across the cosmic time, and therefore wet mergers (i.e., mergers between low-mass galaxies with significant gas fractions) can enhance the gas-phase metal content, while low-mass galaxies evolve, due to the production of stars that pollute the ISM with processed material in the last stages of their evolution (Tissera et al. 2005). (ii) Outflow processes are stronger at higher redshift, making it easier for gas particles to get expelled in shallower potential wells (Troncoso et al. 2014). Hence, the evolution in redshift of the MZR in the low-mass range is stronger than in the high-mass range, where deeper potential wells will produce a more efficient star formation. The different mass range of the fMZR and bMZR is a natural consequence of seeing galaxies in a range of distances with a fixed exposure time, as is the case with the DECaLS survey and AC114.

We ignored evolution effects (galaxies at higher redshift are in different stages of their evolution) because we focused our analysis on studying differences in the shape and the evolution of the MZR. The evolution effects mainly affect the construction of the bMZR because it represents galaxies at redshift 0.34 < z < 0.7. If we take into account the evolution effects, the scatter of the bMZR will decrease. However, we would not be able to produce a significant bMZR because there are not enough galaxies in our sample.

The right panel of Fig. 9 presents the MZR constructed at z ∼ 0.3 with SF galaxies of the AC114 galaxy cluster (hereafter cMZR). The first clear feature in the cMZR is that the slope of the linear fit (0.241  ±  0.008) is slightly shallower than both the fMZR and bMZR slopes (shown with the black line), with a larger scatter of σ = 0.17 and a Spearman rank of r = 0.65 (p-value = 0.009). These values suggest that physical phenomena are suppressing the chemical evolution of high-mass galaxies, increasing the observed scatter and weakening the correlation between mass and metallicity. The second important feature seen is that the high scatter of the cMZR and the effect in the high-mass range of the shallower slope with respect to the bMZR indicates that SF members of AC114 are either similar to or metal-poorer than the background galaxies; the opposite of the expected evolution of the MZR.

Classically, the MZR and its evolution in redshift reported in the literature is constructed with (i) field SF galaxies and (ii) statistically complete samples of SF galaxies at a given redshift, where the role of galaxy cluster environments is negligible. It is well known that at a given mass, gas-phase metallicities are more metal-poor at higher redshift (or at a given oxygen abundance, galaxies are less massive at higher redshift). Several physical phenomena, such as SNe and galactic winds, the infall of pristine gas, outflow processes, IMF variations, and galactic downsizing, together with the formation of galaxies through dry or wet mergers, robustly explain the shape of the MZR. However, all of these are internal physical processes. In dense environments, such as massive galaxy clusters, galaxies can suffer harassment, interaction with closer neighbours (galaxy-galaxy interaction), and/or interaction with the ICM. These external physical processes can change the gas-mass fraction, modifying the gas-phase metal content. Therefore, these effects can change the MZR and may explain the cMZR in Fig. 9.

4.2. Comparison of the fMZR, bMZR, and cMZR with the literature

As a sanity check of the MZRs found above, we proceeded to compare them with those reported in the literature. The fMZR was compared with Maiolino et al. (2008) and Cresci et al. (2012). The bMZR was compared with Savaglio et al. (2005), Maiolino et al. (2008), and Lamareille et al. (2009). Finally, the cMZR was compared with Maiolino et al. (2008) and Vaughan et al. (2020).

Savaglio et al. (2005), Cresci et al. (2012), and Maiolino et al. (2008) use the R23 parameter to estimate metallicities, whereas Lamareille et al. (2009) and Vaughan et al. (2020) use the N2 parameter. Therefore, we followed Kewley & Ellison (2008) to convert R23 and N2 to O3N2.

In Maiolino et al. (2008), the authors parameterize the MZR with the expression 12 + log(O/H) = −0.0864(log M*/M0)2 + K0 at redshifts 0.07, 0.7, 2.2, and 3.5, with M0 and K0 free variables. We interpolated both variables as a function of the redshift to get their parameterization at redshifts 0.25, 0.3, and 0.7.

The left panel of Fig. 10 shows the fMZR compared with field MZRs from (i) Cresci et al. (2012) at z ∼ 0.25, (ii) the Maiolino parameterization at z ∼ 0.25, and (iii) the local Universe MZR of Tremonti et al. (2004). As we can see, there is an agreement between the fMZR and the literature. The galaxies that are not in agreement with the MZRs from the literature are, in fact, those which are beyond the 1σ scatter. Those foreground galaxies that are metal-poor with respect to the MZRs in the literature seem to follow a trend. We discarded that the feature seen is an evolutionary effect because the galaxies that fall in the MZR have redshifts up to 0.28, covering the same range as the outlier galaxies. Three of these four outlier foreground galaxies are at the same redshift: Q4-7, Q4-8, and Q4-35, at 0.169, 0.169, and 0.170, respectively. Here, Q4-7 and Q4-8 are closer neighbors, sharing the same region of the FoV. Q4-35 is part of a crowded region in the FoV. On the other hand, Q3-11 at redshift 0.21 has a closer companion in the FoV but without available spectroscopic data. The fact that these galaxies are closer to other ones may indicate that galaxy-galaxy interactions are responsible for the metal-poor values with respect to the fMZR. This will be discussed in Sect. 5.

thumbnail Fig. 10.

Comparison between the MZRs constructed in this work and the literature. Left panel: comparison between the fMZR, Cresci et al. (2012) MZR, and the Maiolino parameterization, both at redshift z ∼ 0.25. Middle panel: comparison between the bMZR, Savaglio et al. (2005), Lamareille et al. (2009), and the Maiolino parameterization at z ∼ 0.7. Right panel: comparison between the cMZR, Vaughan et al. (2020), and the Maiolino parameterization at z ∼ 0.3. Data points with open diamonds show the galaxies classified as composite in the BPT diagram of Fig. 8.

The middle panel of Fig. 10 shows the bMZR compared with field MZRs from (i) Savaglio et al. (2005), (ii) Lamareille et al. (2009), and (iii) the Maiolino parameterization at z ∼ 0.7. In this panel, our measurements are consistent with those reported in the literature in the high-mass range of the MZRs. On the other hand, the low-mass galaxies in our sample tend to follow the Savaglio et al. (2005) MZR. Again, the background galaxies that fall below the MZRs compared are those that are beyond the 1σ scatter. From these two panels, it is clear that the fMZR and bMZR follow the field MZRs at their respective redshift ranges.

In the right panel of Fig. 10, we compare the AC114-MZR (cMZR) with the Maiolino parameterization and the MZR reported in Vaughan et al. (2020, hereafter V20), with the dotted and dash-dotted green lines, respectively. V20 study a sample of four galaxy clusters at 0.3 < z < 0.6 and compares their cluster members with respect to field galaxies in the mass-metallicity plane at the same redshift. They did not find any difference between those MZRs (see Fig. 7 in V20).

To compare the V20 data with this work, two of the four galaxy clusters were found at a similar redshift of AC114, MS2137, and MACS1931 at z = 0.313 and z = 0.352, respectively. We find that just seven and two galaxies have oxygen abundances estimated, belonging to MS2137 and MACS1931, respectively. The field galaxies in V20 at the redshift of AC114 (0.28 < z < 0.35) span masses from ∼1010 to ∼1011M, whereas our mass estimates in AC114 cover two orders of magnitude starting from ∼109M. Therefore, that sample is not appropriate to compare with our results.

In the right panel of Fig. 10, the linear fit of the AC114 galaxy sample presents a higher downshift than the Maiolino parameterization and V20 with respect to the local Universe MZR. The cMZR is downshifted 0.19 dex with respect to the local Universe MZR, and the scatter of the cMZR is 0.17 dex.

The fact that both the fMZR and bMZR are in agreement with field MZRs at their respective redshifts, together with the cMZR showing a slightly shallower slope, a downshift similar to the bMZR, and a high scatter of 0.17 dex, indicates that there are physical phenomena behind those features, which are likely related to environmental and/or dynamical effects.

4.3. Structure of AC114 and consequences for the chemical evolution of its members

We explored the spatial distribution to gain insights into the changes seen in the cMZR with respect to the fMZR and the bMZR. We first converted the projected RA, Dec and redshift coordinates of the cluster members to physical coordinates. Wright (2006) was used to estimate the respective arcsecond per kiloparsec scale conversion as a function of the redshift for a few select values. A second-degree polynomial was fitted, parameterizing the expression: Kpc ( z ) = 9.800 ( ± 0.300 ) z 2 + 16.700 ( ± 0.200 ) z + 0.335 ( ± 0.007 ) $ \frac{\mathrm{Kpc}}{{{\prime\prime}}}(z) = -9.800( \pm 0.300)\, z^{2} + 16.700( \pm 0.200)\ z + 0.335( \pm 0.007) $, with z ∈ [0.1, 0.7]. At the cluster distance DL = 1670  ±  17 Mpc, the corresponding scale is ∼4.6 Kpc/″, and at the cluster range 0.28 < z < 0.34 the scale spans values from ∼4.3 Kpc/″ to ∼4.9 Kpc/″.

AC114 has a nonspherical shape, being elongated in the southeast-northwest direction (Couch et al. 2001; De Filippis et al. 2004). The dynamical study of P15 shows also that AC114 is a very elongated radial filament, spanning 12 000 km s−1 with a broad velocity dispersion of 1893 82 73 km s 1 $ {\sim} 1893_{-82}^{73} \, \mathrm{km \, s^{-1}} $. These dimensions correspond to ∼330 Mpc in the radial direction and ∼9 Mpc in the perpendicular component, the latter being a lower limit, as the calculations are restricted to the FoV of the VIMOS observations.

The three-dimensional representation of the cluster members is shown in Fig. 11 for SF and passive galaxies with blue and red dots, respectively. Figure 11 shows that the center of the cluster (at 0, 0, 0 in X, Y, Z) is populated mostly by passive galaxies; SF galaxies tend to be found at the outskirts.

thumbnail Fig. 11.

Three-dimensional projection of the AC114 members, expressed in Mpc with respect to the cluster center. The x, y, and z axes correspond to the RA,Dec coordinates and to the redshift, respectively. Blue and red points are SF and passive galaxies, respectively. Colored crosses are the projected distribution of the cluster sample in the xz plane. The black cross indicates the center of the galaxy cluster.

The mass and metallicity profiles are presented in Fig. 12. The top panel shows the mass profile of AC114. Just as in Fig. 11, we observe a high concentration of galaxies toward the cluster center, which is dominated by passive galaxies. A negative radial mass gradient is detected. Black diamonds are the mean stellar mass of the galaxies in bins of ∼23 Mpc. The slope of the mass profile of −0.005 ± 0.002 (Spearman rank is −0.57, p-value = 0.005) shows a weak but statistically significant trend to find massive galaxies in the inner regions of the cluster. This is consistent with AC114 being in an early stage of radial relaxation, as is also suggested by the large velocity dispersion estimated in P15. Furthermore, passive galaxies present a slightly flatter slope (−0.004 ± 0.003, red linear fit) compared with SF galaxies (−0.007 ± 0.003, blue linear fit), which is indicative of ongoing dynamical activity in the galaxy cluster. The bottom panel of Fig. 12 presents a weak, non-statistically significant (Spearman rank of −0.38 with a p-value = 0.11) metallicity gradient of the SF galaxies, with a slope of −0.002  ±  0.002.

thumbnail Fig. 12.

Radial profiles of AC114. Top panel: stellar masses of galaxies as a function of the clustercentric radius. Blue and red crosses are SF and passive galaxies. Black diamonds represent mean mass bins of ∼23 h−1 Mpc. Black, blue, and red lines are the linear fits that show a negative mass-gradient as a function of the clustercentric radius for the complete sample, SF, and passive galaxies, respectively. Bottom panel: metallicity profile of the SF galaxies with available mean metallicity estimations. Black diamonds are mean metallicities in bins of ∼33 h−1 Mpc, and the black line is the linear fit that shows a weak negative gradient.

We applied a kernel density estimator (KDE) to the redshift distribution of the galaxy cluster to search for possible substructures and construct the MZR of the sub-components. We made use of the machine-learning KernelDensity algorithm from the Scikit-learn package6 and used a linear kernel model. Three substructures were identified in the redshift distribution with this procedure. The red lines in the top panel of Fig. 13 show the redshift limits of the substructures, namely SS1 (0.2850 ≤ z < 0.3029) being the front part of the galaxy cluster with three SF galaxies, SS2 (0.3029 ≤ z < 0.3243) for the main component of AC114 with ten SF galaxies, and SS3 (0.3243 ≤ z < 0.3380) being the back part of the galaxy cluster with six SF galaxies.

thumbnail Fig. 13.

Top panel: kernel density estimator in the redshift distribution of the galaxy cluster region: 0.28 ≤ z ≤ 0.34. Red lines show the redshift cuts of the three substructures identified: SS1 in 0.2850 ≤ z < 0.3029, SS2 with 0.3029 ≤ z < 0.3243, and SS3 covering 0.3243 ≤ z < 0.3380. Vertical black segments are the clustercentric radial distance of the galaxies. Bottom panel: galaxies of the three substructures in the mass-metallicity plane. Sky-blue triangles, black dots, and red dots represent SS1, SS2, and SS3, respectively. The black line is the cMZR.

Once the substructures were identified, we went ahead to construct the respective MZRs. The bottom panel of Fig. 13 shows the SF galaxies of the substructures in the mass-metallicity plane. Sky blue, black, and red are colors for SS1, SS2, and SS3, respectively. From SS1, only three SF galaxies were identified, grouped in the same region of the mass-metallicity plane. A statistically significant positive correlation is detected for SS3 (slope of 0.26 ± 0.14, Spearman rank r = 0.860, and p-value = 0.019, black dots). Furthermore, most of the galaxies in the sub-sample lie below the cMZR, suggesting that outskirt galaxies of AC114 are slowly evolving due to external processes that are suppressing the chemical evolution. On the other hand, the main component of AC114 (SS2) shows more scatter (0.28 and 0.33 for SS3 and SS2, respectively) in the mass-metallicity plane (red dots). The slope of the fit of SS2 is 0.19 ± 0.13, and no statistically significant trend is seen (Spearman rank of = 0.54, p-value = 0.13). The central region of AC114 consists of ten SF galaxies and 42 passive galaxies. Therefore, SF galaxies closer to the central region of AC114 are suffering not only from interaction with the ICM, but also more likely to be suffering from interaction with their neighbors.

In the previous subsection we were able to estimate clustercentric radial distances, so we repeated the same KDE procedure but applying it to the clustercentric radial distribution of AC114 galaxies. We divided the galaxies according to the following radial ranges (blue lines in the top panel of Fig. 14): the central region at r ≤ 32.5 h−1Mpc with four SF galaxies, the intermediate region between 32.5 ≤ r < 57.4 Mpc with seven SF galaxies, and the outskirts of AC114 in r > 57.4 Mpc with eight SF galaxies.

thumbnail Fig. 14.

Same as Fig. 13, but in this instance the KDE was applied in the radial distribution. Top panel: kernel density estimator in the radial distribution of the galaxy cluster. Blue lines show radial cuts at r = 32.5, r = 57.4, and r = 140 Mpc. The vertical black segments are the clustercentric radial distances of the galaxies. Bottom panel: galaxies in the mass-metallicity plane of the three radial regions. Green triangles, black dots, and blue dots represent the central, intermediate, and outskirt regions, respectively.

The bottom panel of Fig. 14 shows the SF galaxies of AC114 in the mass-metallicity plane colored with green, black, and blue for the central, intermediate, and external regions, respectively. Despite the central region only having four galaxies, two of them are the metal-richer ones of the galaxy cluster. On the other hand, the intermediate and outskirt regions of AC114 show an increase in oxygen abundances with the increase in the stellar mass, with a slope of 0.35 ± 0.08, Spearman rank r = 0.86, and p-value = 0.014, and a slope of 0.162 ± 0.09, Spearman rank r = 0.40, and p-value = 0.319, respectively.

5. Discussion

P15 and S23 studied the dynamical state of AC114 and the stellar populations of the cluster members. The radial profiles in Fig. 12 show a negative mass gradient as a function of the clustercentric radius, consistent with the cluster being in an early stage of radial relaxation, as is also suggested by P15. The steeper mass gradient up to r ∼ 130 Mpc of SF members suggests that these galaxies are more sensitive to the ongoing dynamical processes of AC114 because of their lower masses with respect to the passive ones. Both the cMZR and the bMZR shown in Fig. 9 are downshifted with respect to the local Universe MZR, showing that AC114 SF members are more metal-poor than the field galaxies at the same redshift, and similar in terms of metallicity to background galaxies at a given mass.

S23, studying stellar populations of SF and passive galaxies, have shown that the stellar MZR of AC114 is steeper than for galaxies in the local Universe. However, the gas-phase cMZR in Fig. 9 has a shallower slope than galaxies in the local Universe. This could be due to galaxies quenching by (i) strangulation (Larson et al. 1980; Peng et al. 2015), where the gas is consumed faster when there is no fresh supply, and (ii) the ICM is boosting the metal-poor inflows. By adding our analysis to that reported in S23, stellar metallicities are lower than gas-phase metallicities because those metals were produced in the past, while the gas-phase metallicities are present-day ones. Therefore, the SF galaxies should be in an early stage of the quenching.

When cluster members were separated into substructures, important imprints of environmental effects were revealed. Figures 13 and 14 shows that galaxies located outside the central region set the slope of the cMZR (red and black dots in Fig. 13, and blue and black dots in Fig. 14).

The MZRs of SF galaxies in massive clusters (> 1014M) appear to follow those of field galaxies at their corresponding redshift. For example, Vaughan et al. (2020), studying four CLASH clusters, namely MACS2129, MACS1311, MACS1931, and MS2131 (M200 ∼ 1014 − 1015M; Tiley et al. 2020), in the redshift range z ∼ 0.2 − 0.6, indicate that in these clusters the MZRs do not show any changes with respect to that constructed with field SF galaxies. Additionally, Lara-López et al. (2022) and Gupta et al. (2016), by using SF members of the Fornax cluster (Mvir ∼ 7 × 1013M; Drinkwater et al. 2001) and the RX J1532+30 (z ∼ 0.35, Mvir ∼ 6.4 × 1014M) galaxy cluster, respectively, are also consistent, with no variations with respect to the field MZR. However, these works report flatter galactic metallicity gradients, negative cluster metallicity gradients, and a negative trend in the residuals of the MZRs as a function of the clustercentric radius. Several physical processes are discussed, but the main physical phenomena that drive those trends are strangulation and RPS.

It is well known that RPS is the physical result of the galaxy-ICM interaction that affects the outskirts of galaxy disks: in the presence of negative metallicity gradients, the stripped-out gas is metal-poorer than their nuclei, which leads to increased oxygen abundances.

Strangulation is produced when the pristine gas inflow ends and the gas reservoir is consumed faster, increasing the SFR until the quenching begins. Therefore, the Hα disks become smaller than the stellar disks, and oxygen abundances are more metal-rich than galaxies that do not experience galaxy-ICM interactions at the same mass.

On the other hand, galaxy clusters such as Coma, Hercules (A2151), MACSJ1115+01, and those in the SDSS DR4 survey, show enhanced metallicities with respect to the field MZR, up to ∼0.05 dex in galaxies affected by RPS and/or galaxy-galaxy interactions inside the virial radius of their host dark matter halos (Ellison et al. 2009; Petropoulou et al. 2011, 2012; Gupta et al. 2016). This indicates that extremely dense environments, such as galaxy clusters with Mvir > 1015M, may show variations in the cluster MZRs with respect to the field ones.

The MZR of AC114 also shows this feature for at least two cluster members when it is compared with the field MZR at z ∼ 0.3, as is indicated in Fig. 10. In this case, we used the Maiolino parameterization at redshift 0.3 to compare with the AC114-MZR. However, the two galaxies are more metal-rich, ∼0.10 dex, which is more than that reported by Ellison et al. (2009). Additionally, the other galaxies (below the linear fit of the cMZR with the red line) drive the shallower slope and have lower metallicities than the field SF galaxies at the same redshift up to ∼0.22 dex. The striking feature is that these galaxies are located outside the central region of AC114 (see Fig. 14). Several reasons can explain this feature: (i) the AGN feedback, which is important at low z, can expel the metal-rich gas component in the central region of SF members with negative metallicity gradients (De Rossi et al. 2017; Wang et al. 2023; ii) the presence of strong metal-poor inflows (Finlator & Davé 2008; Ceverino et al. 2016), where galaxies accreting metal-poor gas decrease their average metallicities; and (iii) mergers and tidal interactions, where the interaction with closer companions can dilute the gas-phase metallicities.

In order to understand the main physical phenomena behind the metal-poor content of the SF galaxies of AC114 with respect to the field MZR at the same redshift, the AGN feedback can be discarded because our cluster sample has just two composite galaxies. One of them is over the local Universe MZR, and the other shows the largest deviation below the MZR in the cluster sample. According to Kauffmann et al. (2003), the location of these two composite galaxies in the mass-metallicity plane can be explained as the former being a LINER-composite galaxy and the latter a Seyfert-composite galaxy. Additionally, Seyfert galaxies have stronger [OIII]λ5007 emission lines than SF galaxies. The O3N2 index anticorrelates linearly with the oxygen abundances, so it may explain the location of this galaxy in the mass-metallicity plane. The same applies to the LINER galaxy in our cluster sample, but in this case, LINERS have strong [NII]λ6584 emission lines. According to observation and simulations from Zahid et al. (2014) and Wang et al. (2023), despite the AGN feedback being important at low redshift, for satellite galaxies (non-cDs in galaxy clusters) this effect is negligible.

On the other hand, a strong metal-poor inflow and dynamical events can explain the observed feature. Montuori et al. (2010), simulating mergers and flybys between disk galaxies, showed that galaxy-galaxy interactions can dilute gas-phase metallicities and also increase the SFR, due to the accretion of metal-poor gas from galaxies with negative metallicity gradients. The estimated dilution agrees with Kewley et al. (2006), which reports that galaxy pair interaction dilutes gas-phase metallicities up to ∼0.2 − 0.3 dex. Those values are also in line with the downshifts of the outlier foreground galaxies with respect to the fMZR, which are ∼0.18 on average. Additionally, Montuori et al. (2010) also explain that galaxy-galaxy interaction increases the gas-phase metallicities after the dilution due to the increase in the SFRs. Therefore, this mechanism increases the scatter of the MZR, as is also seen in Rupke et al. (2008).

This scenario is also supported by ongoing dynamical activity (P15) and by the fact that this galaxy cluster has an extended soft X-ray tail, together with the hardness ratio map, which suggests that this cluster has an ongoing merger activity with a smaller structure (De Filippis et al. 2004).

6. Summary and conclusions

We explored the role of galaxy clusters in the MZR through SF galaxies of the AC114 galaxy cluster. We analyzed optical, intermediate-resolution VIMOS/VLT spectra for a sample of 182 galaxies in the FOV of AC114.

The search in photometric catalogues allowed us to get optical and NIR luminosities for 76% and 66% of the galaxies of the sample from the DECaLS DR10 and VIKING DR4 ESO surveys, respectively. Synthetic photometry was performed using the VIMOS spectra when optical photometry was missing. Optical and NIR-based stellar masses of the galaxies were estimated via theoretical stellar mass-to-light ratios computed from the E-MILES models by considering the age dependence for both sdss−i and Ks band passes. A comparison between optical and NIR-based mass estimations shows agreement between them for SF galaxies, but there is a slight deviation in massive, passive galaxies of 0.23 dex. However, the offset is within the propagated uncertainties of σ = 0.29 mag, so it does not change our results.

We classified the galaxies into ELGs and passive galaxies, based on the presence or absence of [OII]λ3727, [OIII]λλ4959, 5007, Hβ, and Hα emission lines, respectively. Using the BPT diagram, together with the Kewley et al. (2001) and Kauffmann et al. (2003) demarcation limits, we found that the galaxy sample has seven composite galaxies, where two of them (a Seyfert composite and a LINER composite) are cluster members. The remaining galaxies in the sample are, in fact, pure SF galaxies.

Gas-phase metallicities were estimated from the O3N2, which are calibrated on Te-based metallicities.

Our main conclusions can be summarized as follows:

  • Foreground and background MZRs span masses from 108.7 to 1010.3M and 109.4 to 1010.7M, respectively, with both the fMZR and the bMZR showing a scatter of ∼0.12 dex. The fMZR was constructed with less-massive galaxies than bMZR due to an observational bias. However, the downshifts reported show a coherent evolution in redshift of the MZR, where foreground galaxies are 0.09 dex metal-poorer, and background galaxies are 0.20 dex metal-poorer on average with respect to the local Universe galaxies.

  • The AC114-MZR presents different features with respect to field galaxies at similar redshifts. The cMZR span masses from 109.2 to 1010.8M, with 0.17 dex of scatter. A downshift of 0.19 dex on average and a slightly shallower slope is seen with respect to both the fMZR and the bMZR. This shows that, in general, AC114 SF galaxies have lower metallicities at the same mass than field galaxies at the same redshift. The metallicity measurements of the SF cluster members are comparable to the background galaxies in our galaxy sample. Besides, star-forming galaxies within the main component of AC114 and the central region do not show a correlation between mass and gas-phase metallicity, with a large scatter. This gives important clues about the environmental effects that are suppressing the chemical evolution of these galaxies.

  • The mass and metallicity profiles as a function of the clustercentric radius show that AC114 started its radial relaxation phase, in agreement with Proust et al. (2015). The negative mass gradient suggests that SF galaxies are the most sensitive to the dynamical state of the galaxy cluster. Additionally, the negative (although not statistically significant) metallicity gradient can be a sign that massive galaxies located near the central region of AC114 are slightly more metal-rich than those found in external regions.

  • The substructures identified with the KDE show that for SF galaxies in the outskirts of the galaxy cluster the metallicity is dependent on the stellar mass. Moreover, these galaxies are responsible for the shallower slope of the cMZR (compared with the field MZR at z ∼ 0.3), which means that these galaxies are also the more metal-poor ones.

  • When the cMZR is compared with the field MZR at z ∼ 0.3, two galaxies are located above the field MZR by ∼0.10 dex. This feature is consistent with the literature, which suggests that galaxies are suffering from RPS, strangulation, and/or galaxy-galaxy interactions.

  • Galaxies falling below the linear fit of the cMZR are the drivers of the shallower slope. These galaxies are also found outside the central region of AC114. This means that outskirt galaxies are evolving slowly in terms of their chemical content, showing metallicities up to 0.22 dex lower than the local Universe galaxies and on average 0.17 dex lower than the field MZR at z ∼ 0.3. Strong metal-poor inflows produced by dynamical events such as galaxy-galaxy interactions may explain the strikingly low metallicities of the SF members of AC114. This scenario is also supported by the dynamical state of the cluster and the ongoing merger event between the main component of AC114 and a smaller structure, observed in the X-rays.

  • We emphasize the role of the mass of a galaxy cluster. Environmental and dynamical effects drive the chemical evolution of members in a unique way with respect to field galaxies. However, only galaxy clusters with M > 1015M appear to show observable changes in the MZR.


3

Bagpipes is a state of the art Python code for modelling galaxy spectra and fitting spectroscopic and photometric observations. https://bagpipes.readthedocs.io/en/latest/

4

IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation.

Acknowledgments

We thank the anonymous referee for their useful comments which improved the quality of the paper. A special thanks goes to the Dr. Piercarlo Bonifacio for the feedback presented in the early stage of this research. Additionally, we want to thanks Felipe Vivanco and Gustavo Worner, for the opportunity to discuss our results and for the valuable contributions to the development of the Python scripts used in this work. This publication is based on observations collected at the European Southern Observatory under ESO programme ID 083.A-0566. Also, this program has made use of data from the VIKING survey from VISTA at the ESO Paranal Observatory, programme ID 179.A-2004. Data processing has been contributed by the VISTA Data Flow System at CASU, Cambridge and WFAU, Edinburgh.

References

  1. Alonso, M. S., Lambas, D. G., Tissera, P., et al. 2007, MNRAS, 375, 1017 [NASA ADS] [CrossRef] [Google Scholar]
  2. Aller, L. H. 1984, Astrophysics and Space Science Library (Dordrecht: Reidel), 1984 [CrossRef] [Google Scholar]
  3. Astropy-Specutils Development Team 2019, Astrophysics Source Code Library [record ascl:1902.012] [Google Scholar]
  4. Baldwin, J. A., Phillips, M. M., & Terlevich, R. 1981, PASP, 93, 5 [Google Scholar]
  5. Barton, E. J., Geller, M. J., & Kenyon, S. J. 2000, ApJ, 530, 660 [NASA ADS] [CrossRef] [Google Scholar]
  6. Bell, E. F., & de Jong, R. S. 2001, ApJ, 550, 212 [Google Scholar]
  7. Butcher, H., & Oemler, A. 1978, ApJ, 219, 18 [NASA ADS] [CrossRef] [Google Scholar]
  8. Cameron, A. J., Katz, H., & Rey, M. P. 2023, MNRAS, 522, L89 [NASA ADS] [CrossRef] [Google Scholar]
  9. Campusano, L. E., Pelló, R., Kneib, J.-P., et al. 2001, A&A, 378, 394 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245 [Google Scholar]
  11. Carnall, A. C., McLure, R. J., Dunlop, J. S., et al. 2019, MNRAS, 490, 417 [Google Scholar]
  12. Ceverino, D., Sánchez Almeida, J., Muñoz Tuñón, C., et al. 2016, MNRAS, 457, 2605 [NASA ADS] [CrossRef] [Google Scholar]
  13. Couch, W. J., & Sharples, R. M. 1987, MNRAS, 229, 423 [NASA ADS] [CrossRef] [Google Scholar]
  14. Couch, W. J., Barger, A. J., Smail, I., et al. 1998, ApJ, 497, 188 [NASA ADS] [CrossRef] [Google Scholar]
  15. Couch, W. J., Balogh, M. L., Bower, R. G., et al. 2001, ApJ, 549, 820 [NASA ADS] [CrossRef] [Google Scholar]
  16. Cowie, L. L., Songaila, A., Hu, E. M., et al. 1996, AJ, 112, 839 [NASA ADS] [CrossRef] [Google Scholar]
  17. Cresci, G., Mannucci, F., Sommariva, V., et al. 2012, MNRAS, 421, 262 [NASA ADS] [Google Scholar]
  18. De Filippis, E., Bautz, M. W., Sereno, M., et al. 2004, ApJ, 611, 164 [NASA ADS] [CrossRef] [Google Scholar]
  19. de Rossi, M. E., Tissera, P. B., & Scannapieco, C. 2007, MNRAS, 374, 323 [NASA ADS] [CrossRef] [Google Scholar]
  20. De Rossi, M. E., Bower, R. G., Font, A. S., et al. 2017, MNRAS, 472, 3354 [NASA ADS] [CrossRef] [Google Scholar]
  21. Dey, A., Schlegel, D. J., Lang, D., et al. 2019, AJ, 157, 168 [Google Scholar]
  22. Drinkwater, M. J., Gregg, M. D., & Colless, M. 2001, ApJ, 548, L139 [Google Scholar]
  23. Edmunds, M. G., & Pagel, B. E. J. 1984, MNRAS, 211, 507 [NASA ADS] [CrossRef] [Google Scholar]
  24. Ellison, S. L., Patton, D. R., Simard, L., et al. 2008, AJ, 135, 1877 [NASA ADS] [CrossRef] [Google Scholar]
  25. Ellison, S. L., Simard, L., Cowan, N. B., et al. 2009, MNRAS, 396, 1257 [NASA ADS] [CrossRef] [Google Scholar]
  26. Ellison, S. L., Viswanathan, A., Patton, D. R., et al. 2019, MNRAS, 487, 2491 [NASA ADS] [CrossRef] [Google Scholar]
  27. Erb, D. K., Shapley, A. E., Pettini, M., et al. 2006, ApJ, 644, 813 [NASA ADS] [CrossRef] [Google Scholar]
  28. Ferland, G. J., Korista, K. T., Verner, D. A., et al. 1998, PASP, 110, 761 [Google Scholar]
  29. Finlator, K., & Davé, R. 2008, MNRAS, 385, 2181 [NASA ADS] [CrossRef] [Google Scholar]
  30. Fouesneau, M. 2022, https://doi.org/10.5281/zenodo.7016775 [Google Scholar]
  31. Freudling, W., Romaniello, M., Bramich, D. M., et al. 2013, A&A, 559, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Fujita, Y. 2004, PASJ, 56, 29 [NASA ADS] [Google Scholar]
  33. Gao, Y., Bao, M., Yuan, Q., et al. 2018, ApJ, 869, 15 [NASA ADS] [CrossRef] [Google Scholar]
  34. Girardi, L., Bressan, A., Bertelli, G., et al. 2000, A&AS, 141, 371 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Gunn, J. E., & Gott, J. R. 1972, ApJ, 176, 1 [NASA ADS] [CrossRef] [Google Scholar]
  36. Gupta, A., Yuan, T., Tran, K.-V. H., et al. 2016, ApJ, 831, 104 [NASA ADS] [CrossRef] [Google Scholar]
  37. Hamuy, M., Walker, A. R., Suntzeff, N. B., et al. 1992, PASP, 104, 533 [NASA ADS] [CrossRef] [Google Scholar]
  38. Hamuy, M., Suntzeff, N. B., Heathcote, S. R., et al. 1994, PASP, 106, 566 [NASA ADS] [CrossRef] [Google Scholar]
  39. Hopkins, A. M., Miller, C. J., Nichol, R. C., et al. 2003, ApJ, 599, 971 [Google Scholar]
  40. Hummer, D. G., & Storey, P. J. 1987, MNRAS, 224, 801 [NASA ADS] [CrossRef] [Google Scholar]
  41. Kauffmann, G., Heckman, T. M., Tremonti, C., et al. 2003, MNRAS, 346, 1055 [Google Scholar]
  42. Kauffmann, G., White, S. D. M., Heckman, T. M., et al. 2004, MNRAS, 353, 713 [Google Scholar]
  43. Kewley, L. J., & Ellison, S. L. 2008, ApJ, 681, 1183 [Google Scholar]
  44. Kewley, L. J., Dopita, M. A., Sutherland, R. S., et al. 2001, ApJ, 556, 121 [Google Scholar]
  45. Kewley, L. J., Geller, M. J., & Barton, E. J. 2006, AJ, 131, 2004 [NASA ADS] [CrossRef] [Google Scholar]
  46. Kim, K. J., Bayliss, M. B., Noble, A. G., et al. 2023, ApJ, 955, 32 [CrossRef] [Google Scholar]
  47. Kobayashi, C., Springel, V., & White, S. D. M. 2007, MNRAS, 376, 1465 [Google Scholar]
  48. Kobulnicky, H. A., & Kewley, L. J. 2004, ApJ, 617, 240 [CrossRef] [Google Scholar]
  49. Kobulnicky, H. A., & Skillman, E. D. 1996, ApJ, 471, 211 [NASA ADS] [CrossRef] [Google Scholar]
  50. Lamareille, F., Brinchmann, J., Contini, T., et al. 2009, A&A, 495, 53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Lara-López, M. A., Galán-de Anta, P. M., Sarzi, M., et al. 2022, A&A, 660, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. Larson, R. B. 1974, MNRAS, 169, 229 [NASA ADS] [CrossRef] [Google Scholar]
  53. Larson, R. B., Tinsley, B. M., & Caldwell, C. N. 1980, ApJ, 237, 692 [Google Scholar]
  54. Le Fèvre, O., Saisse, M., Mancini, D., et al. 2003, Proc. SPIE, 4841, 1670 [NASA ADS] [CrossRef] [Google Scholar]
  55. Lequeux, J., Peimbert, M., Rayo, J. F., et al. 1979, A&A, 80, 155 [NASA ADS] [Google Scholar]
  56. Maiolino, R., & Mannucci, F. 2019, A&A Rev., 27, 3 [NASA ADS] [CrossRef] [Google Scholar]
  57. Maiolino, R., Nagao, T., Grazian, A., et al. 2008, A&A, 488, 463 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  58. Montuori, M., Di Matteo, P., Lehnert, M. D., et al. 2010, A&A, 518, A56 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  59. Moustakas, J., Zaritsky, D., Brown, M., et al. 2011, ArXiv e-prints [arXiv:1112.3300] [Google Scholar]
  60. Natarajan, P., Kneib, J.-P., Smail, I., et al. 1998, ApJ, 499, 600 [NASA ADS] [CrossRef] [Google Scholar]
  61. Neistein, E., van den Bosch, F. C., & Dekel, A. 2006, MNRAS, 372, 933 [NASA ADS] [CrossRef] [Google Scholar]
  62. O’Donnell, J. E. 1994, ApJ, 422, 158 [Google Scholar]
  63. Peimbert, M. 1967, ApJ, 150, 825 [NASA ADS] [CrossRef] [Google Scholar]
  64. Peng, Y., Maiolino, R., & Cochrane, R. 2015, Nature, 521, 192 [Google Scholar]
  65. Pérez-Montero, E. 2015, Highlights of Spanish Astrophysics VIII, 286 [Google Scholar]
  66. Pettini, M., & Pagel, B. E. J. 2004, MNRAS, 348, L59 [NASA ADS] [CrossRef] [Google Scholar]
  67. Petropoulou, V., Vílchez, J., Iglesias-Páramo, J., et al. 2011, ApJ, 734, 32 [NASA ADS] [CrossRef] [Google Scholar]
  68. Petropoulou, V., Vílchez, J., & Iglesias-Páramo, J. 2012, ApJ, 749, 133 [NASA ADS] [CrossRef] [Google Scholar]
  69. Pilyugin, L. S. 2001, A&A, 369, 594 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  70. Pilyugin, L. S., & Grebel, E. K. 2016, MNRAS, 457, 3678 [NASA ADS] [CrossRef] [Google Scholar]
  71. Press, W. H., & Schechter, P. 1974, ApJ, 187, 425 [Google Scholar]
  72. Proust, D., Yegorova, I., Saviane, I., et al. 2015, MNRAS, 452, 3304 [NASA ADS] [CrossRef] [Google Scholar]
  73. Ricciardelli, E., Vazdekis, A., Cenarro, A. J., et al. 2012, MNRAS, 424, 172 [NASA ADS] [CrossRef] [Google Scholar]
  74. Rupke, D. S. N., Veilleux, S., & Baker, A. J. 2008, ApJ, 674, 172 [NASA ADS] [CrossRef] [Google Scholar]
  75. Sánchez, S. F., Rosales-Ortega, F. F., Jungwiert, B., et al. 2013, A&A, 554, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  76. Savaglio, S., Glazebrook, K., Le Borgne, D., et al. 2005, ApJ, 635, 260 [NASA ADS] [CrossRef] [Google Scholar]
  77. Saviane, I., Yegorova, I., & Proust, D. 2023, MNRAS, 526, 2458 [NASA ADS] [CrossRef] [Google Scholar]
  78. Scannapieco, C., Tissera, P. B., White, S. D. M., et al. 2008, MNRAS, 389, 1137 [CrossRef] [Google Scholar]
  79. Sereno, M., Lubini, M., & Jetzer, P. 2010, A&A, 518, A55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  80. Smail, I., Couch, W. J., Ellis, R. S., et al. 1995, ApJ, 440, 501 [NASA ADS] [CrossRef] [Google Scholar]
  81. Spitoni, E., Calura, F., Mignoli, M., et al. 2020, A&A, 642, A113 [EDP Sciences] [Google Scholar]
  82. Springel, V., White, S. D. M., Jenkins, A., et al. 2005, Nature, 435, 629 [Google Scholar]
  83. Tiley, A. L., Vaughan, S. P., Stott, J. P., et al. 2020, MNRAS, 496, 649 [NASA ADS] [CrossRef] [Google Scholar]
  84. Tissera, P. B., De Rossi, M. E., & Scannapieco, C. 2005, MNRAS, 364, L38 [NASA ADS] [CrossRef] [Google Scholar]
  85. Tonry, J., & Davis, M. 1979, AJ, 84, 1511 [Google Scholar]
  86. Tremonti, C. A., Heckman, T. M., Kauffmann, G., et al. 2004, ApJ, 613, 898 [Google Scholar]
  87. Troncoso, P., Maiolino, R., Sommariva, V., et al. 2014, A&A, 563, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  88. Vaughan, S. P., Tiley, A. L., Davies, R. L., et al. 2020, MNRAS, 496, 3841 [NASA ADS] [CrossRef] [Google Scholar]
  89. Vazdekis, A., Ricciardelli, E., Cenarro, A. J., et al. 2012, MNRAS, 424, 157 [Google Scholar]
  90. Vazdekis, A., Koleva, M., Ricciardelli, E., et al. 2016, MNRAS, 463, 3409 [NASA ADS] [CrossRef] [Google Scholar]
  91. Wang, K., Wang, X., & Chen, Y. 2023, ApJ, 951, 66 [NASA ADS] [CrossRef] [Google Scholar]
  92. Wright, E. L. 2006, PASP, 118, 1711 [NASA ADS] [CrossRef] [Google Scholar]
  93. Zahid, H. J., Dima, G. I., Kudritzki, R.-P., et al. 2014, ApJ, 791, 130 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: Tables:

Tables A.2, A.3, and A.4 are only available at the CDS.

Table A.1.

Compilation of all statistical parameters of the fits performed in Figures 9, 12, 13, and 14. First column: Corresponding figure. Second column: Fitted structure (e.g., the cluster MZR, or the SF mass profile). Third column: Slope of the linear fit. Fourth column: 1σ scatter of the distribution. Fifth column: Spearman correlation test. Sixth column: p−value of the Spearman correlation test.

All Tables

Table A.1.

Compilation of all statistical parameters of the fits performed in Figures 9, 12, 13, and 14. First column: Corresponding figure. Second column: Fitted structure (e.g., the cluster MZR, or the SF mass profile). Third column: Slope of the linear fit. Fourth column: 1σ scatter of the distribution. Fifth column: Spearman correlation test. Sixth column: p−value of the Spearman correlation test.

All Figures

thumbnail Fig. 1.

Observed galaxies (colored circles) on top of an image taken from the DECaLS DR10 viewer for the AC114 cluster. Quadrants Q1, Q2, Q3, and Q4 are counterclockwise, with Q4 in the bottom right clump.

In the text
thumbnail Fig. 2.

Reduced spectra of a passive galaxy (top panel) and an ELG (bottom panel). Red lines indicate absorption features, whereas blue lines indicate emission features.

In the text
thumbnail Fig. 3.

Top panel: distribution of photometric data collected from DECaLS DR10 survey. The blue line indicates the distribution of g-band magnitudes and the orange line indicates the distribution of i-band magnitudes. Bottom panel: distribution of photometric data collected by the VIKING DR4 survey. The blue line indicates the distribution of H-band magnitudes and the orange line indicates the distribution of Ks-band magnitudes.

In the text
thumbnail Fig. 4.

Synthetic spectra reproduced with the Bagpipes models. The blue spectrum represents a passive galaxy. The orange spectrum represents an SF galaxy.

In the text
thumbnail Fig. 5.

Redshift distribution of the galaxy sample. Upper panel: distribution of passive galaxies (red line) and ELGs (blue line). Lower panel: redshift distribution of the complete sample.

In the text
thumbnail Fig. 6.

Comparison between DECaLS DR10 (on the x axis) and those from artificial magnitudes (on the y axis). Colors are the same as in Fig. 5. The gray line is an identity line. Top panel: comparison between synthetic and DECaLS magnitudes in the i band. Bottom panel: comparison between synthetic and DECaLS magnitudes in the i band with the applied correction.

In the text
thumbnail Fig. 7.

Central panel: comparison between i and Ks-based masses. Blue and red dots represent ELGs and passive galaxies, respectively. The gray line corresponds to the 1:1 line ratio. Upper panel: distribution of optical luminosity based stellar masses. The black line shows the mass distribution of the complete sample, whereas the blue and red areas indicate the distribution of ELGs and passive galaxies, respectively. Right panel: distribution of near-infrared luminosity based stellar masses.

In the text
thumbnail Fig. 8.

Baldwin-Philips-Terlevich diagram of the sample. The dashed black line is the KL at z = 0 (Kewley et al. 2001), which excludes pure AGNs. The solid black line is the Kauffmann et al. (2003) limit, which excludes composite galaxies (ionization from star formation and AGNs). The red points are galaxies at the cluster redshift. Black diamonds correspond to foreground and background galaxies. According to both demarcation limits, all ELGs are SF galaxies. Data points with blue open diamonds are galaxies classified as composite galaxies, with two of them being composite galaxies at the cluster redshift and five of them found outside the cluster redshift.

In the text
thumbnail Fig. 9.

Mass-metallicity relations constructed with the SF galaxies of the sample compared with the local Universe MZR from Tremonti et al. (2004). The left panel shows the foreground MZR (fMZR) with blue dots and background MZR (bMZR) with red dots. Lines correspond to linear fits. The shaded area corresponds to the 1σ scatter for the fMZR and bMZR. Right panel: MZR with the SF members of AC114 (cMZR). The black line and shaded area are the linear fit and the respective 1σ scatter. The dashed blue and red lines are the fMZR and bMZR linear fits. The open circles (left) and open circles (right) show galaxies whose stellar mass was estimated using synthetic photometry.

In the text
thumbnail Fig. 10.

Comparison between the MZRs constructed in this work and the literature. Left panel: comparison between the fMZR, Cresci et al. (2012) MZR, and the Maiolino parameterization, both at redshift z ∼ 0.25. Middle panel: comparison between the bMZR, Savaglio et al. (2005), Lamareille et al. (2009), and the Maiolino parameterization at z ∼ 0.7. Right panel: comparison between the cMZR, Vaughan et al. (2020), and the Maiolino parameterization at z ∼ 0.3. Data points with open diamonds show the galaxies classified as composite in the BPT diagram of Fig. 8.

In the text
thumbnail Fig. 11.

Three-dimensional projection of the AC114 members, expressed in Mpc with respect to the cluster center. The x, y, and z axes correspond to the RA,Dec coordinates and to the redshift, respectively. Blue and red points are SF and passive galaxies, respectively. Colored crosses are the projected distribution of the cluster sample in the xz plane. The black cross indicates the center of the galaxy cluster.

In the text
thumbnail Fig. 12.

Radial profiles of AC114. Top panel: stellar masses of galaxies as a function of the clustercentric radius. Blue and red crosses are SF and passive galaxies. Black diamonds represent mean mass bins of ∼23 h−1 Mpc. Black, blue, and red lines are the linear fits that show a negative mass-gradient as a function of the clustercentric radius for the complete sample, SF, and passive galaxies, respectively. Bottom panel: metallicity profile of the SF galaxies with available mean metallicity estimations. Black diamonds are mean metallicities in bins of ∼33 h−1 Mpc, and the black line is the linear fit that shows a weak negative gradient.

In the text
thumbnail Fig. 13.

Top panel: kernel density estimator in the redshift distribution of the galaxy cluster region: 0.28 ≤ z ≤ 0.34. Red lines show the redshift cuts of the three substructures identified: SS1 in 0.2850 ≤ z < 0.3029, SS2 with 0.3029 ≤ z < 0.3243, and SS3 covering 0.3243 ≤ z < 0.3380. Vertical black segments are the clustercentric radial distance of the galaxies. Bottom panel: galaxies of the three substructures in the mass-metallicity plane. Sky-blue triangles, black dots, and red dots represent SS1, SS2, and SS3, respectively. The black line is the cMZR.

In the text
thumbnail Fig. 14.

Same as Fig. 13, but in this instance the KDE was applied in the radial distribution. Top panel: kernel density estimator in the radial distribution of the galaxy cluster. Blue lines show radial cuts at r = 32.5, r = 57.4, and r = 140 Mpc. The vertical black segments are the clustercentric radial distances of the galaxies. Bottom panel: galaxies in the mass-metallicity plane of the three radial regions. Green triangles, black dots, and blue dots represent the central, intermediate, and outskirt regions, respectively.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.