Tangential discontinuities known as cold fronts (CFs) are abundant in groups and clusters of galaxies (GCs). The relaxed, spiral-type CFs were initially thought to be isobaric, but a significant, $10\%$--$20\%$ jump in the thermal pressure $P_t$ was reported when deprojected CFs were stacked, interpreted as missing $P_t$ below the CFs (i.e. at smaller radii $r$) due to a locally-enhanced nonthermal pressure $P_{nt}$. We report a significant ($\sim4.3\sigma$) deprojected jump in $P_t$ across a single sharp CF in the Centaurus cluster. Additional seven CFs are deprojected in the GCs A2029, A2142, A2204, and Centaurus, all found to be consistent (stacked: $\sim1.9\sigma$) with similar pressure jumps. Combining our sample with high quality deprojected CFs from the literature indicates pressure jumps at significance levels ranging between $2.7\sigma$ and $5.0\sigma$, depending on assumptions. Our nominal results are consistent with $P_{nt}\simeq (0.1\mbox{--}0.3)P_t$ just below the CF. We test different deprojection and analysis methods to confirm that our results are robust, and show that without careful deprojection, an opposite pressure trend may incorrectly be inferred. Analysing all available deprojected data, we also find: (i) a nearly constant CF contrast $q$ of density and temperature within each GC, monotonically increasing with the GC mass $M_{200}$ as $q\propto M_{200}^{0.23\pm0.04}$; (ii) hydrostatic mass discontinuities indicating fast bulk tangential flows below all deprojected CFs, with a mean Mach number $\sim0.76$; and (iii) the newly deprojected CFs are consistent (stacked: $\sim2.9\sigma$) with a $1.25^{+0.09}_{-0.08}$ metallicity drop across the CF. These findings suggest that GCs quite generally harbor extended spiral flows.

Virial shocks around galaxy clusters are expected to show a cutoff in the thermal Sunyaev-Zel'dovich (SZ) signal, coincident with a leptonic ring. However, until now, leptonic virial signals were reported only in Coma and in stacked Fermi-LAT clusters, and an SZ virial shock signal was reported only in A2319. We point out that a few clusters --- presently Coma, A2319, and A2142 --- already show a sharp drop in Planck SZ pressure near the virial radius, coincident with a weak LAT $\gamma$-ray excess. These signatures are shown to trace the virial shocks of the clusters, at joint medium to high confidence levels. The electron acceleration rates inferred from $\gamma$-rays are consistent with previous measurements. The combined signal allows a separate measurement of the $\sim0.5\%$ acceleration efficiency and of the accretion rate. Lower limits of order a few are imposed on the shock Mach numbers.

Galaxy clusters are thought to grow by accreting mass through large-scale, strong, yet elusive, virial shocks. Such a shock is expected to accelerate relativistic electrons, thus generating a spectrally-flat leptonic virial-ring. However, until now, only the nearby Coma cluster has shown evidence for a $\gamma$-ray virial ring. We stack Fermi-LAT data for the 112 most massive, high latitude, extended clusters, enhancing the ring sensitivity by rescaling clusters to their virial radii and utilizing the expected flat energy spectrum. In addition to a central unresolved, hard signal (detected at the $\sim 5.8\sigma$ confidence level), probably dominated by AGN, we identify (at the $5.8\sigma$ confidence level) a bright, spectrally-flat $\gamma$-ray ring at the expected virial shock position. The ring signal implies that the shock deposits $\sim 0.6\%$ (with an interpretation uncertainty factor $\sim2$) of the thermal energy in relativistic electrons over a Hubble time. This result, consistent with the Coma signal, validates and calibrates the virial shock model, and indicates that the cumulative emission from such shocks significantly contributes to the diffuse extragalactic $\gamma$-ray and low-frequency radio backgrounds.

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.

Following evidence for an east--west elongated virial ring around the Coma galaxy cluster in a $\sim220$ GeV VERITAS mosaic, we search for corresponding signatures in $>$GeV $\gamma$-rays from Fermi-LAT, and in soft, $\sim0.1$ keV X-rays from ROSAT. For the ring elongation and orientation inferred from VERITAS, we find a $3.4\sigma$ LAT excess, and detect ($>5\sigma$) the expected signature in ROSAT bands R1 and R1+R2. The significances of both LAT and ROSAT signals are maximal near the VERITAS ring parameters. The intensities of the ROSAT, Fermi, and VERITAS signals are consistent with the virial shock depositing $\sim0.3\%$ (with an uncertainty factor of $\sim3$) of its energy over a Hubble time in a nearly flat, $p\equiv - d\ln N_e/d\ln E\simeq 2.0$--$2.2$ spectrum of cosmic-ray electrons. The sharp radial profiles of the LAT and ROSAT signals suggest preferential accretion in the plane of the sky, as indicated by the distribution of neighboring large-scale structure. The X-ray signal gauges the compression of cosmic-rays as they are advected deeper into the cluster.

Large scale structures such as groups and clusters of galaxies show a universal, nearly linear entropy radial profile $K(r)$. Using deprojected 16 clusters and 12 groups from the literature, we find that $K\propto r^{0.96\pm0.01}$, consistent with the mean power-law index $(0.9\mbox{--}1.1)$ of previous studies. A similarly good fit is given by a $\tau\propto r^{0.72\pm0.01}$ ratio between cooling and free-fall times. Both profiles slightly flatten at small radii, as $\tau$ becomes of order unity. The entropy profile is usually attributed to self-similar shock accretion (shown to be inconsistent with the data), to non-standard heat conduction, or to turbulent heating. We argue that a dynamical mechanism is needed to sustain such a universal profile, oblivious to the temperature peak at the edge of the core and to the virial shock at the outskirts, and robust to the presence of ongoing cooling, merger, and active galactic nucleus (AGN) activity. In particular, we show that such a profile can be naturally obtained in a spiral flow, which is likely to underlie most galaxy aggregates according to the ubiquitous spiral patterns and cold fronts observed. Generalizing a two-phase spiral flow model out to the virial radius surprisingly reproduces the thermal profile. A generalized Schwarzschild criterion indicates that observed spiral patterns must involve a convective layer, which may regulate the thermal profile.

Tangential discontinuities, seen as X-ray edges known as cold fronts (CFs), are ubiquitous in cool-core galaxy clusters. We analyze all 17 deprojected CF thermal profiles found in the literature, including three new CFs we tentatively identify (in clusters A2204 and 2A0335). We discover small but significant thermal pressure drops below all nonmerger CFs, and argue that they arise from strong magnetic fields below and parallel to the discontinuity, carrying 10%-20% of the pressure. Such magnetization can stabilize the CFs, and explain the CF-radio minihalo connection.