Free Access
Issue
A&A
Volume 638, June 2020
Article Number A126
Number of page(s) 14
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/202037803
Published online 23 June 2020

© ESO 2020

1. Introduction

The physics driving the empirically established dependence of most galaxy properties on their environment is still poorly constrained (Boselli & Gavazzi 2006, 2014; Park & Hwang 2009). Ram-pressure stripping (Gunn & Gott 1972; Quilis et al. 2000; Vollmer et al. 2001; Bekki & Couch 2003) is often invoked to explain the observations of gas-deficient galaxies or galaxies with irregular appearance in clusters (Kenney et al. 2004; Chung et al. 2007; Poggianti et al. 2017; Jachym et al. 2019; Ebeling & Kalita 2019). It has also been successfully shown to explain the distribution, in projected phase-space, of passive galaxies in clusters (Jaffé et al. 2016, 2018; Yoon et al. 2017). However, because its strength depends on galaxies moving at high relative velocities through a relatively dense gas, ram-pressure stripping is only expected to be efficient in massive clusters of galaxies. In reality the distribution of galaxy properties throughout the large-scale structure of the universe requires a complex combination of multiple physical mechanisms. Some of the galaxies are quenched well before they find their way to the interiors of clusters (Bianconi et al. 2018). At low redshifts, approximately 40% of the galaxies residing outside clusters show no ongoing star formation (Balogh et al. 1999; Verdugo et al. 2008), a fact which cannot be explained with ram-pressure stripping. A collective term, pre-processing, is often used to explain the presence of this passive population and includes a number of physical mechanisms acting on galaxies residing in smaller haloes or filaments. Most likely starvation is the mechanism responsible for this (Larson et al. 1980; Balogh et al. 2000; Bekki et al. 2002; Goto et al. 2003) as it is capable of affecting a large fraction of all the galaxies, the ones that are satellites in a larger halo but not necessarily in a cluster of galaxies. Starvation is the absence of accretion of fresh gas onto galaxies that can be used for star formation. It has been shown that the star formation histories of the passive galaxies in the local universe are consistent with starvation being the main mechanism responsible for quenching their star formation (Peng et al. 2015). Einasto et al. (2014) showed that the fraction of passive galaxies in groups increases with group richness, and that the environment in which the groups reside modulates this trend. There are also other mechanisms that can alter the star formation properties of galaxies in lower density environments (Cowie & Songaila 1977; Fujita 2004). Thus, galaxies that have been pre-processed undergo little change when they reach the interiors of massive clusters. Other galaxies, likely individual, only stop forming stars when they cross the virial radius of a cluster and are observed as jellyfish galaxies (Poggianti et al. 2017). In addition, high-density environments, be it groups, clusters, or non-virialised structures like filaments, are associated not only with quenching, but sometimes also with an increase in star formation (Fujita & Nagashima 1999; Hwang & Lee 2009; Stroe et al. 2015; Deshev et al. 2017; Vulcani et al. 2019).

Abell 963 (hereafter A963) is one of the best studied clusters outside the local universe. It has been targeted with XMM-Newton, Galex, Spitzer, and the LMT (Cybulski et al. 2016). The cluster has been observed in the optical with the Isaac Newton Telescope, and virtually all of its members brighter than R ≤ 19 have been targeted for optical spectroscopy Jaffé et al. (2016, hereafter J16). A963 is one of only two clusters outside the local universe for which we have 21 cm data, a result of a blind interferometric survey with the WSRT called BUDHIES (Verheijen et al. 2007; Deshev et al. 2009; Gogate et al. 2020). The cluster and its constituent galaxies have been extensively studied (Jaffé et al. 2013, J16). The projected phase-space analysis of A963 presented in J16 explains to a remarkable extent the distribution of the galaxies with different properties residing in this cluster, and confirms the important role of ram-pressure stripping as the reason for the high passive fraction of cluster galaxies. However, the projected phase-space analysis, particularly when applied to individual clusters, ignores part of the information contained in the data, like possible deviations from spherical symmetry of the cluster and anisotropic accretion of galaxies. Although some clusters might indeed be close to spherical, this is certainly not the case for clusters experiencing recent mergers or significant accretion. The accretion of galaxies onto clusters is likely not isotropic both in terms of the number of galaxies falling from different directions and the level of pre-processing of these galaxies (Martínez et al. 2016; Kuutma et al. 2017; Crone Odekon et al. 2018). Observations and simulations show that 35%–40% of galaxies, by luminosity, are found in filaments (Tempel et al. 2014). In addition, the galaxies accreted along filaments are subjected to pre-processing (Kuutma et al. 2017; Vulcani et al. 2019), so they have different properties from those accreted individually.

In this article we concentrate on the distribution of galaxies with different properties on the sky, complementing the analysis presented in J16. While this falls short of the desired full, observed, phase-space analysis, we demonstrate that it is very informative about the fate of galaxies infalling onto large clusters, and the extent to which some of them are pre-processed. We concentrate on the galaxy members of A963_1, which is the main cluster (Jaffé et al. 2013), and we refer to it as A963 for the remainder of this article.

Throughout, we use a standard ΛCDM cosmology with H0 = 70 km s−1 Mpc−1, Ωm = 0.3, and ΩΛ = 0.7. A963 is observed at z = 0.205 (∼1.0 Gpc), thus 1″corresponds to 3.36 kpc.

2. Data

We combine a multi-fibre optical spectroscopy of galaxies, with rest-frame 21 cm imaging (WSRT), and XMM-Newton X-ray imaging. All the data used for the analysis presented in this article have been used in previous publications or have been publicly released. Here we briefly describe the data while presenting the measurements of the individual galaxy properties. The lead data set for this study is a combination of two separate optical spectroscopy runs of galaxies in the field of A963, both executed with Hectospec on the MMT1. Together they contain optical spectra of 1712 galaxies within the square degree surrounding the cluster centre. Approximately 3/4 of these data come from the ACReS survey (Haines et al. 2013), which is the spectroscopic part of the LoCuSS survey (Smith et al. 2010). The rest are from a survey published by Hwang et al. (2014). There are many galaxies targeted by both surveys. In this study we select all the targets from ACReS and add the unique data from Hwang et al. (2014). These data formed the majority of the data set (> 93%) analysed by J16.

2.1. Cluster membership

Since most of our data come from the ACReS survey (∼75%), we adopted their selection of cluster members (Haines et al. 2015) based on the distribution of galaxies in projected phase space. Because of the shape of the gravitational potential of the cluster, the member galaxies form a well-defined trumpet shape. While adding the additional data we overplotted them onto the ACReS cluster members, and selected the galaxies overlapping with the already defined cluster region. This process is shown in Fig. 1, where the 481 cluster members are shown with blue squares, while those from the ACReS survey are outlined with black squares. The distribution of the remaining galaxies with available spectroscopy (orange dots) clearly shows the presence of other structures around the main cluster, as discussed in Jaffé et al. (2011). The central redshift of the cluster, (z = 0.205) measured by Deshev et al. (2017) is shown with horizontal grey line on the figure. We use the position of the brightest cluster galaxy as a cluster centre (RA = 10 h 17 m 03 . s 60 $ 10^{\mathrm{h}} 17^{\mathrm{m}} 03{{\overset{\text{ s}}{.}}}60 $, Dec = + 39 ° 02 49 . 920 $ +39^{\circ}02^{\prime}49{{\overset{\prime\prime}{.}}}920 $), measured from the R-band imaging of the cluster with INT, La Palma, collected as part of the BUDHIES survey (Gogate et al. 2020). The X-ray emission does not show any significant displacement from this point on all scales. The final number of cluster members used in this analysis is 386, after applying a stellar mass cut and removing galaxies with active galactic nuclei (Sect. 2.3). The brightest cluster galaxy is not part of this spectroscopic survey.

thumbnail Fig. 1.

Selection of the cluster members (blue squares) of A963 in projected phase space, redshift vs. cluster-centric distance. The cluster members from the ACReS survey have black squares around their symbols. The orange dots show all the galaxies with spectroscopy in the field. The redshift of the cluster (z = 0.205) is shown with the grey horizontal line.

Figure 2 shows the K-band luminosity density based on all the cluster members. The colour map, together with the logarithmically spaced yellow contours, show the main overdensity of A963 together with a number of other structures clustered predominantly along the north–south axis of the cluster. The red contour is drawn at 1.8 K Mpc−2. This level was chosen because it encompasses all the cluster members and is used, together with a circle with a radius of 0.52°, to define the radial extent of the maps of the cluster presented in the following sections.

thumbnail Fig. 2.

K-band luminosity density of the A963 cluster members. Solid squares show cluster members with log(M) ≥ 10.0 M. Small points show cluster members with lower stellar mass. The yellow contour levels are indicated by the tick marks on the colour bar. The green circle has a radius equal to R200 of the cluster. The red contour, at 1.8 K Mpc−2, is used to define the radial extent of the rest of the maps in the article.

2.2. Voronoi binning

Binning of data with variable bin size is a common way of achieving constant signal-to-noise ratio across the bins. One way of calculating adaptive bins on a 2D plane is Voronoi binning, often used with integral-field unit (IFU) observations of galaxies, where the exponential drop in surface brightness requires binning away from the galactic centre. Mapping the properties of galaxies residing in clusters presents similar problems due to the drop in number density away from the cluster centre. We use Voronoi binning, as implemented by Cappellari & Copin (2003) in the python procedure Vorbin. The program is designed to work with IFU data, and as such it requires the data to be on a regular grid with well-defined pixels. We start by putting a square sky area centered on A963 on a regular grid with 90 × 90 cells and with 1.1° on the side. This number of cells is large enough so that the most populated one contains only five galaxies. At the same time it is small enough to guarantee the fast working of Vorbin, as the computational time required scales linearly with the number of pixels to bin. In addition, we changed the Vorbin subroutine calculating the signal-to-noise ratio in the bins to a simple sum in order to achieve a semi-constant number of galaxies per Voronoi bin. We required the Voronoi bins to contain about ten galaxies each in order to get meaningful statistics over each bin. In reality, because the geometry of the final bins is also taken into account (Cappellari & Copin 2003) the final number of galaxies in each bin varies between 9 and 15 with a median of 12. The 386 galaxies are assigned to 33 Voronoi bins. A map of the Voronoi bins with their number, referred to in the rest of the article, and showing the number of galaxies in each bin is shown in Fig. 3. The green circle at the centre of the figure has a radius equal to R200 of the cluster as calculated by J16 (R200 = 1.85 Mpc), based largely on the same data set. The map covers the area where the K-band luminosity density is higher than 1.8 K Mpc−2.

thumbnail Fig. 3.

Map of the individual Voronoi bins. Each bin is coloured according to the total number of galaxies within it. All the bins are numbered, starting with zero. The green circle has a radius equal to R200 of the cluster.

2.3. Stellar mass, spectroscopic completeness, and galaxy classification

The stellar mass and star formation rate of all the cluster members was calculated with Multi-wavelength Analysis of Galaxy Physical properties (MAGPHYS, da Cunha et al. 2008). The fit was based on 18 bands covering wavelengths from far-infrared to far-ultraviolet, assuming Chabrier IMF (see Bianconi et al. 2020, for details). The distribution of stellar mass of A963 cluster members is shown in panel a of Fig. 4. Our completeness limit of log(M) ≥ 10.0 M is shown with the vertical red, dashed line. The stellar mass complete data set comprises 393 galaxies, members of A963. The middle panel “b” of Fig. 4 shows the distribution of K-band luminosities of the cluster members. The red dashed line shows the magnitude limit for target selection of the ACReS survey K = 16.964, equal to K + 2.0. The data used in this work is 89% complete for galaxies brighter than this limit.

thumbnail Fig. 4.

Panel a: stellar mass distribution of the cluster members with the completeness limit indicated by the red, dashed line. Panel b: completeness of our data set with respect to K-band luminosity. The red dashed line is at K = K+2.0- the magnitude limit of ACReS survey. Panel c: completeness distribution on the sky. The green circle has a radius equal to R200 of the cluster.

The presence of AGN at the centre of a galaxy could produce a spectral signature similar to that of ongoing star formation. We use the so-called BPT diagram (Baldwin et al. 1981) to detect emission line galaxies with ionisation state of their interstellar matter (ISM) consistent with ionisation from AGN. The method uses a combination of two line ratios, [O III] 5007 Å/Hβ and [N II] 6583 Å/Hα. A line separating the AGNs from star forming galaxies from Kewley et al. (2001) was used. In total we found seven AGNs among the cluster members. These galaxies were not considered in the following analysis.

In this work we are interested in the properties of galaxies in different parts of the cluster, as projected on the sky. In panel c of Fig. 4 we show the completeness on the sky of the data set used in this work. The spectroscopic completeness is shown as a fraction of all the galaxies with K ≤ 16.964 and JK colour consistent with the main sequence of the cluster. The lowest spectroscopic completeness of our data set is 58%. More than 84% of the sky area covered by this survey has spectroscopic completeness greater than 0.8.

Figure 5 shows the specific star formation rate (SSFR, defined as the star formation rate per unit stellar mass) of the cluster members as a function of their stellar mass. The dashed, black line at constant log(SFR) = −0.5 is used to divide the galaxies into star forming and passive, shown with stars and circles, respectively. Each symbol is colour-coded according to that galaxy’s Dn4000 value (see Sect. 2.3 for details), which is a good indicator of the age of the light-dominating stellar population. The error bar in the corner shows the median 1σ confidence limits on the stellar mass and specific star formation rate from the χ2 fitting with MAGPHYS.

thumbnail Fig. 5.

Specific star formation rate of the A963 member galaxies as a function of their stellar mass. The black, dashed line at constant log(SFR) = −0.5 is used to split the galaxies into star forming, shown with stars, and passive, shown with circles. Symbols are coloured according to the strength of the 4000 Å drop in luminosity. The symbol in the upper right corner shows median uncertainties.

2.3.1. Dn4000: Recently quenched and red star forming galaxies

Passive galaxies exhibit a characteristic break in their continuum emission blueward of 4000 Å. This is due to accumulation of metal absorption lines coming from the atmospheres of low-mass stars and is known to correlate well with the age of the light-dominating stellar population (Bruzual 1983; Poggianti & Barbaro 1997; Balogh et al. 1999). We use the narrow definition of the index Dn4000 (Balogh et al. 1999) for its robustness against reddening. The mean uncertainty of the individual Dn4000 measurements is 0.02. This number is confirmed by comparing the measurements of Dn4000 of galaxies for which more than one observation is available. In these cases the final value of the index is a weighted average of the separate measurements. The Dn4000 index was successfully measured for 381 of the 386 cluster members. The remaining five galaxies only yielded a lower limit to Dn4000 due to weak or no signal in the blue window of the index, and we did not use them in the following analysis.

At redshift of A963 (and in the local Universe), galaxies form two well-defined sequences of Dn4000 versus stellar mass (Geller et al. 2014, 2016; Haines et al. 2017; Deshev et al. 2017), one characterised by old stellar populations and a second containing star forming galaxies whose blue light is dominated by young stars. We used the entire ACReS data set (∼ 27000 galaxies), divided into passive and star forming galaxies, according to the absence or presence of emission lines in their spectra, to estimate the mean Dn4000 for passive galaxies as a function of stellar mass. Figure 6 shows the distribution of offsets from this relation for all A963 member galaxies split into passive and star forming as described in Sect. 2.3, and Fig. 5. The two samples show two distinct distributions with the minimum between them shown by the grey dotted line. We use this information to define two classes of galaxies which will be analysed in Sect. 3. All the passive galaxies to the left of the grey dotted line in Fig. 6 (higher negative offsets) have Dn4000 more typical of star forming galaxies and must have had their star formation quenched within the last 0.5 Gyr, ∼1/3 of the cluster crossing time (assuming σcl = 970 km s−1 and R200 = 1.85 Mpc, J16). We refer to these galaxies as recently quenched (RQ). All the star forming galaxies that are to the right of the grey line have Dn4000 more typical of passive galaxies, indicating that the evolved stars dominate their spectral energy distribution (SED). We refer to these galaxies as red star forming (RSF). In total we find 28 RQ and 35 RSF galaxies, constituting 11% and 27% of the passive and star forming galaxies, respectively.

thumbnail Fig. 6.

Distribution of Dn4000 values of A963 member galaxies. Plotted is the offset from the mean relation between stellar mass and Dn4000 for all the passive galaxies in the ACReS survey.

The value of Dn4000 is sensitive to the presence of even a relatively small amount of stars hotter than A-type, which defines the characteristic timescale on which the index reacts to changes in the star formation properties of the galaxy. We used GALEV evolutionary synthesis models (Kotulla et al. 2009) and Bruzual & Charlot (2003) spectral evolution models to simulate the evolution of the SED of galaxies depending on their ongoing and past star formation. Passive galaxies with low Dn4000 values (RQ) are formed only when the truncation of star formation is very rapid or accompanied by a burst of star formation. In the case of high stellar mass galaxies the latter is often required. The modest specific star formation rates typical for galaxies with log(M) ≥ 10 do not usually produce enough massive stars and their integrated light is dominated by the old stellar population. As a consequence their Dn4000 remains high in the absence of a strong star burst. This issue is further strengthened by the fact that the optical spectroscopy data used here is collected through a single fibre per galaxy. These fibres have limited spatial coverage and are more likely to miss the outskirts of the target galaxies which is often the location with the strongest star formation. A burst producing as little as 1% of the final stellar mass of the galaxy is sufficient to lower its Dn4000 to levels more typical of star forming galaxies for a period of a few ×108 yrs. The rapid quenching and the quenching accompanied by a star burst are expected from galaxies subjected to ram-pressure stripping, where the star bursts occurring in their discs are more centrally concentrated (Jachym et al. 2019). In contrast, galaxies subjected to starvation do not show this decrease in Dn4000. Their Dn4000 quickly reaches values typical of passive galaxies even when they still host some residual star formation (like the RSF class). The lifetime of the RQ phase, defined this way, is much longer than the timescales over which other signatures of star formation, such as optical emission lines, disappear from their spectra (Kennicutt & Evans 2012). This allows us to detect galaxies that are likely subjected to ram-pressure stripping and separate them from the ones subjected to starvation. The typical timescales associated with the latter are ∼4 Gyr (Peng et al. 2015). The strength of the ram-pressure stripping has strong dependence on the orbit of the galaxy, its orientation, and total mass, among other factors, which define a relatively wide range of timescales on which star formation is expected to be quenched. Tonnesen et al. (2007) estimate between 0.5 and 1.5 Gyr, all significantly faster than starvation, and broadly consistent with the timescales over which Dn4000 reacts to changes in the star formation rate.

2.4. 21 cm measurements

The measurements of the 21 cm emission from atomic hydrogen from galaxies in the field of A963 is based on the data from the BUDHIES survey (Verheijen et al. 2007, 2010; Deshev et al. 2009; Jaffé et al. 2013, J16, Deshev 2019; Gogate et al. 2020), acquired with the Westerbork Synthesis Radio Telescope (WSRT2). The measurements for A963 were collected between 2005 and 2008 and amount to a total of ∼1400 hours of integration. This was the first of only a handful of blind interferometric H I surveys at such high redshift, and as such represents a unique data set which allows us to estimate the H I content of individual A963 members, and in stacks of spectra over separate parts of the cluster. This estimation was published by J16 as a function of cluster centric distance. In Sect. 3.3 we present it as a 2D map on the sky.

In order to make the measurements, we extract an H I spectrum at the position and redshift of each star forming cluster member galaxy and stack them within individual Voronoi bins. The extracted spectra spans ±850 km s−1, from a data cube smoothed in frequency to a resolution of 312.5 kHz (∼80 km s−1). The flux S, integrated over the central ∼500 km s−1 of the stacked spectra is converted to H I mass as MH I = 2.36 × 105× d2 × Sdv, where d is the distance and dv is the average channel separation, in this case dv = 40 km s−1.

3. Results

3.1. Distribution of star formation in A963

In Fig. 7 we show a map of the distribution of local fraction of star forming galaxies fSF calculated in Voronoi bins on the sky. The positions of individual galaxies are also indicated with different symbols. Circles with red borders indicate passive galaxies, green for the recently quenched ones (RQ) and white for the rest. Stars with blue borders indicate the star forming galaxies, magenta for the anemic ones (RSF) and cyan for the rest. The green circle in the middle has a radius equal to R200 of the cluster. The black contours show the extended X-ray emission from the cluster and other nearby structures, some of which are considered to be groups of galaxies infalling for the first time onto A963 (Haines et al. 2018). The orientation of the image and the azimuthal angle α, to which we refer in the following sections, are indicated along the large grey circle.

thumbnail Fig. 7.

Distribution of A963 member galaxies and their properties on the plane of the sky. The passive galaxies are shown as red bordered circles, with the recently quenched ones (RQ) in green, the rest in white. Star forming galaxies are shown as blue bordered stars, with the anemic ones (RSF) in magenta and the others in cyan. The local fraction of star forming galaxies, fSF, mapped on the Voronoi bins is colour-coded according to the colour bar. The green circle in the centre has a radius equal to R200 of the cluster. The black contours show the extended X-ray emission. The orientation of the image is shown with the large grey circle, showing the azimuth α as well as the directions towards north and east.

In order to quantify the deviation from spherical symmetry of the cluster and its infall region, we employ a model published by Haines et al. (2015), based on “observations” of the 75 most massive clusters from the Millennium simulations (Springel et al. 2005). The model is shown in Fig. 8. Galaxies falling onto the average of these 75 clusters have had their star formation quenched instantly at a certain time after their first passing through R200 of the cluster. This time delay is shown at the top of the figure, and gives rise to the sequence of thin lines with corresponding colours. The solid black points with error bars show the fraction of star forming galaxies in A963 as a function of cluster centric distance. The same fraction, but measured in Voronoi bins, is shown as empty squares representing the bins with their K-band luminosity weighted centroids. The two sets of points are consistent with each other. The dashed line is a fit of a third-order polynomial to the binned data. The normalisation of the model curves is set so their outermost points have the same fSF as the binned observations at this radius. At r >  1.5 R200 the data is not in agreement with the models even within the Poissonian error bars. Inside that radius the data does not follow a trend similar to that of all the models. This is a strong indication that this cluster is not spherically symmetrical in terms of fraction of star forming galaxies. It should be noted that A963 shows a radial trend that is different from the average ACReS cluster, which shows a steady increase in fSF with distance from the centre (Haines et al. 2015).

thumbnail Fig. 8.

Fraction of star forming galaxies as a function of cluster-centric distance in simulated clusters and in A963. Thin lines of different colour represent galaxies quenched with a delay after their crossing of R200 of the parent cluster, according to the colour bar at the top of the figure. The data for A963 are shown in black (symbols and line). The fSF values calculated for individual galaxies are shown as black points, and the mean fSF values within each Voronoi bin are shown as empty squares. The dashed line shows a third-order polynomial fit to the Voronoi bins. For clarity, a small horizontal offset has been added to the black points.

By minimising the difference between the range of models and the fit to the measurements in Voronoi bins we find that a delay of 2.17 Gyr is closest to the data, in agreement with Haines et al. (2015). If only Voronoi bins with r ≤ R200 were used in the fit a higher delay would be preferred to match the unusually high fraction of star forming galaxies in the centre of this cluster. We use the model with a delay of 2.17 Gyr to populate the Voronoi bins with star forming galaxies, according to the cluster-centric distance of the K-band luminosity weighted centroids of each bin, resulting in the left panel of Fig. 9. The right panel of the figure shows the ratio between the data for A963, presented in Fig. 7, and this spherically symmetrical model. Redder colours represent a lack of star forming galaxies, bluer ones an excess. The model, presented in the left panel of Fig. 9, contains fSF = 0.19 for bins with r ≤ R200, in perfect agreement with Butcher & Oemler (1984), who were the first to estimate the unusually high fraction of star forming galaxies hosted by this cluster. The data deviates significantly from the model containing on average, for bins with centroids at r ≤ R200, 45% higher fSF. However, larger deviations are visible in individual bins. The central bin contains 6 times the expected fSF from the model. To the north and south of it there are bins with 3 and 2.2 times the expected fSF. Outside R200 the deviations from the model showing lack of star forming galaxies are usually associated with overdensities shown with the black contours. Bin number 16, directly east of the cluster centre, contains 70% more star forming galaxies than expected from the model.

thumbnail Fig. 9.

Left: model from Fig. 8, which minimises the difference with the binned data, with a delay of 2.17 Gyr, converted to 2D distribution and used to populate the Voronoi bins with star forming galaxies, according to their K-band luminosity weighted centroids (black points). Right: map of fSF, shown in Fig. 7, divided by the spherically symmetrical model shown on the left. Blue bins show an excess of star forming galaxies, red ones show a lack. The black contours show the K-band luminosity density. The lowest contour is at 40 K Mpc−2. Each following contour is at twice the level of the previous one. The green circle has a radius equal to R200 of the cluster.

There is more than one way of performing adaptive binning of a 2D distribution of points, depending on the distribution of points, the size of the initial grid, and the required population of the Voronoi bins (signal-to-noise ratio). In order to estimate how robust these results are against variations in the sampling of the cluster members and their consecutive binning, we performed 1000 jackknife realisations and present the results in Fig. 10. For each run we randomly select 90% of the 386 cluster members and perform the Voronoi binning on this subsample. Then we calculate the fraction of star forming galaxies per bin. The mean and standard deviation of those 1000 maps is shown in the left and right panels of Fig. 10, respectively. The mean shows qualitatively the same picture as Fig. 7. Because the bins are recalculated every time and because their borders move, the large bin-to-bin variations visible in Fig. 7 are smoothed out. This effect is particularly strong closer to the centre of the cluster where the bins are smaller, due to the increased density of galaxies. This presentation of fSF makes the asymmetrical distribution of galaxy properties in A963 more pronounced, and the association of the concentrations of passive galaxies outside R200 with overall galaxy concentrations is more clear. Overall, the mapped region outside R200 is split into two distinct parts. To the north and west of the cluster, at azimuth 45 < α < 250, the fraction of star forming galaxies is fSF≃ 0.34. On the opposite side fSF≃ 0.55 with separate regions at 0.6. Inside R200 there is similar bimodality with star forming galaxies dominating the south-eastern part of the cluster as close to the cluster centre as 0.5 R200.

thumbnail Fig. 10.

Left: map showing the mean fraction of star forming galaxies in Voronoi bins for 1000 jackknife runs on randomly selected 90% of the data set. Right: standard deviation of fSF in these 1000 runs. The black contours show the K-band luminosity density, as in Fig. 9. The green circle has a radius equal to R200 of the cluster.

We performed tests by varying the size of the initial grid cells and the required S/N in each Voronoi bin and found the results of the jackknife resampling to be very robust against such variations. The right side of Fig. 10 shows the standard deviation of the 1000 jackknife runs. The median standard deviation is 0.1. The real uncertainty is higher as this number does not take into account the uncertainty in estimating star formation rate (median 0.35 dex). We estimate that this variation in the star formation rate would change the overall fraction of star forming galaxies by ≃0.1. The uncertainty on the stellar mass (median 0.07 dex) does not play a role since we are separating galaxies using a constant SFR threshold.

We mention one caveat associated with the varying bin size over the image. This makes the resolution of the image also not constant, even though clear variations are visible among neighbouring pixels. Particularly in the outskirts, where individual Voronoi bins change relatively little between the separate jackknife runs, the resolution of the map is much lower than the pixel size.

Figure 11 shows the distribution on the sky of the specific star formation rate in A963. All 386 cluster members were used to produce this map. Qualitatively this map is very similar to the fSF map presented in Fig. 10 but with a few key differences. The same bimodal behaviour is observed outside R200 of the cluster, but within it there is a much more pronounced region of reduced SSFR, running north–south through the cluster and extending out to the edges of the map.

thumbnail Fig. 11.

Local mean (left) and standard deviation (right) of the specific star formation rate in A963 cluster members. The black contours show the K-band luminosity density as in Fig. 9. The green circle has a radius equal to R200 of the cluster.

3.2. Recently quenched and red star forming galaxies

The left panel of Fig. 12 shows the distribution on the sky of RQ galaxies as a local fraction of the passive galaxies. As in Figs. 10 and 11 the mean fRQ from 1000 jackknife runs is shown. The right panel shows the RSF galaxies as a local fraction of the star forming galaxies. The black contours show the luminosity density of cluster members. The two distributions are very different. While the RSF galaxies tend to cluster around the high-luminosity density regions, the RQ galaxies appear to avoid them. The right panel also shows that the concentrations of star forming galaxies inside R200 seen in Figs. 9 and 10 are in fact dominated by RSF, in agreement with Fig. 11. Outside R200 the RSF distribution is very similar to the distribution of passive galaxies. There is only one region with elevated luminosity density containing a high fraction of RQ galaxies. The east side of the map is richer in RQ galaxies, and the area directly east of the cluster is where RQ galaxies reach most closely to the cluster centre, as do the star forming galaxies shown in Figs. 10 and 11.

thumbnail Fig. 12.

Recently quenched galaxies as a local fraction of the passive galaxies (left) and red, star forming galaxies as a local fraction of the star forming galaxies (right). The black contours show the K-band luminosity density as in Fig. 9. The green circle has a radius equal to R200 of the cluster.

3.3. H I deficiency

The left part of Fig. 13 shows a map of the H I deficiency of the star forming galaxies, members of A963. The map is a result of 1000 Voronoi binnings of randomly selected 90% of the mass complete sample. For every new calculation of the Voronoi bins, we extracted a spectrum from the H I data cube at the location and redshift of all the star forming galaxies in a given bin. These were then brought to the same redshift and averaged. The flux within ±250 km s−1 of the optical redshift of the galaxy was summed and converted to M H i observed $ M_{Hi}^{{\rm{observed}}} $ as explained in Sect. 2.4. In order to estimate how H I deficient each galaxy is, this M H i observed $ M_{Hi}^{{\rm{observed}}} $ was then compared with the expected H I mass- M H i expected $ M_{Hi}^{{\rm{expected}}} $, calculated with a recipe by Haynes & Giovanelli (1984), based on the optical size of the galaxy, assuming that all the star forming galaxies are of late-type morphology. The optical sizes of A963 galaxies were measured with SExtractor (Bertin & Arnouts 1996) as the radius containing 90% of their light in R band, from the INT mosaic imaging (Sect. 2.1). Haynes & Giovanelli (1984) also provide a way of estimating the expected amount of gas for early-type galaxies. We decided not to include the passive galaxies in this analysis because they contain very small amounts of gas and adding them has little effect on the map presented in Fig. 13. The map shows clearly the central part of A963 as well as the regions to the north and south of the cluster as dominated by galaxies highly deficient in H I. The luminosity density contours trace the regions with high deficiency. From almost all directions, except directly south, less H I deficient galaxies can be found well inside R200. Outside this radius the least H I deficient regions are directly east and west of the cluster centre.

thumbnail Fig. 13.

Left: map of the local H I deficiency of star forming galaxies. Lighter colours indicate high deficiency. This is derived from all star forming galaxies. Overlaid are symbols indicating the location of galaxies that are individually detected at 21 cm. The sizes and colours of the symbols are proportional to the offset between the systemic velocities of the gaseous and stellar components. The black contours show the K-band luminosity density, as in Fig. 9. The green circle has a radius equal to R200 of the cluster. Middle: total H I maps of three selected sources overplotted on the R-band image. The contours are drawn at 20, 40, 80, 160, and 320 μJy beam−1. The cluster members with optical spectroscopy are shown in red, and their redshift is indicated. The non-member galaxies with optical spectroscopy are shown in purple. The central redshift of the H I signal is indicated in the top left corner. The beam size is shown in the lowest panel. Right: renzograms of the selected sources. Contours at a constant significance level of 3σ are drawn for every channel containing signal from the source. The noise is estimated from the surrounding area of the sources as a function of frequency. The difference between the central velocity of the channels and the central velocity of the entire H I signal is colour-coded. The proposed optical counterparts are indicated with red arrows.

Overlaid on the H I deficiency map are symbols showing the positions of the galaxies detected individually at 21 cm. In total the BUDHIES survey of A963 detected 127 galaxies (Gogate et al. 2020). Of these 50 are within the velocity range of the cluster and their proposed optical counterparts are within the stellar mass limits of the sample analysed here. The sizes and colours of individual symbols in Fig. 13 are proportional to the offset in systemic velocity between the H I gas and the stellar component. The combined uncertainty of the two redshift estimates is < 50 km s−1 (in red) and ∼71% of the I detected galaxies are below this value. The remaining high offsets could be the result of confusion, although some show characteristics consistent with the expectation for ram-pressure stripping (discussed in the next section). The confusion is driven by the large size of the synthesised beam of WSRT (shown in the bottom panel of the middle column of Fig. 13), and by the limited availability of optical spectroscopic data. This prevents us from uniquely identifying the optical counterparts of the H I detections. The plotted offsets are between the central velocity of the H I profile and the optical redshift of the nearest galaxy with available optical redshift. In addition it is very likely that some individual H I detections have multiple optical counterparts.

3.3.1. Ram-pressure stripping

A few galaxies were found to have significant offsets in their HI and optical redshift measurements. This difference in velocity may be indicative of ram-pressure stripping. Simulations by Tonnesen & Bryan (2010) find that in strong ram-pressure stripping, the bulk of the stripped gas is accelerated up to ∼60% of the wind velocity over a period of 250–500 Myrs. Given the expected infall velocity of ⪆1000 km s−1 (Haines et al. 2015; Barsanti et al. 2016), this would seem to be a natural explanation for the high velocity-offsets.

The six panels on the right of Fig. 13 shows three galaxies which are among the most plausible candidates for ram-pressure stripping in the A963 cluster (we leave a fuller analysis of the rest of the sample to a future work). As well as the strong velocity offsets, all three show non-circular contours (the renzograms show contours at a constant 3σ significance level, colour-coded by velocity channel, following Rupen et al. (1999) and Kregel et al. (2004). That the contours are non-circular and are found well away from the candidate parent galaxy is strong evidence that the H I is at least partially extended into extragalactic space and not directly associated with optical emission. Simulations by Taylor et al. (2017) have shown that displacing large amounts of H I through tidal encounters, without causing easily visible disturbances in the stellar disc, is extraordinarily difficult, whereas the nature of ram-pressure has no difficulty in removing gas while leaving the stellar disc unaffected. Taken together, this is strongly suggestive of a ram-pressure origin. If so, this would be the first direct imaging of ram-pressure stripping at HI wavelengths at this redshift.

There are two significant caveats to this interpretation. First, due to the large physical size of the WSRT beam at this redshift, we cannot rule out contributions from other galaxies within the beam, which could create the illusion of non-circular contours. Further optical spectroscopic redshift measurements are necessary to fully rule this out, and this scenario is especially plausible for the galaxy in the lower section of Fig. 13. However, this is much less likely for the other two cases. Given the sensitivity of the H I observations we would expect to only detect the brightest galaxies visible in each image, and while there may be some contribution from the smaller galaxies visible in the upper section, the middle galaxy shows no obvious possible source of confusion.

The second caveat is that, even though it is an indicator of a disturbance, the magnitude of the velocity offset between the H I and optical is much larger than might be expected. Local galaxies have been shown (Oosterloo & van Gorkom 2005; Chung et al. 2007; Koopmann et al. 2008; Leisman et al. 2016; Hallenbeck et al. 2017; Minchin et al. 2019; Taylor et al. 2020) to usually have only a small fraction of their H I present in their tails, and retain the bulk of their gas in close proximity to their optical discs. Here, the velocity offsets (200–400 km s−1) are high enough that the gas may have been entirely displaced from the parent galaxy. While this is highly unusual, it is not completely unknown for local galaxies. For example VCC 1249 in the Virgo cluster has an H I component that is entirely displaced from the stellar emission and a velocity difference of ∼200 km s−1 (Sancisi et al. 1987; Henning et al. 1993; McNamara et al. 1994; Arrigoni Battaia et al. 2012). The galaxies described here are much more massive and we therefore expect ram-pressure to be stronger in order to remove their gas, likely resulting in different effects than more typical cases that show only faint H I tails. One example of these differences is shown by tuned simulations (Vollmer et al. 2001), which indicate that in massive galaxies the stripped gas can fall back on the stellar disc after the stripping event. Nonetheless, detailed modelling and additional optical spectroscopy will be necessary to properly confirm the nature of these objects.

4. Discussion

4.1. Connection to the large-scale structure

The results outlined in the previous section suggest a strong dependence of galaxy properties on the direction from which they are infalling onto A963. Naturally we want to trace the infalling galaxies to large cluster-centric distances. However, in the absence of new, more complete spectroscopic data, than used by J16, it is difficult to expand on their conclusions regarding the large-scale structure surrounding A963 and its connectivity.

Jaffé et al. (2016) proposed the existence of a network of three filaments of the large-scale structure (LSS) crossing at the location of A963 (their Fig. 7). Two of them, NE and SE, are relatively tight and well defined; another, W, somewhat more disperse on the sky, although equally coherent in velocity as the first two. This is based on galaxies with available spectroscopic redshift from SDSS. At the redshift of A963, SDSS is complete only for the most massive, usually passive galaxies which, as pointed out by J16, makes the mapping of these filaments tentative at best. The existence and location of these filaments finds a strong support in the entropy map shown in Fig. 14, based on XMM-Newton observations. The figure shows a pseudo-entropy map, obtained using a calibrated hardness ratio map that traces the temperature of the cluster’s ICM and the soft band map that traces its density. The arrows show the approximate directions of the proposed three filaments with the W filament bracketed by the two thin arrows. The inflow of low-entropy gas along filaments is observed, at the expected azimuthal angles (Finoguenov et al. 2005). Finoguenov et al. (2006) also observed similar low-entropy structure in A3266 likely related to recent accretion of groups and low-mass clusters.

thumbnail Fig. 14.

Entropy map of the hot gas of A963, based on XMM-Newton observations. The arrows show the approximate direction of inflow expected from the filaments identified by J16 from the large-scale galaxy distribution around the cluster. The × symbol in the centre marks the location of the BCG. The green circle has a radius equal to R200 of the cluster.

4.2. Non-spherical cluster

Clusters of galaxies do not have well-defined shapes as the density distribution of their constituent matter drops progressively as a function of distance to the centre without a well-defined edge (Einasto 1965; Navarro et al. 1997). The distribution of the gravitational mass within the clusters is traceable by strong and weak lensing, but only close to the cluster centre and with considerable uncertainties, those associated with the distribution of background objects, the distribution of their intrinsic orientation, and the presence of other nearby structures along the line of sight (Hwang et al. 2014). Estimating shapes of clusters of galaxies is further complicated by the irregular distribution of the galaxies. The X-ray emitting gas present in the interiors of the clusters can be used as an indicator, but it has a relatively limited radial extent and becomes unobservable at relatively small distances from the centre (Chambers et al. 2002). The distribution of intra-cluster light is known to follow a similar distribution to that of the intra-cluster gas and the overall mass distribution (Montes & Trujillo 2019), but it is also only observable close to the cluster centre. Despite these difficulties the non-sphericity of clusters has long been established observationally (Binggeli 1982; Plionis 1994; Kasun & Evrard 2005; Donahue et al. 2016; Einasto et al. 2018a), and in simulations (Okabe et al. 2018).

We focus on mapping the parts of the clusters where galaxies are being transformed, and so in this section we discuss the shape of A963, as indicated by the varying properties of its member galaxies. This is not necessarily directly related to the density of galaxies (Fig. 9), nor to the distribution of the X-ray emitting gas (Fig. 7). Inside R200 the strongly asymmetrical distribution of the fraction of star forming galaxies is clearly visible in Figs. 7 and 10, and is also supported by the distribution of specific SFR (Fig. 11). Directly east of the cluster centre, between R500 and R200, ∼50 ± 10% of the galaxies are star forming, unlike in any other direction. Within this radial range the fraction of star forming galaxies varies between 10% and 50% with a typical uncertainty of ±10%. Even though the binning does smooth the distribution, the jackknife resampling shows that this result is significant. Apart from the cluster centre, the region to the south-west of it inside R200 is the most passive one with only ∼14% of the galaxies there showing signs of ongoing star formation. This is likely associated with the infalling group discovered by its X-ray emission (named Grp10 by Haines et al. (2018) and group C by J16) and also visible on the luminosity density plots.

Studies have shown a strong correlation between the shape and orientation of the cluster and those of the brightest cluster galaxy (Fuller et al. 1999; Kasun & Evrard 2005). In addition, the orientation of the cluster is often aligned with the large-scale structure surrounding it (Binney & Silk 1979; Kasun & Evrard 2005; Faltenbacher et al. 2005) and other nearby structures up to 30 Mpc away (Binggeli 1982). In cases where the studied cluster is embedded in a massive supercluster its radio and X-ray emission, BCG and infalling substructures are not only aligned with the cluster itself but also with the supercluster axis (Einasto et al. 2018a,b). We used GALFIT (Peng et al. 2002) to model the BCG galaxy of A963 with a double Sersic profile (Sersic 1968) finding an axis ratio of the outer isophotes close to 0.5 with a position angle pointing north (α ≃ 0). The exact value of α changes with surface brightness from 20 to 340 but remains consistent with the alignment of ∼20° projected on the sky, expected from analysing the orientations of dark matter haloes and their central galaxies in Horizon-AGN cosmological hydrodynamical simulation (Okabe et al. 2020). On the other hand, the X-ray emission shows very small ellipticity.

4.3. Anisotropic accretion

Before discussing the galaxy properties in the infall region of A963 we note that because of the decreasing galaxy density away from the cluster centre beyond ∼1.5 R200, the binning method we use only provides information as a function of azimuth, with all radial trends masked out by the binning. The only part of the region outside ∼1.5 R200 where we resolve the cluster in the radial direction is at azimuth 170 < α < 225, where the members of the infalling groups visible in Fig. 2 are raising the galaxy density.

In Fig. 15 we summarise the properties of the galaxies in the infall region of A963 as a function of azimuth. All curves in the figure, except in the middle panel, show the average values over the said radial extent as a function of azimuthal angle, and are normalised to their peak value. The grey shaded regions show the approximate location and azimuth span of the three filaments suggested by J16. The top panel shows the density of galaxies in the infall region of the cluster (1.5 R200 <  r≲ 3 R200) with a solid black line, and at 10 R200 < r <  20 R200 with a lime green dashed line. The former is derived by the K-band luminosity density shown in Fig. 2; the latter is a projected number density of SDSS DR12 spectroscopic targets within the velocity range of the cluster.

thumbnail Fig. 15.

Properties of the galaxy population in the vicinity of A963 and its infall region as a function of the azimuthal angle. Top panel: galaxy number density normalised to its maximum. The black curve shows the average K-band luminosity density in the infall region (between 1.5 R200 and the edge of the field at ≲3 R200). The lime green dashed line shows the mean galaxy number density between 10 and 20 R200. Middle panel: H I deficiency in the infall region. Bottom panel: specific star formation rate (solid purple line), fraction of RQ (dashed green line), and fraction of RSF (red dotted line), expressed as a fraction of their maximum. The grey shaded areas represent the location and relative angular span of the three filaments suggested by J16.

The overall variations in density shown by the lime green line are approximately a factor of two from the peak value with three distinct peaks. One coinciding with the NE filament and the other two located close but not exactly overlapping with the SE and W filaments. The black curve shows much greater variation in density which deviates more from the proposed picture with three filaments.

The middle panel shows the H I deficiency in the infall region of the cluster. An artificial saturation value in deficiency is introduced at 1.4 as any values above that are beyond the sensitivity of the survey and are driven by noise. This curve bears very little resemblance to any of the density curves in the top panel and the location of the filaments. The entire northern half of the cluster is devoid of H I. While the southern overdensity is also very deficient, it does contain galaxies with reduced deficiency. The two biggest patches of reduced deficiency while mostly associated with low-density regions also contain higher density regions.

The bottom panel shows the fractions of recently quenched and red star forming galaxies as well as the SSFR of all the galaxies. All the curves are normalised to their peak values. The SSFR peaks, as expected, at the location of the lowest H I deficiency. However, the W filaments which also seem to show decreased H I deficiency does not bring a corresponding increase in SSFR. The overall variation in SSFR is a factor of four. The green dashed curve showing fRQ follows closely the fSF curve, although with somewhat larger variation probably driven by the relatively small number of RQ galaxies. The fRSF curve is almost exactly a negative of the fRQ one, as noted in the discussion of Fig. 12. The northern part of the cluster seems to offer a very consistent picture of overdensity, traceable out to ∼20 R200 and populated by gas-poor passive galaxies, few of which have been quenched recently; the galaxies that are not passive host a very low level of star formation. The rest of the infall region is much more difficult to interpret. In particular, the southern overdensity, which at least partially consists of multiple groups (see Fig. 2) which contain both passive and star forming galaxies showing a range in H I deficiency, and which have probably arrived at the southern doorstep of the cluster from different directions, considering the presence of three filaments none of which points directly south.

We can try to understand this by comparing with results from simulations, published by Haines et al. (2015). By observing the 75 most massive clusters in the Millennium simulations, they estimate the fraction of backsplash galaxies in the infall region of these clusters. When defined as the galaxies that have already been inside R200 of the cluster, it is expected that ∼30% of the galaxies located at 1.5 R200 projected from the cluster centre are backsplash galaxies. This fraction drops to 0% at ∼3 R200. The varying intensity and azimuth of the accretion of matter onto clusters is also a function of time. A similar picture is observed in the ten most massive clusters from Millennium simulations, shown in Fig. 9 of Haines et al. (2015). When considering a single cluster the accretion is likely much more strongly varying in rate and direction with time. In A963 the main accretion happens along the north–south axis. We can expect that the east–west axis would contain a lower fraction of backsplash galaxies, which can partially explain the overdensity of unprocessed galaxies there. Part of the overdensities along the north–south axis is likely due to backsplash galaxies, although many of the individual groups in the south of the cluster are detected in X-ray and by Dressler-Shectman test (Dressler & Shectman 1988) presented in J16, which is unlikely for backsplash groups as those would be stripped of their intra-cluster gas and would likely be dispersed during their pericentric passage. It is possible that an accretion of a number of groups at an earlier epoch, along the north–south direction has given a long lasting asymmetric shape of the outskirts of A963. Prolonged anisotropic accretion could mimic the effects of pre-processing, when considering individual clusters, because it will likely lead to strongly asymmetric distribution of backsplash galaxies in the infall region. However, the deviations from the spherically symmetrical model shown in Fig. 9 are too large to be entirely due to backsplash galaxies. The northern overdensity is also traceable out to at least 20 R200. The presence of red star forming galaxies in this region (Fig. 12) is a sign that pre-processing is responsible for the quenching of these galaxies, because backsplash galaxies are not expected to have any ongoing star formation.

While the origin of the passive galaxies situated north of the cluster centre is relatively uncertain, there is little doubt about the gas-rich population located east of the cluster. A963 is an X-ray bright cluster and, as demonstrated by J16, it is unlikely that a galaxy would still host any detectable amount of H I after passing through the pericentre of its orbit. Hence these galaxies must be on their first infall towards the cluster. This is supported by the increase in the fraction of RQ galaxies observed at the same location. The fact that we observe an increased fraction of star forming galaxies well inside R200 from the east does not mean they do not experience ram-pressure stripping. We note that the synthesised beam of the WSRT, used for the measurements of the 21 cm emission, is ∼ 50 × 70 kpc at the distance of A963, while the Hectospec fibres used for optical spectroscopy cover the central ∼5 kpc of the galaxy. This means that in most cases the gas disc of the infalling galaxies can be offset significantly without leaving the WSRT beam. The optical measurements on the other hand only accept light from the central region of the galaxies, which is where the star formation continues the longest before quenching due to ram-pressure stripping. For such galaxies the stripping process could be well under way while we still detect them as unaffected. The two galaxies showing strong velocity offset between the stars and gas in this region support the idea of them being subjected to ram-pressure stripping. We should also add that the initial phase of ram-pressure stripping is likely accompanied by an increase in star formation, particularly in the central parts of the disc (Vulcani et al. 2018).

5. Conclusions

We used the optical spectroscopy data from Haines et al. (2013) and Hwang et al. (2014), together with K- and R-band imaging, X-ray imaging, and the deepest 21 cm blind survey (Gogate et al. 2020) to try and understand the origin of the galaxies responsible for the strong Butcher-Oemler effect observed in A963, as well as the role played by anisotropic accretion in shaping its galaxy population. We present a series of maps showing the distribution on the sky of galaxy properties in and around A963. The cluster shows a strongly asymmetrical structure both in its interior (r ≤  R200), and in its infall region (1.5 R200 <  r <  ∼3 R200). These deviations from spherical symmetry cannot be explained by cluster related processes which are expected to show simpler dependence on cluster-centric distance, given the very symmetrical X-ray emission whose peak coincides with the location of the BCG. The reason thus must be related to the large-scale structure surrounding A963 and anisotropic accretion of galaxies both in terms of their number and how pre-processed they are. We find an influx of gas rich, star forming galaxies preferentially from east and west of the cluster centre. They reach deep into the cluster interior, particularly from the east side, and are responsible for the strong Butcher-Oemler effect A963 shows. We also observe signatures of ram-pressure stripping in those areas, and the few passive galaxies there have likely had their star formation quenched recently. The accretion perpendicular to that, from the north and south, is dominated by mostly passive galaxies, often clustered in groups. The few star forming galaxies in these regions show very low specific star formation rates and the passive ones do not show signs of any recent star formation. The connection with the large-scale structure surrounding A963 is, however, inconclusive, being data limited. At cluster-centric distances r <  3 R200 the galaxy properties suggest strong accretion along the north–south axis of pre-processed galaxies and a weaker accretion of unprocessed galaxies along the east–west axis. At distances > 3 R200 three filaments are visible in the overall galaxy distribution, as suggested by J16. We present strong support for the presence of the three filaments in the form of an entropy map, but the properties of the galaxies they bring vary with the western filament bringing into the cluster less processed galaxies than the NE and SE ones. Alongside this we also present the first, tentative, detection at 21 cm of ongoing ram-pressure stripping at z = 0.2.

While this analysis, together with the previously published projected phase-space analysis of A963 (Jaffé et al. 2013, J16), brings us closer to understanding the processes that define the make-up of a massive, intermediate redshift cluster, many questions still remain unanswered. What is the exact connection of A963 with the large-scale structure surrounding it? Why do areas with similar density in the infall region of the cluster seem to bring in galaxy populations with vastly different star formation properties and gas content? The data of the galaxy distribution on larger scales suggests a fragmentation of the filaments around this cluster; is this real or a limitation of the data? Lastly we would like to stress the importance of considering the full dimensionality and all the information available from astronomical observations, when analysing clusters of galaxies.


1

Observations reported here were obtained at the MMT Observatory, a joint facility of the University of Arizona and the Smithsonian Institution.

2

The WSRT is operated by the Netherlands Foundation for Research in Astronomy(NFRA/ASTRON), with financial support by the Netherlands Organisation of Scientific research (NWO).

Acknowledgments

We thank the anonymous referee for the numerous suggestions and the swift processing of this article. We thank Pavel Jáchym for numerous fruitful discussions during this work. This work was supported by the Czech Science Foundation grant 19-18647S and the institutional project RVO 67985815. This work was also supported by institutional research funding IUT26-2 and IUT40-2 of the Estonian Ministry of Education and Research, by the Centre of Excellence “Dark side of the Universe” (TK133) financed by the European Union through the European Regional Development Fund. This research made use of Astropy, (http://www.astropy.org) a community-developed core Python package for Astronomy (Astropy Collaboration 2013, 2018), and the Kapteyn Package for Python (Terlouw & Vogelaar 2015). This research made use of NASAs Astrophysics Data System Bibliographic Services.

References

  1. Arrigoni Battaia, F., Gavazzi, G., Fumagalli, M., et al. 2012, A&A, 543, A112 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Astropy Collaboration (Robitaille, T. P., et al.) 2013, A&A, 558, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Astropy Collaboration (Price-Whelan, A. M., et al.) 2018, AJ, 156, 123 [Google Scholar]
  4. Baldwin, J. A., Phillips, M. M., & Terlevich, R. 1981, PASP, 93, 5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Balogh, M. L., Navarro, J. F., & Morris, S. L. 2000, ApJ, 540, 113 [NASA ADS] [CrossRef] [Google Scholar]
  6. Balogh, M. L., Morris, S. L., Yee, H. K. C., Carlberg, R. G., & Ellingson, E. 1999, ApJ, 527, 54 [NASA ADS] [CrossRef] [Google Scholar]
  7. Barsanti, S., Girardi, M., Biviano, A., et al. 2016, A&A, 595, A73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Bekki, K., & Couch, W. J. 2003, ApJ, 596, L13 [NASA ADS] [CrossRef] [Google Scholar]
  9. Bekki, K., Couch, W. J., & Shioya, Y. 2002, ApJ, 577, 651 [NASA ADS] [CrossRef] [Google Scholar]
  10. Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Bianconi, M., Smith, G. P., Haines, C. P., et al. 2018, MNRAS, 473, L79 [NASA ADS] [CrossRef] [Google Scholar]
  12. Bianconi, M., Smith, G. P., Haines, C. P., et al. 2020, MNRAS, 492, 4599 [CrossRef] [Google Scholar]
  13. Binggeli, B. 1982, A&A, 107, 338 [NASA ADS] [Google Scholar]
  14. Binney, J., & Silk, J. 1979, MNRAS, 188, 273 [NASA ADS] [CrossRef] [Google Scholar]
  15. Boselli, A., & Gavazzi, G. 2006, PASP, 118, 517 [NASA ADS] [CrossRef] [Google Scholar]
  16. Boselli, A., & Gavazzi, G. 2014, A&ARv, 22, 74 [NASA ADS] [CrossRef] [Google Scholar]
  17. Bruzual, A. G. 1983, ApJ, 273, 105 [NASA ADS] [CrossRef] [Google Scholar]
  18. Bruzual, G., & Charlot, S. 2003, MNRAS, 344, 1000 [NASA ADS] [CrossRef] [Google Scholar]
  19. Butcher, H., & Oemler, A., Jr 1984, ApJ, 285, 426 [NASA ADS] [CrossRef] [Google Scholar]
  20. Cappellari, M., & Copin, Y. 2003, MNRAS, 342, 345 [NASA ADS] [CrossRef] [Google Scholar]
  21. Chambers, S. W., Melott, A. L., & Miller, C. J. 2002, ApJ, 565, 849 [NASA ADS] [CrossRef] [Google Scholar]
  22. Chung, A., van Gorkom, J. H., Kenney, J. D. P., & Vollmer, B. 2007, ApJ, 659, L115 [NASA ADS] [CrossRef] [Google Scholar]
  23. Cowie, L. L., & Songaila, A. 1977, Nature, 266, 501 [Google Scholar]
  24. Crone Odekon, M., Hallenbeck, G., Haynes, M. P., et al. 2018, ApJ, 852, 142 [NASA ADS] [CrossRef] [Google Scholar]
  25. Cybulski, R., Yun, M. S., Erickson, N., et al. 2016, MNRAS, 459, 3287 [NASA ADS] [CrossRef] [Google Scholar]
  26. da Cunha, E., Charlot, S., & Elbaz, D. 2008, MNRAS, 388, 1595 [NASA ADS] [CrossRef] [Google Scholar]
  27. Deshev, B. 2019, PhD thesis, Tartu Observatory (Estonia: University of Tartu) [Google Scholar]
  28. Deshev, B., Verheijen, M., van Gorkom, J., et al. 2009, Panoramic Radio Astronomy: Wide-field 1–2 GHz Research on Galaxy Evolution, 24 [Google Scholar]
  29. Deshev, B., Finoguenov, A., Verdugo, M., et al. 2017, A&A, 607, A131 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Donahue, M., Ettori, S., Rasia, E., et al. 2016, ApJ, 819, 36 [Google Scholar]
  31. Dressler, A., & Shectman, S. A. 1988, AJ, 95, 985 [NASA ADS] [CrossRef] [Google Scholar]
  32. Ebeling, H., & Kalita, B. S. 2019, ApJ, 882, 127 [CrossRef] [Google Scholar]
  33. Einasto, J. 1965, Trudy Astrofizicheskogo Instituta Alma-Ata, 5, 87 [NASA ADS] [Google Scholar]
  34. Einasto, M., Lietzen, H., Tempel, E., et al. 2014, A&A, 562, A87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Einasto, M., Deshev, B., Lietzen, H., et al. 2018a, A&A, 610, A82 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Einasto, M., Gramann, M., Park, C., et al. 2018b, A&A, 620, A149 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Faltenbacher, A., Allgood, B., Gottlöber, S., Yepes, G., & Hoffman, Y. 2005, MNRAS, 362, 1099 [NASA ADS] [CrossRef] [Google Scholar]
  38. Finoguenov, A., Böhringer, H., & Zhang, Y. Y. 2005, A&A, 442, 827 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Finoguenov, A., Henriksen, M. J., Miniati, F., Briel, U. G., & Jones, C. 2006, ApJ, 643, 790 [NASA ADS] [CrossRef] [Google Scholar]
  40. Fujita, Y. 2004, PASJ, 56, 29 [NASA ADS] [Google Scholar]
  41. Fujita, Y., & Nagashima, M. 1999, ApJ, 516, 619 [NASA ADS] [CrossRef] [Google Scholar]
  42. Fuller, T. M., West, M. J., & Bridges, T. J. 1999, ApJ, 519, 22 [NASA ADS] [CrossRef] [Google Scholar]
  43. Geller, M. J., Hwang, H. S., Fabricant, D. G., et al. 2014, ApJS, 213, 35 [NASA ADS] [CrossRef] [Google Scholar]
  44. Geller, M. J., Hwang, H. S., Dell’Antonio, I. P., et al. 2016, ApJS, 224, 11 [NASA ADS] [CrossRef] [Google Scholar]
  45. Gogate, et al. 2020, MNRAS, submitted [Google Scholar]
  46. Goto, T., Okamura, S., Sekiguchi, M., et al. 2003, PASJ, 55, 757 [NASA ADS] [CrossRef] [Google Scholar]
  47. Gunn, J. E., & Gott, J. R., III 1972, ApJ, 176, 1 [NASA ADS] [CrossRef] [Google Scholar]
  48. Haines, C. P., Pereira, M. J., Smith, G. P., et al. 2013, ApJ, 775, 126 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  49. Haines, C. P., Pereira, M. J., Smith, G. P., et al. 2015, ApJ, 806, 101 [NASA ADS] [CrossRef] [Google Scholar]
  50. Haines, C. P., Iovino, A., Krywult, J., et al. 2017, A&A, 605, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Haines, C. P., Finoguenov, A., Smith, G. P., et al. 2018, MNRAS, 477, 4931 [NASA ADS] [CrossRef] [Google Scholar]
  52. Hallenbeck, G., Koopmann, R., Giovanelli, R., et al. 2017, AJ, 154, 58 [CrossRef] [Google Scholar]
  53. Haynes, M. P., & Giovanelli, R. 1984, AJ, 89, 758 [NASA ADS] [CrossRef] [Google Scholar]
  54. Henning, P. A., Sancisi, R., & McNamara, B. R. 1993, A&A, 268, 536 [NASA ADS] [Google Scholar]
  55. Hwang, H. S., & Lee, M. G. 2009, MNRAS, 397, 2111 [NASA ADS] [CrossRef] [Google Scholar]
  56. Hwang, H. S., Geller, M. J., Diaferio, A., Rines, K. J., & Zahid, H. J. 2014, ApJ, 797, 106 [CrossRef] [Google Scholar]
  57. Jachym, P., Kenney, J. D. P., Sun, M., et al. 2019, ApJ, 883, 145 [NASA ADS] [CrossRef] [Google Scholar]
  58. Jaffé, Y. L., Aragón-Salamanca, A., Kuntschner, H., et al. 2011, MNRAS, 417, 1996 [NASA ADS] [CrossRef] [Google Scholar]
  59. Jaffé, Y. L., Poggianti, B. M., Verheijen, M. A. W., Deshev, B. Z., & van Gorkom, J. H. 2013, MNRAS, 431, 2111 [CrossRef] [Google Scholar]
  60. Jaffé, Y. L., Verheijen, M. A. W., Haines, C. P., et al. 2016, MNRAS, 461, 1202 [NASA ADS] [CrossRef] [Google Scholar]
  61. Jaffé, Y. L., Poggianti, B. M., Moretti, A., et al. 2018, MNRAS, 476, 4753 [NASA ADS] [CrossRef] [Google Scholar]
  62. Kasun, S. F., & Evrard, A. E. 2005, ApJ, 629, 781 [NASA ADS] [CrossRef] [Google Scholar]
  63. Kenney, J. D. P., van Gorkom, J. H., & Vollmer, B. 2004, AJ, 127, 3361 [NASA ADS] [CrossRef] [Google Scholar]
  64. Kennicutt, R. C., & Evans, N. J. 2012, ARA&A, 50, 531 [NASA ADS] [CrossRef] [Google Scholar]
  65. Kewley, L. J., Dopita, M. A., Sutherland, R. S., Heisler, C. A., & Trevena, J. 2001, ApJ, 556, 121 [Google Scholar]
  66. Koopmann, R. A., Giovanelli, R., Haynes, M. P., et al. 2008, ApJ, 682, L85 [NASA ADS] [CrossRef] [Google Scholar]
  67. Kotulla, R., Fritze, U., Weilbacher, P., & Anders, P. 2009, MNRAS, 396, 462 [NASA ADS] [CrossRef] [Google Scholar]
  68. Kregel, M., van der Kruit, P. C., & de Blok, W. J. G. 2004, MNRAS, 352, 768 [NASA ADS] [CrossRef] [Google Scholar]
  69. Kuutma, T., Tamm, A., & Tempel, E. 2017, A&A, 600, L6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  70. Larson, R. B., Tinsley, B. M., & Caldwell, C. N. 1980, ApJ, 237, 692 [NASA ADS] [CrossRef] [Google Scholar]
  71. Leisman, L., Haynes, M. P., Giovanelli, R., et al. 2016, MNRAS, 463, 1692 [NASA ADS] [CrossRef] [Google Scholar]
  72. Martínez, H. J., Muriel, H., & Coenda, V. 2016, MNRAS, 455, 127 [NASA ADS] [CrossRef] [Google Scholar]
  73. McNamara, B. R., Sancisi, R., Henning, P. A., & Junor, W. 1994, AJ, 108, 844 [NASA ADS] [CrossRef] [Google Scholar]
  74. Minchin, R. F., Taylor, R., Köppen, J., et al. 2019, AJ, 158, 121 [CrossRef] [Google Scholar]
  75. Montes, M., & Trujillo, I. 2019, MNRAS, 482, 2838 [Google Scholar]
  76. Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 [NASA ADS] [CrossRef] [Google Scholar]
  77. Okabe, T., Nishimichi, T., Oguri, M., et al. 2018, MNRAS, 478, 1141 [NASA ADS] [CrossRef] [Google Scholar]
  78. Okabe, T., Nishimichi, T., Oguri, M., et al. 2020, MNRAS, 491, 2268 [Google Scholar]
  79. Oosterloo, T., & van Gorkom, J. 2005, A&A, 437, L19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  80. Park, C., & Hwang, H. S. 2009, ApJ, 699, 1595 [NASA ADS] [CrossRef] [Google Scholar]
  81. Peng, C. Y., Ho, L. C., Impey, C. D., & Rix, H.-W. 2002, AJ, 124, 266 [NASA ADS] [CrossRef] [Google Scholar]
  82. Peng, Y., Maiolino, R., & Cochrane, R. 2015, Nature, 521, 192 [NASA ADS] [CrossRef] [Google Scholar]
  83. Plionis, M. 1994, ApJS, 95, 401 [NASA ADS] [CrossRef] [Google Scholar]
  84. Poggianti, B. M., & Barbaro, G. 1997, A&A, 325, 1025 [NASA ADS] [Google Scholar]
  85. Poggianti, B. M., Moretti, A., Gullieuszik, M., et al. 2017, ApJ, 844, 48 [NASA ADS] [CrossRef] [Google Scholar]
  86. Quilis, V., Moore, B., & Bower, R. 2000, Science, 288, 1617 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  87. Rupen, M. P. 1999, in Spectral Line Observing II: Calibration and Analysis, eds. G. B. Taylor, C. L. Carilli, & R. A. Perley, ASP Conf. Ser., 180, 229 [Google Scholar]
  88. Sancisi, R., Thonnard, N., & Ekers, R. D. 1987, ApJ, 315, L39 [NASA ADS] [CrossRef] [Google Scholar]
  89. Sersic, J. L. 1968, Atlas de Galaxias Australes (Cordoba, Argentina: Observatorio Astronomico) [Google Scholar]
  90. Smith, G. P., Khosroshahi, H. G., Dariush, A., et al. 2010, MNRAS, 409, 169 [NASA ADS] [CrossRef] [Google Scholar]
  91. Springel, V., White, S. D. M., Jenkins, A., et al. 2005, Nature, 435, 629 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  92. Stroe, A., Oosterloo, T., Röttgering, H. J. A., et al. 2015, MNRAS, 452, 2731 [NASA ADS] [CrossRef] [Google Scholar]
  93. Taylor, R., Davies, J. I., Jáchym, P., et al. 2017, MNRAS, 467, 3648 [Google Scholar]
  94. Taylor, R., Köppen, J., Jáchym, P., et al. 2020, AJ, 159, 218 [CrossRef] [Google Scholar]
  95. Tempel, E., Stoica, R. S., Martínez, V. J., et al. 2014, MNRAS, 438, 3465 [NASA ADS] [CrossRef] [Google Scholar]
  96. Terlouw, J. P., & Vogelaar, M. G. R. 2015, Kapteyn Package, version 2.3, (Groningen: Kapteyn Astronomical Institute). Available from http://www.astro.rug.nl/software/kapteyn [Google Scholar]
  97. Tonnesen, S., & Bryan, G. L. 2010, ApJ, 709, 1203 [NASA ADS] [CrossRef] [Google Scholar]
  98. Tonnesen, S., Bryan, G. L., & van Gorkom, J. H. 2007, ApJ, 671, 1434 [NASA ADS] [CrossRef] [Google Scholar]
  99. Verdugo, M., Ziegler, B. L., & Gerken, B. 2008, A&A, 486, 9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  100. Verheijen, M., van Gorkom, J. H., Szomoru, A., et al. 2007, ApJ, 668, L9 [Google Scholar]
  101. Verheijen, M., Deshev, B., van Gorkom, J., et al. 2010, ArXiv e-prints [arXiv:1009.0279] [Google Scholar]
  102. Vollmer, B., Cayatte, V., Balkowski, C., & Duschl, W. J. 2001, ApJ, 561, 708 [NASA ADS] [CrossRef] [Google Scholar]
  103. Vulcani, B., Poggianti, B. M., Gullieuszik, M., et al. 2018, ApJ, 866, L25 [NASA ADS] [CrossRef] [Google Scholar]
  104. Vulcani, B., Poggianti, B. M., Moretti, A., et al. 2019, MNRAS, 487, 2278 [CrossRef] [Google Scholar]
  105. Yoon, H., Chung, A., Smith, R., & Jaffé, Y. L. 2017, ApJ, 838, 81 [NASA ADS] [CrossRef] [Google Scholar]

All Figures

thumbnail Fig. 1.

Selection of the cluster members (blue squares) of A963 in projected phase space, redshift vs. cluster-centric distance. The cluster members from the ACReS survey have black squares around their symbols. The orange dots show all the galaxies with spectroscopy in the field. The redshift of the cluster (z = 0.205) is shown with the grey horizontal line.

In the text
thumbnail Fig. 2.

K-band luminosity density of the A963 cluster members. Solid squares show cluster members with log(M) ≥ 10.0 M. Small points show cluster members with lower stellar mass. The yellow contour levels are indicated by the tick marks on the colour bar. The green circle has a radius equal to R200 of the cluster. The red contour, at 1.8 K Mpc−2, is used to define the radial extent of the rest of the maps in the article.

In the text
thumbnail Fig. 3.

Map of the individual Voronoi bins. Each bin is coloured according to the total number of galaxies within it. All the bins are numbered, starting with zero. The green circle has a radius equal to R200 of the cluster.

In the text
thumbnail Fig. 4.

Panel a: stellar mass distribution of the cluster members with the completeness limit indicated by the red, dashed line. Panel b: completeness of our data set with respect to K-band luminosity. The red dashed line is at K = K+2.0- the magnitude limit of ACReS survey. Panel c: completeness distribution on the sky. The green circle has a radius equal to R200 of the cluster.

In the text
thumbnail Fig. 5.

Specific star formation rate of the A963 member galaxies as a function of their stellar mass. The black, dashed line at constant log(SFR) = −0.5 is used to split the galaxies into star forming, shown with stars, and passive, shown with circles. Symbols are coloured according to the strength of the 4000 Å drop in luminosity. The symbol in the upper right corner shows median uncertainties.

In the text
thumbnail Fig. 6.

Distribution of Dn4000 values of A963 member galaxies. Plotted is the offset from the mean relation between stellar mass and Dn4000 for all the passive galaxies in the ACReS survey.

In the text
thumbnail Fig. 7.

Distribution of A963 member galaxies and their properties on the plane of the sky. The passive galaxies are shown as red bordered circles, with the recently quenched ones (RQ) in green, the rest in white. Star forming galaxies are shown as blue bordered stars, with the anemic ones (RSF) in magenta and the others in cyan. The local fraction of star forming galaxies, fSF, mapped on the Voronoi bins is colour-coded according to the colour bar. The green circle in the centre has a radius equal to R200 of the cluster. The black contours show the extended X-ray emission. The orientation of the image is shown with the large grey circle, showing the azimuth α as well as the directions towards north and east.

In the text
thumbnail Fig. 8.

Fraction of star forming galaxies as a function of cluster-centric distance in simulated clusters and in A963. Thin lines of different colour represent galaxies quenched with a delay after their crossing of R200 of the parent cluster, according to the colour bar at the top of the figure. The data for A963 are shown in black (symbols and line). The fSF values calculated for individual galaxies are shown as black points, and the mean fSF values within each Voronoi bin are shown as empty squares. The dashed line shows a third-order polynomial fit to the Voronoi bins. For clarity, a small horizontal offset has been added to the black points.

In the text
thumbnail Fig. 9.

Left: model from Fig. 8, which minimises the difference with the binned data, with a delay of 2.17 Gyr, converted to 2D distribution and used to populate the Voronoi bins with star forming galaxies, according to their K-band luminosity weighted centroids (black points). Right: map of fSF, shown in Fig. 7, divided by the spherically symmetrical model shown on the left. Blue bins show an excess of star forming galaxies, red ones show a lack. The black contours show the K-band luminosity density. The lowest contour is at 40 K Mpc−2. Each following contour is at twice the level of the previous one. The green circle has a radius equal to R200 of the cluster.

In the text
thumbnail Fig. 10.

Left: map showing the mean fraction of star forming galaxies in Voronoi bins for 1000 jackknife runs on randomly selected 90% of the data set. Right: standard deviation of fSF in these 1000 runs. The black contours show the K-band luminosity density, as in Fig. 9. The green circle has a radius equal to R200 of the cluster.

In the text
thumbnail Fig. 11.

Local mean (left) and standard deviation (right) of the specific star formation rate in A963 cluster members. The black contours show the K-band luminosity density as in Fig. 9. The green circle has a radius equal to R200 of the cluster.

In the text
thumbnail Fig. 12.

Recently quenched galaxies as a local fraction of the passive galaxies (left) and red, star forming galaxies as a local fraction of the star forming galaxies (right). The black contours show the K-band luminosity density as in Fig. 9. The green circle has a radius equal to R200 of the cluster.

In the text
thumbnail Fig. 13.

Left: map of the local H I deficiency of star forming galaxies. Lighter colours indicate high deficiency. This is derived from all star forming galaxies. Overlaid are symbols indicating the location of galaxies that are individually detected at 21 cm. The sizes and colours of the symbols are proportional to the offset between the systemic velocities of the gaseous and stellar components. The black contours show the K-band luminosity density, as in Fig. 9. The green circle has a radius equal to R200 of the cluster. Middle: total H I maps of three selected sources overplotted on the R-band image. The contours are drawn at 20, 40, 80, 160, and 320 μJy beam−1. The cluster members with optical spectroscopy are shown in red, and their redshift is indicated. The non-member galaxies with optical spectroscopy are shown in purple. The central redshift of the H I signal is indicated in the top left corner. The beam size is shown in the lowest panel. Right: renzograms of the selected sources. Contours at a constant significance level of 3σ are drawn for every channel containing signal from the source. The noise is estimated from the surrounding area of the sources as a function of frequency. The difference between the central velocity of the channels and the central velocity of the entire H I signal is colour-coded. The proposed optical counterparts are indicated with red arrows.

In the text
thumbnail Fig. 14.

Entropy map of the hot gas of A963, based on XMM-Newton observations. The arrows show the approximate direction of inflow expected from the filaments identified by J16 from the large-scale galaxy distribution around the cluster. The × symbol in the centre marks the location of the BCG. The green circle has a radius equal to R200 of the cluster.

In the text
thumbnail Fig. 15.

Properties of the galaxy population in the vicinity of A963 and its infall region as a function of the azimuthal angle. Top panel: galaxy number density normalised to its maximum. The black curve shows the average K-band luminosity density in the infall region (between 1.5 R200 and the edge of the field at ≲3 R200). The lime green dashed line shows the mean galaxy number density between 10 and 20 R200. Middle panel: H I deficiency in the infall region. Bottom panel: specific star formation rate (solid purple line), fraction of RQ (dashed green line), and fraction of RSF (red dotted line), expressed as a fraction of their maximum. The grey shaded areas represent the location and relative angular span of the three filaments suggested by J16.

In the text

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

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

Initial download of the metrics may take a while.