We study the disruption process of hierarchical 3-body systems with bodies of comparable mass. Such systems have long survival times that vary by orders of magnitude depending on the initial conditions. By comparing with 3-body numerical integrations, we show that the evolution and disruption of such systems can be statistically described as a simple random-walk process in the outer-orbit's energy, where the energy-exchange per pericenter passage (step-size) is calculated from the initial conditions. In our derivation of the step-size, we use previous analytic results for parabolic encounters, and average over the (Kozai-Lidov) oscillations in orbital parameters, which are faster then the energy diffusion timescale. While similar random-walk models were studied before, this work differs in two manners: (a) this is the first time that the Kozai-Lidov averaged step-size is derived from first principles and demonstrated to reproduce the statistical evolution of numerical ensembles without fitting parameters, and (b) it provides a characteristic life-time, instead of answering the binary question (stable/unstable), set by case-specific criteria.

In the hierarchical paradigm of structure formation, galaxy clusters are the largest objects ever to virialize. They are thought to grow by accreting mass through large scale, strong virial shocks. Such a collisionless shock is expected to accelerate relativistic electrons, thus generating a spectrally flat leptonic virial ring. However attempts to detect virial rings have all failed, leaving the shock paradigm unconfirmed. Here we identify a virial $\gamma$-ray signal by stacking Fermi-LAT data for 112 clusters, enhancing the ring sensitivity by rescaling clusters to their virial radii and utilizing the anticipated spectrum. In addition to a central unresolved, hard signal (detected at the nominal $5.8\sigma$ confidence level), probably dominated by active galactic nuclei, we identify ($5.9\sigma$) a bright, spectrally flat $\gamma$-ray ring at the expected shock position. It corresponds to $\sim 0.6\%$ (with an uncertainty factor $\sim2$) thermal energy deposition in relativistic electrons over a Hubble time. This result validates the shock paradigm, calibrates its parameters, and indicates that the cumulative emission from such shocks significantly contributes to the diffuse extragalactic $\gamma$-ray and radio backgrounds.