## Measuring the ICM velocity structure within the A3266 galaxy cluster

E. Gatuzz 1 , J. Sanders 1 , A. Liu 1 , A. Fabian 2 , C. Pinto 3 , H. Russell , D. Eckert 4 5 , S. Walker 6 J. ZuHone 7 and R. Mohapatra 8

- 1 Max-Planck-Institut für extraterrestrische Physik, Gießenbachstraße 1, 85748 Garching, Germany e-mail: egatuzz@mpe.mpg.de
- 2 Institute of Astronomy, Madingley Road, Cambridge CB3 0HA, UK
- 3 INAF - IASF Palermo, Via U. La Malfa 153, I-90146 Palermo, Italy
- 4 School of Physics &amp; Astronomy, University of Nottingham, University Park, Nottingham NG7 2RD, UK
- 5 Department of Astronomy, University of Geneva, Ch. d'Ecogia 16, CH-1290 Versoix, Switzerland
- 6 Department of Physics and Astronomy, University of Alabama in Huntsville, Huntsville, AL 35899, USA
- 7 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA, 02138, USA
- 8 Department of Astrophysical Sciences, Princeton University, NJ 08544, USA

Received XXX; accepted YYY

## ABSTRACT

Wepresent a detailed analysis of the velocity structure of the hot intracluster medium (ICM) within the A3266 galaxy cluster, including new observations taken between June and November 2023. Firstly, morphological structures within the galaxy cluster were examined using a Gaussian Gradient Magnitude (GGM) and adaptively smoothed GGM filter applied to the EPIC-pn X-ray image. Then, we applied a novel XMM-Newton EPIC-pn energy scale calibration, which uses instrumental Cu K α as reference for the line emission, to measure line-of-sight velocities of the hot gas within the system. This approach enabled us to create two-dimensional projected maps for velocity, temperature, and metallicity, showing that the hot gas displays a redshifted systemic velocity relative to the cluster redshift across all fields of view. Further analysis of the velocity distribution through non-overlapping circular regions demonstrated consistent redshifted velocities extending up to 1125 kpc from the cluster core. Additionally, the velocity distribution was assessed along regions following surface brightness discontinuities, where we observed redshifted velocities in all regions, with the largest velocities reaching 768 ± 284 km s. Moreover, we computed the velocity Probability Density Function (PDF) from the velocity map. / We applied a normality test, finding that the PDF adheres to an unimodal normal distribution consistent with theoretical predictions. Lastly, we computed a velocity structure function (VSF) for this system using the measured line-of-sight velocities. These insights advance our understanding of the dynamic processes within the A3266 galaxy cluster and contribute to our broader knowledge of ICM behavior in merging galaxy clusters.

Key words. X-rays: galaxies: clusters - Galaxies: clusters: intracluster medium - Galaxies: clusters: individual: A3266 -

## 1. Introduction

The intracluster medium (ICM) is a hot, ionized gas that fills the space between galaxies in a galaxy cluster, emitting X-rays due to its high temperature of several million kelvins. Measuring the velocities of the intracluster medium (ICM) is crucial for several reasons. Turbulent motions provide additional pressure support, particularly at large radii, a ff ecting cluster mass estimates and calculations of hydrostatic equilibrium (Lau et al. 2009). Simulations predict a close connection between entropy, temperature, density, and velocity power spectra, which should be tested (Gaspari et al. 2014). Velocity measurements can help constrain AGN feedback models, as the distribution of energy from AGN feedback within the cluster depends on the balance between turbulence and sound waves or shocks. Additionally, these motions contribute to the transport of metals within the ICM through processes such as the uplift and sloshing caused by AGN (e.g., Simionescu et al. 2008). The microphysics of the ICM, including viscosity, can also be probed by measuring velocities (ZuHone et al. 2016, 2018). Furthermore, velocity measurements can directly detect the sloshing of gas in cold fronts, which can persist for several Gyr (Ascasibar &amp; Markevitch 2006).

The Hitomi observatory directly measured random and bulk motions in the ICM using the Fe-K emission lines, thanks to its high spectral resolution microcalorimeter SXS X-ray detector. In the core of the Perseus cluster, it measured a bulk flow gradient of 150 km s across 60 kpc of the cluster core and a line-of-sight / velocity dispersion of 164 ± 10 km s between radii of 30 / -60 kpc (Hitomi Collaboration et al. 2016). These results, obtained for a very limited spatial region, showed that the Perseus core is not strongly turbulent despite the apparent impact of the AGN and its jets on the surrounding ICM and the merger-driven sloshing motions. Such a low level of turbulence argues against it providing distributed heating, although sound waves or shocks could plausibly do this as well as heating from merging activity (Fabian et al. 2017). Unfortunately, no further measurements were made in other clusters or di ff erent regions of Perseus due to the loss of Hitomi , although the new XRISM mission is expected to extend the Fe K velocity measurements in more clusters.

Suzaku also obtained velocity measurements in several systems using the Fe-K line, thanks to its accurate calibration. Tamura et al. (2014) placed upper limits on relative velocities of 300 km s over scales of 400 kpc in Perseus. Ota &amp; Yoshida /

(2016) examined several clusters with Suzaku , though systematic errors from the calibration were likely around 300 km s, / and its point spread function (PSF) was large. Other methods of measuring velocities include using XMM-Newton RGS spectra to measure line widths (e.g., Pinto et al. 2015) and examining resonant scattering (e.g., Ahoranta et al. 2016), both limited to the cluster core. Indirect measurements of velocity structure in clusters include analyzing the power spectrum of density fluctuations and linking this via simulations to the velocity spectrum (e.g., Zhuravleva et al. 2014), or examining the magnitude of thermodynamic perturbations (Hofmann et al. 2016). However, these methods are model-dependent.

Sanders et al. (2020, hereafter SAN 20) present a novel + technique using background X-ray lines observed in the spectra of the XMM-Newton EPIC-pn detector to calibrate the absolute energy scale of the detector to better than 150 km s at Fe-K. / Using this technique, SAN + 20 mapped the bulk velocity distribution of the ICM over a large fraction of the Perseus and Coma clusters. For the Perseus cluster, they detected evidence of sloshing associated with a cold front. In the Coma cluster, SAN + 20 found that the velocity of the gas closely matches the optical velocities of the two central galaxies, NGC 4874 and NGC 4889. However, the velocity structure near the cluster center was not studied for either source due to the lack of o ff set observations covering the central region.

In recent years, four large XMM-Newton observation programs were approved to observe the Virgo, Centaurus, and Ophiuchus clusters and apply this method. Using this technique, we have obtained accurate velocity measurements with uncertainties down to ∆ v ∼ 100 km s. The analysis of the Virgo cluster / shows signatures of both AGN outflows and gas sloshing in the velocity field of the cluster (Gatuzz et al. 2022b). For the Centaurus cluster, we have found (a) an overall gradient in velocities at large radii, (b) a complex velocity profile near the cluster center, and (c) changes in velocities along the cold fronts (Gatuzz et al. 2022a). In the Ophiuchus cluster, we have identified a large redshifted-blueshifted interface located approximately 150 kpc east of the cluster core (Gatuzz et al. 2023c). This interface, which shows signs of a lack of metal mixing, is traced by discontinuities in the X-ray surface brightness. These large observational campaigns also allow us to measure the velocity structure functions (VSF) within the systems (Gatuzz et al. 2023a). VSFs constitute a useful diagnostic tool to study the turbulence in a medium, since they represent the variation of velocity with scale (Federrath et al. 2010, 2021; Seta et al. 2023). Furthermore, a detailed analysis of the metallicity distribution in the context of the velocity structure has been carried out for these sources (Gatuzz et al. 2023b,d,e).

The A3266 galaxy cluster is an excellent candidate for studying the ICM velocity distribution. It is one of the most brilliant clusters in the sky (Edge et al. 1990) and one of the brightest merging clusters with temperatures exceeding 7 keV. ASCA observations of the cluster provide evidence of a merger, detecting a temperature variation along the merger axis (Henriksen et al. 2000). Chandra observations identified a cooler filamentary region centered on the central cD galaxy and aligned along the merger axis (Henriksen &amp; Tittley 2002). Finoguenov et al. (2006) estimated the mass of this cooler, denser material to be 1 3 . × 10 13 M with high metallicity using ⊙ XMM-Newton observations, suggesting it results from a merger with a mass ratio of 1:10. Sauvageot et al. (2005) proposed a similar scenario, attributing the low-entropy material to a merger in a direction close to the plane of the sky.

Fig. 1: Top panel: RGB image of A3266. The bands used are 0.2-0.5 keV (R), 0.5-1.0 keV (G) and 1.0-2.0 keV (B). We applied a σ = 1 Gaussian smoothing to the exposure corrected and background subtracted image. Bottom panel: X-ray surface brightness exposure corrected and background subtracted of the A3266 cluster in the 0.5 to 9.25 keV energy range. The point-like sources were excluded in the analysis. The 10-level contours of the X-ray image are shown in white.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00837v1-figure-802b3e50-6eab-41e8-bd0f-5ca429ee4542.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00837v1-figure-1c3c2bfc-f90e-45a1-8a9d-40ad467a2217.jpg)

Numerical simulations computed by Sauvageot et al. (2005) present two possible scenarios for the merger: the subcluster may have entered from the southwest, passing the core 0.8 Gyr ago and now nearing turnaround, or it could have merged from the northeast and is now exiting to the southwest after passing the core 0.15-0.20 Gyr ago. Ettori et al. (2019) analyzed XMMNewton data out to the cluster's virial radius, finding masses of M 500 = 8 8 . × 10 14 M and ⊙ M 200 = 15 × 10 14 M , with radii ⊙ of R 500 = 1 43 Mpc and . R 200 = 2 33 Mpc, assuming an NFW . mass model (Navarro et al. 1996). Dehghan et al. (2017) measured spectroscopic redshifts of 1300 sources within the cluster,

Table 1: XMM-Newton observations of the A3266 cluster.

|      ObsID | RA          | DEC         | Date Start-time   |   Exposure (ks) |
|------------|-------------|-------------|-------------------|-----------------|
| 0105262001 | 04:31:37.00 | -61:34:02.0 | 2000-09-23        |             9.3 |
| 0105261101 | 04:32:35.49 | -61:29:01.0 | 2000-09-25        |            15   |
| 0105262101 | 04:32:35.49 | -61:29:01.0 | 2000-09-25        |             9.7 |
| 0105261001 | 04:30:13.30 | -61:21:02.0 | 2000-09-27        |            15.4 |
| 0105262201 | 04:30:13.30 | -61:21:02.0 | 2000-09-27        |             9.2 |
| 0105260701 | 04:32:26.70 | -61:18:46.0 | 2000-10-01        |            22.1 |
| 0105262301 | 04:31:15.80 | -61:15:33.0 | 2000-10-03        |            34.1 |
| 0105260901 | 04:30:23.60 | -61:30:32.0 | 2000-10-09        |            26   |
| 0105260801 | 04:31:15.80 | -61:15:33.0 | 2000-10-11        |            21.8 |
| 0105262501 | 04:31:37.00 | -61:34:02.0 | 2003-03-15        |            13.6 |
| 0921880301 | 04:30:02.04 | -61:22:25.5 | 2023-06-13        |           132.8 |
| 0921880401 | 04:30:52.45 | -61:33:54.6 | 2023-06-15        |           132.9 |
| 0921880501 | 04:32:01.22 | -61:17:33.4 | 2023-06-25        |            78   |
| 0921880601 | 04:32:32.01 | -61:30:22.5 | 2023-06-27        |            80   |
| 0921880801 | 04:30:52.45 | -61:33:54.6 | 2023-06-29        |            76.9 |
| 0921880101 | 04:32:01.22 | -61:17:33.4 | 2023-10-17        |           117.9 |
| 0921880201 | 04:32:32.01 | -61:30:22.5 | 2023-11-22        |           122.6 |
| 0921880701 | 04:30:02.04 | -61:22:25.5 | 2023-11-24        |            85.9 |

identifying that the cluster core, which can be split into two components, decomposes into six groups and filaments to the north. They concluded that a range of continuous dynamical interactions are taking place, rather than a simple NE-SW merger.

A3266 was the target of a calibration observation for the new eROSITA X-ray telescope (Predehl et al. 2021) during its calibration and performance verification (CalPV) phase (Dennerl et al. 2020). In their data analysis, Sanders et al. (2022a) confirmed that the source is not a simple merging cluster but rather involves three di ff erent systems merging with the main body: the NE structure, the NW structure, and the W structure. These systems are seen as low-entropy material associated with highermetallicity gas. However, the data did not conclusively reveal whether there is a merger shock on the west side of the central core. Pressure jumps in the southeast and west directions were identified towards the outskirts of the cluster. The presence of several subclusters that are merging or will soon merge makes A3266 an ideal laboratory for detailed study of the ICM velocity structure within such a system.

We present an analysis of the ICM velocity structure within the A3266 galaxy cluster using XMM-Newton observations. The outline of this paper is as follows. We describe the data reduction process in Section 2 and analyze the X-ray images in Section 3. The fitting procedure and the results are shown in Section 4 while a discussion of the results is included in Section 5. Finally, Section 6 presents the conclusions and summary. Throughout this paper we assume the distance of A3266 to be z = 0 0594 . (Dehghan et al. 2017) and a concordance Λ CDM cosmology with Ω m = 0 3, . ΩΛ = 0 7, and . H 0 = 70 km s -1 Mpc -1 .

## 2. A3266 XMM-Newton observations

## 2.1. EPIC-pn spectra

We analyzed eight A3266 observations taken between June and November 2023 (PI: E. Gatuzz, ID:092188). Additionally, we included ten observations that are available in the archive. Table 1 lists information on the observations, including IDs, coordinates, dates, and clean exposure times. All XMM-Newton European Photon Imaging Camera (EPIC-pn, Strüder et al. 2001) spectra were reduced with the Science Analysis System (SAS 1 , version 19.1.0). We processed the observations with the epchain SAS tool, including only single-pixel events (PATTERN == 0) and filtering the data with FLAG == 0 to avoid bad pixels and regions

1 https://www.cosmos.esa.int/web/xmm-newton/sas

close to CCD edges. A 1.0 cts s rate threshold was applied to / filter out bad time intervals from flares.

To measure the velocity structure within the A3266 cluster, we applied the technique described in Sanders et al. (2020) and Gatuzz et al. (2022b) to calibrate the absolute energy scale of the EPIC-pn detector using the instrumental background X-ray lines identified in the detector spectra. This calibration method includes corrections for (1) the average gain of the detector during the observation, (2) the spatial gain variation across the detector over time, and (3) the energy scale as a function of detector position and time. We created X-ray images and identified point sources using the SAS task edetect\_chain , with a likelihood parameter det\_ml &gt; 10. In the following analysis we have excluded the point sources.

## 3. X-ray images

The top panel of Figure 1 shows an RGB image of the cluster, with energy bands corresponding to 0.2-0.5 keV (R), 0.5-1.0 keV (G), and 1.0-2.0 keV (B). The exposure-corrected images were smoothed with a Gaussian of σ = 1. The bottom panel of Figure 1 displays the X-ray exposure-corrected image in the 0.5-9.25 keV energy band. A Gaussian smoothing of σ = 4 pixels was applied, and the 10-level contour of the X-ray image is shown in white.

The top panel of Figure 2 shows the total exposure time in the 4.0-9.25 keV energy range. The middle panel displays the normalized number of counts in each 1.59 arcsec pixel within the Fe-K complex (6.50-6.90 keV), after subtracting the neighboring continuum images (6.06-6.43 keV and 6.94-7.16 keV, rest frame). The instrumental Cu ring is included to indicate in which regions the redshift can be measured. This image demonstrates a good distribution of counts in the cluster center region, allowing for detailed analysis of the Fe-K complex. The bottom panel shows the normalized number of counts in each 1.59 arcsec pixel around the Cu background line. The plot reveals the Cu hole described in Sanders et al. (2020) for each exposure. Since velocity measurements can only be obtained using the Cu background line, it is crucial to have a large number of counts in regions where the Fe-K complex redshift will be measured.

To study the structures within the galaxy cluster in detail, we filtered the data using a Gaussian Gradient Magnitude (GGM) filter (GGM, Sanders et al. 2016a,b). Figure 3 shows the gradient magnitude of the X-ray image convolved with Gaussians of scales 4 pixels (8 arcsec, top panel) and 8 pixels (16 arcsec, bottom panel). Point sources were cosmetically filled using the values of neighboring pixels before filtering. A 700-800 kpc elliptical surface brightness edge surrounding the core is clearly shown. The image does not accurately depict the relative magnitudes of features but highlights these structures.

We also applied the adaptively smoothed GGM filter developed in Sanders et al. (2022a). To apply this filter, the image and exposure map were convolved by a Gaussian with a σ corresponding to a radius containing 1024 counts. The log 10 value of this map was then taken, and the gradient was calculated by taking the di ff erence between neighboring pixels along the two axes and adding them in quadrature. Figure 4 shows the resulting image, clearly revealing the sharp edge near the cluster core.

Figure 5 shows the fractional di ff erence in the 0.5 to 2 keV surface brightness from the average at each radius. The black regions correspond to point-like sources whose sizes have been artificially increased due to image smoothing. The X-ray image was smoothed by a 2 σ Gaussian. Multiple clear sharp discontinuities are visible, which have also been identified by Sanders

et al. (2022a). We analyze these structures in detail in subsection 4.3.

## 4. Spectral fitting procedure

The spectra were fitted in the 4.0-9.25 keV energy band. For each spatial region analyzed, we combined the spectra from di ff erent observations, using the number of counts in the fitting spectral region as the weighting factor. We modeled the ICM X-ray emission with the apec thermal emission model version 3.0.9 (Foster et al. 2019), which includes corrections derived from the Hitomi observations (Hitomi Collaboration et al. 2018). A tbabs component (Wilms et al. 2000) was also included to account for the Galactic X-ray absorption, corresponding to 2 26 . × 10 20 cm -2 (Kalberla et al. 2005). However, we note that for the energy range analyzed, the absorbing component has no e ff ect on the modeling (Gatuzz et al. 2022b,a, 2023c). The free parameters in the model are the temperature, metallicity (i.e., metal abundances with He fixed at cosmic levels), redshift, and normalization.

For the X-ray spectral fits, we used the xspec data fitting package (version 12.14.0 ). We assumed 2 cash statistics (Cash 1979). Errors are quoted at the 1σ confidence level unless otherwise stated. Abundances are given relative to Lodders et al. (2009). As background components, we included CuK α , NiK α , ZnK α , and CuK β instrumental emission lines, as well as a power-law component with its photon index fixed at 0.136 (the average value obtained from the archival observations analyzed in Sanders et al. 2020).

## 4.1. Spectral maps

We created maps of the cluster spectral properties using a method developed in Sanders et al. (2020), based on fitting spectra from dynamically sized ellipses. To create these maps, we used elliptical regions with a 2:1 axis ratio, rotated so that the longest axis lies tangentially to the vector pointing to the central core. The radii of the ellipses adaptively change in a grid with a spacing of 0.25 arcmin to achieve approximately 1000 counts in the Fe-K complex after continuum subtraction, thereby reducing the velocity uncertainties. The semi-major axis range from 6 1 arcmin to 2 0 arcmin. We combined the spectra obtained per . . region for all observations and created weighted-average ancillary response files (ARFs) while using the same response matrix for all the pixels. The spectra were then fitted as described in Section 4.

Figure 6 shows the velocity map (left panel) and the velocity uncertainties (right panel) obtained. We use a blue-red color scale for the velocity map to indicate whether the gas is moving towards the observer (blue) or away (red) to improve clarity. We found a mean velocity of ∼ 312 km s for the entire map. For / illustrative purposes, Figure 7 shows the velocity map relative to the A3266 cluster, after subtracting the mean value 312 km s / obtained from the spectral map. The mean optical velocity for galaxies from Dehghan et al. (2017) for all regions within the entire map is ∼ 1727 5 km s, larger than the ICM overall ve-. / locity. Figure 8 shows the spectral maps for temperature (top panel), metallicity (middle panel), and emission measure (bottom panel). We found that the metallicity is not uniform within the cluster, decreasing towards the northwest, corresponding to the low fractional di ff erence shown in Figure 5. The metallicities we found are lower than those measured by Sanders et al.

2 https://heasarc.gsfc.nasa.gov/xanadu/xspec/

Table 2: A3266 cluster best-fit parameters for annular regions.

| Region   | apec model               | apec model      | apec model                | apec model                | apec model   |
|----------|--------------------------|-----------------|---------------------------|---------------------------|--------------|
|          | kT                       | Z               | z ( × 10 - 3 )            | norm ( × 10 - 3 )         | cstat / dof  |
| 1        | 7 . 19 + 0 . 64 - 0 45   | 0 . 19 ± 0 . 03 | 12 . 69 + 0 . 41 - 0 45   | 59 . 94 + 2 . 36 - 2 25   | 1470 / 1199  |
| 2        | . 7 . 00 + 0 . 48 - 0 26 | 0 . 21 ± 0 . 02 | . 17 . 74 + 0 . 37 - 0 49 | . 61 . 01 ± 1 . 63        | 1498 / 1199  |
| 3        | . 7 . 46 ± 0 . 44        | 0 . 22 ± 0 . 02 | . 19 . 30 ± 0 . 49        | 61 . 65 ± 1 . 57          | 1494 / 1199  |
| 4        | 7 . 45 ± 0 . 46          | 0 . 21 ± 0 . 02 | 19 . 80 ± 0 . 52          | 60 . 22 ± 1 . 48          | 1484 / 1199  |
| 5        | 7 . 34 ± 0 . 49          | 0 . 18 ± 0 . 02 | 19 . 60 ± 0 . 55          | 62 . 25 ± 1 . 90          | 1529 / 1199  |
| 6        | 7 . 34 ± 0 . 51          | 0 . 20 ± 0 . 02 | 18 . 38 ± 0 . 54          | 59 . 98 + 1 . 88 - 1 81   | 1490 / 1199  |
| 7        | 7 . 40 ± 0 . 51          | 0 . 19 ± 0 . 02 | 19 . 07 + 0 . 57 - 0 54   | . 58 . 43 + 1 . 93 - 1 84 | 1499 / 1199  |
| 8        | 7 . 65 ± 0 . 33          | 0 . 21 ± 0 . 02 | . 34 . 37 ± 0 . 63        | . 61 . 13 ± 1 . 27        | 1405 / 1199  |

(2022a). However, it is important to note that our analysis only includes the high-energy band ( &gt; 4 keV), so the absolute abundances are not entirely reliable. The temperature map tends to be more uniform, but as indicated before, the actual temperature distribution cannot be accurately determined since we are modeling only the high-energy band.

## 4.2. Velocity radial profiles

We proceeded to analyze annular regions to study the radial velocity structure. The thickness of these rings increases as the square root of the distance from the center of A3266. Figure 9 shows the exact spatial regions analyzed. Table 2 presents the best-fit results obtained for each region, numbered from the center outward, while Figure 10 displays the resulting velocity structure. We achieved velocity measurements with an accuracy of up to ∆ v ∼ 190 km s (for ring 4). The largest redshift blueshift rel-/ / ative to A3266 is (805 ± 284) km s for ring 5. As observed in / the spectral map (see subsection 4.1), we obtained a redshifted velocity distribution for almost all regions. For distances less than 240 kpc, the velocities tend to increase with radius. We found a mean velocity of ∼ 283 km s for all the regions. We / computed a χ 2 test based on the null hypothesis that the velocities are consistent with a constant velocity model and obtained a p -value = 0 33, . thus the observed variations in velocity are not due to random chance. Figure 11 shows the temperature (top panel) and metallicity (bottom panel) distribution. We found a uniform distribution for both parameters, similar to previous profiles computed by Ghirardini et al. (2019); Sanders et al. (2022b). However, since our spectral analysis does not include the soft energy band, these values could be underestimated. Furthermore, the apparent correlation between the emission measure and metallicity is most likely an artifact as the lower the metallicity the larger should be the emission measure to obtain lines of comparable fluxes.

## 4.3. Discontinuities in the surface-brightness

We manually selected multiple regions based on the fractional di ff erence in 0.5-2 keV X-ray surface brightness to study the velocity structure. The top panel of Figure 12 shows the regions used for the analysis, which are similar to those identified in Sanders et al. (2022a). The regions can be identified as follows: the central core (region 1), cavities (regions 2, 3, 4), the W structure (region 5), the NE structures (regions 6, 7), the bridge (region 8), and the NW structure (region 9). Table 3 shows the bestfit results. As shown in the spectral map (see subsection 4.1), we found redshifted velocities for all regions, with the largest velocities corresponding to cavity 2 (626 ± 284 km s) and the NE struc-/ ture (region 6, 596 ± 274 km s). The velocities for the remaining /

Table 3: A3266 cluster best-fit parameters for manually selected regions following the A3266 merging structure.

| Region   | apec model               | apec model      | apec model                  | apec model                  | apec model   |
|----------|--------------------------|-----------------|-----------------------------|-----------------------------|--------------|
|          | kT                       | Z               | z ( × 10 - 3 )              | norm ( × 10 - 3 )           | cstat / dof  |
| 1        | 7 . 26 + 0 . 72 - 0 51   | 0 . 17 ± 0 . 03 | 19 . 03 + 0 . 68 - 0 75     | 59 . 09 + 2 . 81 - 2 75     | 1455 / 1199  |
| 2        | . 7 . 26 + 0 . 53 - 0 45 | 0 . 22 ± 0 . 03 | . 18 . 87 + 0 . 57 - 0 56   | . 61 . 62 + 1 . 78 - 1 77   | 1470 / 1199  |
| 3        | . 7 . 58 ± 0 . 67        | 0 . 17 ± 0 . 03 | . 17 . 84 + 0 . 68 - 0 65   | . 59 . 96 + 2 . 72 - 2 62   | 1478 / 1199  |
| 4        | 7 . 36 + 0 . 63 - 0 56   | 0 . 18 ± 0 . 03 | . 15 . 62 + 0 . 55 - 0 53   | . 59 . 97 + 2 . 16 - 2 07   | 1494 / 1199  |
| 5        | . 7 . 75 ± 0 . 48        | 0 . 19 ± 0 . 02 | . 24 . 54 + 0 . 68 - 0 65   | . 60 . 27 + 1 . 93 - 1 91   | 1489 / 1199  |
| 6        | 7 . 04 + 0 . 57 - 0 35   | 0 . 19 ± 0 . 03 | . 15 . 69 + 0 . 42 - 0 51   | . 61 . 51 + 2 . 10 - 2 07   | 1461 / 1199  |
| 7        | . 7 . 04 + 0 . 49 - 0 32 | 0 . 21 ± 0 . 02 | . 18 . 16 + 0 . 43 - 0 50   | . 60 . 51 + 1 . 66 - 1 71   | 1536 / 1199  |
| 8        | . 7 . 53 ± 0 . 70        | 0 . 18 ± 0 . 03 | . 17 . 78 + 0 . 72 - 0 68   | . 60 . 39 + 2 . 77 - 2 72   | 1466 / 1199  |
| 9        | 6 . 89 + 0 . 77 - 0 . 25 | 0 . 17 ± 0 . 03 | . 15 . 54 + 0 . 38 - 0 . 65 | . 59 . 56 + 2 . 86 - 2 . 77 | 1473 / 1199  |

structures were similar, considering the uncertainties. We have also included the mean optical redshift for galaxies from Dehghan et al. (2017) within each one of the regions obtained. The mean galactic velocities are within the ICM velocities except for regions 1, 4, 7, and 8. The mean ICM velocity for all regions is ∼ 275 km s, and the mean galactic velocities are / ∼ 525 km s. / However, the regions following the surface brightness discontinuities do not correspond to the sub-structures identified in Dehghan et al. (2017). Figure 13 shows the temperature (top panel) and the metallicity (bottom panel) for these regions. Both parameters tend to be uniform across the di ff erent regions, with a hint of increasing metallicity in region 2.

## 4.4. Probability distribution function

Figure 14 shows the velocity probability distribution function (PDF) computed from the velocity map, weighted by area. We computed normality tests, including the Shapiro-Wilk (Shapiro &amp; Wilk 1965, SHA + 65) and D'Agostino and Pearson's (D'Agostino &amp; Pearson 1973, DAG + 73), to determine if a Gaussian well models the data. We found that the PDF follows a normal distribution as the p -value is larger than the α = 0 05 level . for both tests (0.63 for SHA + 65 and 0.71 for DAG + 73). The overall redshift obtained from the Gaussian fit is 524 ± 21 km s / with velocity dispersion of σ v = 68 ± 12 km s. /

## 4.5. Velocity structure function

Following Gatuzz et al. (2023a), we computed the first-order velocity structure function by taking the weighted average of the di ff erence between the line-of-sight velocities ( v ) of two points separated by r . Mathematically, we define it as:

$$\delta v ( r ) & = \frac { \sum _ { x } w ( x + \mathbf e _ { 1 } r, x ) | v ( x + \mathbf e _ { 1 } r ) - v ( x ) | } { \sum _ { x } w ( x + \mathbf e _ { 1 } r, x ) } \quad \begin{pmatrix} \text{ent regio} \\ \text{lected ex} \\ \text{measure} \\ \text{$\cdot$} \end{pmatrix}$$

where x denotes the position of any point in the dataset and e1 denotes a unit-vector in any direction. We bin δ v into logarithmically-spaced bins of separation r . Figure 15 shows the VSF weighted by the velocity uncertainties. The uncertainties are large, as the velocity di ff erences are smaller than the noise level.

## 5. Discussion

We have found that the overall gas moves away from the observer throughout the entire system (Figure 6), relative to the cluster redshift, with velocities increasing from ∼ 300 -800

km s. After subtracting the average redshift, the velocity spec-/ tral map does show an opposite blueshifted-redshifted structure around the cluster center (Figure 7), with the blueshifted region displaying a lower metallicity. However, there is no clear spiral pattern associated with gas sloshing (see, for example, Figure 18 in Sanders et al. 2020). Dehghan et al. (2017) found that the A3266 core has a velocity dispersion of ∼ 1460 km s and that / the dispersion of the whole cluster is ∼ 1340 km s in their analy-/ sis of spectroscopic redshifts of sources within the cluster. They concluded that there is a range of continuous dynamical interactions taking place within this system.

We did not find hints for a multimodal PDF, compared, for example, with the Ophiuchus cluster (Gatuzz et al. 2023a). Simulations predict a σ v value in the ∼ 200 -400 km s range / with indications of isotropy for the hot ICM phase (Ayromlou et al. 2023). For the Perseus cluster, Hitomi measured up to σ v ∼ 220 km s, while / XMM-Newton measurements by (Gatuzz et al. 2023a) show a velocity dispersion of σ v = 344 ± 18 km s / (Virgo), σ v = 371 ± 12 km s (Centaurus) and / σ v = 847 ± 92 km s / (Ophiuchus). Whereas the Perseus cluster velocity dispersion measured by Hitomi was due to line broadening, the XMMNewton measurements provide dispersion velocities from the range of shifts. The VSF results, on the other hand, can be interpreted as upper limits &lt; 600 km s up to / ∼ 3000 kpc scales, in good agreement with the VSFs computed for the Virgo, Centaurus and Ophiuchus clusters (see Figure 3 in Gatuzz et al. 2023a).

Dehghan et al. (2017) measured an overall optical redshift of the entire cluster and related structures to be ∼ 200 ± 50 km s / with respect to the cluster core with velocity dispersion of 1337 ± 67 km s, such high dispersion does not allow a direct / comparison with the hot ICM velocity distribution. However, we note that they cover more considerable distances than our analysis. They indicated that such high-velocity dispersion is not due to virialized motions but rather due to relative motions of the cluster components. For example, the western core component substructure has a redshift of ∼ 600 km s with respect to the / cluster core. Other redshifted sub-structures in their analysis include the subclusters 2 and 3, with velocities of ∼ 559 km s and / ∼ 393 km s, respectively (Dehghan et al. 2017, see their Table / 3). Large velocities can be observed due to subgroups, strong merger shocks, and or substructures (Liu et al. 2015, 2016; Ota / &amp;Yoshida 2016; Liu et al. 2018).

The most exciting feature observed in our analysis is the mean velocity of the hot gas, redshifted with respect to the A3266 cluster, measured at ∼ 300 km s (see Sections 4.1, 4.2 / and 4.3). Moreover, the mean optical velocity values tend to be &gt; 2 times larger than the mean ICM velocities for the di ff erent regions analyzed (i.e., the spectral map and the manually selected extraction regions). As indicated by the optical redshifts measured by Dehghan et al. (2017), the o ff set between the ICM and the cluster is due mainly to substructures in the cluster associated with the merging.

## 6. Conclusions and summary

We have analyzed the velocity structure of the hot ICM within the A3266 galaxy cluster. We measured line-of-sight velocities of the gas using the technique developed by Sanders et al. (2020); Gatuzz et al. (2022b,a), which allows a better calibration of the absolute energy scale of the XMM-Newton EPIC-pn detector. Our main findings and conclusions are:

- 1. We studied the morphological structures within the galaxy cluster by applying a GGM and adaptatively smoothed GGM filter to the EPIC-pn X-ray image.
- 2. We made 2D projected maps for the velocity, temperature, and metallicity. We have found that the hot gas within the system displays a redshifted velocity with respect to the cluster redshift along all the fields of view. After subtracting the average redshift, the velocity spectral map does show an opposite blueshifted-redshifted structure, with the blueshifted region displaying a lower metallicity.
- 3. We studied the velocity distribution by creating nonoverlapping circular regions. Our results are consistent with redshifted velocities as far as 1125 kpc from the cluster core.
- 4. We analyzed the velocity distribution along regions following the discontinuities in the surface brightness. We found redshifted velocities for all regions with the largest velocities to be 805 ± 284 km s. /
- 5. We compute the velocity PDF from the velocity map. By applying a normality test, we found that the PDF follows an unimodal normal distribution, as predicted by simulations.
- 6. We obtained upper limits for the VSF of &lt; 600 km s up to / ∼ 3000 kpc scales.

Finally, the present work will be followed by a detailed study of the temperature and abundance distribution by modeling the soft X-ray spectra of the A3266 galaxy cluster.

Acknowledgements. This work was supported by the Deutsche Zentrum für Luftund Raumfahrt (DLR) under the Verbundforschung programme (Messung von Schwapp-, Verschmelzungs- und Rückkopplungsgeschwindigkeiten in Galaxienhaufen.). This work is based on observations obtained with XMMNewton, an ESA science mission with instruments and contributions directly funded by ESA Member States and NASA. This research was carried out on the High Performance Computing resources of the raven cluster at the Max Planck Computing and Data Facility (MPCDF) in Garching operated by the Max Planck Society (MPG).

## References

Ahoranta, J., Finoguenov, A., Pinto, C., et al. 2016, A&amp;A, 592, A145 Ascasibar, Y. &amp; Markevitch, M. 2006, ApJ, 650, 102

Ayromlou, M., Nelson, D., Pillepich, A., et al. 2023, arXiv e-prints, arXiv:2311.06339 Cash, W. 1979, ApJ, 228, 939 D'Agostino, R. &amp; Pearson, E. S. 1973, 60, 613, see correction ? . Dehghan, S., Johnston-Hollitt, M., Colless, M., &amp; Miller, R. 2017, MNRAS, 468, 2645 Dennerl, K., Andritschke, R., Bräuninger, H., et al. 2020, in Society of PhotoOptical Instrumentation Engineers (SPIE) Conference Series, Vol. 11444, Space Telescopes and Instrumentation 2020: Ultraviolet to Gamma Ray, ed. J.-W. A. den Herder, S. Nikzad, &amp; K. Nakazawa, 114444Q Edge, A. C., Stewart, G. C., Fabian, A. C., &amp; Arnaud, K. A. 1990, MNRAS, 245, 559 Ettori, S., Ghirardini, V ., Eckert, D., et al. 2019, A&amp;A, 621, A39 Fabian, A. C., Walker, S. A., Russell, H. R., et al. 2017, MNRAS, 464, L1 Federrath, C., Klessen, R. S., Iapichino, L., &amp; Beattie, J. R. 2021, Nature Astronomy, 5, 365 Federrath, C., Roman-Duval, J., Klessen, R. S., Schmidt, W., &amp; Mac Low, M. M. 2010, A&amp;A, 512, A81 Finoguenov, A., Henriksen, M. J., Miniati, F., Briel, U. G., &amp; Jones, C. 2006, ApJ, 643, 790 Foster, A., Smith, R., &amp; Brickhouse, N. S. 2019, in American Astronomical Society Meeting Abstracts, Vol. 233, American Astronomical Society Meeting Abstracts #233, 251.05 Gaspari, M., Churazov, E., Nagai, D., Lau, E. T., &amp; Zhuravleva, I. 2014, A&amp;A, 569, A67 Gatuzz, E., Mohapatra, R., Federrath, C., et al. 2023a, MNRAS, 524, 2945 Gatuzz, E., Sanders, J. S., Canning, R., et al. 2022a, MNRAS, 513, 1932 Gatuzz, E., Sanders, J. S., Dennerl, K., et al. 2023b, MNRAS, 520, 4793 Gatuzz, E., Sanders, J. S., Dennerl, K., et al. 2023c, arXiv e-prints, arXiv:2303.17556 Gatuzz, E., Sanders, J. S., Dennerl, K., et al. 2023d, MNRAS, 525, 6394

- Gatuzz, E., Sanders, J. S., Dennerl, K., et al. 2023e, MNRAS, 526, 396

Gatuzz, E., Sanders, J. S., Dennerl, K., et al. 2022b, MNRAS, 511, 4511 Genel, S., Vogelsberger, M., Springel, V., et al. 2014, MNRAS, 445, 175 Ghirardini, V., Eckert, D., Ettori, S., et al. 2019, A&amp;A, 621, A41 Henriksen, M., Donnelly, R. H., &amp; Davis, D. S. 2000, ApJ, 529, 692 Henriksen, M. J. &amp; Tittley, E. R. 2002, ApJ, 577, 701 Hitomi Collaboration, Aharonian, F., Akamatsu, H., et al. 2016, Nature, 535, 117 Hitomi Collaboration, Aharonian, F., Akamatsu, H., et al. 2018, PASJ, 70, 12 Hofmann, F., Sanders, J. S., Nandra, K., Clerc, N., &amp; Gaspari, M. 2016, A&amp;A, 585, A130

- Kalberla, P. M. W., Burton, W. B., Hartmann, D., et al. 2005, A&amp;A, 440, 775
- Lau, E. T., Kravtsov, A. V., &amp; Nagai, D. 2009, ApJ, 705, 1129
- Liu, A., Yu, H., Diaferio, A., et al. 2018, ApJ, 863, 102
- Liu, A., Yu, H., Tozzi, P., &amp; Zhu, Z.-H. 2015, ApJ, 809, 27
- Liu, A., Yu, H., Tozzi, P., &amp; Zhu, Z.-H. 2016, ApJ, 821, 29
- Lodders, K., Palme, H., &amp; Gail, H.-P. 2009, Landolt Börnstein [ arXiv:0901.1149 ]
- Marinacci, F., Vogelsberger, M., Pakmor, R., et al. 2018, MNRAS, 480, 5113
- Naiman, J. P., Pillepich, A., Springel, V ., et al. 2018, MNRAS, 477, 1206
- Navarro, J. F., Frenk, C. S., &amp; White, S. D. M. 1996, ApJ, 462, 563
- Nelson, D., Pillepich, A., Genel, S., et al. 2015, Astronomy and Computing, 13, 12
- Nelson, D., Pillepich, A., Springel, V., et al. 2018, MNRAS, 475, 624 Ota, N. &amp; Yoshida, H. 2016, PASJ, 68, S19
- Pillepich, A., Nelson, D., Hernquist, L., et al. 2018a, MNRAS, 475, 648

Pillepich, A., Springel, V ., Nelson, D., et al. 2018b, MNRAS, 473, 4077

- Pinto, C., Sanders, J. S., Werner, N., et al. 2015, A&amp;A, 575, A38
- Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2016, A&amp;A, 594, A13
- Predehl, P., Andritschke, R., Arefiev, V., et al. 2021, A&amp;A, 647, A1
- Sanders, J. S., Bi ffi , V., Brüggen, M., et al. 2022a, A&amp;A, 661, A36
- Sanders, J. S., Bi ffi , V., Brüggen, M., et al. 2022b, A&amp;A, 661, A36
- Sanders, J. S., Dennerl, K., Russell, H. R., et al. 2020, A&amp;A, 633, A42
- Sanders, J. S., Fabian, A. C., Russell, H. R., Walker, S. A., &amp; Blundell, K. M. 2016a, MNRAS, 460, 1898
- Sanders, J. S., Fabian, A. C., Taylor, G. B., et al. 2016b, MNRAS, 457, 82 Sauvageot, J. L., Belsole, E., &amp; Pratt, G. W. 2005, A&amp;A, 444, 673
- Seta, A., Federrath, C., Livingston, J. D., &amp; McClure-Gri ffi ths, N. M. 2023, MNRAS, 518, 919
- Shapiro, S. S. &amp; Wilk, M. B. 1965, Biometrika, 52, 591
- Simionescu, A., Werner, N., Finoguenov, A., Böhringer, H., &amp; Brüggen, M. 2008, A&amp;A, 482, 97
- Springel, V., Pakmor, R., Pillepich, A., et al. 2018, MNRAS, 475, 676
- Struble, M. F. &amp; Rood, H. J. 1999, ApJS, 125, 35
- Strüder, L., Briel, U., Dennerl, K., et al. 2001, A&amp;A, 365, L18
- Tamura, T., Yamasaki, N. Y., Iizuka, R., et al. 2014, 782, 38
- Torrey, P., V ogelsberger, M., Genel, S., et al. 2014, MNRAS, 438, 1985
- Vogelsberger, M., Genel, S., Sijacki, D., et al. 2013, MNRAS, 436, 3031
- Vogelsberger, M., Genel, S., Springel, V., et al. 2014a, Nature, 509, 177
- Vogelsberger, M., Genel, S., Springel, V., et al. 2014b, MNRAS, 444, 1518

Weinberger, R., Springel, V., Hernquist, L., et al. 2017, MNRAS, 465, 3291

Wilms, J., Allen, A., &amp; McCray, R. 2000, ApJ, 542, 914

Zhuravleva, I., Churazov, E., Schekochihin, A. A., et al. 2014, Nature, 515, 85 ZuHone, J. A., Miller, E. D., Bulbul, E., &amp; Zhuravleva, I. 2018, ApJ, 853, 180

- ZuHone, J. A., Miller, E. D., Simionescu, A., &amp; Bautz, M. W. 2016, ApJ, 821, 6

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00837v1-figure-22b0d864-cbff-4d79-b35b-e608fc008a76.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00837v1-figure-a0ee6de2-e763-4a0c-8fd4-32235053bc4d.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00837v1-figure-c05328de-0dac-40f7-852f-ead3cf9d61b3.jpg)

Fig. 2: Top panel: total exposure time (s) in the 0.5-9.25 keV energy range. Middle panel: Normalized Fe-K count map, showing the number of counts in each 1.59 arcsec pixel in the Fe-K complex with the instrumental Cu ring included to indicate in which regions the redshift can be measured. Bottom panel: Normalized background Cu count map, showing the number of counts in each 1.59 arcsec pixel in the Cu K α line energy range.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00837v1-figure-1d4693a5-84ba-470f-893d-03bbeeaeb29e.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00837v1-figure-1df2432f-b841-4a74-af5e-0792ed478597.jpg)

-61.29

Fig. 3: Edge-filtered X-ray images using a Gaussian Gradient Magnitude (GGM) filter on scales of 4 pixels ( top panel ) and 8 pixels ( bottom panel ).

Fig. 4: Edge-filtered X-ray images using an adaptive Gaussian Gradient Magnitude (GGM) filter from Sanders et al. (2022b) using a signal to noise ratio of 32, where the scale shows the gradient as the log(10) change per 4 arcsec pixel.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00837v1-figure-9980ed27-13b6-4cb2-baf4-736705128b60.jpg)

Fig. 5: Subtracted fractional di ff erence from the average at each radius for the 0.5-2 keV X-ray surface brightness. The pointlike sources were excluded in the analysis. The X-ray image was smoothed by a 2 σ Gaussian.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00837v1-figure-13b4ee32-335f-4d5e-ba92-9b2255a17cfc.jpg)

Fig. 6: Left panel: velocity map (km s) relative to the A3266 cluster ( / z = 0 0594). Maps were created by moving 2:1 elliptical regions . (i.e., rotated to lie tangentially to the nucleus) containing ∼ 750 counts in the Fe-K region. Right panel: 1 σ velocity uncertainties. White circles correspond to point sources which were excluded from the analysis. The spectral fitting is done in the 4 -9 25 keV . energy range.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00837v1-figure-f09b3568-6d58-4db5-a45d-422c94f1b325.jpg)

Fig. 7: Velocity map (km s) relative to the A3266 cluster, after / subtracting the mean value 312 km s obtained from the spectral / map described in Section 6.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00837v1-figure-0c8dad1f-73cb-4457-acde-9def06c8f7fc.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00837v1-figure-b093d5cb-a54b-4d6d-b84e-68f50b835d70.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00837v1-figure-dbfa9a1e-a420-469c-95a3-2e856f90a183.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00837v1-figure-aa5c458b-3105-4ead-ba28-e84d851cc923.jpg)

Fig. 8: Top panel: temperature map in units of keV. Middle panel: Abundance map relative to abundances from Lodders et al. (2009). Bottom panel: emission measure map. The spectral fitting is done in the 4 -9 25 keV energy range. .

Fig. 9: A3266 annular regions analyzed in subsection 4.2. Green circles correspond to point-like sources excluded from the analysis.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00837v1-figure-1b2c318f-cd09-47f1-8ee5-d64abe74079a.jpg)

Fig. 10: Velocities obtained for each annular region described in subsection 4.2 (numbered from the center to the outside). The A3266 redshift is indicated with a horizontal line. The spectral fitting is done in the 4 -9 25 keV energy range. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00837v1-figure-f7b46cd4-599d-4ee6-85cc-e7ece60960bf.jpg)

Fig. 11: Temperature (top panel) and metallicity (bottom panel) profiles obtained from the best fit results for each annular region described in subsection 4.2. The spectral fitting is done in the 4 -9 25 keV energy range. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00837v1-figure-71b5ebfb-4a94-4a49-8892-7fcfabc75421.jpg)

Fig. 12: Top panel: Manually selected regions following the fractional di ff erence in 0.5-2 keV X-ray surface brightness. Green circles correspond to point-like sources excluded from the analysis. Bottom panel: Velocities obtained for manually selected regions following the A3266 merging structure. The A3266 redshift is indicated with a horizontal line. The spectral fitting is done in the 4 -9 25 keV energy range. We have also included . the mean optical redshift for galaxies from Dehghan et al. (2017) within each one of the regions obtained.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00837v1-figure-2cfd4ca2-4b7e-42f0-9f85-6608317a541d.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00837v1-figure-37842c07-f7e2-4c67-a2d5-9dcf521e1f76.jpg)

Fig. 13: Temperature (top panel) and metallicity (bottom panel) profiles obtained from the best fit results for manually selected regions following the A3266 merging structure. The spectral fitting is done in the 4 -9 25 keV energy range. .

Fig. 14: Top panel Velocity PDF of the hot ICM obtained for the A3266 velocity map. The blue line corresponds to the best-fit using a Gaussian model.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00837v1-figure-00420ccd-e9c3-4a5d-9426-b22a0fa129d7.jpg)

Fig. 15: VSFs of the hot ICM within A3266. Power laws with slope ∼ 1 3 (i.e., Kolmogorov) and / ∼ 2 3 are also included. The / shaded region shows the 1 -σ interval around the mean VSF. The vertical 'dashdot' line indicates the PSF scale.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00837v1-figure-6cd3606c-7a94-4dee-9c25-44785af0f72b.jpg)
E. Gatuzz J. Sanders A. Liu A. Fabian C. Pinto H. Russell D. Eckert S. Walker J. ZuHone R. Mohapatra
Measuring the ICM velocity structure within the A3266 galaxy cluster
We present a detailed analysis of the velocity structure of the hot
intracluster medium (ICM) within the A3266 galaxy cluster, including new
observations taken between June and November 2023. Firstly, morphological
structures within the galaxy cluster were examined using a Gaussian Gradient
Magnitude (GGM) and adaptively smoothed GGM filter applied to the EPIC-pn X-ray
image. Then, we applied a novel {\it XMM-Newton} EPIC-pn energy scale
calibration, which uses instrumental Cu K$\alpha$ as reference for the line
emission, to measure line-of-sight velocities of the hot gas within the system.
This approach enabled us to create two-dimensional projected maps for velocity,
temperature, and metallicity, showing that the hot gas displays a redshifted
systemic velocity relative to the cluster redshift across all fields of view.
Further analysis of the velocity distribution through non-overlapping circular
regions demonstrated consistent redshifted velocities extending up to 1125 kpc
from the cluster core. Additionally, the velocity distribution was assessed
along regions following surface brightness discontinuities, where we observed
redshifted velocities in all regions, with the largest velocities reaching $768
\pm 284$ km/s. Moreover, we computed the velocity Probability Density Function
(PDF) from the velocity map. We applied a normality test, finding that the PDF
adheres to an unimodal normal distribution consistent with theoretical
predictions. Lastly, we computed a velocity structure function (VSF) for this
system using the measured line-of-sight velocities. These insights advance our
understanding of the dynamic processes within the A3266 galaxy cluster and
contribute to our broader knowledge of ICM behavior in merging galaxy clusters.