The following article is Open access

Exploring the AGN–Ram Pressure Stripping Connection in Local Clusters

, , , , , , , , , , and

Published 2022 March 9 © 2022. The Author(s). Published by the American Astronomical Society.
, , Citation Giorgia Peluso et al 2022 ApJ 927 130 DOI 10.3847/1538-4357/ac4225

Download Article PDF
DownloadArticle ePub

You need an eReader or compatible software to experience the benefits of the ePub3 file format.

0004-637X/927/1/130

Abstract

Ram pressure stripping (RPS) by the intracluster medium is one of the most advocated mechanisms that affect the properties of cluster galaxies. A recent study based on a small sample has found that many galaxies showing strong signatures of RPS also possess an active galactic nucleus (AGN), suggesting a possible correlation between the two phenomena. This result has not been confirmed by a subsequent study. Building upon previous findings, here we combine MUSE observations conducted within the GASP program and a general survey of the literature to robustly measure the AGN fraction in ram-pressure-stripped cluster galaxies using Baldwin–Phillips–Terlevich emission line diagrams. Considering a sample of 115 ram-pressure-stripped galaxies with stellar masses ≥ 109M, we find an AGN fraction of ∼27%. This fraction strongly depends on stellar mass: it raises to 51% when only ram-pressure-stripped galaxies of masses M* ≥ 1010 M are considered. We then investigate whether the AGN incidence is in excess in ram-pressure-stripped galaxies compared to nonstripped galaxies using as a comparison a sample of noncluster galaxies observed by the MaNGA survey. Considering mass-matched samples, we find that the incidence of AGN activity is significantly higher (at a confidence level >99.95%) when RPS is in the act, supporting the hypothesis of an AGN–ram pressure connection.

Export citation and abstract BibTeX RIS

Original content from this work may be used under the terms of the Creative Commons Attribution 4.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.

1. Introduction

Both theoretical and observational studies concur that there is a strong connection between the presence of an active galactic nucleus (AGN) and the host galaxy properties (see Kormendy et al. 2013, and references therein), suggesting that internal processes might regulate the AGN activity and, conversely, the AGN activity might be relevant for shaping galaxy properties. AGNs are preferentially found in more massive galaxies (M* > 109 M; see, e.g., Kauffmann et al. 2003; Decarli et al. 2007; Juneau et al. 2011; Pimbblet et al. 2013; Sabater et al. 2013; Lopes et al. 2017; Rodríguez del Pino et al. 2017; Sánchez et al. 2018) and the mass of the host galaxy is the main parameter driving the level of AGN activity (Magliocchetti et al. 2020). However, it is still debated if other factors, such as dense galaxy environment or galaxy clusters (e.g., Pimbblet et al. 2013) have an impact on the presence of AGNs in galaxies.

Despite the vast literature on this topic (e.g., Kauffmann et al. 2004; Best et al. 2007; Silverman et al. 2009; von der Linden et al. 2010; Hwang et al. 2012; Martini et al. 2013; Sabater et al. 2013; Ehlert et al. 2014; Silverman & David 2015; Coldwell et al. 2017; Lopes et al. 2017; Marziani et al. 2017; Argudo-Fernández et al. 2018; Gordon et al. 2018; Koulouridis et al. 2018; Magliocchetti et al. 2018), different studies have reached quite opposite results, most likely due to the different techniques adopted to identify AGNs, select the samples, and characterize the environment. Using a spectroscopic sample, Dressler et al. (1985) first suggested that the fraction of AGNs in clusters (∼1%) is significantly lower than in the field (∼5%). Similarly, Lopes et al. (2017), identifying AGNs in the Sloan Digital Sky Survey (SDSS) using optical emission lines and Baldwin–Phillips–Terlevich (BPT; Baldwin et al. 1981) diagrams, found that AGNs favor environments typical of the field, low-mass groups, or cluster outskirts. Using the same data set but a different cluster sample, von der Linden et al. (2010) showed instead that the AGN fraction does not change as a function of environment, nor of clustercentric distance (see also Miller et al. 2003). Similar conclusions were obtained by Martini et al. (2007), Lehmer et al. (2007), Sivakoff et al. (2008), and Arnold et al. (2009), exploiting X-ray data. Yet, the radio AGN fraction seems to be much higher in clusters than in the field (Best et al. 2007; Sabater et al. 2013).

Considering local density as a proxy for environment, Kauffmann et al. (2004) found that AGN host galaxies with strong [O iii] emission are twice as frequent in low-density regions than in high-density regions (see also Miller et al. 2003; Montero-Dorta et al. 2009). In contrast, Amiri et al. (2019) did not find any effect of the galaxy density on nuclear activity. Sabater et al. (2013) found that (at fixed mass) the prevalence of optical AGNs is a factor of 2–3 lower in the densest environments (see also Man et al. 2019), but increases by a factor of ∼2 in the presence of strong one-on-one interactions. Gilmour et al. (2007) showed that X-ray-selected AGNs lie predominantly in moderate dense regions.

The expected connection between AGN incidence and properties with the environment has roots in the fact that the characteristics of AGNs are strongly linked to the conditions of the available gas, which in turn can be affected by the galaxy environment. So any environmentally specific physical mechanism that has the potential to affect the galaxy gas can impact the AGN activity. For example, mergers—which most frequently happen in the field—have frequently been cited as a method to fuel AGNs (e.g., Sanders et al. 1988) and a number of morphological studies claim an excess of post-merger systems in their AGN samples (Bahcall et al. 1997; Canalizo & Stockton 2001; Urrutia et al. 2008; Letawe et al. 2010; Smirnova et al. 2010).

Another process able to affect the gas supply in galaxies is ram pressure stripping (RPS; Gunn et al. 1972). This mechanism happens most efficiently in clusters and massive groups (Hester 2006) and it is due to the pressure exerted by the intracluster medium (ICM) on the galaxy interstellar medium (ISM). This interaction can produce many visible effects on the galaxy, such as altering its less-bound gas, giving rise to wakes of stripped material departing from the main galaxy body (Kenney et al. 2004; van Gorkom 2004; Fumagalli et al. 2014; Poggianti et al. 2017a) and inducing a quenching of star formation (Vollmer et al. 2001; Tonnesen et al. 2007; Vulcani et al. 2020). Prior to complete gas removal, it has been observed that ram pressure can also increase the star formation rate in galaxies (Crowl & Kenney 2006; Merluzzi et al. 2013; Vulcani et al. 2018) and simulations support this finding (Kronberger et al. 2008; Kapferer et al. 2009; Tonnesen et al. 2009; Bekki 2014), suggesting that the increased pressure initially helps compress the gas and triggers increased star formation. The same mechanism that initially promotes star formation can also fuel the AGN during the RPS process: gas can be funneled toward the galaxy centers due to gravitational instabilities and the spiraling toward the center of clumps that lose angular momentum (Schulz & Struck 2001; Tonnesen et al. 2009; Ramos-Martínez et al. 2018). The funneling of gas toward the galaxy center can also ignite the central supermassive black hole. Theoretical models (Tonnesen et al. 2009) have indeed demonstrated that gas inflows can fuel the central AGN in ram-pressure-stripped galaxies, possibly due to the presence of magnetic fields (Ramos-Martínez et al. 2018). The enhanced accretion onto the black hole can then produce heating and outflows due to AGN feedback (Ricarte et al.2020).

Therefore RPS might be simultaneously responsible for an enhanced AGN activity and the appearance of tails of stripped material. This scenario has been first proposed by Poggianti et al. (2017b, hereafter P17b) to explain the very high incidence (6/7) of AGN detected in a sample of galaxies strongly affected by RPS, also called jellyfish galaxies (see also Maier et al. 2022). That analysis is based on integral-field spectroscopic data coming from the GAs Stripping Phenomena (GASP) in galaxies (Poggianti et al. 2017a). Subsequent GASP studies have led to the identification of AGN-driven outflows (Radovich et al. 2019) and a compelling case for AGN feedback in action (George et al. 2019).

The P17b analysis is based on a small sample and importantly is only composed of jellyfish galaxies with very striking tails and all massive galaxies. Thus, it leaves open the possibility that the AGN activity could be only related to the (rather short) peak phase of stripping and/or only to the galaxy mass regardless of RPS.

A subsequent study that did not find a high incidence of AGNs (Roman-Oliveira et al. 2018) 6 analyzed a sample of ram-pressure-stripped galaxies in a supercluster at z ∼ 0.2 and found only 5/70 AGN, according to optical line diagnostics. Their sample spans a wide galaxy stellar mass range (from 109 to 1011.5 M) and is based on visual identification of the candidates. At odds with GASP, none of these candidates have IFS data to confirm they are indeed affected by RPS. It also includes RPS candidates with different degrees of stripping, while as said above the P17b study includes only very dramatic cases.

In this paper, we build on previous results and aim at estimating the incidence of AGN on the largest possible sample of ram-pressure-stripped galaxies to date.

Throughout the paper, we adopt a Chabrier (2003) initial mass function (IMF) in the mass range 0.1–100 M. The cosmological constants assumed are Ωm = 0.3, ΩΛ = 0.7, and H0 = 70 km s−1 Mpc−1.

2. Data sets and Galaxy Samples

In this paper we aim at characterizing the incidence of AGN activity among ram-pressure-stripped cluster galaxies. We start our analysis by considering a sample of galaxies drawn from the GASP (Poggianti et al. 2017a) survey, then, to increase the statistics and more robustly support the results, we also gather a literature sample of all ram-pressure-stripped galaxies identified by different authors in the last four decades. Finally, we exploit the fifteenth data release of the MaNGA survey (DR15; Aguado et al. 2019) to build a control sample of galaxies with characteristics similar to the ram-pressure-stripped galaxies, but considering only those not in clusters and therefore presumably not affected by strong RPS.

2.1. The GASP Sample and Data

GASP is a project aimed at studying gas removal processes, mainly due to ICM–ISM interaction, using a sample of 114 galaxies. More specifically, it comprises both RPS candidates and undisturbed galaxies located in clusters, groups, and field, spanning a range in stellar masses from 109 to 3.2 × 1011 M and a redshift range of 0.04 < z < 0.07. All galaxies were selected on the basis of B-band imaging coming from three different surveys: WINGS (Fasano et al. 2006), OMEGAWINGS (Gullieuszik et al. 2015), and PM2GC (Calvi et al.2011).

GASP is based on an ESO Large Program carried out with the integral-field spectrograph MUSE, mounted at the VLT, whose large field of view ($1^{\prime} \times 1^{\prime} $) and high, but seeing-limited, spatial resolution (0farcs2 pixel−1, seeing of 1''), allow us to cover the galaxy outskirts and possible tails of gas departing from the main body of the galaxies up to ten times the galaxy effective radius (i.e., ∼10 Re ) with a resolution of ∼1 kpc at the galaxy redshifts.

For our analysis, we select only cluster members that have been confirmed to be ram pressure stripped based on the MUSE data; in fact, they all have extraplanar Hα emission in various stages of stripping (B. Poggianti et al. 2021, in preparation), from weak/initial stripping (JStage = 0.5) to significant tails (JStage = 1) to extreme tails longer than the stellar disk diameter (JStage = 2, so-called jellyfish galaxies) to truncated disks corresponding to a late stage of RPS (Jstage = 3), for a total of 51 galaxies. From now on we will call this sample GASP-RPS. All of these are morphologically late-type and star-forming galaxies.

We make use of the fluxes of the emission-only component of the lines Hα, Hβ, [O iii]λ 5007 Å and [N ii]λ 6583Å measured with the KUBEVIZ code (Fossati et al. 2016) from the continuum-subtracted MUSE cubes corrected for both Galactic and intrinsic extinction, as described in detail in Poggianti et al. (2017a). Stellar masses are taken from Vulcani et al. (2018) and are computed using the SINOPSIS spectrophotometric code (Fritz et al. 2017) by summing up the masses of all the spaxels within the galaxy disk (Gullieuszik et al. 2020).

To characterize the ionization mechanism acting on the gas and therefore identify galaxies with AGNs, we inspect the BPT diagnostic diagram [N ii]/Hα versus [O iii]/Hβ ratios (BPT-N ii; Baldwin et al. 1981). We consider only spaxels with a signal-to-noise (S/N) ratio greater than 3 for all the lines used. We use the relation from Kauffmann et al. (2003) (K03) to separate star-forming from composite regions, the Kewley et al. (2001) (K01) line to identify AGNs, 7 and the Sharp et al. (2010) (SB10) relation to discriminate between Seyferts and LINERs. We classify a galaxy as an AGN host if in its central (3'') region there are at least 20 spaxels 8 that have a Seyfert or LINER classification, otherwise we flag it as star-forming. In the galaxies with high values of extinction (AV , as measured by the Balmer decrement) in the central cores, that might prevent us from identifying an AGN, we further inspect the LINER classified spaxels: a biconical shape of their distribution suggests extended ionized regions and therefore indicates the presence of the AGN.

For those galaxies with a central LINER/AGN classification, we have carefully checked the emission line fits, in particular the Hβ line given that an underestimate of Hβ flux would lead to an overestimate of the [O iii]/Hβ ratio, mimicking line ratios typical of AGNs.

2.2. Ram Pressure Stripping Candidates from the Literature

We have performed a systematic literature search of all the ram-pressure-stripped galaxies identified by 2020 December. These galaxies were studied exploiting a wide variety of observational techniques, including radio (e.g., Gavazzi et al. 1995), submillimeter (e.g., Scott et al. 2013; Jáchym et al. 2014, 2019), infrared (e.g., Sivanandam et al. 2010, 2014), optical (e.g., Gavazzi et al. 1995, 2001; Sun et al. 2007; Sivanandam et al. 2010; Sun et al. 2010; Yagi et al. 2010; Fumagalli et al. 2014; Fossati et al. 2016; Gavazzi et al. 2017; Roberts et al. 2020), UV (e.g., Smith et al. 2010), and X-ray (e.g., Sun et al. 2006, 2010). The assembled sample is therefore greatly heterogeneous and while for some galaxies it has been confirmed that RPS is the only acting mechanism, in some other cases galaxies are most likely undergoing both RPS and tidal interactions. As our aim is to include all ram-pressure-stripped galaxies and collect a sample as large as possible, we consider also the latter cases. However, we remove cases where a merger or tidal interaction is the main cause of the galaxy transformation. 9

We narrow down our search to galaxies for which we retrieve information at any wavelength on the ionization mechanism of the central emission, 10 obtaining a total sample of 80 galaxies (from now on LIT-RPS sample). All these turn out to have some active star formation (in addition to the eventual AGN activity) in the available literature.

The LIT-RPS sample is located in the redshift range 0.001 ≤ z ≤ 0.34, plus a galaxy at z = 0.73, and covers a stellar mass range of 1.3 × 108 < M*/M < 2.0 × 1011. Stellar masses have been collected from the literature and homogenized to the same Chabrier (2003) IMF (as in GASP). When a stellar mass estimate was not available (4/80), we computed it using the available photometric data following the Bell et al. (2001) approach, as described in the Appendix.

To assess the strength of RPS signatures and compare them with the GASP JStage classification, four of us (GP, BMP, BV, AM) visually inspected the available images in the literature for all the galaxies. Following the scheme described in Section 2.1, we assign a flag indicating the extent of the tail (JStage) based on the Hα emission (if available) and also a general JStage based on any wavelength observed (JStagegen). In the case of multiple images with different resolutions or at different wavelengths showing a different extent of the tail, we always consider the wavelength with the longest visible tail to assign the JStagegen. The classifiers agreed in most of the cases. In the discrepant cases, each galaxy was inspected together by the classifiers to ensure homogeneity and to find a consensus. This visual inspection also confirmed that the LIT-RPS sample is composed of morphologically late-type galaxies (spirals or irregulars).

2.3. The MaNGA Sample

Mapping Nearby Galaxies at Apache Point Observatory (MaNGA; Bundy 2015) is an integral-field spectroscopic survey observing galaxies at 0.01 ≤ z ≤ 0.15 using the BOSS Spectrograph (Smee et al. 2013) mounted at the 2.5 m SDSS telescope (Gunn et al. 2006), which covers a spectral range from 3600 to 10300 Å, with a resolution of R ∼ 1400 at 4000 Å and R ∼ 2600 at 9000 Å.

We exploit the MaNGA DR15 release and use the outputs of the Pipe3D pipeline (Sánchez et al. 2016a, 2016b, 2018). More specifically, we use the Pipe3D-v2_4_3 11 catalog, which contains integrated properties, characteristics, and gradients of different quantities for 4656 galaxies. Of interest for our work are integrated stellar masses and star formation rates obtained from the Hα emission line that we convert to our adopted Chabrier (2003) IMF.

We first exclude from the sample 75 duplicate galaxies and then select only galaxies with a specific star formation rate (sSFR) > 10−11yr−1, for a total of 2509 galaxies. The latter selection allows us to consider only star-forming galaxies, as are galaxies in both the GASP-RPS and LIT-RPS samples.

As we aim at assembling a sample not affected by RPS, we crossmatch our sample with the environmental catalog by Tempel et al. (2014), who provide halo mass estimates based on Navarro et al. (1997) profiles. Using a searching radius of 5'', we obtain a match for 2061 galaxies, 861 of which are located in structures with halo masses $\mathrm{log}({M}_{{\rm{h}}}/{\text{}}{M}_{\odot })\lt 13.0$, 12 therefore are most likely isolated (Yang et al. 2007).

Finally, to reduce the effect of a different morphological mix among the different samples, we use the visual morphological classification from the MaNGA Value Added Catalogs 13 that is based on inspection of image mosaics using a new reprocessing of SDSS and Dark Energy Legacy Survey (DESI) images, following the methods from Hernández-Toledo et al. (2010) and exclude 70 early-type (Ellipticals, S0s, and S0as) and two unclassified galaxies. 14

To assemble the final MaNGA sample, we consider only galaxies that in a circular aperture of 3'' diameter (i.e., the SDSS fiber size) centered on the galaxy have at least 20 spaxels with S/N > 3 for all lines that will be used to detect the presence of an AGN. 782 galaxies pass this selection and constitute our reference sample, called MaNGA-Ref. This sample covers a redshift range 0.0024 < z < 0.1439. We note that the MaNGA fiber core diameter (2'' ) is similar to the typical seeing value (2farcs5). At the median redshift of our MaNGA-Ref sample (z = 0.0317), the MaNGA spatial resolution of 2'' corresponds to 1.27 kpc and at the 75% redshift percentile (z = 0.043) to 1.7 kpc. This is only slightly worse than the GASP spatial resolution of ∼1 kpc (median z = 0.05), which is dominated by the seeing (∼1''). 15

To identify AGN, we inspect the BPT-N ii maps provided by the online tool MARVIN 16 and use the same classification criteria as for the RPS samples. We count the number of spaxels classified as AGN (i.e., Seyfert + LINER), star-forming (SF), or composite. If the number of spaxels classified as Seyfert or LINER is larger than 20 in a circular aperture of ${3}^{{\prime} ^{\prime} }$ diameter, we classify the galaxy as AGN, otherwise as SF.

3. Results I: The Incidence of AGN among Ram-pressure-stripped Galaxies

In this section we present the sample of AGN hosts in the GASP-RPS and in the LIT-RPS samples separately and quantify the incidence of AGNs among ram-pressure-stripped galaxies. We also investigate if these fractions depend on the properties of the ram-pressure-stripped galaxies, such as stellar mass and Jstage. In the following section we will quantitatively compare these fractions controlling for the different mass distribution to those of the MaNGA-Ref sample.

3.1. GASP-RPS

In the GASP-RPS sample, seven galaxies are already known to host an AGN; six of them were presented in P17b and one (JO36), in Fritz et al. (2017). The latter is an edge-on disk hosting an obscured AGN which is not directly identified using BPT diagrams due to strong dust absorption. However, evidence for the AGN presence comes from extra-nuclear LINER-like emission with a cone morphology and the AGN is detected by Chandra as a point-like X-ray source (Fritz et al.2017).

Among the P17b candidates, JO194 was classified as a LINER and its combined line ratios are better reproduced by an AGN model (Radovich et al. 2019), while JO201, JO206, JO204, JW100, and JO135 are classified as Seyfert galaxies according to BPT diagrams and are all Seyfert 2. JO204 and JO135 also have extended emission line regions ionized by the AGN. Four of these galaxies display AGN outflows (Radovich et al. 2019).

Having inspected all other GASP cluster member ram-pressure-stripped galaxies, we find other five AGN candidates that are presented in Figure 1. The stripping characteristics of these galaxies were discussed in previous works and only summarized here, but the analysis of their central ionization mechanism is shown here for the first time. JO49 has unwinding tails due to RPS (Bellhouse et al. 2021). It presents a central LINER-like region surrounded by a thin composite-like ring, which in the BPT diagram corresponds to a long finger of points encompassing the composite region extending well beyond the K01 line. We note that JO49 hosts also an X-ray source of luminosity 1.2 × 1041 erg s−1, detected by XMM-Newton (Webb et al. 2020) 17 .

Figure 1.
Standard image High-resolution image
Figure 1.

Figure 1. Left: BPT-NII diagnostic diagram for all spaxels with S/N > 3. The red dotted and continuous lines are defined as in Kewley et al. (2001) and Kauffmann et al. (2003), respectively. The green lines are taken from Sharp et al. (2010). Right: galaxy map color-coded according to the BPT-NII classification; red lines are the stellar emission isocontours corresponding to the galactic disk edges.

Standard image High-resolution image

JO85, another unwinding ram-pressure-stripped galaxy (Bellhouse et al. 2021), has fewer LINER-like points than JO49 embedded in a Composite-like region, but it is highly obscured by dust (AV ∼ 2.7 in the central region as measured by the Balmer decrement) and has a central Chandra point source with a luminosity of 5.0 × 1040 erg s−1 (Evans et al. 2020).

JO147 (first described by Merluzzi et al. 2013; see also Poggianti et al. 2019) is an inclined, highly extincted disk and is stripped in the northwest direction. We find that it has LINER-like opposite cones embedded in wider composite cones. Its luminosity in the X-ray band observed by XMM-Newton is 2.4 × 1041 erg s−1 (Webb et al. 2020).

JO171 is a Hoag-like ring galaxy with long tails stripped in the north direction (Moretti et al. 2018). It has central AGN-powered spaxels (Seyfert 2) in the inner kpc.

Finally, JW39 has long tails originating from unwinding spiral arms (Bellhouse et al. 2021). It has a LINER-like circular central region surrounded by a larger circular area with composite emission.

The latter two galaxies have no available central X-ray counterparts, from neither XMM-Newton nor Chandra.

To summarize, with respect to the sample of AGNs described in P17b and Radovich et al. (2019), we find an additional Seyfert 2 and four LINER-like galaxies, yielding a total sample of 12 AGN hosts in the GASP-RPS sample. Their main properties are summarized in Table 1.

Table 1. AGN Candidates in the GASP Sample

IDR.A.Decl. z Cluster $\mathrm{log}{M}_{* }/{M}_{\odot }$ JstageAGN FlagReferences
JO85351.1306816.868150.0355A258910.713This paper
JO3618.24758315.5914880.0407A16010.834Fritz et al. (2017)
JO194359.25284−34.6805880.042A405911.223 P17b
JO204153.44513−0.9141820.0424A95710.622 P17b
JO20110.376208−9.262750.0446A8510.822 P17b
JO4918.6827090.2861360.0451A16810.723This paper
JO147201.70721−31.3959750.0506A355811.023This paper
JO206318.447542.4762180.0511IIZW10811.022 P17b
JO171302.56125−56.6418230.0521A366710.622This paper
JO135194.26791−30.3750880.0544A353211.022 P17b
JW100354.1044321.1507020.0619A262611.422 P17b
JW39196.0321219.2106910.0663A166811.223This paper

Note. Columns: (1) GASP ID; (2–3) coordinates of the optical center; (4) galaxy redshift; (5) host cluster; (6) galaxy stellar masses (Vulcani et al. 2018); (7) Jstage (B. Poggianti et al. 2021, in preparation); (8) AGN classification; (9) works in which the source is presented. The adopted AGN flag for both GASP-RPS and LIT-RPS galaxies ranges from 0 to 6: 0 means that star formation is the dominant ionization process at the galaxy center according to BPT-N ii classification; 1, 2, 3 if the galaxy hosts a Seyfert 1, Seyfert 2 or LINER-like nucleus, respectively, again according to the BPT diagram; 4 if the AGN has been detected through the X-ray signal, but not in the optical; 5 when the galaxy is classified as a radio galaxy; 6 when the source is classified as an AGN, without any specification on the type.

Download table as:  ASCIITypeset image

The left panel of Figure 2 shows the mass distribution of galaxies hosting an AGN compared to the entire GASP-RPS sample. While RPS galaxies cover a mass range of 8.7 $\leqslant \mathrm{log}({M}_{* }/{M}_{\odot })\leqslant 11.5$, AGN hosts are among the most massive galaxies in the sample, having all $\mathrm{log}({M}_{* }/{M}_{\odot })\geqslant 10.5$.

Figure 2.

Figure 2. Stellar mass distributions for all galaxies (black histogram) and for galaxies hosting an AGN (red histogram). From left to right: the GASP-RPS, LIT-RPS, and MaNGA-Ref samples.

Standard image High-resolution image

We are now in the position of computing the fraction of AGNs (fAGN) over the total (AGN+SF) number of galaxies, considering different subsamples, as summarized in Table 2. The AGN fraction in the total GASP-RPS sample is ${0.24}_{-0.05}^{+0.06}$, with uncertainties computed as binomial errors. Restricting the sample to $\mathrm{log}({M}_{* }/{M}_{\odot })\geqslant $ 10.5, this fraction becomes ${0.71}_{-0.12}^{+0.10}$.

Table 2. AGN Fractions in the GASP-RPS Sample

NAGN/NTOT fAGN Jstage $\mathrm{log}({M}_{* }/{M}_{\odot })$
12/51 ${0.24}_{-0.05}^{+0.06}$ ≥0.5All
12/17 ${0.71}_{-0.12}^{+0.10}$ ≥0.5 ≥10.5
0/16 ${0.0}_{-0.0}^{+0.06}$ =0.5All
1/13 ${0.08}_{-0.05}^{+0.11}$ =1All
10/18 ${0.56}_{-0.12}^{+0.11}$ =2All
1/4 ${0.25}_{-0.15}^{+0.25}$ =3All

Note. Considering galaxies of different mass, and characterized by different Jstages. Errors are binomial.

Download table as:  ASCIITypeset image

Considering the various stages of stripping (Figure 3), the frequency of AGNs increases with the strength of RPS signatures; no galaxies with AGN activity have Jstage = 0.5, while the AGN incidence increases among moderate stripping galaxies (Jstage = 1, 8%) and is particularly high among Jstage = 2 galaxies, where it reaches 56%. Still one out of four galaxies in a late stage of RPS (truncated disks, Jstage = 3) has an AGN. Interestingly, using the same methods of the current analysis, only two AGNs are found in the GASP non-RPS sample of star-forming galaxies, which consists of 49 galaxies (Vulcani et al. 2021; B. Poggianti et al. 2021, in preparation).

Figure 3.

Figure 3. Stacked histograms for galaxies of different Jstages and divided among centrally star-forming galaxies (light blue histogram) and AGNs (dark blue histogram) according to the BPT-N ii classification for GASP-RPS (left) and LIT-RPS (right). Percentages are AGN fractions in the corresponding bin of Jstage and Jstagegen.

Standard image High-resolution image

In the GASP-RPS sample, out of the 17 galaxies more massive than $\mathrm{log}({M}_{* }/{M}_{\odot })\geqslant 10.5$, 10 (∼58%) have Jstage = 2 and, vice-versa, ∼55% (10/18) of the Jstage = 2 galaxies are more massive than $\mathrm{log}({M}_{* }/{M}_{\odot })\geqslant 10.5$. All of them host an AGN. It is significant that none of the massive galaxies have a Jstage = 0.5. This result suggests a tight correlation between stellar mass and Jstage. The correlation is probably linked with the higher capability of massive galaxies to retain gas. While low-mass galaxies are already completely stripped when they approach the densest regions in clusters, high-mass galaxies more easily hold onto their gas (Jaffé et al. 2018; Luber et al., submitted to ApJ) and experience RPS in these dense regions, where the gas removal is the most intense. Since AGNs are preferentially located in the most stripped and massive galaxies, we cannot state which of these two parameters is more connected to the presence of an AGN.

3.2. LIT-RPS

The catalog of the 82 literature ram-pressure-stripped galaxies with SF/AGN information is presented in Table 3. Stellar masses and AGN classification, along with the source for those values are given in Table 4.

Table 3. Portion of the Catalog of Ram-pressure-stripped Galaxies

NameR.A.Decl. z ClusterJstageJstagegen Alternative Names
MIP001417–3023033.5693−30.38430.2955A2744...2.0F1228
HLS001427–3023443.61065−30.395810.3033A2744...2.0F0083
NGC 156665.00175−54.937810.005Dorado0.00.0...
ID345149.91912.52810.727CGr322.02.0...
LEDA36382175.7352319.966210.02427A13672.02.0CGCG97073

Note. Columns: (1) galaxy name; (2) and (3) equatorial coordinates of the galaxy center from SIMBAD; (4) galaxy redshift; (5) host cluster; (6) and (7) Jstage and general Jstage, defined in Section 3.2; (8) alternative names. Stellar masses flagged with the asterisk (*) are computed by means of photometric data as described in the Appendix. This table is published in its entirety in a machine-readable format. A portion is shown here for guidance regarding its form and content.

Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table is available.

Download table as:  DataTypeset image

Table 4. Stellar Masses and AGN Classifications

Name $\mathrm{log}({M}_{* }/{M}_{\odot })$ AGN References
MIP001417–3023039.6 Rawle et al. (2014)0 Owers et al. (2012)Owers et al. (2012); Rawle et al. (2014)
HLS001427–30234410.9 Rawle et al. (2014)1 Owers et al. (2012)Owers et al. (2012); Rawle et al. (2014)
NGC 156610.8 Elagali et al. (2019)1 Véron-Cetty & Véron (2006)Elagali et al. (2019)
ID34510.3 Boselli et al. (2019)1 Boselli et al. (2019)Boselli et al. (2019)
LEDA363829.5 Mendel et al. (2014)0 SDSS (2004)Gavazzi et al. (1995, 2001)
   Sivanandam et al. (2014); Boselli et al. (2018); Yagi et al. (2017)

Notes. Columns: (1) galaxy name; (2) logarithm of the stellar masses, adopting Chabrier (2003) IMF (references from which the values have been taken are reported below); (3) AGN classification (see Table 1) and related reference; (4) references which present a characterization of the galaxy as an RPS candidate. The four galaxies with classifications equal to 4 and 5 (e.g., where the AGN is spotted observing them in X and radio) resulted to be star-forming in the optics. In the text we analyze the consequence to change their AGN flag to 0. This table is published in its entirety in a machine-readable format. A portion is shown here for guidance regarding its form and content.

References. Abramson & Kenney (2014); Abramson et al. (2016); Best et al. (2012); Birchall et al. (2020); Boselli et al. (2005, 2006, 2015, 2016, 2018, 2019); Chemin et al. (2005); Chen et al. (2020); Chung et al. (2007, 2009); Chyży et al. (2007); Consolandi et al. (2017); Cortese et al. (2007); Cramer et al. (2020, 2019); Damas-Segovia et al. (2016); Davies et al. (2020); Decarli et al. (2007); Ebeling et al. (2019, 2014); Elagali et al. (2019); Fossati et al. (2012, 2016); Fruscione et al. (1990); Fumagalli et al. (2011, 2014); Gavazzi (1989); Gavazzi et al. (1984, 1995, 2001, 2017, 2018); Gu et al. (2013); Gullieuszik et al. (2020); Hester et al. (2010); Ho et al. (1997); Jáchym et al. (2014, 2017, 2019); Kalita et al. (2019); Kantharia et al. (2008); Kenney et al. (1995, 2004, 2008, 2014, 2015); Lee et al. (2017, 2018); Mahajan et al. (2010); McPartland et al. (2016); Mendel et al. (2014); Merluzzi et al. (2010, 2013, 2016); Minchin et al. (2019); Nucita et al. (2017); Oosterloo et al. (2005); Owen et al. (2006); Owers et al. (2012); Rahman et al. (2011); Rakshit et al. (2017); Rawle et al. (2014); Roberts et al. (2020); Ruszkowski et al. (2014); Salim et al. (2018, 2016); Scott et al. (2010); SDSS (2004); Sivanandam et al. (2010, 2014); Smith et al. (2010); Sobral et al. (2015); Stein et al. (2017); Stroe et al. (2020); Sun et al. (2005, 2006, 2007, 2010); Tschöke et al. (2001); Véron-Cetty & Véron (2003, 2006); Vollmer & Vollmer (2003); Vollmer et al. (2004, 2009); Wang et al. (2020); Weżgowiec et al. (2012); Yagi et al. (2010, 2013, 2017); Yoshida et al. (2002, 2004, 2012); Zhang et al. (2013).

Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table is available.

Download table as:  DataTypeset image

One of these galaxies has broad optical lines typical of Seyfert 1. For all the other galaxies the AGN classification is based on the BPT-N ii diagnostic. For ∼68% (55/81) of them the classification in published results is based on spectroscopic observations published in dedicated papers, either from the integral-field unit (Merluzzi et al. 2013; Fossati et al. 2016; Merluzzi et al. 2016; Consolandi et al. 2017; Boselli et al. 2019; Stroe et al. 2020), long-slit or fiber spectra (Véron-Cetty & Véron 2003; Owen et al. 2006; Cortese et al. 2007; Mahajan et al. 2010 18 ; Owers et al. 2012; Ebeling et al. 2019). For the other 32% (26/81) of the galaxies, we instead use the online AGN classification based on the analysis of emission line ratios extracted from integrated spectra of the central circular aperture (r ∼ 3'') observed with the SDSS fiber (Data Release 8, from now on DR8; Aihara et al. 2011) as analyzed by Brinchmann et al. (2004), Kauffmann et al. 2003, and Tremonti et al. (2004) in the Value Added Catalog MPA/JHU. For 24 galaxies, we had both the DR8 automatic classification and information about the central source from individual publications in the literature. In these cases we favored the latter.

We note that four of the 81 galaxies also have information coming from either X-ray or radio data (Winkler & Winkler 1992; Best et al. 2012; Owers et al. 2012; Ebeling et al. 2019; Kalita et al. 2019; Caglar et al. 2020). While their position on a BPT-N ii diagram suggests they are star-forming, the additional data instead classify them as AGNs. In what follows we will therefore discuss how results change if we include or exclude these four objects.

Overall, 24/82 galaxies host an AGN (∼30%). If we disregard the AGN classification based on X-ray or radio data and consistently consider only the BPT-N ii classification, the fraction above becomes 20/82 (∼24%).

The central panel in Figure 2 shows the mass distribution of the galaxies with and without an AGN. The entire sample spans a mass range 8.1 $\lt \mathrm{log}({M}_{* }/{M}_{\odot })\lt $ 11.4. Similar to what was found for GASP, most of the AGNs are massive galaxies, even though in this sample there are also a few less massive AGN hosts. Above the GASP AGN mass limit ($\mathrm{log}({M}_{* }/{M}_{\odot })\gt 10.5$) the AGN fraction becomes ${0.80}_{-0.12}^{+0.08}$.

Table 5 reports the AGN fraction for the different subsamples considered, including that for galaxies of different Jstagegen. The trend of the AGN fraction with Jstagegen is weaker than in GASP-RPS (see also Figure 3), with the percentages ranging between 17% and 27% but being consistent within the large errors in all Jstagesgen. 19

Table 5. AGN Fractions for the LIT-RPS Sample

NAGN/NTOT fAGN Jstagegen $\mathrm{log}({M}_{* }/{M}_{\odot })$
20/82 (24/82) ${0.24}_{-0.04}^{+0.05}$ (${0.29}_{-0.05}^{+0.05}$)AllAll
12/15 (12/15) ${0.80}_{-0.12}^{+0.08}$ (${0.80}_{-0.12}^{+0.08}$)All≥10.5
4/15 (5/15) ${0.27}_{-0.10}^{+0.13}$ (${0.33}_{-0.11}^{+0.13}$)=0All
2/10 (2/10) ${0.20}_{-0.10}^{+0.15}$ (${0.20}_{-0.10}^{+0.15}$)=0.5All
3/18 (4/18) ${0.17}_{-0.07}^{+0.11}$ (${0.22}_{-0.08}^{+0.11}$)=1All
8/30 (10/30) ${0.27}_{-0.07}^{+0.09}$ (${0.33}_{-0.08}^{+0.09}$)=2All

Note. Considering galaxies of different mass ranges and characterized by different Jstagegen. Errors on fractions are binomial. Values outside/in brackets are the fractions computed ignoring/considering the four galaxies classified as AGNs based on radio and X data.

Download table as:  ASCIITypeset image

We remind the reader that while the AGN classification and mass estimates among GASP galaxies are homogeneous, for the LIT-RPS sample we based the former on a number of different data and indicators. In addition, stellar masses have been computed following many different approaches and so, even though homogenized to the same IMF, there could be some systematics among the different galaxies. Finally, we recall that the Jstagegen flag is based on a very heterogeneous set of images in terms of wavelengths, depth, quality, and therefore results must be taken with caution.

4. Results II: Is the AGN Fraction among Ram-pressure-stripped Galaxies Higher than in Nonstripped Galaxies?

In the previous section we have quantified the incidence of AGNs in ram-pressure-stripped galaxies. We have seen that they represent 24% of the overall, both for GASP-RPS and LIT-RPS samples. In the following we will always exclude masses < 109 M, in all samples. Table 6 presents the AGN fractions in GASP-RPS and LIT-RPS separately for stellar masses ≥109 M and ≥1010 M. 20 The fractions in the LIT-RPS sample are always higher than in GASP-RPS (0.29 versus 0.24 and 0.55 versus 0.46, respectively for the two mass bins), but are compatible within the binomial errors. Including also the four X-ray/radio AGNs in the literature, fractions are slightly higher. We note that also the stellar mass distributions of the two samples are similar (Figure 4), and indeed a Kolmogorov–Smirnov (KS) test cannot exclude that they are drawn from the same parent distribution.

Figure 4.

Figure 4. Normalized stellar mass distributions of the GASP (black histogram), MaNGA (purple-dotted histogram), and literature sample (light red histogram). For the Montecarlo, we have selected galaxies above the vertical-dotted line, i.e., with masses $\mathrm{log}({M}_{* }/{M}_{\odot })\geqslant 9$.

Standard image High-resolution image

Table 6. AGN Fractions and Binomial Errorbars for the GASP-RPS, LIT-RPS, and ALL-RPS Samples in Two Different Mass Bins

 GASP-RPS LIT-RPS ALL-RPS
$\mathrm{log}({M}_{* }/{M}_{\odot })$ NAGN/NTOT fAGN   NAGN/NTOT fAGN   NAGN/NTOT fAGN
≥9.012/50 ${0.24}_{-0.06}^{+0.07}$  19/65(23/65) ${0.29}_{-0.05}^{+0.06}({0.35}_{-0.06}^{+0.06})$  31/115(35/115) ${0.27}_{-0.04}^{+0.04}$ (${0.30}_{-0.04}^{+0.04}$)
≥10.012/25 ${0.46}_{-0.09}^{+0.10}$  17/31(19/31) ${0.55}_{-0.09}^{+0.09}$ (${0.61}_{-0.09}^{+0.08}$) 29/57(31/57) ${0.51}_{-0.07}^{+0.07}({0.54}_{-0.07}^{+0.07})$

Note. For the LIT-RPS and ALL-RPS samples, values in parentheses are obtained considering also galaxies identified as AGNs on the basis of X-ray or radio data.

Download table as:  ASCIITypeset image

It is therefore appropriate to join the two RPS samples to obtain the largest possible statistics 21 and derive the total AGN fractions: ${0.27}_{-0.04}^{+0.04}$ at masses ≥109 M and ${0.51}_{-0.07}^{+0.07}$ for masses ≥1010 M. These fractions are high but less extreme than the fraction that would have been inferred from the P17b results, where 6/7 galaxies were AGNs with a corresponding fraction of ${0.86}_{-0.09}^{+0.18}$. This is due to the fact that the 2017 sample was composed of all massive Jstage = 2 galaxies, and as we have seen in the previous sections these are the most favorable conditions for AGN activity in ram-pressure-stripped galaxies.

We now aim at establishing whether the AGN frequency is connected to RPS and therefore we compare the measured fractions to those obtained exploiting the MaNGA-Ref sample, used as representative of non-ram-pressure-stripped field galaxies.

As for the other samples, also in MaNGA-Ref, AGNs are located preferentially among the most massive galaxies (right panel of Figure 2). The AGN fraction is ${0.15}_{-0.01}^{+0.01}$ 22 above $\mathrm{log}({M}_{* }/{M}_{\odot })\geqslant 9.0$ and ${0.28}_{-0.02}^{+0.02}$ for $\mathrm{log}({M}_{* }/{M}_{\odot })\geqslant 10.0$.

A KS test excludes that the MaNGA-Ref mass distribution is drawn from the same parent distribution of the GASP+LIT sample (Figure 4). Since the probability to find an AGN increases with galaxy mass, to properly compare the fractions obtained from MaNGA-Ref and ALL-RPS we need to control for the different mass distributions. We perform a bootstrap random extraction of the MaNGA sample to create 10,000 subsamples with the same mass distribution of the ALL-RPS sample matching the number of ALL-RPS galaxies in bins of 0.3 dex in stellar mass. For each of the extracted samples we compute the AGN fraction fAGN. We repeat the random extraction considering separately two stellar mass ranges, M* ≥ 109 M and M* ≥ 1010 M. Violin plots with the fAGN distributions for the two mass ranges, their medians, and the 25th and 75th percentiles are shown in Figure 5. We find that the median fAGN of the 10,000 realizations of mass-matched MaNGA galaxies are fAGN = 0.18 for M ≥ 109 M and 0.35 for M ≥ 1010 M. These values are lower than the corresponding values in the ALL-RPS sample, which are 0.27 and 0.51, respectively. In order to assess the significance of the difference between the RPS and non-RPS samples, we compute the pivotal confidence intervals of the bootstrap distribution and find that the mass-matched MaNGA fractions are lower than the ALL-RPS fractions at the 99.99% confidence level for galaxies with M ≥ 109 M (>99.99% if we include the four radio/X-ray AGNs), and at the 99.96% level for M ≥ 1010 M (>99.99% if the four radio/X-ray AGNs are included) (see Figure 5).

Figure 5.

Figure 5. Comparison of the AGN fraction in the different samples. Red and orange lines refer to the ALL-RPS sample: the AGN fraction for galaxies with M* ≥ 109 M is shown by the thick orange line, that for galaxies with M* ≥ 1010 M by the thick red line. The matched shaded areas indicate by how much fractions change if we consider also the AGN classified on the basis of X-ray or radio data (see text for details). Blue and light blue violin plots refer to the MaNGA-Ref sample, for the two mass bins as indicated in the labels. They show the probability density of the bootstrap random extractions mass-matched to the ALL-RPS sample, at different AGN fraction values, smoothed by a kernel density estimator. Gray horizontal dashed and dotted lines represent median values and 25% and 75% percentiles of the AGN fraction, respectively. Values of the pivotal confidence intervals of the bootstrap distribution are also reported: the mass-matched MaNGA fractions are lower than the ALL-RPS fractions at the 99.99% confidence level for galaxies with M ≥ 109 M and at the 99.96% level for M ≥ 1010 M.

Standard image High-resolution image

Since the three samples considered span slightly different redshift ranges, we performed the bootstrap random extractions also limiting all samples to z ≤ 0.075 (the GASP redshift limit). Results remained unchanged as fractions are affected only at the 1% level at most. Finally, we also tried comparing mass-matched MaNGA samples separately with GASP-RPS and LIT-RPS. Though clearly the statistics decrease, we still find high probabilities that the mass-matched MaNGA sample has lower AGN fractions than the RPS samples (81.4% and 95.4% for GASP-RPS in the two mass ranges, and >99.99% for both LIT-RPS samples).

From our analysis the incidence of AGN activity among ram-pressure-stripped galaxies is significantly higher than that in the MaNGA field control sample. A ram-pressure-stripped galaxy has a 1.5 times higher probability to host an AGN than a similar non-ram-pressure-stripped galaxy. This effect is not driven by different stellar mass distributions and points to a connection between RPS and AGN activity.

A larger (of the order of hundreds), homogeneous sample of ram-pressure-stripped galaxies with integral-field spectroscopy would be needed to place these results on more solid ground. Since this is currently unavailable, the analysis presented here collects the best available data sets for addressing the question of the AGN–RPS connection. There are however several caveats worth stressing.

First of all, the AGN fraction depends strongly on the criteria adopted when using the BPT diagram. In this paper, we are including both LINER-like and Seyfert AGN, in order to capture also low luminosity AGN. This is done in all samples considered in a similar manner, so it should not affect the relative incidence and the main conclusions of this work, but the pure AGN fractions will strongly depend on the initial choice.

Second, although great care has been taken to ensure the most homogeneous analysis possible, the data sets are clearly nonhomogeneous. Even GASP and MaNGA, which are both based on integral-field data for every galaxy, have been observed with different instruments, thus have different resolutions, spaxel size, etc., and span a slightly different redshift range (see above for invariance of AGN fractions with the redshift interval adopted). The literature sample, obviously, is in itself very heterogeneous, with the spectroscopic information coming from many different sources. The results shown in this paper should therefore be taken with caution and revisited once large homogeneous samples will become available.

Third, in principle, it is possible that the high AGN fraction we observe in ram-pressure-stripped galaxies is not a consequence of RPS itself. If the AGN incidence in cluster star-forming galaxies was higher in general than in similar galaxies in the field, the differences with respect to MaNGA would go in the same direction as what we observe. However, as mentioned also above, the AGN fraction in the GASP non-RPS sample is small (2/49) (Vulcani et al. 2021). This sample is composed both of cluster and field undisturbed galaxies. If we consider only the GASP cluster control sample (star-forming and late-types), there is no AGN (B. Poggianti et al. 2021, in preparation). So, this caveat is unlikely to be responsible for our results.

Finally, we note that we are not studying the global AGN fraction in clusters, but the occurrence of AGN activity in a very specific class of cluster galaxies: those with clear signs of RPS, which are all star-forming and late-type galaxies and thus represent a small fraction of the total cluster galaxy population that is dominated by early-type galaxies. Therefore, our results cannot be used to infer the total AGN fraction in clusters and not necessarily show similar trends.

5. Summary

In this paper we have investigated the occurrence of AGN activity in ram-pressure-stripped galaxies in local clusters, comparing it with the AGN frequency in a control sample of field galaxies. In all cases, we rely on BPT diagnostic diagrams based on the [N ii] line. All the galaxies analyzed in this paper are star-forming and morphologically late-type galaxies.

First, we assembled two samples of ram-pressure-stripped galaxies. We have used the MUSE data of 51 galaxies observed in the context of the GASP survey (GASP-RPS) finding a Seyfert 2 and 4 LINER-like AGN hosts previously unknown, in addition to the seven galaxies already discussed in P17b and Fritz et al. (2017). We have then conducted a search in the literature assembling a sample of 82 ram-pressure-stripped galaxies for which it was possible to retrieve information on their nuclear activity (either from IFU or slit/fiber) (LIT-RPS).

We find similar fractions of AGNs in GASP and in literature ram-pressure-stripped galaxies, with the AGN incidence being slightly higher in the literature than in GASP, but consistent within the uncertainties. Overall, the AGN fraction in the total GASP-RPS+LIT-RPS sample is ${0.27}_{-0.04}^{+0.04}$ at masses M* ≥109 M and ${0.51}_{-0.07}^{+0.07}$ at M* ≥ 1010 M. Thus, more than half of the ≥1010 M ram-pressure-stripped galaxies show AGN activity.

We then compare these findings with those for a sample of galaxies drawn from the MaNGA survey and inhabiting dark matter haloes with masses ≤1013 M. With this halo mass cut we ensure that rich groups and clusters are excluded, hence these galaxies are not undergoing significant RPS and this can serve as a control field sample. We perform a bootstrap random extraction from the MaNGA sample to create 10,000 realizations with the same stellar mass distribution of the RPS sample.

Our two main results can be summarized as follows:

  • 1.  
    The great majority of galaxies hosting an AGN, in all three samples considered, are high-mass galaxies. There are just very few galaxies with an AGN at a mass below 1010 M (none below 1010.5 M in GASP). As a consequence, the AGN fractions are higher above these limits, and very low below. Another factor that could be playing a role is the ram-pressure strength or phase (Jstage): the highest AGN fractions are observed among the most strongly ram-pressure-stripped galaxies with the longest tails. However, with the current samples it is hard to disentangle between mass and Jstage effects.
  • 2.  
    Even after matching the galaxy mass distributions, the AGN incidence in the field MaNGA sample is lower than in the RPS sample at the ≥99.96% confidence level. Overall, a ram-pressure-stripped galaxy has a 1.5 times higher probability to host an AGN than a similar non-ram-pressure-stripped galaxy. This supports the hypothesis that ram pressure can trigger the AGN activity.

We thank the Referee for the useful suggestions that improved the presentation of the work. We warmly thank Jong-Ho Shinn from the Korea Astronomy and Space Science Institute for useful discussion regarding the statistical analysis. G.P. thanks A. Werle for the helpful discussion. Based on observations collected at the European Organization for Astronomical Research in the Southern Hemisphere under ESO program 196.B-0578. This project has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation program (grant agreement No. 833824). We acknowledge the financial contribution from the grant PRIN MIUR 2017 n.20173ML3WW_001 (PI: Cimatti), from the INAF main-stream funding program (PI: Vulcani), and from the agreement ASI-INAF n.2017-14- H.0 (PI: A. Moretti). Y.J. acknowledges financial support from CONICYT PAI (Concurso Nacional de Inserción en la Academia 2017) grant No. 79170132 and FONDECYT Iniciación 2018 grant No. 11180558. J.F. acknowledges financial support from the UNAM- DGAPA-PAPIIT IN111620 grant, México.

Appendix: Mass Estimates

To compute stellar masses for those galaxies for which they are missing in the literature. we use the Bell et al. (2001) relation between the mass-to-light ratio of a galaxy and its color:

Equation (A1)

where Lλ is the luminosity in a band, indicated with λ, COL is a photometric color and aλ and bλ are coefficients depending on both λ and COL. For our calculations, we used the Bell et al. (2001) tables for a solar metallicity Z = 0.02 and a Bruzual & Charlot (2003) SSP model, converting from a Salpeter (1955) to a Chabrier (2003) IMF subtracting a factor ‒0.24. For one galaxy, 235144-260358 (Cortese et al. 2007), in order to use the formula (A1) we first converted HST magnitudes to a UBV photometric system with the use of calibration equations for the Advanced Camera for Surveys (ACS) presented in Sirianni et al. (2005). We assume a typical 0.3 dex uncertainty on the computed stellar masses, which we take as bin size of the stellar mass distribution.

Footnotes

  • 6  

    See also Boselli et al. (2021) which appeared unrefereed on astroph on September 28 2021.

  • 7  

    The choice of the K01 demarcation line to identify AGNs is a conservative choice that minimizes the AGN spaxels. Note that recent works (see Law et al. 2021 and references therein) find a demarcation line in the N ii-BPT diagram closer to the K03 separation, but we choose the most conservative one to minimize the contamination from star-forming regions.

  • 8  

    We adopted this number upon visual inspection of the maps. This choice allows us to have enough spaxels to identify possible AGNs with high confidence, still focusing on the central part of the galaxy.

  • 9  

    The candidate merging systems removed from the sample are F0237 (Owers et al. 2012), NGC 4294, NGC 4299, and NGC 4302 (Pappalardo et al. 2012; Vollmer et al. 2013).

  • 10  

    A list of all ram-pressure-stripped galaxies known in the literature, regardless of AGN information, will be published in J. Crossett et al. 2021, in preparation.

  • 11  
  • 12  

    We verified that results are insensitive to the exact choice of this threshold, exploring log halo masses up to 13.6. We decided to use a conservative cut (13.0) to avoid the possibility that ram-pressure-stripped galaxies in groups contaminate the sample.

  • 13  
  • 14  

    For consistency with the other samples, we have applied the morphological cut, but all the results remain unchanged if no morphological criterion is applied.

  • 15  

    This is why in GASP the data cubes have been filtered with a 5 × 5 kernel (Poggianti et al. 2017a).

  • 16  
  • 17  

    The 4XMM-DR10 catalog contains source detections covering an energy interval from 0.2 keV to 12 keV. On the other hand, the Chandra energy range goes from 0.5 to 7 keV.

  • 18  

    We note that for two galaxies, GMP3618 and D100, Mahajan et al. (2010) give different results with respect to the classification reported in the DR7 (Abazajian et al. 2009) and DR8 (Aihara et al. 2011) analysis even though they used DR7 data to build up BPT-N ii.

  • 19  

    The subsample of LIT-RPS galaxies with Jstage is too small to study trends with the length of the Hα tails.

  • 20  

    From now on we exclude from this analysis ID345 at z = 0.73, considered a redshift outlier.

  • 21  

    Note that the galaxy JO147 appears in both samples, from now on we will just consider it once.

  • 22  

    Previous MaNGA works (e.g., Sánchez et al. 2018) have found a significantly lower AGN incidence. However previous analysis has not applied any cut in sSFR as we do, and have adopted much more stringent definitions (emission line ratios above the Kewley demarcation lines considering all the BPT diagrams simultaneously and Hα equivalent width >1.5 Å in the central regions) for AGNs, therefore results are not directly comparable.

Please wait… references are loading.
10.3847/1538-4357/ac4225