A Lack of Evidence for Global Ram-pressure Induced Star Formation in the Merging Cluster Abell 3266

We have selected Abell 3266 to search for ram-pressure induced star formation as a global property of a merging cluster. Abell 3266 (z = 0.0594) is a high mass cluster that features a high velocity dispersion, an infalling subcluster near to the line of sight, and a strong shock front. These phenomena should all contribute to making Abell 3266 an optimum cluster to see the global effects of RPS induced star formation. Using archival X-ray observations and published optical data, we cross-correlate optical spectral properties ([OII, H$\beta$]), indicative of starburst and post starburst, respectively with ram-pressure, $\rho$v$^{2}$, calculated from the X-ray and optical data. We find that post-starburst galaxies, classified as E+A, occur at a higher frequency in this merging cluster than in the Coma cluster and at a comparable rate to intermediate redshift clusters. This is consistent with increased star formation due to the merger. However, both starburst and post-starburst galaxies are equally likely to be in a low or high ram pressure environment. From this result we infer that the duration of the starburst phase must be very brief so that: (1) at any time only a small fraction of the galaxies in a high ram pressure environment show this effect, and (2) most post-starburst galaxies are in an environment of low ram pressure due too their continued orbital motion in the cluster.


INTRODUCTION
An estimated 20% of all galaxies are in clusters of galaxies (Bahcall 1988). These large structures are typically ∼3 million light years or more in radius and are roughly circular in projection. A hot ( 10 7−8 K), X-ray emitting intracluster medium (ICM) permeates all galaxy clusters. The X-ray spectra from the plasma has a thermal bremsstrahlung continuum with emission lines from high ionization states of elements such as O, Fe, Si, and S. These elements are formed in stars and their presence in the ICM implies a connection to star formation, which takes place predominately in galaxies (Jiménez-Teja et al. 2019). The elements have an average global abundance of 0.3 Solar for rich clusters with kT ∼5 keV (de Plaa et al. 2007). However, cluster radial metallicity profiles show a central peak of Fe as well as other elements (Simionescu et al. 2009). While studies of sedimentation of metals to the center of the cluster could account, at least in part, for the excess metals (Medved et al. 2014), inclusion of the intracluster magnetic field gives a prohibitively long time for this process (Sarazin 1988). With the expected slow settling time via diffusive models and inhibition from a likely non-radial magnetic field, sedimentation would be inefficient. Therefore, much of the metal enhancement must have originated from local galaxy evolutionary processes that are more efficient in the cluster center. (Gunn & Gott 1972) originally suggested that as galaxies move through the ICM, there is a possibility that their interstellar medium (ISM) could be removed through RPS. The high velocity of the galaxy through the ICM and the local density of the ICM, which increases toward the center, results in a ram pressure that is higher in the inner 500 kpc than the outer parts of the cluster (Smith et al. 2010). The increased metals in a region of generally higher gas density (and therefore higher ram-pressure) is consistent with RPS being associated with the metal enrichment. Though there are competing galaxy evolutionary processes that favor the cluster center (Henriksen & Byrd 1996), images of gas removal in galaxies by RPS are evident in the jelly-fish galaxies.
Modeling ram-pressure stripping of the hot coronae of galaxies in clusters is consistent with the trend of smaller cluster galaxy coronae compared to the hot X-ray halos of field galaxies. This indicates that at least some gas removal occurs within cluster environments for recently accreted galaxies (Vijayaraghavan & Ricker 2017). In fact, evidence of ram-pressure stripping of M86 in the Virgo cluster was seen in the X-ray soon after it was predicted (Forman et al. 1979). Later observations mapped the structure of its 150 kpc tidal tail in the X-ray with high resolution (Randall et al. 2008) providing details of the stripping process. More recently, an ample number of observations showing that RPS occurs in galaxy clusters (Crowl & Kenney 2006); (Merluzzi et al. 2013); (Kenney et al. 2014) have become available. Simulations and observations generally agree that the outer disk (> 20 kpc) is stripped first and the gas forms a tail in the wake of the galaxy (Gullieuszik et al. 2017); (Roediger & Bruggen 2007); (Tonnesen & Bryan 2010); (Sun et al. 2007); (Fumagalli et al. 2014); (Cortese et al. 2006). Simulations are in general agreement that ram pressure increases the galactic star formation rate ;  and there are examples of increased star formation observed during RPS (Ramatsoku et al. 2019); (Owers et al. 2012); (Ebling et al. 2014); (Poggianti et al. 2015); (Koopmann & Kenney 2004). (Poggianti et al. 2015) report a significant enhancement in star formation due to RPS in a large sample of clusters using optical spectroscopy. However, there are counter examples of RPS where enhanced star formation is not seen; for example, in the Virgo cluster (Crowl & Kenney 2008).
A general feature of simulations that are characterized by enhanced star formation is that ram pressure causes a compression of the gas within the disk of the galaxy. This transfer of gas from the outer regions of the galaxy to the inner, increases the gas density, resulting in an increase in the star formation rate. However, other simulations do not see enhanced star formation in the disk (Tonnesen & Bryan 2012) or only a slight enhancement for galaxies with select physical parameters and orientation (Troncoso Iribarren 2016); (Bekki 2014); (Steinhauser et al. 2016). Simulations of RPS during a galaxy cluster merger ;  show increased star formation due to local enhancements in gas density from the cluster wide shocks that typically accompany mergers and the broadened velocity dispersion. This phase of enhanced star formation in the disk and gas halo stripping (Bekki & Couch 2003) lead to an overall reduction of the gas available to the galaxy, resulting in a subsequent decrease in the overall star formation rate. Thus, there is an optimum time to observe the increased star formation rate as it should only be enhanced for a relatively short time period after ram pressure stripping begins . The length of time depends on the strength of the ram pressure acting on the galaxy and the velocity of the galaxy, since this would affect the length of the RPS phase. As the galaxy encounters denser regions of the ICM, and hence stronger ram pressure, higher star formation rates should occur. This is directly observed in the Jellyfish galaxies found in merging clusters (Vulcani et al. 2018); (Vijayaraghavan & Ricker 2013); (McPartland et al. 2016); (Ebling et al. 2014). The increase in star formation should be accompanied by an observed enhancement in the equivalent width of emission or absorption lines. The stripping timescale in the center of the cluster is estimated at 200 -500 Myr (Steinhauser et al. 2016); (Gullieuszik et al. 2017); (Smith et al. 2010) so that if a galaxy cluster is observed after this time, any enhancement in these lines may no longer be observable. (Poggianti et al. 2015) has reported increased star formation via RPS in clusters of various morphology using a spectroscopic data set. In this paper, we will specifically address the connection with cluster merger. We have chosen a cluster with an ongoing merger to look for a global effect of increased star formation since there is some disagreement about the importance of this phase in the simulations. Abell 3266 is a well studied cluster with an ongoing merger (Finoguenov et al. 2006); (Henriksen & Tittley 2002); (Henriksen et al. 2000); (Dehghan et al. 2017). The cluster has a high velocity dispersion in the central region, 1367 km s −1 , while the outer region of the cluster exhibits a velocity dispersion similar to other rich clusters, 1000 km s −1 (Henriksen & Tittley 2002). These velocities are consistent with the trend in values derived by (Dehghan et al. 2017) though smaller as these authors use a less restrictive velocity criterion for cluster membership and therefore model the larger cluster environment. The merger process creates more high velocity galaxies in the cluster as the kinetic energy of the infalling subcluster is turned into internal cluster energy via dynamical friction thus inflating the velocity dispersion in the central region. The cluster velocity distribution also has a high velocity component identified with a subcluster, which contains 30 galaxies (Henriksen et al. 2000). Further evidence of a merger is an asymmetric gas distribution and regions of shocked gas (Henriksen & Tittley 2002). There should be higher ram pressure in these local regions of increased galaxy velocity and post-shock, higher density gas. Depending on the timescale of stripping, as discussed above, we expect to see a correlation between high local ram pressure and the equivalent width of [OII] or Hβ due to an increase in star formation. If no correlation is seen between ram pressure and these lines, it could be that this phase of star formation has already ceased. However, a cluster with an ongoing merger is the best case for catching the relatively short duration process of induced star formation.
Given the large number of groups in the cluster found optically that also appear to have enhanced X-ray emission (Dehghan et al. 2017), we calculate the local gas density for each galaxy. The local gas density is obtained from an annulus with inner radius of 25 kpc and outer radius of 50 kpc. These values are nominally chosen to both exclude emission from a disk similar to the Milky Way and restrict gas density to the local value encountered by the galaxy. The angular scale for Abell 3266 is 70.4 kpc arcmin −1 so the annulus dimensions are sub arc minute. The point-spread function for the PSPC is described by a Gaussian with a width that varies with distance from the detector center from 15 arc sec to 50 arc sec at 1 keV (Hasinger et al. 1992). Thus the regions are folded through the point spread function (see Table 2). The energy range for each region was restricted to 0.5-2.4 keV, which maximizes source counts while minimizing background contamination. Only galaxies with an X-ray signal-to-noise count rate greater than 3 were used to obtain a gas density. The count rate is converted to a flux using an APEC model with specified values for temperature, abundance, and column density. The fluxes and normalizations are shown in Table 3. The ion density is calculated from the normalization, which contains the emission integral, using the sampled volume and an average electron per ion corresponding to the assumed abundances. The ion density is given in Table 4. The abundances assumed for the count rate to flux conversion as well as the density calculation is 90% H and 10% He, with metals ignored. The emission weighted PSPC temperature, 5.8 keV (Henriksen et al. 2000) is used throughout the cluster as the PSPC is insensitive to temperature changes in the X-ray emitting gas. The column density is set to 2×10 20 cm −2 (Bekhti et al. 2016). The ram pressure depends on the velocity of each galaxy in addition to the ion density and mass. (Henriksen et al. 2000) found that the range of galaxies in the Abell 3266, following a clipping procedure described in detail there, includes galaxies with heliocentric velocity between 15,000 and 21,000 km s −1 . Figure 1 shows the velocity histogram for the galaxies used in our analysis.
Of the original 263 galaxies in the WINGS survey, 51 are outside of the velocity range and are likely foreground or background. Another 35 galaxies are eliminated due to low signal-to-noise X-ray data or non-physical [OII] values. The heliocentric velocity of the main cluster core found by (Henriksen et al. 2000) is 17,804 km s −1 giving a redshift of 0.0594, consistent with (Dehghan et al. 2017). We use 0.0594 as the redshift throughout the paper. The velocity of each galaxy is found by subtracting the galaxy velocity from the heliocentric cluster velocity. The galaxy velocities are given in Table 4. The distribution of galaxies, as seen in Figure 1, is asymmetric and the peak is offset from zero, further evidence of the merger.
In general, the gas density is higher toward the center, and lower in the outskirts. However, in a cluster merger, there is an asymmetry in the X-ray surface brightness, as can be seen in Figure 2. That is why a simple radial analysis, which assumes azimuthal symmetry is insufficient, justifying the analysis here that uses local density, to reflect gas distribution asymmetries from the merger.

RESULTS
The ram pressure for each galaxy, ρv 2 , was calculated using the gas densities and velocities in Table 4. Figure 3 shows the histogram of ram pressure for the 177 galaxies in this study. As can be seen by observing the histogram, the galaxies experience a large range of ram pressure. We adopted the condition for high ram pressure used by (Kapferer et al. 2009) of 5× 10 −11 dynes cm −2 . Of the 177 galaxies used, 61 galaxies were found to have ram pressures that can be considered "high" (i.e. ≥ 5 × 10 −11 dyn cm −2 . The remaining 116 galaxies were found to have "low" ram pressures. The calculated ram pressures for each of the 177 galaxies are listed in Table 4. We would expect to see galaxies with high ram pressures in regions where the density is high, or in galaxies with high velocities with respect to the cluster. In Figure 2, the distribution of these galaxies can be seen overlaid on the x-ray image of the cluster. Galaxies with low and high ram pressure appear throughout the cluster though high ram pressure galaxies tend to be closer to the center. Figure 2. Positions of galaxies overlaid on the xray image of Abell 3266 color coded to show ram pressure strength. Red depicts a galaxy with ram pressure > 5 × 10 −11 dyn cm −2 and blue a galaxy with ram pressure ≤ 5× 10 −11 dyn cm −2 . Figure 3. Histogram of the ram pressure for each galaxy. The first two bins contain all galaxies with low ram pressure Star formation in galaxies is enhanced by the ram pressure they experience. This is seen in both simulations and observations. To assess whether a galaxy is currently experiencing or recently experienced an enhancement in star formation, we inspected the equivalent width of Hβ and [OII] for each galaxy in the study.   It is believed that the surrounding environment can impact how the galaxy evolves (e.g. (Abadi et al. 1999); (Miller & Owen 2003); (Roediger & Bruggen 2007); (Fritz et al. 2011); (Fritz et al. 2014). If there is a correlation between ram pressure and star formation in Abell 3266, it should be apparent when comparing a histogram of equivalent widths for galaxies with high ram pressures verses galaxies with low ram pressures. Figures 6 and 7 show the histograms of the equivalent width of Hβ and [OII], respectively, plotted based on whether the galaxy has high or low ram pressure. Note that If the equivalent width of Hβ is positive (negative), this corresponds to absorption (emission). Perhaps there is a weak correlation that is buried in the large, inclusive samples. Figure 6. Histogram of Hβ equivalent width based on the ram pressure of the galaxy. The solid histogram corresponds to high ram pressure (≥ 5 x 10 −11 dyn cm −2 ). The black dashed histogram corresponds to low ram pressure (< 5 x 10 11 dyn cm −2 ). Each bin has a width of 0.5Å. The low and high ram pressure galaxies span the same range of absorption line strength though the high absorption line galaxies have a higher fraction of high ram pressure galaxies. Low absorption line galaxies are preferentially low ram pressure.   . Plot of equivalent width of Hβ verses ram pressure for sample A, the 34 galaxies with 3 sigma significance. There is no visible preference for galaxies with Hβ in emission to be in a region of high ram pressure.
To evaluate this hypothesis, we addressed the fact that the equivalent widths of Hβ and [OII] obtained from WINGS do not all have a signal to noise of 3σ or greater. For this reason we have created two subsamples which will be referred to as Sample A and Sample B. Sample A consists of all galaxies with line emission greater than 3σ significance. Sample A consists of a total of 34 galaxies for Hβ and 24 for O [II]. Figures 8 and 9 show a plot of O[II] verses ram pressure and Hβ verses ram pressure, respectively, for the Sample A galaxies. Sample B is a subgroup of Sample A, consisting of all galaxies with a signal to noise greater than 3σ for both Hβ and [OII]. Sample B consists of a total of 12 galaxies.   4.1. E+A Galaxies (Fisher et al. 1998) defines an E+A galaxy as one with Balmer absorptions lines that are (Hδ+Hγ+Hβ)/3 > 4Å, and [OII] emission lines are < 5Å. For our criteria, we come as close as we can to this with our database by requiring an equivalent width of either Hβ or Hδ > 4Å, and equivalent width of [OII] < -5Å. There are 6 galaxies out of the 71 (that have spectroscopic coverage) or 8.5% that meet these classification requirements. Of those 6, 3 have "high" ram pressures > 5 x 10 −11 , while the other 3 have "low" ram pressures. Therefore, the data do not support a correlation between the E+A phenomenon and ram pressure. However, the frequency of E+A galaxies, 8.5% is substantially higher than the fraction found in a non-merging cluster, Coma (Caldwell et al. 1993) and this indicates a galaxy evolutionary effect associated with the merger.

DISCUSSION
The relationship between ram pressure and star formation during a cluster merger was investigated for Abell 3266. This galaxy cluster was chosen because it shows evidence of a merger in progress: very high velocity dispersion in the center, an asymmetric density distribution, and regions of shock heated gas. Thus, it is the best opportunity to catch an increase in star formation due to RPS before the galaxies are quenched. Spectral type O and B type stars emit continuum radiation that ionizes Hydrogen in the surrounding medium. The Hydrogen ions then recombine with free electrons producing recombination emission lines. These emission lines are believed to correspond to HII regions surrounding O and B type stars. Other stellar classes, such as A and F type stars, do not produce strong enough radiation to ionize the surrounding gas and thus are lacking Hydrogen emission lines. Because of the short lifetime of the high mass stars, this emission line should track starburst galaxies. We found that the galaxies exhibiting Hβ in emission do not correlate with high ram pressure in our study. We do not see evidence of a widespread RPS induced starburst phase associated with the merger. There may be evidence of an increased level of recent star formation due to the merger, but we don't find a direct link to ram pressure. We then turn our attention to the possibility of post-starburst galaxies and their correlation with ram pressure.
An episode of rapid star formation would consume a large percentage of the gas available in the galaxy. This consumption, coupled with ram pressure stripping, would lead to a termination of star formation. While the brief timescale of enhanced star formation via RPS is critical to whether or not a correlation is observed it is possible to detect the post-starburst galaxies as E+A galaxies. An E+A galaxy has a characteristic spectrum, displaying strong Balmer absorption and little or no [OII] emission ( (Tran et al. 2003); (Fritz et al. 2014); (Caldwell et al. 1993). The star formation observed in an E+A galaxy is understood to have occurred within the last 1.5 Gyr (Couch & Sharples 1987(@). (Tran et al. 2003) observed that the number of E+A galaxies found in intermediate redshift clusters is greater than the number found in the field. This implies that environment is important in their formation. Such events as ram pressure stripping could play a role in the increased number of post starburst galaxies found in clusters verses the field. (Tran et al. 2003) and (Dressler et al. 1999) both observe that E+A's are generally located outside of the cluster core, which may be a result of their continued motion after RPS induced starburst.
The [OII] emission comes from the medium surrounding hot, young stars, namely O and B type stars. The optimum conditions for producing [OII] emission can be found in HII regions (Hogg et al. 1998). This makes the [OII] emission line a good indication of young stellar populations. If the [OII] emission line is very weak, or absent, it is likely that there is no current star formation occurring in the galaxy. Pairing this information with strong Balmer absorption lines allows us to classify galaxies as post starburst. If there was a recent starburst in a galaxy, there will still be some O and B type stars present which have not reached the end of their lifetime ( 1-10 Myr); however, there will be very little gas remaining in the former HII regions, resulting in a lack of [OII] emission. The A and F type stars, with lifetimes on the order of 1 Gyr or longer, will still provide strong Balmer absorption lines in the spectra of those galaxies. We find here that 8.5% or 6 out of 71 of the galaxies that have both OII and Hβ or Hδ available are E+A galaxies. (Tran et al. 2003) find E+A's make up a non-negligible component of the cluster population (7% -13%) at intermediate redshifts. These authors point out that compared to the low E+A fraction in Coma (3%) (Caldwell et al. 1993) their results show the E+A fraction evolves strongly with redshift. We find a comparable fraction, 8.5% in Abell 3266, a low redshift merging cluster to the intermediate redshift clusters. Within our limited sample we can say the E+A fraction is higher in this merging cluster compared to the nearby Coma cluster. The 6 E+A galaxies are split evenly between high and low ram pressure. While this E+A analysis does not show a correlation between post starburst galaxies and high ram pressure, the lack of immediacy in this phenomena means that it is likely the E+A galaxies are observed in a different gas environment than the one they were formed in.
The lack of a strong global correlation between ram-pressure and optical star formation indicators may also be complicated by the multiplicity of cluster environmental effects on galaxy evolution. For example, radio observations reveal an asymmetric, truncated distribution of HI in the spiral galaxy, UGC 07049, in the NGC 4065 group of galaxies (Freeland et al. 2010). In addition, radio continuum emission reveals strong star formation in the galaxy. The authors conclude that a combination of tidal and ram-pressure stripping are most likely necessary to produce the HI deficiency and disk asymmetry. In a study of three Virgo Cluster galaxies, in one of the galaxies, NGC 4264, X-ray observations at the position of magnetic field enhancements do not show compressional heating as expected from a strong adiabatic shock due to ram pressure. This favors tidal interactions rather than ram pressure stripping. Conversely, NGC 4569's radio polarized ridge shows a higher temperature, indicative of ram-pressure shocked gas. The third galaxy, NGC 2276, shows no clear indications of a shock, which is consistent with the observed distortions coming from tidal interaction (Wezgowiec et al. 2012).
Finally, though there is evidence of an ongoing merger in Abell 3266, it may have previously had a major merger at an earlier epoch that resulted in cluster wide quenching of star formation diminishing the effect of the present merger. Support for this idea comes from studies at intermediate redshift. (Vaughan et al. 2020) find that star formation takes place throughout the stellar disk in redshift, 0.3 -0.6 galaxies though cluster galaxies have a significantly smaller Hα disk than field galaxies. The authors concluded that ram-pressure stripping and disk strangulation can account for this result. The high level of E+A galaxies we find is consistent with a larger merger within the last ∼billion years.

CONCLUSIONS
Through this work the properties of ram pressure and Hβ equivalent width are examined for 177 galaxies located in the merging galaxy cluster Abell 3266. There are cases of ram pressure induced star formation within galaxies, the so called Jellyfish galaxies that appear to be associated with cluster merger. Although Abell 3266 is an ongoing merger the expected relationship between ram pressure and Hβ equivalent width was not observed. This shows a lack of evidence for ongoing global star formation due to ram pressure at this point of the merger. To test whether star formation had recently ended, we identified galaxies whose spectra were characteristic of E+A galaxies, or post starburst galaxies. There did not appear to be a correlation with E+A galaxies and ram pressure either. However, there is an enhancement in the number of E+A galaxies in this merging cluster compared to a quiescent cluster or the field, which is consistent with a major merger in the past billion years. The evidence contained in this work does not provide clear support that star formation is currently enhanced globally by ram pressure globally. But because simulations show that the enhanced star formation phase is short and depends sensitively on the environment conditions the galaxy encounters, we conclude that this phase is episodic, involving only a small fraction of the cluster galaxies at any epoch. The lack of correlation is also likely the result of the galaxy's continued motion after the RPS induced star formation phase.