X-Ray Analysis of an Off-Axis Merger Stage Binary Galaxy Cluster: PSZ2 G279.79+39.09

TL;DR

This study presents a detailed X-ray analysis of the bimodal galaxy cluster PSZ2 G279.79+39.09 at $z=0.29$ using archival XMM-Newton and Chandra data to diagnose its merger state. By fitting surface-brightness profiles with multi-component $eta$ models, deriving hydrostatic masses, and constructing temperature, pressure, and entropy maps, the authors characterize an off-axis merger with a hot, high-entropy bridge and a residual cool core in the eastern subcluster. The results yield a mass ratio of $ ext{M}_{500,E}/ ext{M}_{500,W} oughly1:1.7$, a weak merger shock with Mach number $M\,\approx\,1.41$, and a projected core separation of about $1.35$ Mpc, suggesting a pre- or near-core-passage configuration. The work highlights the role of thermodynamic bridges and morphological indicators in constraining merger timelines and motivates follow-up optical spectroscopy and tailored hydrodynamic simulations to refine the dynamical scenario.

Abstract

We present an X-ray analysis of the merging galaxy cluster system PSZ2 G279.79+39.09 ($z=0.29$) using archival XMM-Newton and Chandra observations. The surface brightness image is bimodal, elongated east-west with a projected core separation of $\sim 1.35$ Mpc. We measure gas temperatures of 5.36 keV for the eastern subcluster (PSZ-E) and 5.44 keV for the western component (PSZ-W). Assuming isothermal intracluster gas, the hydrostatic masses are $\log(M_{500}/M_\odot)=14.76$ for PSZ-E and 14.54 for PSZ-W, implying a mass ratio of $\sim 1:1.7$. PSZ-E shows X-ray concentration indices of $c_{40}/c_{400}=0.124$ and $c_{100}/c_{500}=0.278$, together with a centroid shift of $w=0.016$, indicating a disturbed halo that still hosts a compact cool core; PSZ-W is comparably disturbed even in its core. Both subclusters exhibit ICM asymmetries consistent with ram-pressure stripping, and PSZ-W displays an X-ray tail extending nearly to the outskirts of PSZ-E. The orientation and length of this tail support an off-axis merger geometry. Thermodynamic maps reveal a hot ($\sim 7.3$ keV), high-pressure, high-entropy bridge between the cores. From the Rankine-Hugoniot temperature jump, we infer a Mach number $M=1.41^{+0.33}_{-0.30}$, consistent with a weak merger shock propagating at $1620^{+500}_{-420}$ km s$^{-1}$. These results indicate a merger with a non-zero impact parameter, likely observed near core passage ($\lesssim 0.5$ Gyr before or after), with the pre-pericenter scenario slightly preferred based on the projected separation and thermodynamic structure.

Figures (5)

• Figure 1: Left: The adaptively-smoothed vignetting-corrected background subtracted XMM-Newton image at 0.5 - 2.0 keV band. Right: The 4.5$\sigma$ smoothed background subtracted Chandra image at 0.5 - 2.0 keV band. The dashed circles represent each sub-cluster's ICM, which will later be used to obtain the mean cluster spectral properties. The apparent differences in morphology between the XMM-Newton and Chandra images are due to instrumental effects. The Chandra image (ACIS-I) has finer angular resolution and smaller smoothing, revealing compact features and apparent discontinuities where CCD gaps intersect diffuse emission, while the broader PSF and adaptive smoothing in XMM-Newton produce a smoother, more symmetric appearance.
• Figure 2: Left: The adaptively-smoothed vignetting-corrected background subtracted XMM-Newton image at 0.5 - 2.0 keV band demonstrating the selected region to generate surface brightness profile. Right: PSF-convolved X-ray surface-brightness profile along the merger axis. The model (black) combines a double-$\beta$ eastern core, single-$\beta$ western core, and flat-topped mesa bridge ($n=8$).
• Figure 3: Top Left: X-ray temperature map obtained using adaptive circular binning and generating spectral templates via APEC model in XSPEC. Top Middle: The pressure map. Top Right: The entropy map. Bottom Left: The temperature error map. Bottom Middle: The pressure error map. Bottom Right: The entropy error map. The selected cold and hot regions for spectral fitting are shown for visual aid. Due to low count statistics in the bridge region, some spaxels yield significantly high values with substantial errors.
• Figure 4: The projection angle ($\alpha$) as a function of the relative projected radial velocity difference ($V_r$) of our sub-clusters assuming a zero angular‐momentum head-on merger case. The black curve represents the limit of bound solutions due to the Newtonian criterion. The red and the blue dashed lines represent the bound and the unbound solutions, respectively. Three magenta points represent the possible two BI and one UO solutions. The red dotted line corresponds to 68% confidence ranges. Using the photometric redshifts from 2024Kluge, the radial velocity difference between the two BCGs is estimated to be $\sim780\pm2160$ km s$^{-1}$, indicating that the redshift difference is not statistically significant. Nevertheless, to perform a simple Newtonian binding criterion test, we assume a nominal 10% uncertainty on the radial velocity difference (see the vertical gray solid and dashed lines) to explore the potential dynamical configuration of the system. A detailed spectroscopic follow-up campaign will be required to accurately determine the true velocity difference and dynamical state of the two subclusters.
• Figure 5: Posterior probability distributions of the three-dimensional infall velocities (left), three-dimensional separations (middle), and orbital phase parameters (right) for all dynamically allowed off-axis orbital families. The orbital phase variable $f_{t}$ represents the dimensionless time along the orbit: $f_{t,\mathrm{in}}$ measures the fractional time from apocenter to pericenter for the bound--incoming (first-infall) solutions; $f_{t,\mathrm{out}}$ measures the fractional time from pericenter to apocenter for the bound--outgoing (post-pericenter) solutions; and $f_{t,\mathrm{unbound}}$ is shown via the absolute hyperbolic anomaly $|F_{\rm hyp}|$ for unbound trajectories. Top row: bound--incoming (BI) solutions. Middle row: bound--outgoing (BO) solutions. Bottom row: unbound hyperbolic solutions. Vertical black, purple, and blue lines mark the weighted mode, median, and mean of each distribution, respectively.