2025ApJ...978..122L

In a binary merger of two subclusters with comparable masses, a pair of merger shocks are typically generated, often manifesting as double radio relics. Using cosmological hydrodynamic simulations, we identify major merger events with mass ratio $1\lt {{ \mathcal M }}_{1}/{{ \mathcal M }}_{2}\lesssim 4$ and impact parameter b/r_vir,1 ≲ 1, where r_vir,1 is the virial radius of the heavier subcluster. We analyze merger shock surfaces approximately 1 Gyr after the pericenter passage, focusing on their morphology and the distribution of the Mach number, M_s, of their constituent shock zones. The shock surfaces exhibit an elongated shape with a minor-to-major axis ratio of ∼0.6-0.9 and cover the area of ∼5%-20% of the enclosed sphere. The area ratio of the two shock surfaces roughly scales as ${A}_{{\rm{ss}},{\rm{2}}}/{A}_{{\rm{ss}},{\rm{1}}}\propto {{ \mathcal M }}_{1}/{{ \mathcal M }}_{2}$, typically positioning the larger shock ahead of the lighter subcluster. The axis connecting the two subclusters generally does not pass through the centers of the shock surfaces, due to the nonzero impact parameter and ambient turbulent flows. The distribution of M_s of shock zones on each surface can be approximated by a log-normal function, peaking at M_s,peak ≈ 2-4.5 and extending up to ∼10. The surface-area-weighted and X-ray- emissivity-weighted average Mach numbers are comparable, with ${\langle {M}_{s}\rangle }_{\mathrm{area}}\approx 2.3-4.4$ and ${\langle {M}_{s}\rangle }_{X}\approx 2-4$. In contrast, the cosmic-ray energy flux-weighted average Mach numbers are higher with ${\langle {M}_{s}\rangle }_{\mathrm{CR}}\approx 3-5$. This discrepancy aligns with the differences between Mach numbers derived from X-ray and radio observations of radio relic shocks. On the other hand, we find that mostly ${\langle {M}_{s}\rangle }_{X}\gtrsim 2$ for simulated merger shocks, although shocks with M_X–ray ≲ 2 are often reported in observations.