# Environmental Effects on the Metal Enrichment of Low Mass Galaxies in Nearby Clusters

Petropoulou V.11affiliationmark: ,Vílchez J.11affiliationmark: , Iglesias-Páramo J.11affiliationmark: 22affiliationmark: 11affiliationmark: Instituto de Astrofísica de Andalucía- C.S.I.C., Glorieta de la Astronomía, 18008 Granada, Spain;
22affiliationmark: Centro Astronómico Hispano Alemán, C/ Jesús Durbán Remón 2-2, 04004 Almería, Spain;
###### Abstract

In this paper we study the chemical history of low-mass star-forming (SF) galaxies in the local Universe clusters Coma, A1367, A779, and A634. The aim of this work is to search for the imprint of the environment on the chemical evolution of these galaxies. Galaxy chemical evolution is linked to the star formation history (SFH), as well as to the gas interchange with the environment, and low-mass galaxies are well known to be vulnerable systems to environmental processes affecting both these parameters. For our study we have used spectra from the SDSS-III DR8. We have examined the spectroscopic properties of SF galaxies of stellar masses , located from the core to the cluster’s outskirts. The gas-phase O/H and N/O chemical abundances have been derived using the latest empirical calibrations. We have examined the mass-metallicity relation of cluster galaxies finding well defined sequences. The slope of these sequences, for galaxies in low-mass clusters and galaxies at large cluster-centric distances, follows the predictions of recent hydrodynamic models. A flattening of this slope has been observed for galaxies located in the core of the two more massive clusters of the sample, principally in Coma, suggesting that the imprint of the cluster environment on the chemical evolution of SF galaxies should be sensitive to both the galaxy mass and the host cluster mass. The HI gas content of Coma and A1367 galaxies indicate that low-mass SF galaxies, located at the core of these clusters, have been severely affected by ram-pressure stripping. The observed mass-dependent enhancement of the metal content of low-mass galaxies in dense environments seems plausible, according to hydrodynamic simulations. This enhanced metal enrichment could be produced by the combination of effects such as wind reaccretion, due to pressure cofinement by the intra-cluster medium (ICM), and the truncation of gas infall, as a result of the ram-pressure stripping. Thus, the properties of the ICM should play an important role in the chemical evolution of low-mass galaxies in clusters.

galaxy clusters: general — galaxy clusters: individual(Coma, Abell 1656, Abell 1367, Abell 779, Abell 634)
slugcomment: Not to appear in Nonlearned J., 45.

## 1 Introduction

The chemical evolution of a galaxy is linked to the SFH, as well as to the gas interchange with the environment via inflows or outflows. In this sense, galaxy metallicity, could be an observable parameter providing information on the impact of the environment on the galaxy SFH and/or the galaxy gas content.

SFH has been found to be related to the environment in which a galaxy inhabits: dense regions such as galaxy clusters, are mainly populated by red passive spheroids, whereas SF galaxies dominate in the field (Dressler, 1980; Balogh et al., 1999; Poggianti et al., 1999; Treu et al., 2003; Balogh et al., 2004; Rines et al., 2005; Finn et al., 2005; Haines et al., 2007; Poggianti et al., 2008). Though, intense discussions on the underlying cause of this observational trend are still ongoing. One interpretation has been that this trend could be the mere result of the fact that high density regions have favored the formation mostly of massive galaxies, which in turn convert their gas in stars faster than dwarf galaxies (Kennicutt, 1998) and they have already finished forming stars earlier in the past (Merlin & Chiosi, 2006). Recent observational evidences support the idea that the stellar mass function can be associated to the environment (Vulcani et al., 2011; Bolzonella et al., 2010).

Another reason could be that the cluster environment causes galaxies to transform their properties, as they move from low density regions into the cluster core. A number of plausible mechanisms have been proposed, including interactions with the ICM, interactions with the cluster gravitational potential and small scale galaxy-galaxy interactions (see Treu et al., 2003, and references therein). Each one of these mechanisms is expected to be effective in different regions of clusters and their outskirts, and to affect star formation in different ways and time-scales (see e.g. the review by Boselli & Gavazzi, 2006).

For low-mass galaxies, however, the picture appears much clearer. Due to their low mass surface density and rotation velocity, dwarf galaxies are the most vulnerable systems to environment-related processes that quench star formation, such as the stripping of their halo and/or disk gas. Haines et al. (2007), based on a large sample of dwarf galaxies, have found that passive dwarfs are found only in very dense environments, typical of a cluster virial region, or as satellites of a more massive companion. Studies of three of the richest clusters in the local Universe, A2199 (Haines et al., 2006), Coma and A1367 (Mahajan et al., 2010), suggest that star formation in dwarf galaxies is quenched only in the center of the clusters, in contrast with massive galaxies that can become passive in all environments. All these evidences seem to imply that there exist fundamental differences in the evolution of giant and dwarf galaxies in clusters: while various physical mechanisms could be co-responsible for the evolution of massive galaxies, the evolution of dwarf galaxies seems to be primarily driven by the environment in which they are found.

In accordance to the current cosmological paradigm (e.g Springel et al., 2005), clusters at have been found to accreate late-type galaxies along the filamentary structures that compose the cosmic web (e.g. Smith et al., 2011). Mahajan et al. (2011) have found that post-starburst dwarfs are preferentially located in the infall regions of Coma and A1367, suggesting that these galaxies experience a sudden quenching of star formation due to the interaction with the ICM. Moreover, before the gas gets totally stripped off, dwarf galaxies could experience an enhancement of star formation, either during their infall into the cluster along the filaments (Mahajan et al., 2010; Porter et al., 2008), or in the first stages of their encounter with the hot ICM (pressure triggered star formation, Abramson et al., 2011; Sun et al., 2007; Gavazzi et al., 2001). Thus SF dwarf galaxies are excellent probes to test the influence of the environment on the process of star formation and galaxy evolution.

The gas interchange of a galaxy with the environment is observationally well confirmed. For example, the ram-pressure stripping (RPS, Gunn & Gott, 1972) can remove the interstellar medium (ISM) of a galaxy as it moves through the hot intracluster gas. Evidences of ongoing gas stripping, both in low and high luminosity galaxies, have turned out to be frequent in local Universe clusters, such as the Virgo cluster (Fumagalli et al., 2011; Abramson et al., 2011; Kenney & Koopmann, 1999), A1367 (Gavazzi et al., 2001), and A3627 (Sun et al., 2007). In the Coma cluster core several galaxies have been observed showing clear signs of gas stripping, such as H or ultra-violet tails (Yagi et al., 2010; Smith et al., 2010). Additional observational evidences exist confirming the gas stripping in clusters, such as truncated ionized gas disks (Koopmann & Kenney, 2004; Koopmann et al., 2006; Cedrés et al., 2009; Jaffé et al., 2011) and disturbed gas kinematics (Jaffé et al., 2011). Moreover, the atomic gas properties of cluster galaxies seem to tell the same story (Haynes et al., 1984; Solanes et al., 2001; Chung et al., 2009; Levy et al., 2007).

Simulations (e.g. Abadi et al., 1999; Mori & Burkert, 2000; Tonnesen et al., 2007) suggest that RPS is efficient even up to the cluster virial radius and can remove completely the gas content in time-scales of yr, comparable to the cluster crossing time. Additionally, RPS has been found to be a multi-stage process (Roediger & Hensler, 2005) and recent simulations indicate that both varying ICM density (Brüggen & De Lucia, 2008; Tecce et al., 2010) and inhomogeneities in the ISM (Tonnesen & Bryan, 2009) are expected to play an important role in the gas stripping effect. Tecce et al. (2011) report that the gradient of ram pressure becomes steeper with increasing cluster mass; in the massive clusters (), the ram pressure at the core is times higher than at 111 is the radius from the cluster center which defines a sphere with interior mean density 200 times the critical density and is used as an approximation of virial radius., whereas in galaxy group-sized haloes the ram pressure at is of the central value. Bekki (2009) simulations suggest that even moderately strong ram pressure, e.g. in clusters of or even in groups, could strip the hot gas halos of galaxies, an effect known as strangulation. In these simulations, the stripping of galactic halo gas seems to be more efficient than that of disc gas, and the efficiency increases with decreasing galaxy mass.

Gas infall has also been proposed to take place in both low-mass irregular and high-mass spiral galaxies, in order to explain broad-band colors, gas fractions, SFRs, and metallicities at low and high redshifts (see Dalcanton, 2007, and references therein). Additionally, interacting galaxies can undergo nuclear metal dilution due to gas inflows (Kewley et al., 2010; Montuori et al., 2010; Rupke et al., 2010; Michel-Dansac et al., 2008), resulting in the flattening of the gas-phase metallicity gradients and altering their positions on the mass-metallicity relation (MZR).

All these environmental effects described, on both the SFH and the gas content of a galaxy, would be expected to leave their imprint on galaxy metallicity. Even the most recent events of star formation can enrich the gas with metals, so gas metallicities should be very sensitive to trace ongoing changes on the SFH and gas content (Ellison et al., 2009). SF galaxies, through their ionized gas emission, provide this observable information on nebular metallicity. Thus, studying the chemical abundances of SF galaxies in clusters can help appreciate the effect of the environment on the chemical evolution of cluster galaxies.

Previous works (Mouhcine et al., 2007; Cooper et al., 2008; Ellison et al., 2009) on the metallicity of SF galaxies as a function of the environment have not provided conclusive results, as small variations of dex have been derived. These works however, have been limited to the high stellar mass range, where there are not conclusive evidences whether the cluster environment could change significantly the SFH. Additionally, they have focused on the effect of local galaxy density, which possibly is not the only relevant parameter. Highlighting the cluster environment impact, HI deficient spirals in Virgo and Pegasus I clusters, have been found to be more chemically enriched than normal spirals (Skillman et al., 1996; Robertson et al., 2011), while this correlation between HI deficiency and metallicity has not been observed for field galaxies (Robertson et al., 2011).

The gas metallicity of SF dwarf galaxies in local Universe clusters has also been addressed. In Virgo dIrr, the gas metallicity has not been found to show a clear trend with the environment (Vilchez, 1995; Lee et al., 2003; Vaduvescu et al., 2007); however, Lee et al. (2003) have found that five of these dIrr are gas deficient with respect to field dIrr at comparable oxygen abundances, and this gas deficiency correlates with the X-Ray surface brightness of the ICM. Some SF dwarfs in Hydra cluster (Duc et al., 2001) have been found metal-rich for their luminosities and Vaduvescu et al. (2011), comparing the MZR of SF dwarfs in Hydra, Fornax and Virgo, have suggested that differences in the MZR seem to exist for galaxies in such different environments. Going to a more massive cluster (A2151), Petropoulou et al. (2011) have been able to observe dwarf galaxies, located in the cluster core, showing higher metallicities for their mass.

In the present study we extend these previous works, investigating the chemical history of SF dwarf galaxies in four clusters in the local Universe: Coma (A1656), A1367, A779, and A634. Coma and A1367 are among the nearest very rich galaxy clusters, and both belong to the large structure called the Coma supercluster. Recent insightful evidences on the assembly history and the SFH of the low-mass galaxy population in Coma supercluster (Smith et al., 2011; Mahajan et al., 2011, 2010) seem to indicate that it provides exemplary conditions for searching the potential imprints of the cluster environment on galaxy chemical enrichment. To investigate whether the mass of the host cluster could play a significant role in the chemical evolution of SF dwarf galaxies, we include to our sample two lower-mass clusters: A779 and A634. In total, these clusters are the four Abell clusters with SDSS spectroscopic data, that fulfill the criteria: i) to be visible from the northern hemisphere ( -25 deg), and ii) to be located at the same distance at Mpc (). Thus, the sample of clusters studied in the present work belong to a semi-spheric shell of the local Universe222In this semi-spheric shell are found three additional Abell clusters: A400, A539 and A2666, but they do not have SDSS spectroscopic data. Multifiber spectroscopic observations for these clusters will be presented in a forthcoming work (Petropoulou et al.,in prep)., and span a wide range of halo masses from to .

We have used the latest spectroscopic release SDSS-III DR8, where new emission line measures have been provided, after correcting the spectra for the underlying stellar population. This is an important issue when deriving gas-phase metallicities, reducing biases that previous works could have been suffering. Additionally, three out of four of the present sample clusters were not included in previous releases than DR7. Anyway, they would have been excluded by the previous studies on the gas-phase metallicity of SDSS cluster galaxies (Ellison et al., 2009; Cooper et al., 2008; Mouhcine et al., 2007), due to the redshift cut-off implemented in these works. Moreover, our analysis include, apart from oxygen abundance, the N/O ratio, a relevant observable to appreciate the effects on cluster galaxies (see also Petropoulou et al., 2011).

This paper is organized as follows: In §2 we describe the selection of the sample of SF dwarf galaxies in four nearby clusters, in §3 we derive the spectroscopic properties of our galaxy sample and in §4 their gas-phase metal abundances. In §5 we discuss the MZR of our sample clusters and in §6 the chemical enrichment of the SF dwarf galaxies relative to their HI mass content. In §7 we discuss the environmental effects that could affect the chemical evolution of cluster galaxies and finally in §8 we briefly summarize the findings of the present work. In this work we adopt the cosmological parameters =73 km s Mpc, , .

## 2 The Galaxy Sample

In this paper we study the chemical properties of SF dwarf galaxies in the central region (up to ) of four nearby clusters (Coma, A1367, A779, and A634), using spectroscopic data of SDSS-III DR8. In Table 1 we give details on the cluster properties and in Table 2 we give details on the regions studied in this work (see the captions of both Tables).

Coma is a very massive cluster, with very high velocity dispersion and X-ray luminosity (Ledlow et al., 2003). There are two central cD galaxies NGC 4874 and 4889, and a subcluster, projected to Mpc to the south-west of the core, and centered on NGC 4839, appears to be merging with the main cluster (Briel et al., 1992; Neumann et al., 2003). A1367 reveals a complex dynamical state (Cortese et al., 2004), with an elongated X-ray emission, showing multiple clumps, supporting a multiple merger scenario (Donnelly et al., 1998; Sun & Murray, 2002). Additionally, substructures falling into the cluster core (Cortese et al., 2006; Gavazzi et al., 2003) suggest an early stage of formation of the cluster. Cortese et al. (2008a), based on UV photometry, have found that star formation in Coma is substantially suppressed compared to that in the field, while A1367 has an abundance of bright SF galaxies (see Iglesias-Páramo et al., 2002). The optical luminosity function of galaxies has a much steeper faint-end slope in Coma than in A1367 (Iglesias-Páramo et al., 2003). On the other hand, A779 and A634 are two clusters of very low velocity dispersion. A634 X-ray luminosity is at the ROSAT detection limits, while A779 presents a nearly circular X-ray emission around the cD galaxy NGC 2832.

Cluster galaxies were selected using spectroscopic redshift information. In Figure 1 we plot the velocity histograms of all galaxies with SDSS spectroscopic data in the areas considered. We perform a gauss fit to the velocity distribution of each cluster. The central velocity obtained by the fit is in good agreement with the mean cluster velocity given by NED, indicated in Figure 1. We consider as cluster galaxies all the galaxies with velocities , where is the mean cluster velocity and is the dispersion given by the fit. In Table 2 we give the velocity range considered for each cluster, and the total number of galaxies with SDSS spectroscopic data found within the area considered and within the respective velocity range.

To select the dwarf galaxies, as the four clusters are located at similar distances (m-M35.0), we adopt the criterion used by Mahajan et al. (2010, 2011): M15, which results in a selection of , assuming (from Blanton et al., 2001) and the SDSS upper magnitude limit r = 17.77. This criterion yields to a sample of galaxies of log mass range M.

From this sample we select the SF dwarf galaxies adopting the criteria: SDSS BPT class=1 (classification as star-forming, based on the BPT diagram, Baldwin et al., 1981), and signal to noise (S/N) for the lines [O iii] 5007, [N ii] 6584, H, H, [S ii] 6717, [S ii] 6731. Finally, we have taken special care to include to our final sample of SF dwarf galaxies one spectrum per galaxy (a few galaxies have multiple observations), as well as to exclude spectra that correspond to distinct H ii regions of some parts of galaxies, because these spectra provide substantially underestimated mass values (see §3). The number of SF dwarf galaxies for each cluster is given in Table 2. For these galaxies we calculate the distance from the cluster center given in units of the . In Table 2 we also give the number of SF dwarf galaxies located in the cluster core, at distance , and the total number of dwarf galaxies at for each cluster.

We note that there are more SF dwarf galaxies, showing low S/N emission lines (the galaxies with BPT class=2 in SDSS), in the areas of the clusters of our sample, specially of Coma. These galaxies could be suffering the quenching effect of the cluster environment, illustrated well in studies based on the H EWs of SF galaxies (e.g. Balogh et al., 2004; Rines et al., 2005). In the present work, given that a high S/N emission line spectrum is required to derive properly spectroscopic properties (e.g. reddening coefficients) and gas metallicities, we do not include in our galaxy sample the low S/N SF population (we call them class 2 population). By doing this, we probably consider the more recently accreted galaxies to the cluster environment, and as a consequence the potential suppressing of their star-formation could be observed in an early stage. As we will discuss in §5, the behavior of class 2 population would not change the conclusions of this work.

Poggianti et al. (2006), based on a large sample of clusters, sampling the whole mass range from groups to massive clusters, have been able to show that the fraction of SF galaxies depends on galaxy mass, both for clusters at high-z and low-z. These authors found that about 20% of the galaxies in clusters at with km s are SF. Due to the SN restrictions, and the SDSS spectroscopic data incompleteness (specially in clusters, due to constraints in the fiber placement, see Blanton et al., 2005a, b), our galaxy sample is not a complete sample of the SF galaxies in these clusters. Thus, we should not compare our data with the rates found by Poggianti et al. (2006). However, the goal of this work is to investigate the imprint of the cluster environment on the chemical evolution of SF galaxies in clusters, the general quenching effect in clusters being well established by several previous works.

Figure 2 shows the color g-i histogram of dwarf galaxies (green line), SF dwarf galaxies (blue line), and SF dwarf galaxies to for each cluster (note that logarithmic scale has been used for Coma). It is well illustrated that the dwarf galaxy population in these clusters is composed by two main populations, red and blue, producing two maxima in the color histogram at g-i and g-i, following the general bimodal distribution of the galaxy population (Strateva et al., 2001; Baldry et al., 2004). Coma, has a larger fraction of red-sequence over blue SF dwarfs, in agreement with its low fraction of spiral galaxies (Dressler, 1980). The dynamically younger cluster A1367, in turn, hosts an almost equivalent blue and red dwarf population, as well as a larger number of SF dwarfs inside . These differences are in agreement with previous findings on the color-magnitude diagrams of these clusters (e.g. Mouhcine et al., 2011; Terlevich et al., 2001). We note that in all clusters there are some SF galaxies with red colors and that the SF dwarfs at span the whole color range. We have performed Kolmogorov-Smirnov (KS) tests which have shown that the SF dwarf galaxies to distances and those to distances are statistically indistinguishable with respect to their color for Coma, with a probability of 89%. Smaller probabilities have been found for the three remaining clusters (for A779 and A634 possibly due to the lower galaxy number).

Kewley et al. (2005) argued that using the 3-arcsec diameter aperture of the SDSS spectra to estimate SFRs and metallicities of late type galaxies located at could introduce systematic biases. This is due to the significant radial metallicity gradients that can be found in SF spiral galaxies (e.g Kennicutt et al., 2003; Magrini et al., 2007; Bresolin et al., 2009). However, in this work we focus our attention to dwarf/irregular galaxies, where spatial metallicity gradients are not expected to be important (Kobulnicky & Skillman, 1997; van Zee et al., 2006; van Zee & Haynes, 2006). Thus, taking into account the smaller diameter of dwarfs, we do not expect aperture biases to be important for the galaxies of the present sample.

## 3 Spectroscopic Properties

SDSS-III DR8 provides new emission line measures (Tremonti et al., 2004; Brinchmann et al., 2004), derived after correcting the spectra for the underlying stellar continuum, using the high-resolution population synthesis models of Bruzual & Charlot (2003). The Balmer emission lines can be severely affected by the absorption of the underlying stellar component (Martín-Manjón et al., 2008, 2010) and this would affect the derived values of the reddening coefficient, the line indices and, as a consequence the gas-phase metallicity (see §4). This effect is expected to be negligible here, allowing a reliable study of the galaxies’ gas metallicity. We find that DR8 measures of Balmer emission lines such as H and H are significantly improved as compared to DR7 spectroscopic data.

We use the Balmer emission lines H, H, H (and where the S/N permits also H) to derive the reddening coefficients c(H) for our sample of SF galaxies, using Case B approximation (Osterbrock, 1989), adopting the Cardelli et al. (1989) extinction law, and taking into account the line measurements error. Mahajan et al. (2010) provide Spitzer MIPS 24m flux for a sample of Coma cluster galaxies. From this sample, 25 galaxies belong to our sample of SF dwarfs. For these 25 galaxies we find that the derived reddening coefficient c(H shows a tight correlation with their flux at 24, in the line of Relaño et al. (2010), who suggest that the dust responsible for the Balmer extinction should be emitting at 24 .

The presence of both singly and doubly ionized oxygen line transitions in the optical wavelengths has permitted to develop an efficient metallicity calibration (Pagel et al., 1979) based on the indicator R=([O ii] 3727+[O iii] 4959+[O iii] 5007)/H. Obtaining the [O ii] 3727 line measures from DR8, however, resulted problematic. Although some spectra present high S/N [O ii] 3727 line, no measure is provided by DR8. For the galaxies that SDSS does not provide [O ii] 3726,3729 measures with S/N, but [O ii] 3727 is detectable in the spectrum with S/N, we measure the [O ii] 3727 integrated flux. We have verified that our measures and DR8 give consistent emission line fluxes for several emission lines (e.g. [N ii] /H and [N ii] /[O ii] , when there is [O ii] ). Table 3 presents the number of galaxies for which we have SDSS [O ii] 3726,3729 measures, the number of galaxies for which we measure integrated [O ii] 3727, the total number of galaxies with [O ii] measure and the total number of SF dwarf galaxies for A1656 and A1367.

Despite this effort, we could not obtain [O ii] 3727 for all galaxies belonging to the four clusters of z0.023. This is because [O ii] 3727, redshifted to this velocity, lie at the edge of the wavelength range covered by SDSS spectroscopy. We have performed Kolmogorov-Smirnov (KS) tests and have verified that the subsample of galaxies with [O ii] measured (by DR8 or this work) and the whole sample of SF dwarf galaxies are statistically indistinguishable in their properties: , metallicity (by PMC09, see §4) and mass. We conclude that the lack of [O ii] measures is a random effect, due to the radial velocity of each galaxy for the cluster velocity dispersion, and does not correlate with any of the fundamental galaxy properties.

We have used the reddening corrected line fluxes for our sample of galaxies to compute five standard optical line ratios, that we show in Figure 3, combined into three commonly used diagnostic diagrams. Our sample galaxies show typical line ratios of normal SF dwarf galaxies. Additionally, the galaxies located at span the whole range of values in all diagrams, not showing any peculiarity with respect to the spectroscopic properties of their ionized gas. Figure 4 shows the excitation as measured by [O iii] /H versus EW(H) (that we derived from the stellar absorption corrected H flux and the adjacent continuum as provided in DR8) for the SF dwarf galaxies in Coma (left) and A1367 (right). Filled symbols represent galaxies to distances and again, we see that the galaxies in the cluster core are not different from the rest regarding excitation.

We have also obtained the median estimate of the total stellar mass of the galaxies from the MPA-JHU spectroscopic catalog (Kauffmann et al., 2003). We have found that this galaxy mass estimate is consistent with the mass we derived using the k-correct algorithm (Blanton & Roweis, 2007). We also have found a good agreement between the MPA-JHU mass estimates and Mouhcine et al. (2011) galaxy masses, derived by K band imaging, for 7 galaxies of our sample that belong also to the sample of Mouhcine et al. (2011).

## 4 Oxygen Abundances

Accurate abundance measurements require the determination of the electron temperature, which is usually obtained from the ratios of auroral to nebular line intensities such as [O iii] 4363 over [O iii] 5007. This is often referred to as the “direct” method. When auroral lines are not detected, empirical methods are generally used, based on the ratios of bright forbidden lines to hydrogen recombination lines. Among these, the most widely used abundance indicator is the R333R=([O ii] 3727+[O iii] 4959+[O iii] 5007)/H index (although still carrying significant hazards, e.g. see Pérez-Montero & Díaz, 2005). An other extensively used parameter to derive gas phase metallicities is N2444N2=log([N ii] 6584/H) (Storchi-Bergmann et al., 1994; Pérez-Montero & Contini, 2009, PMC09 from now). The use of this parameter has two important advantages: the relation between N2 and the oxygen abundance is single-valued and the emission lines involved are very close in wavelength, so the N2 parameter is almost free of uncertainties introduced by reddening correction or flux calibration. A third very useful index is O3N2555O3N2=log{([O iii] 5007/H)/([N ii] 6583/H)} (Alloin et al., 1979), which shows a relatively tight and linear relationship with 12+log(O/H) (Pettini & Pagel, 2004, from now on PP04) for O3N2.

The calibration of these parameters can be empirical (e.g. Pettini & Pagel, 2004; Pilyugin & Thuan, 2005; Pilyugin et al., 2010, from now P10), that is using direct measurements of the oxygen abundance of H ii regions in the Local Universe, or theoretical (e.g. McGaugh, 1991; Kewley & Dopita, 2002; Tremonti et al., 2004; Nagao et al., 2006; Dors et al., 2011), that is using photoionization models covering different ranges of physical parameters. Though, it has long been reported (e.g. Kewley & Ellison, 2008) that there exist systematic differences that can reach up to dex between most model and empirical calibrations, specially in the higher metallicity regime (but see Dors et al., 2011).

A strong hint that the direct measures and the empirical calibrations should yield metallicities closer to reality was given by Bresolin et al. (2009). These authors derived the metallicities of 28 H ii regions in the galaxy NGC300 using the direct method, and found that these metallicities agree with the abundances derived from young massive stars, with stellar and nebular abundances giving virtually coincident slopes and intercepts of radial gradients. An important quantity that should be taken into account when comparing stellar and nebular oxygen abundances is the amount of oxygen depleted onto dust grains in ionized nebulae. According to Peimbert & Peimbert (2010) this fraction amounts to about 0.10 dex for low-mass galaxies. Bresolin et al. (2009) discussed this effect and found that, even considering a gas depletion factor of dex, the intercepts of the nebular and stellar abundance gradients would be consistent. Then they compared their direct estimates of the metallicity of the ionized gas with the most widely used model calibrations, and they report that model calibrations yield higher metallicities of about a factor of two in the O/H range derived. On the contrary, the abundances derived using empirical calibrations use to be in good agreement with their direct measurements. Another strong hint was given recently by Dors et al. (2011). These authors produced models that reproduce O/H estimates consistently with the values derived using the direct method, for the upper range of O/H values, where all previous models showed systematic discrepancies.

Simón-Díaz & Stasińska (2011) derived the O/H abundance ratio using 13 B type stars from the Orion star forming region and obtained an excellent agreement with the O/H recombination line abundances derived for the Orion nebula by Esteban et al. (2004). For the Orion nebula Simón-Díaz & Stasińska (2011) considered a correction of 0.12 dex due to the dust grain depletion, that is in excellent agreement with the value estimated by Mesa-Delgado et al. (2009). Since recombination lines are generally very weak they are not observed in our sample of galaxies. In any case, this analysis is based on the relative comparison of galaxy metallicities; thus no relevant effect is expected from the derivation of absolute abundances.

In this work we use the empirical calibrations of P10, N2 calibration of PMC09, and O3N2 calibration of PP04 to derive the oxygen abundances. We note that for our sample of galaxies always O3N2 , lying in the valid range of PP04 calibration. In Figure 5 (left) we compare the oxygen abundances derived using the PMC09 and PP04 calibrations for the galaxies in the two most populated clusters of our sample, and we can see a very good agreement. The use of the P10 empirical calibration requires the measure of the [O ii] 3727 line, that we do not have for all our sample galaxies, see §3. In Figure 5 (middle) we compare the oxygen abundances derived using the P10 and PMC09 calibration, for those galaxies with [O ii] measured. We see that the PMC09 calibration yields slightly higher metallicities (dex) than P10 calibration, in agreement with our previous considerations in Petropoulou et al. (2011). We also derive the oxygen abundance of galaxies with [O ii] measures using the model of Dors et al. (2011), based on the diagnostic diagram [O iii] /[O ii] versus [N ii] /[O ii] 666We here adopt the value 12+log(O/H)=8.69 for the solar oxygen abundance (Asplund et al., 2009).. We find that these models systematically underestimate metallicities for 12+log(O/H)8.2. On the contrary, for 12+log(O/H)8.2, we find a good agreement with the empirical PIL10 estimates (see Figure 5, right).

In order to avoid limiting our sample only to galaxies with [O ii] measures, in the following we use the PMC09 calibration of the galaxy oxygen abundance. We are aware that PMC09 could slightly overestimate (dex) the galaxy O/H abundance, but as discussed later on, we do not expect this to affect our discussion, as we always perform careful comparison of metallicities derived using the same method. We also use the N2S2 calibration of PMC09 to derive N/O ratios for our SF dwarf galaxies as in Petropoulou et al. (2011).

In order to illustrate the properties of our galaxies, in Figure 6 we plot the H flux of the SF region covered by the SDSS fibre, for the SF dwarf galaxies in Coma (triangles) and A1367 (diamonds), versus the EW(H). The H fluxes have been corrected for extinction (from the derived c(H), see §3), and have been transformed to an equivalent number of ionizing photons Q(H), under the assumption that no ionizing photons escape the galaxies (Osterbrock, 1989). The points are color coded to the derived galaxy metallicity. Then we add the loci of Starburst99 models (Leitherer et al., 1999) for instantaneous bursts of total mass (continuous lines) and (dashed lines) in massive stars for metallicities Z=0.02 (in red), Z=0.008 (green), and Z=0.004 (black). According to these models, and assuming single bursts, the SF regions studied would span a mass range of .

## 5 The Mass-Metallicity Relation

In Figure 7 we plot the galaxy metallicity, as measured by 12+log(O/H), versus galaxy stellar mass, for the SF dwarf galaxies in Coma (A1656, triangles), A1367 (stars), A779 (squares), and A634 (circles). Filled symbols correspond to galaxies to a distance from the cluster center. We see that the SF dwarf galaxies follow well defined sequences on these plots. Having a look at the MZR of A1656 and A1367, we observe that at low galaxy masses (), the derived metallicities cover a range of dex, while at higher masses () the relation becomes tighter, appearing to shape a triangle instead of a simple linear correlation. Additionally, we observe that, for the same bin of mass, galaxies at are preferentially located at the upper part of the global sequences. For the two clusters (A779, A634) of lower mass, we do not verify neither of these two features, as we discuss in the following.

We first note that the behavior observed seems independent of the abundance calibration used. When we use the PP04 calibration of O3N2, we obtain the same mass-metallicity sequences and the SF dwarf galaxies in the core of A1656 and A1367 still crowd the upper part of the MZR (we can not do the same using P10 calibration because we do not have enough galaxies with measured [O ii] inside ). As a consequence, the effect observed seems to be significant, despite the intrinsic errors of the N2 and O3N2 calibrations.

Second, in §3 we have verified that our sample galaxies do not show any systematic differences regarding excitation, as a function of their location in the cluster (in or out ). Thus, the higher metallicities derived for galaxies at , with higher [N ii] /H , are not biased by excitation effects (see Berg et al., 2011, for a thorough discussion on this effect).

The third sanity check has been to test whether the position of the SF dwarfs at on the MZR, is the result of some bias in the mass estimate. We have found that galaxy stellar mass, versus other galaxy properties, such as broad-band colors and excitation, do not show any correlated dependance of the cluster-centric distance of the galaxies. After all these tests, we assume that the shift observed toward higher metallicities, for galaxies at , mainly in Coma, should not be the result of some bias in the abundance and/or mass estimates.

In order to explore further the trends that appear in Figure 7, we have performed a bisector linear fit to the MZR of each cluster, considering first the galaxies at , second the galaxies at , and third all the galaxies together. In Table 4 we present the obtained slopes and the corresponding errors from the dispersion around the linear fit. We observe that the MZR fits for galaxies at have slopes for all the clusters of our sample. Lee et al. (2006) have derived the MZR for 25 nearby dwarf irregular galaxies (), extending the well-known SDSS MZR (Tremonti et al., 2004) to the lower mass range. Lee et al. (2006) have found a very tight correlation, over the whole stellar mass range, with slope ; thus the slopes derived here are in very good agreement with their results. As we discuss further in §7 this value of the slope is supported by the semi-analytic models and hydrodymanic simulations of Finlator & Davé (2008) and Davé et al. (2011) where momentum-driven wind scalings are introduced to explain the MZR.

For Coma, and to less extend for A1367, a flattening of the MZR is observed when considering galaxies inside , with the fits showing clearly smaller slopes, within the quoted errors (see Table 4). This is made evident in Figure 7, where we plot with a continuous line the linear fit for galaxies at and with dashed line the fit corresponding to galaxies at .777We note here that the class 2 population refered to in §2, span the higher mass range, where we have seen that the MZR becomes tighter, and had it been included, would not have significantly changed the derived slopes. For A779 and A634, considering the large errors (due to the reduced number of SF dwarf galaxies in these clusters), we do not appreciate any difference in the slope inside and outside the cluster core. In Figure 7 we plot the overall fit to the MZR for A779 and A634.

In Figure 8 we plot the N/O ratio versus galaxy stellar mass, and again we find a good correlation, even tighter than the MZR, for all our clusters. Again, filled symbols correspond to galaxies to a distance from the cluster center. As before, galaxies in the cluster core of Coma and A1367 tend to crowd the upper part of the global sequence. Oxygen, produced in Type II SNe, is typically released after 10 Myr, while the bulk of nitrogen is produced and released over a substantially longer period, 250 Myr. This is an important piece of information to take into account when one attempts to identify the mechanisms relevant for the chemical evolution of cluster galaxies. As we will discuss further in §7, this chemical “clock” seems to imply that if there is an environmental effect driving the observed difference in the abundance ratio, this should be acting since at least yr ago.

Observing the triangular shape of the MZR of A1656 and A1367 in Figure 7, we investigate whether there could be a physical cause for the galaxies showing lower metallicity for the same bin of mass. The MZR has been found to show a second parameter dependence on SFR (Amorín et al., 2010; Mannucci et al., 2010; Lara-López et al., 2010; Cresci et al., 2011). This can be easily understood by the correlated behavior of gas-phase metallicity and SFR after a gas infall event: when a galaxy accreates metal-poor gas, its metallicity is diluted to values typical of lower mass galaxies, while its final mass increases, and concequently the galaxy moves below the MZR. In the same time, the presence of large amount of gas stimulates star-formation, and this is the reason why galaxies with higher SFR have lower metallicities at a given stellar mass. This should be a transient phase, and when the galaxy will consume the gas, producing the new metals, will return to the mean MZR (Davé et al., 2011; Dalcanton, 2007).

Mannucci et al. (2010) claimed that the MZR is the projection in the 2D space of a more fundamental 3D relation between stellar mass, gas metallicity and SFR and proposed the quantity (SFR) which defines a projection of the MZR that minimizes the scatter of local galaxies. We have used the median estimate of the total SFR provided by DR8 (Brinchmann et al., 2004) to explore whether the quantity could decrease the scatter observed in Figure 7. The SDSS SFR estimates have been derived by combining emission line measurements of the fibre spectrum and applying aperture corrections by fitting models to the photometry outside the fibre (as in Gallazzi et al., 2005; Salim et al., 2007). In Figure 9 we plot the oxygen abundance 12+log(O/H) versus for the SF dwarfs in A1656 and A1367, and we see that the scatter drops for the galaxies at lower metallicities for the same bin of mass (the Spearman’s correlation coefficient increases from 0.73 in Figure 7 to 0.77 in Figure 9 for A1656, and from 0.80 to 0.83 for A1367). However, the separation we observe in the MZR between galaxies inside and outside is kept equally.

Moreover, in this metal dilution scenario, the N/O ratio should not be affected (e.g. the “green pea” galaxies show normal N/O ratios for their mass, Amorín et al., 2010), and the relation between N/O versus stellar mass should not show correlated scatter. However, as Figure 8 shows, the N/O ratio segregates galaxies inside , these galaxies also showing higher N/O ratios for the same bin of mass. The effect of infalling gas is expected to be more relevant for galaxies at high redshift (e.g. Cresci et al., 2010; Tacconi et al., 2010; Dekel et al., 2009), thus the fundamental metallicity relation indroduced by Mannucci et al. (2010) considerably accounts for the evolution of the MZR with redshift (Erb et al., 2006; Maiolino et al., 2008; Mannucci et al., 2009). But in the local Universe, where our sample of clusters are located, different mechanisms should be invoked to explain the trends observed.

The question arises now as to whether the galaxies inside have evolved in a different way than the other cluster galaxies, rendering them chemically more enriched than the galaxies at , or alternatively whether less metallic galaxies, for the same bin of mass, are more vulnerable to the quenching effect of the cluster environment, and this could be the reason we do not observe them preferentially in the cluster core. To investigate the second idea we searched whether there is any observable trend with the distance from the cluster center of the morphological type of our sample of SF galaxies.

We have calculated the standard concentration index , where and are the radii enclosing 90 and 50 per cent of the Petrosian r-band luminosity of the galaxy. We have found that our sample galaxies show typical values as related to late-type galaxies (, see e.g Shimasaku et al., 2001; Strateva et al., 2001) and we see no trend of the value in and out . We have additionally checked that galaxies in the same bin of mass do not show any concentration-metallicity correlation. A more detailed morphological study could be of interest here; Penny & Conselice (2011) found remarkably smooth structures of dwarf galaxies in the Perseus cluster, and suggest that dwarfs in cluster cores should be highly dark matter dominated to prevent their tidal disruption by the cluster potential. In the present study we can not verify the hypothesis that SF galaxies inside are more concentrated/compact, and consequently more resistent to the hostile cluster environment. Thus, morphology does not seem to explain their relative position in the upper part of the MZR; this should be the result of some different evolution.

To investigate whether the scatter observed in the MZR could be related to differences in the underlying stellar population, we have obtained for all our cluster galaxies the spectral index (4000) (Balogh et al., 1999), after correction for emission lines, as given by SDSS DR8. All our SF dwarfs show , which corresponds to typical ages Gyr. Additionally, we have seen that the gas-phase metallicity does not show any correlated behavior with the age of the underlying population. In §3 we have shown that SF dwarfs at in Coma do not show any observable difference in the g-i color distribution as compared to the galaxies at . Thus, the more metallic SF dwarf galaxies found in the cluster core, specially of Coma, are neither older nor redder than the rest of SF galaxies.

To quantify the effect observed in the MZR for Coma and A1367, in the upper panel of Figure 10 we plot the mean difference of the derived 12+log(O/H) abundance for each galaxy with respect to the 12+log(O/H) given by the bisector linear fit, as a function of the cluster-centric radial distance, sampled in a bin of 0.5 . The errors correspond to the standard deviation from the mean value of log(O/H)-log(O/H). In the lower panel we plot the same for the log(N/O) ratio. We observe that in Coma, the mean difference in 12+log(O/H), in the closest bin to the cluster core, is positive and above the rms error (i.e. positive for all objects), and can reach above dex. The same is found for the log(N/O) ratio, which can get up to dex difference in the cluster core. In A1367 smaller differences are obtained, of dex in log(O/H) and dex in log(N/O). This divergence, combined to the fact that for the low-mass clusters of our sample, this trend has not been revealed at all, seem to indicate that the cluster mass is a most relevant parameter. As it will be discussed in §7, what seems to drive the disparate evolution of SF dwarf galaxies in Coma appears to be related to the properties of the forceful ICM of this cluster.

We explore the dependence of the trend observed for our cluster galaxies as a function of the local galaxy density, using the density estimator .888Defined as the logarithm of the density obtained to the average of the projected distances to the fourth and fifth nearest neighbors, including all galaxies with SDSS spectra, within the cluster region (see Petropoulou et al., 2011). In the upper panel of Figure 11 we plot the mean difference of the derived 12+log(O/H) for each galaxy with respect to the 12+log(O/H) given by the bisector linear fit, as a function of , in a bin of 0.5 dex. The errors correspond to the standard deviation from the mean value of log(O/H)-log(O/H). In the lower panel we plot the same difference for the log(N/O) ratio. In the highest local density bin we find a mean difference in log(O/H) of dex for both Coma and A1367, in agreement with the previous findings (e.g. Ellison et al., 2009). Log(N/O) presents a measured difference of dex for both clusters.

We note here that estimating the local galaxy density for clusters with large velocity dispersions as in Coma (1000 km s) and A1367, bears significant hazards. Allowing neighbors to have velocity differences of the order of , could result in counting out galaxies in the cluster core, and thus underestimating density. Converselly, allowing velocity differences of 2 or 3 could introduce a severe bias to relatively isolated objects, overestimating their density. In Figure 11 we use derived permitting neigbors to have velocity differences of 2000 km s, but we have checked that the behavior does not change using velocity difference of , , and .

Previous works have pointed out the effect of the environment on the gas-phase metallicity of galaxies (Mouhcine et al., 2007; Cooper et al., 2008; Ellison et al., 2009) and have found that galaxies within dense environments show statistically higher metallicities of dex at the same bin of galaxy mass. These works, however, have related this behavior to local galaxy density. Specially Ellison et al. (2009) have discussed the importance of cluster membership, and although they have found that enhanced metallicities are present at distances from the cluster center, they have concluded that the enhancement observed is driven by local overdensity and does not depend on cluster properties such as , or cluster mass.

There is a general good relation between the cluster-centric distance and the local galaxy density, within the cluster virial radius (e.g. Rines et al., 2005), so the metallicity enhancement of the cluster core galaxies is expected to appear as a function of the local galaxy density as well. However, in the present work we find that the enhancement of the gas-phase metallicity is sensitive to the cluster mass, and appears to be more prominent if we consider the region of a massive cluster such as Coma. This evidence points towards a possible connection of the chemical enrichement of cluster galaxies with their ICM properties.

## 6 Chemical Enrichment vs Hi Mass

In clusters, the gas interchange of a galaxy with its environment is expected to be very relevant. Observable gas tails give evidence of ongoing gas stripping and a significant fraction of spirals have been found to have truncated ionized gas disks compared to their stellar disks (see §1). Recently, Jaffé et al. (2011), based in a large sample of emission line galaxies from the EDisCS999European Southern Observatory Distant Cluster Survey sample at intermediate redshifts (), found that the fraction of kinematically disturbed galaxies increases with cluster velocity dispersion and decreases with distance from the cluster centre, but remains constant with projected galaxy density. In addition, disturbed gas kinematics does not co-occur with morphological distortions as traced by optical (HST) imaging, suggesting the mechanism that affects most the gas of cluster galaxies has to be linked with the ICM.

The atomic gas is expected to be the first to suffer the ram-pressure stripping in the cluster environment. Indeed, cluster spiral galaxies are more HI deficient than similar galaxies in the field, and this deficiency appears to be increasing towards the cluster center (Haynes et al., 1984). Additionally, because of the shallower gravitational potentials, dwarfs can be extremely fragile to RPS (Solanes et al., 2001). In the Virgo cluster core, Chung et al. (2009) found many HI deficient galaxies, and at intermediate distances, at 1Mpc from the center, they found a remarkable number of galaxies with long, one-sided HI tails pointing away from M87. Truncated HI disks were found even in Pegasus I cluster (Levy et al., 2007), a cluster with a low level of X-ray emission.

An interesting method to check for modulations of the HI content of SF galaxies is through comparing the galaxy chemical enrichment with the theoretical values predicted by the “closed-box” model (Edmunds, 1990). According to this model, the mass fraction of metals should be a direct function of the gas mass fraction:

 μ=Mgas/(Mgas+M⋆)

Here, the gas mass is considered to be the mass of HI () with a correction for neutral helium . We neglect the contribution due to molecular hydrogen, as this seems to be small and difficult to evaluate. According to Israel (1997), for Magellanic type irregular galaxies, the mass ratio is 0.2, which would mean an increase of 0.08 dex in the estimate. However, the CO to conversion relation is not well known for low metallicity systems and the fraction might be even higher (e.g. Keres et al., 2003).

Comparing with chemical evolution models is actually possible only when information on the HI mass of a galaxy is available. Two of our sample clusters have recently published HI data. Cortese et al. (2008b) have presented 21 cm HI line observations of the central part of A1367, as part of the Arecibo Galaxy Environment Survey (AGES). This sample covers of the total region covered by our sample and of the region with and their HI lower mass limit, at A1367 distance, is M. They have detected 57 galaxies that belong to A1367 cluster, out of which 17 also belong to our sample of SF dwarf galaxies. Haynes et al. (2011) have recently released a catalog with 21 cm HI line sources, covering of the final ALFALFA101010Arecibo Legacy Fast Arecibo L-band Feed Array (ALFALFA). survey area. This release includes of the area of A1656 covered by the present work, also of the area with , and their low HI mass limit at the Coma cluster distance, is M ( at the Virgo cluster distance, Giovanelli et al., 2005). We have found 31 objects in common with our sample of SF dwarf galaxies.111111We have derived the using the standard formula (Haynes et al., 2011) and adopting the same mean distance for all cluster galaxies.

In Figure 12 (left) we plot the oxygen abundance versus the gas mass fraction (in the form ) for A1656 (triangles) and A1367 (stars) for the galaxies we have HI data (filled symbols mean ). In Figure 12 (right) we plot the gas mass fraction versus the stellar mass. We add (with smaller points) all the SF dwarf galaxies of our sample that do not have HI measurements, but they are located within the regions mapped by AGES and ALFALFA (up to the present release). We assign to these galaxies the HI mass detection limit of each survey, this being an upper limit of the HI mass of these galaxies (the arrows indicate the direction to which these upper limits could be displaced). Again, open and filled small symbols mean outside and inside the cluster respectively.

We observe that A1656 dwarf galaxies located to distances (filled triangles), all except one, have been assigned an upper limit. In contrast, inside of A1367 we have both upper limits and detections. This could be an evidence of the stronger and more effective ram-pressure exerted by the ICM of A1656 than that of A1367. Tecce et al. (2010) have found that at low-z, the mean ram-pressure is dyn cm in clusters with virial masses , while in lower mass clusters of , the mean ram pressure is approximately one order of magnitude lower. At distances (open triangles), there are measured HI masses as well as upper limits for both clusters. Theoretical models have started to investigate the effects of the ICM turbulence, substructure and shocks during groups infall, finding that RPS can be effective even at the cluster outskirts and fillaments, as well as during tidal interactions (e.g. Tonnesen & Bryan, 2008; Kapferer et al., 2008).

We then compare the position of our galaxies in the left plot of Figure 12, with model predictions for the yield. The green continuous line indicates the model value , that is the theoretical yield of oxygen expected for a Salpeter IMF and constant star formation rate, for stars with rotation following Meynet & Maeder (2002) models (van Zee & Haynes, 2006). Both, infall and outflow of well-mixed material will result in effective yields that are less than the true yield, as the enriched material is either diluted or lost from the system. The blue dashed line corresponds to a lower yield , almost 1/4 of the model prediction, which is the effective yield found for a fraction of the sample of isolated dwarf irregular (dI) galaxies of van Zee & Haynes (2006).

Some of the SF dwarf galaxies located at cluster-centric distances in both clusters, and a fraction of the SF dwarfs in A1367 cluster core (), appear to follow the theoretical closed-box model predictions (considering our metallicity typical uncertainty dex, see §4). These galaxies in the core of A1367 could be “newcomers” (see also Petropoulou et al., 2011), observed before the action of RPS lowers significantly their atomic gas content, shifting them towards lower effective yield values.

We also observe a tendency to lower effective yields as metallicity increases, suggesting that these galaxies host lower HI mass than the expected by the closed-box model (a similar behavior was found for A2151 SF galaxies, see Petropoulou et al., 2011, Fig 14). The atomic gas deficiency is confirmed when we compare the HI content of our sample galaxies (measurements or upper limits) with the of field counterparts of the same absolute magnitude, as calculated following Toribio et al. (2011). Two representative examples are the dwarf galaxies CGCG 97-073 and CGCG 97-079 (marked on the plot) that have been investigated before (Iglesias-Páramo et al., 2002; Gavazzi et al., 2001), showing long (H) tails that reveal their recent interaction with the ICM. The HI content of CGCG 97-079 appears to be severely affected: this galaxy has the lower HI mass detected and lays below in the left panel of Figure 12. CGCG 97-073 in turn does not appear to have lost a large amount of HI yet.

We conclude that, to the distance of Coma and A1367, given the detection limit of the HI data, we observe the atomic gas of galaxies that still contain a substantial fraction of their HI mass. The almost absolute lack of detections in the central part of A1656, indicates that galaxies there have been severely affected by the cluster environment, and the ram-pressure stripping has partly or completely removed their atomic gas content, rendering them undetectable in 21 cm HI line. In turn, in the central part of A1367 there are some SF dwarf galaxies still with detectable HI mass (i.e. possible “newcomers”).

## 7 Discussion

It is well established that dense environments, such as cluster cores, present larger fractions of passive galaxies than the field (e.g. Dressler, 1980; Poggianti et al., 1999; Treu et al., 2003; Finn et al., 2005; Poggianti et al., 2006). This confirms the important influence of the environment on galaxy evolution. Although the big debate on nature versus nurture remains open, accumulated evidence seems to indicate that, in the particular case of dwarf galaxies, environment is a fundamental driver of their evolution (Mahajan et al., 2010; Haines et al., 2007, 2006). It has been found (Smith et al., 2011; Porter et al., 2008) that significant infall of low-mass SF galaxies exists along the filamentary structures onto the densest clusters. Reaching cluster cores, dwarf galaxies can experience a star-burst event, induced by their first interaction with the ICM (Gavazzi et al., 2001; Sun et al., 2007; Levy et al., 2007; Petropoulou et al., 2011). Then, RPS is expected to be very efficient (e.g. Tonnesen et al., 2007) and it could strip the gas of a dwarf galaxy in very short timescales (yr) converting a large fraction of dwarf galaxies, first into post-starburst as those observed in Coma (Mahajan et al., 2011; Poggianti et al., 2004), and finally into passive galaxies.

Despite the fateful “switching off” of star formation predicted by the models of RPS, we have found quite a number of SF dwarf galaxies in our clusters. The selection of our galaxy sample to have strong emission lines means that they do still contain ionized gas. However, as we show in §6, a large fraction of these galaxies do not contain the amount of HI mass expected by the closed-box model, indicating that they should be suffering RPS which has first affected their neutral gas content. Tecce et al. (2010) simulations seem to indicate that at , 50% of galaxies inside the virial radius () have experienced important ram-pressure effects. These pressures appear to have a strong effect on the cold gas content: 70% of the simulated galaxies within have been found to be completely depleted of cold gas. Additionally, Tecce et al. (2010) models have suggested that the rate at which the cold gas of a galaxy is stripped, depends on the halo virial mass of the host cluster; thus less massive galaxies within massive clusters are the most affected.

In order to investigate the chemical evolution of the SF dwarf galaxies in the present sample of clusters, we have derived their MZR. The MZR of galaxies is well established (e.g. Tremonti et al., 2004; Lee et al., 2006). There exist a variety of chemical evolution models and hydrodynamical simulations advocating different physical mechanisms to explain this relation. For example, the ejection of metal-enriched gas by galactic outflows, triggered by e.g. (multiple) supernovae explosions, could be more efficient in systems with shallower potential wells (e.g. Larson, 1974; Marconi et al., 1994; De Lucia et al., 2004). A variable efficiency of star formation, increasing with galactic mass (often named “downsizing”), could also be invoked to reproduce the MZR (e.g. Matteucci, 1994; Tissera et al., 2005; de Rossi et al., 2007; Calura et al., 2009). For dwarf galaxies, where supernova winds are expected to be more important, a combination of downsizing and winds could affect the shape and slope of the MZR (Tassis et al., 2008; Spitoni et al., 2010). Other mechanisms, such as dilution caused by infall and variations in the IMF have also been proposed, however they show problems reconciling dwarf SFRs and metallicities (see the review by Tolstoy et al., 2009).

Recently it has been found that hydrodynamic simulations that incorporate momentum-driven wind scalings provide among the most successful overall fits to a wide range of observed galaxy properties (Finlator & Davé, 2008; Oppenheimer et al., 2010; Davé et al., 2011). Finlator & Davé (2008) developed a simple analytic model to understand the MZR, where the gas-phase metallicity of a galaxy is set by a balance between inflow and outflow, plus star formation. In this “equilibrium” model the mass outflow rate (, i.e. the balance between inflow and outflow) scales inversely with circular velocity, and this naturally reproduces a slope of the MZR at . Lee et al. (2006) observationally obtained this slope of the MZR for a sample of nearby dIrr galaxies, with a remarkably small dispersion. This value of the slope is in agreement with the slope given by the linear fit of our sample galaxies, located either in the low-mass clusters A779 and A634, or at the outer regions of Coma and A1367.

Davé et al. (2011) discuss the effect of SFR and of environment on the MZR, and find that the overall trends are consistent with the expectations from the “equilibrium” model. The scatter observed in the MZR for our sample galaxies does not correlate with SFR (§5), thus it does not seem to be related to metal dilution suffered by the galaxies showing lower metallicities at the same bin of mass. Regarding environment, Davé et al. (2011) hydrodynamic simulations reproduce low mass galaxies in dense environments showing higher metallicities for the same bin of mass than galaxies at lower densities. The authors suggest that this should not be the result of galaxies processing more gas into stars, but rather the consequence of the environmental dependence of wind recycling. Much of the material entering into the galaxies’ ISM at z=0 is recycled winds (Oppenheimer et al., 2010), that is the reaccretion of the wind material when the wind velocity does not exceed the galaxy escape velocity (see also Bekki et al., 2009). This mechanism accounts for the fact that metallicity differences seem to disappear for high mass galaxies: at high galaxy mass, wind recycling is so effective that metallicity approaches the theoretical yield regardless of environment. Instead, at low masses, most ejected material could normally escape the galaxy (Oppenheimer et al., 2010), unless the galaxy resides in a dense environment, with a hot gaseous halo that could significantly slow winds.

In §5 we have seen that SF dwarf galaxies in the Coma cluster core, and to less extend in A1367, show higher metallicities for the same bin of mass than galaxies outside . This results in the flattening of the MZR for dwarfs to distances . The amount of this flattening appears to be related to cluster mass, with the less massive clusters of our sample not showing at all this behavior. This observed trend matches well the above wind recycling scenario. Models suggest that at high density environments (i.e. in cluster cores) galactic winds can be suppressed by the high pressure of the ICM (Schindler et al., 2005; Kapferer et al., 2006, 2009). Accordingly this could cause faster recycling and thus preventing galaxies from loosing their metals. The suppression to take place needs ICM pressures dyn cm (Schindler & Diaferio, 2008). Tecce et al. (2011) have found that on the outskirts of clusters as (or more) massive as Coma () the ram-pressure reaches dyn cm, the same order of magnitude as in the core of clusters of virial mass . In turn, the ram-pressure in the core of massive clusters is expected to be times higher. Consequently, as we go to clusters of higher X-ray luminosity, higher would be the pressure exerted and the wind suppression is expected to be more effective. The wind suppression in the cluster cores could approximate the no-wind case depicted at Davé et al. (2011) hydrodynamic simulations, where the slope of the MZR is found to be flatter than when considering wind scalings.

The important role of the density of the ICM can be supported by Poggianti et al. (2009a, b) findings. These authors have found that the fraction of post-starburst galaxies depends on cluster mass, and the fraction of spiral galaxies in clusters in the Local Universe is anticorrelated with (whereas no trend is observed with cluster velocity dispersion). Now it seems to turn up that it is not only to the quenching efficiency, but also to the chemical evolution of low-mass cluster galaxies that the properties of the ICM could play an important role.

Another possibility that emerges from the “equilibrium” model scenario is that galaxies in denser regions could result having higher metallicities as a concequence of curtailed inflow, due to gas stripping and/or strangulation. The relation is the result of the balance between inflow and outflow processes. If gas accretion is truncated by some sort of environmental suppression, then the galaxy consume its gas to form metals along the locus of and will move above the mean MZR (Davé et al., 2011). If accretion does not restarts, the galaxy will end up exhausting all the available gas for star-formation and finally be transformed to a passive galaxy. Accordingly, our SF dwarf galaxies that lie above the MZR could be in this particular phase of environmental quenching, where enhanced metal enrichment precedes the switching off of star-formation. As we show in §6 the majority of galaxies inside seem to have suffered important gas removal, specially in A1656.

The above picture could describe well what is happening to the galaxies in the Coma cluster core, where SF dwarfs show on average dex higher metallicities than the overall MZR. By reaching the outskirts of this massive cluster, the ISM-ICM interaction can produce star-burst events, that accelerate gas depletion in cluster dwarf galaxies, as compared to their isolated counterparts. In the same time, galactic winds should have started to get suppressed, preventing metal lost from the low-mass galaxies, that otherwise would have been expected to be very efficient. As RPS gets shearing the gas content (starting from the halo gas reservoir and the HI disk), the infall of pristine gas gets truncated, and the chemical enrichment follows a steeper path on the mass-metallicity plane. By the end of their star formation, the low-mass galaxies could have experienced a significant metallicity enhancement. Considering that this shutting-off in Coma could take yr (the crossing-time of this cluster), there is enough time also for nitrogen to get delivered to the ISM, yielding the observed trend in the N/O ratio. In Figure 13 we give a schematic representation of this scenario.

Moreover, low-mass galaxies could reach clusters already chemically enhanced up to some degree. Recent works (Mahajan et al., 2010; Porter et al., 2008) have found an increased star formation in galaxies falling into clusters along supercluster-scale filaments. Additionally, simulations (e.g Bekki, 2009) suggest that even moderately strong ram pressure, such as in group-like environments, could strip the hot gas halos of galaxies, with efficiency increasing at the low-mass regimes. Consequently, low-mass galaxies, could suffer gas exhaustion by both gas depletion through star formation and some stripping process such as strangulation, also in their ways towards clusters. Both mechanisms would lead to chemical enhancement (product of the increased star formation as compared to isolated dwarf galaxies in the first case, and of the lulled infall in the second). This scenario could explain the trends observed for the oxygen abundance and N/O ratio even in A1367 (a cluster of ), as well as the elevated metallicities observed for some galaxies at in both Coma and A1367.

Carter et al. (2002) have invoked pressure confinement by the intracluster medium to explain the cluster-centric radial gradient of the stellar metallicity for a sample of passive galaxies in Coma. Smith et al. (2009), however, have argued that the observed gradient should be interpreted as a trend in age rather than metallicity. Younger passively evolving dwarfs could have arisen by the transformation of infalling field late-type galaxies, indicating that the build-up of the passive population is an ongoing process in the outskirts of the cluster, related to environment-driven processes (regarding the infalling galaxy population in Coma, see also Aguerri et al., 2004).

In a subsequent work, Smith et al. (2011) have found that the cluster-centric age gradient for the red-sequence dwarfs in Coma is a global trend, and is not driven by the ongoing merger of the NGC 4839 group to the south west of Coma. These authors have compared their results with the predictions from simulated cluster assembly histories (the Millenium Simulation, Springel et al., 2005) and have argued that in order to reproduce the strength of the age gradient observed for the red-sequence dwarfs in Coma, either a dominant burst or a gradual decline in the star-formation rate has to be invoked; models with very abrupt quenching would lead to shallower age trend. These findings support the scenario proposed here: recent star-formation, within a massive cluster like Coma, could yield the chemical enhancement of dwarf cluster galaxies.

Finally, a mechanism for the chemical enrichment of cluster galaxies to be considered could be the presence of enriched inflows. Historically it has been thought that there was a remarkable uniformity in the metal abundance of the ICM () as a function of cluster mass, ranging from cooling-core to noncooling-core clusters (see Werner et al. (2008) for a review on the chemical enrichment in the ICM, and Schindler & Diaferio (2008) for a review on the processes proposed to explain this enrichment). However, recently, the spatially resolved analysis of the chemical composition of the ICM has revealed that this is not uniformly enriched in metals. ICM abundance gradients are common in clusters, showing a peak in the central region and a decline outwards (e.g. De Grandi et al., 2004; Leccardi & Molendi, 2008; Lovisari et al., 2011). The central metallicity of the ICM can reach even over solar values, associated to the presence of the brightest cluster galaxy (RPS could also play a role in the metallicity enhancement in cluster centers). The question arises as to whether infall of enriched material of the ICM could affect the metallicity of cluster galaxies. However, detailed modeling would be needed, regarding the cooling of the ICM and the accretion mechanisms, which is beyond the scope of the present work.

## 8 Summary and Conclusions

In this work we have studied the chemical history of low-mass SF galaxies in four clusters in the local Universe. The sample clusters belong to a semi-spheric shell of the local Universe ( -25 deg and ) and span a mass range from to . The regions studied cover the clusters’ core up to 3. We have been searching for the potential imprints of the cluster environment on the galaxy chemical enrichment.

We have used the latest SDSS spectroscopic release DR8. SF galaxies have been selected on the basis of their SDSS emission line fluxes, and a limit in magnitud has been applied to select dwarf galaxies. Considering low-mass galaxies, apperture biases are not expected to be important. We note that DR8 spectroscopic data have been corrected for the underlying stellar continuum; this being an important improvement when studing nebular gas properties. We have found that our SF dwarf galaxies show typical line ratios of normal H ii galaxies.

Gas-phase metallicities of the O/H and N/O ratio, have been derived carefully using different empirical and model calibrations. The accurate mass estimates provided by SDSS DR8 have been used, to derive the MZR of the cluster galaxies. Well defined sequences have been found in the MZ plane, and we have observed a decrease of the scatter when the correction for the SFR has been applied. The value derived for the slope of the MZR is in agreement with the predictions of hydrodynamic models, which use momentum-driven winds to reproduce the MZR. Well defined sequences have also been derived in the N/O versus mass plane.

For the more massive clusters of this sample, Coma and A1367, the galaxies located at cluster-centric distances , are preferentially located at the upper part of the global sequences of O/H and N/O versus mass. This increase in metallicity is mass dependent, being higher at the lower mass bins, and in the core of Coma reaches on average dex in O/H. This effect yields the flattening of the MZR for SF dwarf galaxies within the core of these massive clusters.

The metal enhancement of SF dwarfs in the cluster core has been found to be more important when considering the region of the most massive cluster Coma (). Despite the general good relation of local galaxy density with cluster-centric distance, this effect appears diluted in terms of local galaxy density, suggesting that the relevant parameter able to affect the chemical evolution of SF dwarf galaxies should be the presence of a dense ICM.

Finally, we have related the metallicity of our SF dwarf galaxies, with their HI mass content, derived using available 21 cm data. We then have compared with the predictions of the so-called “closed-box“ model and the normal HI content of isolated couterparts, and we have found that SF dwarf galaxies in the cores of A1367 and Coma should be suffering an important ram-pressure stripping.

We discuss that the properties of the ICM could be a key parameter to the chemical evolution of low-mass cluster galaxies. Efficient gas stripping (ram pressure or/and strangulation) and effective metal retention, due to the suppression of galactic winds, could lead to a different chemical enrichment scenario for cluster galaxies.

The present sample consists of four clusters, and although the trends observed could be insightful, meaningful coclusions could be taken out only by improving considerably statistics. In a future work we will investigate further the connection of the chemical enrichment of cluster galaxies with the properties of the ICM, for a larger set of clusters, sampling a wide range of X-ray luminosities.

We thank Dr. O. Dors and Dr. E. Perez-Montero for sharing their models to derive galaxy metallicity. We also thank the anonymous referee for useful suggestions that helped improving the paper. V.P. acknowledge financial support from the Spanish Ministerio de Ciencia e Innovación under grant FPU AP2006-04622. We also acknowledge financial support by the Spanish PNAYA project ESTALLIDOS (AYA2010-21887-C04-01) and CSD2006 00070 â1st Science with GTCâ from the CONSOLIDER 2010 programme of the Spanish MICINN. Funding for SDSS-III has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, and the U.S. Department of Energy Office of Science. The SDSS-III web site is http://www.sdss3.org/. SDSS-III is managed by the Astrophysical Research Consortium for the Participating Institutions of the SDSS-III Collaboration including the University of Arizona, the Brazilian Participation Group, Brookhaven National Laboratory, University of Cambridge, University of Florida, the French Participation Group, the German Participation Group, the Instituto de Astrofisica de Canarias, the Michigan State/Notre Dame/JINA Participation Group, Johns Hopkins University, Lawrence Berkeley National Laboratory, Max Planck Institute for Astrophysics, New Mexico State University, New York University, Ohio State University, Pennsylvania State University, University of Portsmouth, Princeton University, the Spanish Participation Group, University of Tokyo, University of Utah, Vanderbilt University, University of Virginia, University of Washington, and Yale University. This research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration.

## References

• Abadi et al. (1999) Abadi, M. G., Moore, B., & Bower, R. G. 1999, MNRAS, 308, 947
• Abramson et al. (2011) Abramson, A., Kenney, J. D. P., Crowl, H. H., et al. 2011, AJ, 141, 164
• Aguerri et al. (2004) Aguerri, J. A. L., Iglesias-Paramo, J., Vilchez, J. M., & Muñoz-Tuñón, C. 2004, AJ, 127, 1344
• Alloin et al. (1979) Alloin, D., Collin-Souffrin, S., Joly, M., & Vigroux, L. 1979, A&A, 78, 200
• Amorín et al. (2010) Amorín, R. O., Pérez-Montero, E., & Vílchez, J. M. 2010, ApJ, 715, L128
• Asplund et al. (2009) Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481
• Baldry et al. (2004) Baldry, I. K., Glazebrook, K., Brinkmann, J., et al. 2004, ApJ, 600, 681
• Baldwin et al. (1981) Baldwin, J. A., Phillips, M. M., & Terlevich, R. 1981, PASP, 93, 5
• Balogh et al. (1999) Balogh, M. L., Morris, S. L., Yee, H. K. C., Carlberg, R. G., & Ellingson, E. 1999, ApJ, 527, 54
• Balogh et al. (2004) Balogh, M., et al. 2004, MNRAS, 348, 1355
• Bekki et al. (2009) Bekki, K., Tsujimoto, T., & Chiba, M. 2009, ApJ, 692, L24
• Bekki (2009) Bekki, K. 2009, MNRAS, 399, 2221
• Berg et al. (2011) Berg, D. A., Skillman, E. D., & Marble, A. R. 2011, ApJ, 738, 2
• Bird et al. (1995) Bird, C. M., Davis, D. S., & Beers, T. C. 1995, AJ, 109, 920
• Blanton et al. (2001) Blanton, M. R., et al. 2001, AJ, 121, 2358
• Blanton et al. (2005a) Blanton, M. R., Schlegel, D. J., Strauss, M. A., et al. 2005, AJ, 129, 2562
• Blanton et al. (2005b) Blanton, M. R., Lupton, R. H., Schlegel, D. J., et al. 2005, ApJ, 631, 208
• Blanton & Roweis (2007) Blanton, M. R., & Roweis, S. 2007, AJ, 133, 734
• Bolzonella et al. (2010) Bolzonella, M., Kovač, K., Pozzetti, L., et al. 2010, A&A, 524, A76
• Boselli & Gavazzi (2006) Boselli, A., & Gavazzi, G. 2006, PASP, 118, 517
• Bresolin et al. (2009) Bresolin, F., Gieren, W., Kudritzki, R.-P., Pietrzyński, G., Urbaneja, M. A., & Carraro, G. 2009, ApJ, 700, 309
• Briel et al. (1992) Briel, U. G., Henry, J. P., & Boehringer, H. 1992, A&A, 259, L31
• Brinchmann et al. (2004) Brinchmann, J., Charlot, S., White, S. D. M., et al. 2004, MNRAS, 351, 1151
• Brüggen & De Lucia (2008) Brüggen, M., & De Lucia, G. 2008, MNRAS, 383, 1336
• Bruzual & Charlot (2003) Bruzual, G., & Charlot, S. 2003, MNRAS, 344, 1000
• Calura et al. (2009) Calura, F., Pipino, A., Chiappini, C., Matteucci, F., & Maiolino, R. 2009, A&A, 504, 373
• Calzetti et al. (1994) Calzetti, D., Kinney, A. L., & Storchi-Bergmann, T. 1994, ApJ, 429, 582
• Cardelli et al. (1989) Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245
• Carter et al. (2002) Carter, D., Mobasher, B., Bridges, T. J., et al. 2002, ApJ, 567, 772
• Cedrés et al. (2009) Cedrés, B., Iglesias-Páramo, J., Vílchez, J. M., Reverte, D., Petropoulou, V., & Hernández-Fernández, J. 2009, AJ, 138, 873
• Chung et al. (2009) Chung, A., van Gorkom, J. H., Kenney, J. D. P., Crowl, H., & Vollmer, B. 2009, AJ, 138, 1741
• Cooper et al. (2008) Cooper, M. C., Tremonti, C. A., Newman, J. A., & Zabludoff, A. I. 2008, MNRAS, 390, 245
• Cortese et al. (2004) Cortese, L., Gavazzi, G., Boselli, A., Iglesias-Paramo, J., & Carrasco, L. 2004, A&A, 425, 429
• Cortese et al. (2006) Cortese, L., Gavazzi, G., Boselli, A., Franzetti, P., Kennicutt, R. C., O’Neil, K., & Sakai, S. 2006, A&A, 453, 847
• Cortese et al. (2008a) Cortese, L., Gavazzi, G., & Boselli, A. 2008, MNRAS, 390, 1282
• Cortese et al. (2008b) Cortese, L., Minchin, R. F., Auld, R. R., et al. 2008, MNRAS, 383, 1519
• Cresci et al. (2010) Cresci, G., Mannucci, F., Maiolino, R., Marconi, A., Gnerucci, A., & Magrini, L. 2010, Nature, 467, 811
• Cresci et al. (2011) Cresci, G., Mannucci, F., Sommariva, V., et al. 2011, arXiv:1110.4408
• Dalcanton (2007) Dalcanton, J. J. 2007, ApJ, 658, 941
• Davé et al. (2010) Davé, R., Finlator, K., Oppenheimer, B. D., Fardal, M., Katz, N., Kereš, D., & Weinberg, D. H. 2010, MNRAS, 404, 1355
• Davé et al. (2011) Davé, R., Finlator, K., & Oppenheimer, B. D. 2011, MNRAS, 416, 1354
• De Grandi et al. (2004) De Grandi, S., Ettori, S., Longhetti, M., & Molendi, S. 2004, A&A, 419, 7
• Dekel et al. (2009) Dekel, A., Birnboim, Y., Engel, G., et al. 2009, Nature, 457, 451
• De Lucia et al. (2004) De Lucia, G., Kauffmann, G., & White, S. D. M. 2004, MNRAS, 349, 1101
• de Rossi et al. (2007) de Rossi, M. E., Tissera, P. B., & Scannapieco, C. 2007, MNRAS, 374, 323
• Donnelly et al. (1998) Donnelly, R. H., Markevitch, M., Forman, W., et al. 1998, ApJ, 500, 138
• Dors et al. (2011) Dors, O. L., Jr., Krabbe, A., Hägele, G. F., & Pérez-Montero, E. 2011, MNRAS, 415, 3616
• Dressler (1980) Dressler, A. 1980, ApJ, 236, 351
• Duc et al. (2001) Duc, P.-A., Cayatte, V., Balkowski, C., Thuan, T. X., Papaderos, P., & van Driel, W. 2001, A&A, 369, 763
• Edmunds (1990) Edmunds, M. G. 1990, MNRAS, 246, 678
• Ellison et al. (2009) Ellison, S. L., Simard, L., Cowan, N. B., Baldry, I. K., Patton, D. R., & McConnachie, A. W. 2009, MNRAS, 396, 1257
• Erb et al. (2006) Erb, D. K., Shapley, A. E., Pettini, M., Steidel, C. C., Reddy, N. A., & Adelberger, K. L. 2006, ApJ, 644, 813
• Erb (2008) Erb, D. K. 2008, ApJ, 674, 151
• Esteban et al. (2004) Esteban, C., Peimbert, M., García-Rojas, J., et al. 2004, MNRAS, 355, 229
• Finlator & Davé (2008) Finlator, K., & Davé, R. 2008, MNRAS, 385, 2181
• Finn et al. (2005) Finn, R. A., et al. 2005, ApJ, 630, 206
• Fumagalli et al. (2011) Fumagalli, M., Gavazzi, G., Scaramella, R., & Franzetti, P. 2011, A&A, 528, A46
• Gallazzi et al. (2005) Gallazzi, A., Charlot, S., Brinchmann, J., White, S. D. M., & Tremonti, C. A. 2005, MNRAS, 362, 41
• García-Rojas & Esteban (2007) García-Rojas, J., & Esteban, C. 2007, ApJ, 670, 457
• Gavazzi et al. (2001) Gavazzi, G., Boselli, A., Mayer, L., Iglesias-Paramo, J., Vílchez, J. M., & Carrasco, L. 2001, ApJ, 563, L23
• Gavazzi et al. (2003) Gavazzi, G., Cortese, L., Boselli, A., et al. 2003, ApJ, 597, 210
• Gavazzi et al. (2010) Gavazzi, G., Fumagalli, M., Cucciati, O., & Boselli, A. 2010, A&A, 517, A73
• Giovanelli et al. (2005) Giovanelli, R., Haynes, M. P., Kent, B. R., et al. 2005, AJ, 130, 2598
• Gonzalez Delgado et al. (1999) González Delgado, R.M., Leitherer, C., Heckman, T.M. 1999, ApJS, 125, 489
• Gunn & Gott (1972) Gunn, J. E., & Gott, J. R., III 1972, ApJ, 176, 1
• Guseva et al. (2001) Guseva, N.G., Izotov, Y.I., Papaderos, P., Chaffee, C.B., Foltz, C.B., Green, R.F., Fricke, K.J., Noeske, K.G. 2001, A&A, 378, 756
• Haines et al. (2006) Haines, C. P., La Barbera, F., Mercurio, A., Merluzzi, P., & Busarello, G. 2006, ApJ, 647, L21
• Haines et al. (2007) Haines, C. P., Gargiulo, A., La Barbera, F., Mercurio, A., Merluzzi, P., & Busarello, G. 2007, MNRAS, 381, 7
• Haynes et al. (1984) Haynes, M. P., Giovanelli, R., & Chincarini, G. L. 1984, ARA&A, 22, 445
• Haynes et al. (2011) Haynes, M. P., Giovanelli, R., Martin, A. M., et al. 2011, arXiv:1109.0027
• Huang & Sarazin (1996) Huang, Z., & Sarazin, C. L. 1996, ApJ, 461, 622
• Iglesias-Páramo et al. (2002) Iglesias-Páramo, J., Boselli, A., Cortese, L., Vílchez, J. M., & Gavazzi, G. 2002, A&A, 384, 383
• Iglesias-Páramo et al. (2003) Iglesias-Páramo, J., Boselli, A., Gavazzi, G., Cortese, L., & Vílchez, J. M. 2003, A&A, 397, 421
• Israel (1997) Israel, F. P. 1997, A&A, 328, 471
• Jaffé et al. (2011) Jaffé, Y. L., Aragón-Salamanca, A., Kuntschner, H., et al. 2011, MNRAS, 417, 1996
• Kapferer et al. (2006) Kapferer, W., Ferrari, C., Domainko, W., et al. 2006, A&A, 447, 827
• Kapferer et al. (2008) Kapferer, W., Kronberger, T., Ferrari, C., Riser, T., & Schindler, S. 2008, MNRAS, 389, 1405
• Kapferer et al. (2009) Kapferer, W., Kronberger, T., Breitschwerdt, D., et al. 2009, A&A, 504, 719
• Kauffmann et al. (2003) Kauffmann, G., Heckman, T. M., White, S. D. M., et al. 2003, MNRAS, 341, 33
• Kenney & Koopmann (1999) Kenney, J. D. P., & Koopmann, R. A. 1999, AJ, 117, 181
• Kennicutt (1998) Kennicutt, R. C., Jr. 1998, ApJ, 498, 541
• Kennicutt et al. (2003) Kennicutt, R. C., Jr., Bresolin, F., & Garnett, D. R. 2003, ApJ, 591, 801
• Keres et al. (2003) Keres, D., Yun, M. S., & Young, J. S. 2003, ApJ, 582, 659
• Kewley & Dopita (2002) Kewley, L. J., & Dopita, M. A. 2002, ApJS, 142, 35
• Kewley et al. (2005) Kewley, L. J., Jansen, R. A., & Geller, M. J. 2005, PASP, 117, 227
• Kewley & Ellison (2008) Kewley, L. J., & Ellison, S. L. 2008, ApJ, 681, 1183
• Kewley et al. (2010) Kewley, L. J., Rupke, D., Zahid, H. J., Geller, M. J., & Barton, E. J. 2010, ApJ, 721, L48
• Kobulnicky & Skillman (1997) Kobulnicky, H. A., & Skillman, E. D. 1997, ApJ, 489, 636
• Koopmann et al. (2006) Koopmann, R. A., Haynes, M. P., & Catinella, B. 2006, AJ, 131, 716
• Koopmann & Kenney (2004) Koopmann, R. A., & Kenney, J. D. P. 2004, ApJ, 613, 866
• Lara-López et al. (2010) Lara-López, M. A., et al. 2010, arXiv:1005.0509
• Larson (1974) Larson, R. B. 1974, MNRAS, 169, 229
• Leccardi & Molendi (2008) Leccardi, A., & Molendi, S. 2008, A&A, 487, 461
• Ledlow et al. (2003) Ledlow, M. J., Voges, W., Owen, F. N., & Burns, J. O. 2003, AJ, 126, 2740
• Lee et al. (2003) Lee, H., McCall, M. L., & Richer, M. G. 2003, AJ, 125, 2975
• Lee et al. (2006) Lee, H., Skillman, E. D., Cannon, J. M., et al. 2006, ApJ, 647, 970
• Leitherer et al. (1999) Leitherer, C., Schaerer, D., Goldader, J. D., et al. 1999, ApJS, 123, 3
• Levy et al. (2007) Levy, L., Rose, J. A., van Gorkom, J. H., & Chaboyer, B. 2007, AJ, 133, 1104
• Lovisari et al. (2011) Lovisari, L., Schindler, S., & Kapferer, W. 2011, A&A, 528, A60
• Magrini et al. (2007) Magrini, L., Vílchez, J. M., Mampaso, A., Corradi, R. L. M., & Leisy, P. 2007, A&A, 470, 865
• Mahajan et al. (2010) Mahajan, S., Haines, C. P., & Raychaudhury, S. 2010, MNRAS, 404, 1745
• Mahajan et al. (2011) Mahajan, S., Haines, C. P., & Raychaudhury, S. 2011, MNRAS, 412, 1098
• Maiolino et al. (2008) Maiolino, R., et al. 2008, A&A, 488, 463
• Mannucci et al. (2009) Mannucci, F., Cresci, G., Maiolino, R., et al. 2009, MNRAS, 398, 1915
• Mannucci et al. (2010) Mannucci, F., Cresci, G., Maiolino, R., Marconi, A., & Gnerucci, A. 2010, MNRAS, 1314
• Marconi et al. (1994) Marconi, G., Matteucci, F., & Tosi, M. 1994, MNRAS, 270, 35
• Martín-Manjón et al. (2008) Martín-Manjón, M. L., Mollá, M., Díaz, A. I., & Terlevich, R. 2008, MNRAS, 385, 854
• Martín-Manjón et al. (2010) Martín-Manjón, M. L., García-Vargas, M. L., Mollá, M., & Díaz, A. I. 2010, MNRAS, 403, 2012
• Matteucci (1994) Matteucci, F. 1994, A&A, 288, 57
• McCarthy et al. (2008) McCarthy, I. G., Frenk, C. S., Font, A. S., et al. 2008, MNRAS, 383, 593
• Merlin & Chiosi (2006) Merlin, E., & Chiosi, C. 2006, A&A, 457, 437
• Mesa-Delgado et al. (2009) Mesa-Delgado, A., Esteban, C., García-Rojas, J., et al. 2009, MNRAS, 395, 855
• Meynet & Maeder (2002) Meynet, G., & Maeder, A. 2002, A&A, 390, 561
• McGaugh (1991) McGaugh, S. S. 1991, ApJ, 380, 140
• Michel-Dansac et al. (2008) Michel-Dansac, L., Lambas, D. G., Alonso, M. S., & Tissera, P. 2008, MNRAS, 386, L82
• Montuori et al. (2010) Montuori, M., Di Matteo, P., Lehnert, M. D., Combes, F., & Semelin, B. 2010, A&A, 518, A56
• Mori & Burkert (2000) Mori, M., & Burkert, A. 2000, ApJ, 538, 559
• Mouhcine et al. (2007) Mouhcine, M., Baldry, I. K., & Bamford, S. P. 2007, MNRAS, 382, 801
• Mouhcine et al. (2011) Mouhcine, M., Kriwattanawong, W., & James, P. A. 2011, MNRAS, 412, 1295
• Nagao et al. (2006) Nagao, T., Maiolino, R., & Marconi, A. 2006, A&A, 459, 85
• Neumann et al. (2003) Neumann, D. M., Lumb, D. H., Pratt, G. W., & Briel, U. G. 2003, A&A, 400, 811
• Oppenheimer et al. (2010) Oppenheimer, B. D., Davé, R., Kereš, D., et al. 2010, MNRAS, 406, 2325
• Osterbrock (1989) Osterbrock, D. E. 1989, Research supported by the University of California, John Simon Guggenheim Memorial Foundation, University of Minnesota, et al. Mill Valley, CA, University Science Books, 1989, 422 p.,
• Pagel et al. (1979) Pagel, B. E. J., Edmunds, M. G., Blackwell, D. E., Chun, M. S., & Smith, G. 1979, MNRAS, 189, 95
• Peimbert & Peimbert (2010) Peimbert, A., & Peimbert, M. 2010, ApJ, 724, 791
• Penny & Conselice (2011) Penny, S. J., & Conselice, C. J. 2011, EAS Publications Series, 48, 197
• Pérez-Montero & Contini (2009) Pérez-Montero, E., & Contini, T. 2009, MNRAS, 398, 949
• Pérez-Montero & Díaz (2005) Pérez-Montero, E., & Díaz, A. I. 2005, MNRAS, 361, 1063
• Petropoulou et al. (2011) Petropoulou, V., Vílchez, J., Iglesias-Páramo, J., Papaderos, P., Magrini, L., Cedrés, B., & Reverte, D. 2011, ApJ, 734, 32
• Pettini & Pagel (2004) Pettini, M., & Pagel, B. E. J. 2004, MNRAS, 348, L59
• Pilyugin & Thuan (2005) Pilyugin, L. S., & Thuan, T. X. 2005, ApJ, 631, 231
• Pilyugin et al. (2010) Pilyugin, L. S., Vílchez, J. M., & Thuan, T. X. 2010, ApJ, 720, 1738
• Poggianti et al. (1999) Poggianti, B. M., Smail, I., Dressler, A., Couch, W. J., Barger, A. J., Butcher, H., Ellis, R. S., & Oemler, A., Jr. 1999, ApJ, 518, 576
• Poggianti et al. (2004) Poggianti, B. M., Bridges, T. J., Komiyama, Y., et al. 2004, ApJ, 601, 197
• Poggianti et al. (2006) Poggianti, B. M., von der Linden, A., De Lucia, G., et al. 2006, ApJ, 642, 188
• Poggianti et al. (2008) Poggianti, B. M., Desai, V., Finn, R., et al. 2008, ApJ, 684, 888
• Poggianti et al. (2009a) Poggianti, B. M., Aragón-Salamanca, A., Zaritsky, D., et al. 2009, ApJ, 693, 112
• Poggianti et al. (2009b) Poggianti, B. M., et al. 2009, ApJ, 697, L137
• Porter et al. (2008) Porter, S. C., Raychaudhury, S., Pimbblet, K. A., & Drinkwater, M. J. 2008, MNRAS, 388, 1152
• Relaño et al. (2010) Relaño, M., Monreal-Ibero, A., Vílchez, J. M., & Kennicutt, R. C. 2010, MNRAS, 402, 1635
• Rines et al. (2005) Rines, K., Geller, M. J., Kurtz, M. J., & Diaferio, A. 2005, AJ, 130, 1482
• Robertson et al. (2011) Robertson, P., Shields, G. A., & Blanc, G. A. 2011, arXiv:1112.5975
• Roediger & Hensler (2005) Roediger, E., & Hensler, G. 2005, A&A, 433, 875
• Rupke et al. (2010) Rupke, D. S. N., Kewley, L. J., & Barnes, J. E. 2010, ApJ, 710, L156
• Salim et al. (2007) Salim, S., Rich, R. M., Charlot, S., et al. 2007, ApJS, 173, 267
• Schindler et al. (2005) Schindler, S., Kapferer, W., Domainko, W., et al. 2005, A&A, 435, L25
• Schindler & Diaferio (2008) Schindler, S., & Diaferio, A. 2008, Space Sci. Rev., 134, 363
• Simón-Díaz & Stasińska (2011) Simón-Díaz, S., & Stasińska, G. 2011, A&A, 526, A48
• Shimasaku et al. (2001) Shimasaku, K., Fukugita, M., Doi, M., et al. 2001, AJ, 122, 1238
• Smith et al. (2009) Smith, R. J., Lucey, J. R., Hudson, M. J., Allanson, S. P., Bridges, T. J., Hornschemeier, A. E., Marzke, R. O., & Miller, N. A. 2009, MNRAS, 392, 1265
• Smith et al. (2010) Smith, R. J., et al. 2010, MNRAS, 408, 1417
• Smith et al. (2011) Smith, R. J., Lucey, J. R., Price, J., Hudson, M. J., & Phillipps, S. 2011, arXiv:1108.3836
• Skillman et al. (1996) Skillman, E. D., Kennicutt, R. C., Jr., Shields, G. A., & Zaritsky, D. 1996, ApJ, 462, 147
• Solanes et al. (2001) Solanes, J. M., Manrique, A., García-Gómez, C., et al. 2001, ApJ, 548, 97
• Spitoni et al. (2010) Spitoni, E., Calura, F., Matteucci, F., & Recchi, S. 2010, A&A, 514, A73
• Springel et al. (2005) Springel, V., White, S. D. M., Jenkins, A., et al. 2005, Nature, 435, 629
• Stasińska et al. (2007) Stasińska, G., Tenorio-Tagle, G., Rodríguez, M., & Henney, W. J. 2007, A&A, 471, 193
• Storchi-Bergmann et al. (1994) Storchi-Bergmann, T., Calzetti, D., & Kinney, A. L. 1994, ApJ, 429, 572
• Strateva et al. (2001) Strateva, I., Ivezić, Ž., Knapp, G. R., et al. 2001, AJ, 122, 1861
• Struble & Rood (1999) Struble, M. F., & Rood, H. J. 1999, ApJS, 125, 35
• Sun & Murray (2002) Sun, M., & Murray, S. S. 2002, ApJ, 576, 708
• Sun et al. (2007) Sun, M., Donahue, M., & Voit, G. M. 2007, ApJ, 671, 190
• Tacconi et al. (2010) Tacconi, L. J., Genzel, R., Neri, R., et al. 2010, Nature, 463, 781
• Tassis et al. (2008) Tassis, K., Kravtsov, A. V., & Gnedin, N. Y. 2008, ApJ, 672, 888
• Tecce et al. (2010) Tecce, T. E., Cora, S. A., Tissera, P. B., Abadi, M. G., & Lagos, C. D. P. 2010, MNRAS, 408, 2008
• Tecce et al. (2011) Tecce, T. E., Cora, S. A., & Tissera, P. B. 2011, MNRAS, 416, 3170
• Terlevich et al. (2001) Terlevich, A. I., Caldwell, N., & Bower, R. G. 2001, MNRAS, 326, 1547
• Tissera et al. (2005) Tissera, P. B., De Rossi, M. E., & Scannapieco, C. 2005, MNRAS, 364, L38
• Tolstoy et al. (2009) Tolstoy, E., Hill, V., & Tosi, M. 2009, ARA&A, 47, 371
• Tonnesen et al. (2007) Tonnesen, S., Bryan, G. L., & van Gorkom, J. H. 2007, ApJ, 671, 1434
• Tonnesen & Bryan (2008) Tonnesen, S., & Bryan, G. L. 2008, ApJ, 684, L9
• Tonnesen & Bryan (2009) Tonnesen, S., & Bryan, G. L. 2009, ApJ, 694, 789
• Toribio et al. (2011) Toribio, M. C., Solanes, J. M., Giovanelli, R., Haynes, M. P., & Martin, A. M. 2011, ApJ, 732, 93
• Tremonti et al. (2004) Tremonti, C. A., et al. 2004, ApJ, 613, 898
• Treu et al. (2003) Treu, T., Ellis, R. S., Kneib, J.-P., Dressler, A., Smail, I., Czoske, O., Oemler, A., & Natarajan, P. 2003, ApJ, 591, 53
• Vaduvescu et al. (2011) Vaduvescu, O., Kehrig, C., Vilchez, J. M., & Unda-Sanzana, E. 2011, A&A, 533, A65
• Vaduvescu et al. (2007) Vaduvescu, O., McCall, M. L., & Richer, M. G. 2007, AJ, 134, 604
• van Zee & Haynes (2006) van Zee, L., & Haynes, M. P. 2006, ApJ, 636, 214
• van Zee et al. (2006) van Zee, L., Skillman, E. D., & Haynes, M. P. 2006, ApJ, 637, 269
• Vilchez (1995) Vilchez, J. M. 1995, AJ, 110, 1090
• Vulcani et al. (2011) Vulcani, B., Poggianti, B. M., Fasano, G., et al. 2011, arXiv:1111.0832
• Werner et al. (2008) Werner, N., Durret, F., Ohashi, T., Schindler, S., & Wiersma, R. P. C. 2008, Space Sci. Rev., 134, 337
• White et al. (1997) White, D. A., Jones, C., & Forman, W. 1997, MNRAS, 292, 419
• Yagi et al. (2010) Yagi, M., et al. 2010, AJ, 140, 1814
You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters