Collisionless shocks are often studied in two spatial dimensions (2D), to gain insights into the 3D case. We analyze diffusive shock acceleration for an arbitrary number $N\in\mathbb{N}$ of dimensions. For a non-relativistic shock of compression ratio $\mathcal{R}$, the spectral index of the accelerated particles is $s_E=1+N/(\mathcal{R}-1)$; this curiously yields, for any $N$, the familiar $s_E=2$ (i.e., equal energy per logarithmic particle energy bin) for a strong shock in a mono-atomic gas. A precise relation between $s_E$ and the anisotropy along an arbitrary relativistic shock is derived, and is used to obtain an analytic expression for $s_E$ in the case of isotropic angular diffusion, affirming an analogous result in 3D. In particular, this approach yields $s_E = (1+\sqrt{13})/2 \simeq 2.30$ in the ultra-relativistic shock limit for $N=2$, and $s_E(N\to\infty)=2$ for any strong shock. The angular eigenfunctions of the isotropic-diffusion transport equation reduce in 2D to elliptic cosine functions, providing a rigorous solution to the problem; the first function upstream already yields a remarkably accurate approximation. We show how these and additional results can be used to promote the study of shocks in 3D.

The intracluster medium of a galaxy cluster often shows an extended quasi-spiral structure, accentuated by tangential discontinuities known as cold fronts (CFs). These discontinuities are thought to isolate between low-entropy, high-metallicity gas inside (i.e., below) the CF that was lifted from the center of the cluster by some radial factor $f_i$, and high-entropy, low-metallicity gas outside the CF that was pushed inward by a factor $f_o$. We find broad support for such a picture, by comparing the entropy and metallicity discontinuities with the respective azimuthal averages, using newly deprojected thermal profiles in clusters A2029, A2142, A2204 and Centaurus, supplemented by deprojected CFs from the literature. In particular, the mean advection factors $f_K$ and $f_Z$, inferred respectively from entropy and metallicity, strongly correlate ($\mathcal{R}=0.7^{+0.2}_{-0.3}$) with each other, consistent with large-scale advection. However, unlike sloshing simulations, in which the inside/outside phases are an inflow/outflow settling back to equilibrium after a violent perturbation, our results are more consistent with an outflow/inflow, with the fast, Mach $\mathcal{M}_i\sim0.8$ gas inside the CF being a rapidly heated or mixed outflow, probably originating from the cD galaxy, and gas outside the CF being a $\mathcal{M}_o\sim0.03$, slowly-cooling inflow. In particular, entropy indicates an outside advection factor $f_{Ko}\simeq 1.33\pm0.04$ that is approximately constant in all CFs, gauging the distance traversed by inflowing gas within a cooling time. In contrast, $1.1\lesssim f_{Ki}\lesssim 2.5$ and $1\lesssim f_Z\lesssim 17$ vary considerably among clusters, and strongly correlate ($3.1\sigma\text{-}4.2\sigma$) with the virial mass, $f_{Ki}\propto M_{200}^{0.14\pm0.07}$ and $f_Z\propto M_{200}^{1.4\pm0.4}$, suggesting that each cluster sustains a quasi-steady spiral flow.

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.

The spectral index $s$ of particles diffusively accelerated in a relativistic shock depends on the unknown angular diffusion function $\mathcal{D}$, which itself depends on the particle distribution function $f$ if acceleration is efficient. We develop a relaxation code to compute $s$ and $f$ for an arbitrary functional $\mathcal{D}$ that depends on $f$. A local $\mathcal{D}(f)$ dependence is motivated and shown, when rising (falling) upstream, to soften (harden) $s$ with respect to the isotropic case, shift the angular distribution towards upstream (downstream) directions, and strengthen (weaken) the particle confinement to the shock; an opposite effect on $s$ is found downstream. However, variations in $s$ remain modest even when $\mathcal{D}$ is a strong function of $f$, so the standard, isotropic-diffusion results remain approximately applicable unless $\mathcal{D}$ is both highly anisotropic and not a local function of $f$. A mild, $\sim 0.1$ softening of $s$, in both 2D and 3D, when $\mathcal{D}(f)$ rises sufficiently fast, may be indicated by ab-initio simulations.

We point out that particles accelerated in a non-relativistic shock of compression ratio $r$ attain the standard, $p=(r+2)/(r-1)$ spectral index only under certain conditions. Previous derivations of the spectrum, based on the approximations of spatial diffusion or negligible anisotropy, are shown to be unjustified for a general scattering function. We explain and demonstrate numerically that in contrast to previous claims, $p$ can substantially deviate from the standard result for anisotropic scattering. We prove analytically that the standard approach is nevertheless valid in the limit of an isotropic medium. Additional spectral modifications, for example by motions of scattering modes at intermediate optical depths from the shock, are discussed.

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.

Shocks produced by the accretion of infalling gas in the outskirt of galaxy clusters are expected in the hierarchical structure formation scenario, as found in cosmological hydrodynamical simulations. Here, we report the detection of a shock front at a large radius in the pressure profile of the galaxy cluster A2319 at a significance of $8.6\sigma$, using Planck thermal Sunyaev-Zel'dovich data. The shock is located at $(2.93 \pm 0.05) \times R_{500}$ and is not dominated by any preferential radial direction. Using a parametric model of the pressure profile, we derive a lower limit on the Mach number of the infalling gas, $\mathcal{M} > 3.25$ at 95\% confidence level. These results are consistent with expectations derived from hydrodynamical simulations. Finally, we use the shock location to constrain the accretion rate of A2319 to $\dot{M} \simeq (1.4 \pm 0.4) \times 10^{14}$ M$_\odot$ Gyr$^{-1}$, for a total mass, $M_{200} \simeq 10^{15}$ M$_\odot$.

Diffusive shock acceleration (DSA) by relativistic shocks is thought to generate the $dN/dE\propto E^{-p}$ spectra of charged particles in various astronomical relativistic flows. We show that for test particles in one dimension (1D), $p^{-1}=1-\ln\left[\gamma_d(1+\beta_d)\right]/\ln\left[\gamma_u(1+\beta_u)\right]$, where $\beta_u$ ($\beta_d)$ is the upstream (downstream) normalized velocity, and $\gamma$ is the respective Lorentz factor. This analytically captures the main properties of relativistic DSA in higher dimensions, with no assumptions on the diffusion mechanism. Unlike 2D and 3D, here the spectrum is sensitive to the equation of state even in the ultra-relativistic limit, and (for a J{\"u}ttner-Synge equation of state) noticeably hardens with increasing $1<\gamma_u<57$, before logarithmically converging back to $p(\gamma_u\to\infty)=2$. The 1D spectrum is sensitive to drifts, but only in the downstream, and not in the ultra-relativistic limit.

We analyze the IceCube four-year neutrino data in search of a signal from the Fermi bubbles. No signal is found from the bubbles or from their dense shell, even when taking into account the softer background. This imposes a conservative $\xi_i<8\%$ upper limit on the cosmic-ray ion (CRI) acceleration efficiency, and an $\eta\equiv \xi_e/\xi_i \gtrsim0.006$ lower limit on the electron-to-ion ratio of acceleration efficiencies (at the $2\sigma$ confidence level). For typical $\xi_i$, a signal should surface once the number of IceCube neutrinos increases by $\sim$an order of magnitude, unless there is a $<$PeV cutoff on the CRI spectrum.

Galaxy clusters, the largest gravitationally bound objects in the Universe, are thought to grow by accreting mass from their surroundings through large-scale virial shocks. Due to electron acceleration in such a shock, it should appear as a $\gamma$-ray, hard X-ray, and radio ring, elongated towards the large-scale filaments feeding the cluster, coincident with a cutoff in the thermal Sunyaev-Zel'dovich (SZ) signal. However, no such signature was found until now, and the very existence of cluster virial shocks has remained a theory. We find preliminary evidence for a large, $\sim 5$ Mpc minor axis $\gamma$-ray ring around the Coma cluster, elongated towards the large scale filament connecting Coma and Abell 1367, detected at the nominal $2.7\sigma$ confidence level ($5.1\sigma$ using control signal simulations). The $\gamma$-ray ring correlates both with a synchrotron signal and with the SZ cutoff, but not with Galactic tracers. The $\gamma$-ray and radio signatures agree with analytic and numerical predictions, if the shock deposits $\sim 1\%$ of the thermal energy in relativistic electrons over a Hubble time, and $\sim 1\%$ in magnetic fields. The implied inverse-Compton and synchrotron cumulative emission from similar shocks can significantly contribute to the diffuse extragalactic $\gamma$-ray and low frequency radio backgrounds. Our results, if confirmed, reveal the prolate structure of the hot gas in Coma, the feeding pattern of the cluster, and properties of the surrounding large scale voids and filaments. The anticipated detection of such shocks around other clusters would provide a powerful new cosmological probe.

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.

The nature of the bipolar, $\gamma$-ray Fermi bubbles (FB) is still unclear, in part because their faint, high-latitude X-ray counterpart has until now eluded a clear detection. We stack ROSAT data at varying distances from the FB edges, thus boosting the signal and identifying an expanding shell behind the southwest, southeast, and northwest edges, albeit not in the dusty northeast sector near Loop I. A Primakoff-like model for the underlying flow is invoked to show that the signals are consistent with halo gas heated by a strong, forward shock to $\sim$keV temperatures. Assuming ion--electron thermal equilibrium then implies a $\sim10^{56}$ erg event near the Galactic centre $\sim7$ Myr ago. However, the reported high absorption-line velocities suggest a preferential shock-heating of ions, and thus more energetic ($\sim 10^{57}$ erg), younger ($\lesssim 3$ Myr) FBs.

Current measurements of the $\gamma$-ray Fermi bubbles (FB) are based on model-dependent tracers, carry substantial systematic uncertainties, and are at some tension with each other. We show that gradient filters pick out the FB edges, which are found to smoothly connect to the bipolar X-ray structure emanating from the Galactic center, thus supporting the interpretation of the FBs as a Galactic-scale phenomenon. The sharp edges facilitate a direct, model-free measurement of the peripheral FB spectrum. The result is strikingly similar to the full FB-integrated spectrum, softened by a power law of index $\eta\simeq (0.2\mbox{--}0.3)$. This is naturally explained, in both hadronic and leptonic models, if cosmic rays are injected at the edge, and diffuse away preferentially at higher energies $E$. The inferred, averaged diffusion function in the (more plausible) leptonic model, $D(E)\simeq 10^{29.5}(E/10\mbox{ GeV})^{0.48\pm0.02}\mbox{ cm}^2\mbox{ s}^{-1}$, is consistent with estimates for Kraichnan-like turbulence. Our results, in particular the minute spatial variations in $\eta$, indicate that the FB edge is a strong, Mach $\gtrsim5$, forward shock.

Compressible flows around blunt objects have diverse applications, but current analytic treatments are inaccurate and limited to narrow parameter regimes. We show that the gas-dynamic flow in front of an axisymmetric blunt body is accurately derived analytically using a low order expansion of the perpendicular gradients in terms of the parallel velocity. This reproduces both subsonic and supersonic flows measured and simulated for a sphere, including the transonic regime and the bow shock properties. Some astrophysical implications are outlined, in particular for planets in the solar wind and for clumps and bubbles in the intergalactic medium. The bow shock standoff distance normalized by the obstacle curvature is $\sim 2/(3g)$ in the strong shock limit, where $g$ is the compression ratio. For a subsonic Mach number $M$ approaching unity, the thickness $\delta$ of an initially weak, draped magnetic layer is a few times larger than in the incompressible limit, with amplification $\sim ({1+1.3M^{2.6}})/({3\delta})$.

Compressible flows around blunt objects have diverse applications, but present analytic treatments are inaccurate and limited to narrow parameter regimes. We show that the flow in front of an axisymmetric body is accurately derived analytically using a low order expansion of the perpendicular gradients in terms of the parallel velocity. This reproduces both subsonic and supersonic flows measured and simulated for a sphere, including the transonic regime and the bow shock properties.

Following a suggestion that a directed relativistic explosion may have a universal intermediate asymptotic, we derive a self-similar solution for an ultra-relativistic jetted blast wave. The solution involves three distinct regions: an approximately paraboloid head where the Lorentz factor $\gamma$ exceeds $\sim1/2$ of its maximal, nose value; a geometrically self-similar, expanding envelope slightly narrower than a paraboloid; and an axial core in which the (cylindrically, henceforth) radial flow $u$ converges inward towards the axis. Most ($\sim 80\%$) of the energy lies well beyond the leading, head region. Here, a radial cross section shows a maximal $\gamma$ (separating the core and the envelope), a sign reversal in $u$, and a minimal $\gamma$, at respectively $\sim 1/6$, $\sim1/4$, and $\sim3/4$ of the shock radius. The solution is apparently unique, and approximately agrees with previous simulations, of different initial conditions, that resolved the head. This suggests that unlike a spherical relativistic blast wave, our solution is an attractor, and may thus describe directed blast waves such as in the external shock phase of a $\gamma$-ray burst.

We present the results of a very deep (500 ks) Chandra observation, along with tailored numerical simulations, of the nearest, best resolved cluster cold front in the sky, which lies 90 kpc (19 arcmin) to the north-west of M 87. The northern part of the front appears the sharpest, with a width smaller than 2.5 kpc (1.5 Coulomb mean free paths; at 99 per cent confidence). Everywhere along the front, the temperature discontinuity is narrower than 4-8 kpc and the metallicity gradient is narrower than 6 kpc, indicating that diffusion, conduction and mixing are suppressed across the interface. Such transport processes can be naturally suppressed by magnetic fields aligned with the cold front. Interestingly, comparison to magnetohydrodynamic simulations indicates that in order to maintain the observed sharp density and temperature discontinuities, conduction must also be suppressed along the magnetic field lines. However, the northwestern part of the cold front is observed to have a non-zero width. While other explanations are possible, the broadening is consistent with the presence of Kelvin-Helmholtz instabilities (KHI) on length-scales of a few kpc. Based on comparison with simulations, the presence of KHI would imply that the effective viscosity of the intracluster medium is suppressed by more than an order of magnitude with respect to the isotropic Spitzer-like temperature dependent viscosity. Underneath the cold front, we observe quasi-linear features that are ~10 per cent brighter than the surrounding gas and are separated by ~15 kpc from each other in projection. Comparison to tailored numerical simulations suggests that the observed phenomena may be due to the amplification of magnetic fields by gas sloshing in wide layers below the cold front, where the magnetic pressure reaches ~5-10 per cent of the thermal pressure, reducing the gas density between the bright features.

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.

Magnetization in highly conductive plasmas is ubiquitous to astronomical systems. Flows in such media can be described by three path functions $\Lambda_\alpha$, or, for a steady flow, by two stream functions $\lambda_\kappa$ and an additional field such as the mass density $\rho$, velocity $v$, or travel time $\Delta t$. While typical analyses of a frozen magnetic field $\boldsymbol{B}$ are problem-specific and involve nonlocal gradients of the fluid element position $\boldsymbol{x}(t)$, we derive the general, local (in $\Lambda$ or $\lambda$ space) solution $\boldsymbol{B}=(\partial\boldsymbol{x}/\partial\Lambda_\alpha)_t \tilde{B}_\alpha \rho/\tilde{\rho}$, where Lagrangian constants denoted by a tilde are directly fixed at a boundary hypersurface $\tilde{H}$ on which $\boldsymbol{B}$ is known. For a steady flow, $\tilde{\rho}\boldsymbol{B}/\rho=(\partial\boldsymbol{x}/\partial\lambda_\kappa)_{\Delta t}\tilde{B}_\kappa+\boldsymbol{v}\tilde{B}_3/\tilde{v}$; here the electric field $\boldsymbol{E}\sim(\tilde{B}_2\boldsymbol{\nabla}\lambda_1-\tilde{B}_1\boldsymbol{\nabla}\lambda_2)/\tilde{\rho}$ depends only on $\lambda_\kappa$ and the boundary conditions. Illustrative special cases include compressible axisymmetric flows and incompressible flows around a sphere, showing that viscosity and compressibility enhance the magnetization and lead to thicker boundary layers. Our method is especially useful for directly computing electric fields, and for addressing upstream magnetic fields that vary in spacetime. We thus estimate the electric fields above heliospheres and magnetospheres, compute the draping of magnetic substructure around a planetary body, and demonstrate the resulting inverse polarity reversal layer. Our analysis can be incorporated into existing hydrodynamic codes that are based on stream or path functions. (Abridged)

We examine the conditions for the revival of the stalled accretion shock in core-collapse supernovae, in the context of the neutrino heating mechanism. We combine one dimensional simulations of the shock revival process with a derivation of a quasi-stationary approximation, which is both accurate and efficient in predicting the flow. In particular, this approach is used to explore how the evolution of the system depends on the shock radius, $R_S$, and velocity, $V_S$ (in addition to other global properties of the system). We do so through a phase space analysis of the shock acceleration, $a_S$, in the $R_S-V_S$ plane, shown to provide quantitative insights into the initiation of runaway expansion and its nature. In the particular case of an initially stationary ($V_S=0,\;a_S=0$) profile, the prospects for an explosion can be reasonably assessed by the initial signs of the partial derivatives of the shock acceleration, in analogy to a linear damped/anti-damped oscillator. If $\partial a_S/\partial R_S<0$ and $\partial a_S/\partial V_S>0$, runaway expansion will likely occur after several oscillations, while if $\partial a_S/\partial R_S>0$, runaway expansion will commence in a non-oscillatory fashion. These two modes of runaway correspond to low and high mass accretion rates, respectively. We also use the quasi-stationary approximation to assess the advection-to-heating timescale ratio in the gain region, often used as an explosion proxy. Indeed, this ratio does tend to $\sim1$ in conjunction with runaway conditions, but neither this unit value nor the specific choice of the gain region as a point of reference appear to be distinct conditions in this regard.

We review the physics of relativistic shocks, which are often invoked as the sources of non-thermal particles in pulsar wind nebulae (PWNe), gamma-ray bursts (GRBs), and active galactic nuclei (AGN) jets, and as possible sources of ultra-high energy cosmic-rays. We focus on particle acceleration and magnetic field generation, and describe the recent progress in the field driven by theory advances and by the rapid development of particle-in-cell (PIC) simulations. In weakly magnetized or quasi parallel-shocks (where the magnetic field is nearly aligned with the flow), particle acceleration is efficient. The accelerated particles stream ahead of the shock, where they generate strong magnetic waves which in turn scatter the particles back and forth across the shock, mediating their acceleration. In contrast, in strongly magnetized quasi-perpendicular shocks, the efficiencies of both particle acceleration and magnetic field generation are suppressed. Particle acceleration, when efficient, modifies the turbulence around the shock on a long time scale, and the accelerated particles have a characteristic energy spectral index of ~ 2.2 in the ultra-relativistic limit. We discuss how this novel understanding of particle acceleration and magnetic field generation in relativistic shocks can be applied to high-energy astrophysical phenomena, with an emphasis on PWNe and GRB afterglows.

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.

The rotation of a black hole splits its spectrum in two, yet only one sector is known in the highly-damped regime. We find the second, at least partly oblate sector, with quasinormal modes approaching the total reflection frequencies \omega(n>>1) ~ m\Omega - 2\pi iT(n-s), where \Omega and T are the horizon angular velocity and temperature, s is the field spin, and m is an azimuthal eigenvalue. Some physical implications are discussed.

The explosion of a core-collapse supernova can be approximated by the breakdown of steady-state solutions for accretion onto a proto-neutron star (PNS). We analytically show that as the neutrino luminosity exceeds a critical value L_c, the neutrinosphere pressure exceeds the hydrostatic limit even for an optimal shock radius R. This yields L_c \propto M^2 T^2 (with logarithmic corrections) and R \propto M/T, in agreement with numerical results, where M, T are the PNS mass, neutrino temperature. The near-critical flow can be approximated as a ballistic shell on top of an isothermal layer.

We argue that bulk spiral flows are ubiquitous in the cool cores (CCs) of clusters and groups of galaxies. Such flows are gauged by spiral features in the thermal and chemical properties of the intracluster medium, by the multi-phase properties of CCs, and by X-ray edges known as cold fronts. We analytically show that observations of piecewise-spiral fronts impose strong constraints on the CC, implying the presence of a cold, fast flow, which propagates below a hot, slow inflow, separated by a slowly rotating, trailing, quasi spiral, tangential discontinuity surface. This leads to the nearly logarithmic spiral pattern, two-phase plasma, \rho \sim r^{-1} density (or T \sim r^{0.4} temperature) radial profile, and ~100kpc size, characteristic of CCs. By advecting heat and mixing the gas, such flows can eliminate the cooling problem, provided that a feedback mechanism regulates the flow. In particular, we present a quasi-steady-state model for an accretion-quenched, composite flow, in which the fast phase is an outflow, regulated by active galactic nucleus bubbles, reproducing the observed low star formation rates and explaining some features of bubbles such as their R_b \propto r size. The simplest two-component model reproduces several key properties of CCs, so we propose that all such cores harbor a spiral flow. Our results can be tested directly in the next few years, for example by ASTRO-H.

Some galaxy clusters show diffuse radio emission in the form of giant halos (GHs) on Mpc scales or minihalos (MHs) on smaller scales. Comparing VLA and XMM radial profiles of several such clusters, we find a universal linear correlation between radio and X-ray surface brightness, valid in both types of halos. It implies a halo central emissivity nu*j_nu=10^{-31.4+-0.2}*(n/10^-2cm^-3)^2*(T/T_0)^{0.2+-0.5} erg/s/cm^3, where T and T_0 are the local and central temperatures, and n is the electron number density. We argue that the tight correlation and the scaling of j_nu, combined with morphological and spectral evidence, indicate that both GHs and MHs arise from secondary electrons and positrons, produced in cosmic-ray ion (CRI) collisions with a strongly magnetized, B >3 muG intracluster gas. When the magnetic energy density drops below that of the microwave background, the radio emission weakens considerably, producing halos with a clumpy morphology (e.g. RXCJ2003.5-2323 and A2255) or a distinct radial break. We thus measure a magnetic field B=3 muG at a radius r~110 kpc in A2029 and r~50 kpc in Perseus. The spectrum of secondaries, produced from hadronic collisions of ~20 GeV CRIs, reflects the energy dependence of the collision cross section. We use the observed spectra of halos, in particular where they steepen with increasing radius or frequency, to (i) measure B~10(nu/700 MHz) muG, with nu the spectral break frequency; (ii) identify a correlation between the average spectrum and the central magnetic field; and (iii) infer a CRI spectral index s <-2.7 and energy fraction xi_p~10^{-3.6+-0.2} at particle energies above 10 GeV. Our results favor a model where CRIs diffuse away from their sources (which are probably supernovae, according to a preliminary correlation with star formation), whereas the magnetic fields are generated by mergers in GHs and by core sloshing in MHs.

Some galaxy clusters show diffuse radio emission in the form of peripheral relics (so far attributed to primary, shock-(re)accelerated electrons) or central halos. Analysing radio and X-ray data from the literature, we find new connections between halos and relics, such as a universal linear relation between their peak radio brightness and the gas column density. Our results indicate that halos, relics, and halo-relic bridges in a cluster, all arise from the same, homogeneous cosmic ray (CR) ion (CRI) distribution. We analytically derive the signature of synchrotron emission from secondary electrons and positrons (CREs) produced in hadronic CRI collisions, for an arbitrary magnetic field evolution. In our model, flat spectrum halos (both giant and minihalos) arise from steady-state magnetic fields, whereas relics and steep halos reflect recent or irregular magnetic growth. This naturally explains the properties of halos, relics, and the connections between them, without invoking particle (re)acceleration in weak shocks or turbulence. We find CRI energy densities in the range u_p~10^-[12.4,13.3]erg cm^-3, with a spectral index s_p=-2.20+-0.05, and identify a ~0.1 magnetic energy fraction in some halos and behind relics, as far as 2Mpc from the cluster's centre. The CRI homogeneity suggests strong CR diffusion, D(100GeV)>~10^32cm^2s^-1. The strong magnetisation imposes strict upper limits on >10GeV CRE (re)acceleration in weak shocks (efficiency <10^-4) and turbulence; indeed, each weak shock slightly lowers the energy fraction of flat CRs.

Cold fronts (CFs) - density and temperature plasma discontinuities - are ubiquitous in cool cores of galaxy clusters, where they appear as X-ray brightness edges in the intracluster medium, nearly concentric with the cluster center. We analyze the thermodynamic profiles deprojected across core CFs found in the literature. While the pressure appears continuous across these CFs, we find that all of them require significant centripetal acceleration beneath the front. This is naturally explained by a tangential, nearly sonic bulk flow just below the CF, and a tangential shear flow involving a fair fraction of the plasma beneath the front. Such shear should generate near-equipartition magnetic fields on scales ~<50 pc from the front, and could magnetize the entire core. Such fields would explain the apparent stability of cool-core CFs and the recently reported CF-radio minihalo association.

Cold fronts (CFs) are found in most galaxy clusters, as well as in some galaxies and groups of galaxies. We propose that some CFs are relics of collisions between trailing shocks. Such a collision typically results in a spherical, factor ~1.4-2.7 density/temperature discontinuity. These CFs may be found as far as the virial shock, unlike in other CF formation models. As a demonstration of this effect, we use one dimensional simulations where halo reverberations involving periodic collisions between the virial shock and outgoing secondary shocks exist. These collisions yield a distinctive, concentric geometric sequence of CFs which trace the expansion of the virial shock.

Cold fronts (CFs) are found in most galaxy clusters, as well as in some galaxies and groups of galaxies. We propose that some CFs are relics of merging between two shocks propagating in the same direction. Such shock mergers typically result in a quasi-spherical, factor ~1.4-2.7 discontinuity in density and in temperature. These CFs may be found as far out as the virial shock, unlike what is expected in other CF formation models. As a demonstration of this effect, we use one dimensional simulations of clusters and show that shock induced cold fronts form when perturbations such as explosions or mergers occur near the cluster's centre. Perturbations at a cluster's core induce periodic merging between the virial shock and outgoing secondary shocks. These collisions yield a distinctive, concentric, geometric sequence of CFs which trace the expansion of the virial shock.

We analyze the distribution of stars of arbitrary mass function xi(m) around a massive black hole (MBH). Unless xi is strongly dominated by light stars, the steady-state distribution function approaches a power-law in specific energy x=-E/(m*sigma^2)<x_max with index p=m/4M_0, where E is the energy, sigma is the typical velocity dispersion of unbound stars, and M_0 is the mass averaged over m*xi*x_{max}^p. For light-dominated xi, p can grow as large as 3/2 - much steeper than previously thought. A simple prescription for the stellar density profile around MBHs is provided. We illustrate our results by applying them to stars around the MBH in the Milky Way.

We study relativistic unmagnetized collisionless shocks using unprecedentedly large particle-in-cell simulations of two-dimensional pair plasma. High energy particles accelerated by the shock are found to drive magnetic field evolution on a timescale >10^4 plasma times. Progressively stronger magnetic fields are generated on larger scales in a growing region around the shock. Shock-generated magnetic fields and accelerated particles carry >1% and >10% of the downstream energy flux, respectively. Our results suggest limits on the magnetization of relativistic astrophysical flows.

We calculate analytically the transmission and reflection amplitudes for waves incident on a rotating black hole in d=4, analytically continued to asymptotically large, nearly imaginary frequency. These amplitudes determine the asymptotic resonant frequencies of the black hole, including quasinormal modes, total-transmission modes and total-reflection modes. We identify these modes with semiclassical bound states of a one-dimensional Schrodinger equation, localized along contours in the complexified r-plane which connect turning points of corresponding null geodesics. Each family of modes has a characteristic temperature and chemical potential. The relations between them provide hints about the microscopic description of the black hole in this asymptotic regime.

A Bohr-Sommerfeld equation is derived for the highly-damped quasinormal mode frequencies omega(n>>1) of rotating black holes. It may be written as 2\int_C(p_r+ip_0)dr=(n+1/2)h, where p_r is the canonical momentum conjugate to the radial coordinate r along null geodesics of energy hbar*omega and angular momentum hbar*m, p_0=O(omega^0), and the contour C connects two complex turning points of p_r. The solutions are omega(n) = - m*omega_0 - i(phi + n*delta), where {omega_0,delta}>0 are functions of the black-hole parameters alone. Some physical implications are discussed.

Particle acceleration in relativistic shocks is studied analytically in the test-particle, small-angle scattering limit, for an arbitrary velocity-angle diffusion function D. Accurate analytic expressions for the spectral index s are derived using few (2-6) low-order moments of the shock-frame angular distribution. For isotropic diffusion, previous results are reproduced and justified. For anisotropic diffusion, s is shown to be sensitive to D, particularly downstream and at certain angles, and a wide range of s values is attainable. The analysis, confirmed numerically, can be used to test collisionless shock models and to observationally constrain D. For example, strongly forward- or backward-enhanced diffusion downstream is ruled out by GRB afterglow observations.

Observations of gamma-ray burst afterglows suggest that the correlation length of magnetic field fluctuations downstream of relativistic non-magnetized collisionless shocks grows with distance from the shock to scales much larger than the plasma skin depth. We argue that this indicates that the plasma properties are described by a self-similar solution, and derive constraints on the scaling properties of the solution. For example, we find that the scaling of the characteristic magnetic field amplitude with distance from the shock is B \propto D^{s_B} with -1<s_B<=0, that the spectrum of accelerated particles is dn/dE \propto E^{-2/(s_B+1)}, and that the scaling of the magnetic correlation function is <B_i(x)B_j(x+dx)> \propto x^{2s_B} (for x>>D). We show that the plasma may be approximated as a combination of two self-similar components: a kinetic component of energetic particles and an MHD-like component representing "thermal" particles. We argue that the latter may be considered as infinitely conducting, in which case s_B=0 and the scalings are completely determined (e.g. dn/dE \propto E^{-2} and B \propto D^0). Similar claims apply to non- relativistic shocks such as in supernova remnants, if the upstream magnetic field can be neglected. Self-similarity has important implications for any model of particle acceleration and/or field generation. For example, we show that the diffusion function in the angle \mu of momentum p in diffusive shock acceleration models must satisfy D_{\mu\mu}(p,D) = D^{-1}D'_{\mu\mu}(p/D), and that a previously suggested model for the generation of large scale magnetic fields through a hierarchical merger of current-filaments should be generalized. A numerical experiment testing our analysis is outlined (Abridged).

We analytically study diffusive particle acceleration in relativistic, collisionless shocks. We find a simple relation between the spectral index s and the anisotropy of the momentum distribution along the shock front. Based on this relation, we obtain s = (3beta_u - 2beta_u*beta_d^2 + beta_d^3) / (beta_u - beta_d) for isotropic diffusion, where beta_u (beta_d) is the upstream (downstream) fluid velocity normalized to the speed of light. This result is in agreement with previous numerical determinations of s for all (beta_u,beta_d), and yields s=38/9 in the ultra-relativistic limit. The spectrum-anisotropy connection is useful for testing numerical studies and for constraining non-isotropic diffusion results. It implies that the spectrum is highly sensitive to the form of the diffusion function for particles travelling along the shock front.

Strong intergalactic shocks are a natural consequence of structure formation in the universe. They are expected to deposit large fractions of their energy in relativistic electrons (xi_e~0.05 according to SNR observations) and magnetic fields (xi_B~0.01 according to cluster halo observations). We calculate the synchrotron emission from such shocks using an analytical model, calibrated with a hydrodynamical LCDM simulation. The resulting signal composes a large fraction of the extragalactic radio background (ERB) below 500 MHz. The associated angular fluctuations dominate the sky for frequencies nu<10 GHz and angular scales arcmin-deg (after a modest removal of point sources), provided that xi_e*xi_B>3*10^-4. The fluctuating signal is most pronounced for nu<500 MHz, dominating the sky even for xi_e*xi_B=5*10^-5. The signal will be easily observable by next generation radio telescopes such as LOFAR and SKA, and is marginally observable with present telescopes. It may be identified using cross-correlations with tracers of large scale structure, possibly even in existing <10 GHz CMB anisotropy maps and high resolution ~1 GHz radio surveys. Detection of the signal will provide the first identification of intergalactic shocks and of the WHIM, and gauge the unknown intergalactic magnetic field. We show that existing observations of the diffuse <500 MHz radio background are well fit by a simple, double-disk Galactic model, precluding a direct identification of the diffuse ERB. Modelling the frequency-dependent anisotropy pattern observed at very low (1-10 MHz) frequencies can disentangle the distributions of Galactic cosmic-rays, ionized gas and magnetic fields. Space missions such as ALFA will thus provide important insight into the structure and composition of our Galaxy (abridged).

Strong intergalactic shocks are a natural consequence of structure formation in the universe. These shocks are expected to deposit large fractions of their energy in relativistic electrons (xi_e~0.05 of the thermal energy according to supernova remnant observations) and magnetic fields (xi_B~0.01 according to cluster halo observations). We discuss the expected synchrotron emission from such shocks, and the observational consequences for next generation radio telescopes such as the Square Kilometer Array. We present an analytical model, calibrated and verified based on a hydrodynamical LCDM simulation. The resulting signal composes a large fraction (up to a few 10%) of the extragalactic radio background below 500 MHz. The associated angular fluctuations, e.g. delta T_l>260(xi_e*xi_B/5*10^-4)(nu/100 MHz)^-3 K for multipoles 400<l<2000, dominate the radio sky for frequencies<10 GHz and angular scales 1 arcmin-1 deg (after a modest removal of point sources), provided that xi_e*xi_B>3*10^-4. The fluctuating signal is most pronounced for nu<500 MHz, dominating the sky there even for xi_e*xi_B=5*10^-5. We find that the signal is easily observable by next generation telescopes such as the SKA, and marginally observable with present telescopes. The signal may be identified using cross-correlations with tracers of large-scale structure (e.g. gamma-ray emission from intergalactic shocks). Detection of the signal will provide the first identification of intergalactic shocks and of the warm-hot intergalactic medium (believed to contain most of the baryons in the low redshift universe), and gauge the unknown intergalactic magnetic field.

Measurements of the diffuse extragalactic gamma-ray background (EGRB) are complicated by a strong Galactic foreground. Estimates of the EGRB flux and spectrum, obtained by modeling the Galactic emission, have produced a variety of (sometimes conflicting) results. The latest analysis of the EGRET data found an isotropic flux I_x=1.45+-0.05 above 100 MeV, in units of 10^-5 s^-1 cm^-2 sr^-1. We analyze the EGRET data in search for robust constraints on the EGRB flux, finding the gamma-ray sky strongly dominated by Galactic foreground even at high latitudes, with no conclusive evidence for an additional isotropic component. The gamma-ray intensity measured towards the Galactic poles is similar to or lower than previous estimates of I_x. The high latitude profile of the gamma-ray data is disk-like for 40<|b[deg]|<70, and even steeper for |b|>70; overall it exhibits strong Galactic features and is well fit by a simple Galactic model. Based on the |b|>40 data we find that I_x<0.5 at a 99% confidence level, with evidence for a much lower flux. We show that correlations with Galactic tracers, previously used to identify the Galactic foreground and estimate I_x, are not satisfactory; the results depend on the tracers used and on the part of the sky examined, because the Galactic emission is not linear in the Galactic tracers and exhibits spectral variations across the sky. The low EGRB flux favored by our analysis places stringent limits on extragalactic scenarios involving gamma-ray emission, such as radiation from blazars, intergalactic shocks and production of ultra-high energy cosmic rays and neutrinos. We suggest methods by which future gamma-ray missions such as GLAST and AGILE could indirectly identify the EGRB.

Structure formation in the intergalactic medium (IGM) produces large-scale, collisionless shock waves, where electrons can accelerate to highly relativistic energies. Such electrons can Compton scatter cosmic microwave background photons up to gamma-ray energies. We study the radiation emitted in this process using a hydrodynamic cosmological simulation of a LCDM universe. This radiation, extending beyond TeV energies, has roughly constant energy flux per decade in photon energy, in agreement with the predictions of Loeb & Waxman (2000). Assuming that a fraction xi_e=0.05 of the shock thermal energy is transferred to relativistic electrons, as inferred from collisionless non-relativistic shocks in the interstellar medium, we find that the radiation energy flux, e^2(dJ/de)~ 50-160 eV cm^-2 s^-1 sr^-1, constitutes ~10% of the extragalactic gamma-ray background. The associated point-sources are too faint to account for the ~60 unidentified EGRET gamma-ray sources, but GLAST should resolve several sources associated with large-scale IGM structures for xi_e~0.03, and many more sources for larger xi_e. The intergalactic origin of the radiation can be verified through a cross-correlation with, e.g., the galaxy distribution that traces the same structure. Its shock-origin may be tested by a cross-correlation with radio synchrotron radiation, emitted as the same electrons gyrate in post-shock magnetic fields. We predict that GLAST and Cherenkov telescopes such as MAGIC, VERITAS and HESS should resolve gamma-rays from nearby (redshifts z < 0.01) rich galaxy clusters, perhaps in the form of a \~5-10 Mpc diameter ring-like emission tracing the cluster accretion shock, with luminous peaks at its intersections with galaxy filaments detectable even at z~0.025.