2026ApJ...998..201L

Double radio relics, pairs of diffuse radio features located on opposite sides of merging galaxy clusters, are a rare subclass of radio relics that are believed to trace merger shocks and provide valuable constraints on plasma acceleration models and merger history. With the number of known double relics growing in recent and upcoming radio surveys, statistical analyses of their properties are becoming feasible. In this study, we utilize the cosmological magnetohydrodynamics zoomed-in simulations TNG-Cluster, in combination with TNG300-1, to examine the statistical properties of double radio relics. The simulated double-relic pairs exhibit a wide range of luminosity ratios, broadly consistent with the observations. We find that the two relics in a given double system often differ significantly in their shock properties and magnetic field strengths. This diversity implies that the observed brightness asymmetry in the pair cannot be explained by a single factor alone, but instead reflects an interplay of multiple physical parameters. Nevertheless, double radio relics tend to align with the collision axis within ∼30°, and their separation (d_drr) correlates tightly with the time since collision (TSC) as TSC [Gyr] = 0.52d_drr/R_500c - 0.24, allowing it to be inferred with an accuracy of ∼0.2 Gyr. With the statistical samples of simulated radio relics, we predict that low-mass clusters will constitute the dominant population of double radio relic systems detected with upcoming surveys such as the Square Kilometre Array. These results demonstrate that double radio relics can serve as robust probes of merger dynamics and plasma acceleration, and that simulations provide critical guidance for interpreting the large samples expected from next-generation radio surveys.