Showing 1–50 of 61 results for author: Keshet, U

Coma cluster $γ$-ray and radio emission is consistent with a secondary electron origin for the radio halo

Authors: Doron Kushnir, Uri Keshet, Eli Waxman

Abstract: Observations of diffuse, non-thermal radio emission spanning several megaparsecs have been documented in over 100 galaxy clusters. This emission, classified as giant radio halos (GHs), mini halos, and radio relics based mainly on their location and morphology, is interpreted as synchrotron radiation and implies the presence of relativistic electrons and magnetic fields in the intra-cluster medium… ▽ More Observations of diffuse, non-thermal radio emission spanning several megaparsecs have been documented in over 100 galaxy clusters. This emission, classified as giant radio halos (GHs), mini halos, and radio relics based mainly on their location and morphology, is interpreted as synchrotron radiation and implies the presence of relativistic electrons and magnetic fields in the intra-cluster medium (ICM). GHs were initially thought to be generated by secondary electrons resulting from inelastic $p+p\rightarrow X+π^{\pm}$ collisions. However, recent literature has leaned towards primary-electron turbulent (re)acceleration models, partly due to claimed upper limits on the $γ$-ray emission from $π^0$ decay. We demonstrate that the observed GH and $γ$-ray flux in the Coma cluster are consistent with a secondary origin for the GH across a broad range of magnetic field values. Although the constraints on magnetic field configuration are not stringent, they align well with previous estimates for Coma. Within this magnetic field range, the energy density of cosmic-ray protons (CRp) constitutes a few percent to tens of percent of the ICM energy density, as predicted and observed for a sample of radio-emitting galaxy clusters. Notably, we detect a rise in the ratio of CRp to ICM energy densities towards the outer regions of the cluster. This phenomenon was anticipated to arise from either adiabatic compression of CRp accelerated by accretion shocks or, more likely, from strong CRp diffusion. △ Less

Submitted 19 April, 2024; originally announced April 2024.

Comments: 7 pages, 3 figures

Galaxy cluster virial-shock sources in eROSITA catalogs

Authors: Gideon Ilani, Kuan-Chou Hou, Gil Nadler, Uri Keshet

Abstract: Following the recent identification of discrete ROSAT and radio sources associated with the virial shocks of MCXC clusters and groups, we examine if the early eROSITA-DE data release (EDR) shows virial-shock X-ray sources within its $140$ deg$^2$ field. EDR catalog sources are stacked and radially binned around EDR catalog clusters and groups. The properties of the excess virial-shock sources are… ▽ More Following the recent identification of discrete ROSAT and radio sources associated with the virial shocks of MCXC clusters and groups, we examine if the early eROSITA-DE data release (EDR) shows virial-shock X-ray sources within its $140$ deg$^2$ field. EDR catalog sources are stacked and radially binned around EDR catalog clusters and groups. The properties of the excess virial-shock sources are inferred statistically by comparing the virial-shock region to the field. An excess of X-ray sources is found narrowly localized at the $2.0<r/R_{500}<2.25$ normalized radii, just inside the anticipated virial shocks, of the resolved 532 clusters, for samples of both extended ($3σ$ for 534 sources) or bright ($3.5σ$ for 5820 sources; $4σ$ excluding the low cluster-mass quartile) sources. The excess sources are on average extended ($\sim 100$ kpc), luminous ($L_X\simeq 10^{43-44}$ erg s$^{-1}$), and hot ($\sim$keV), consistent with infalling gaseous halos crossing the virial shock. The results agree with the stacked ROSAT-MCXC signal, showing the higher $L_X$ anticipated at EDR redshifts and a possible dependence upon host mass. Localized virial-shock spikes in the distributions of discrete radio, X-ray, and probably also $γ$-ray sources are new powerful probes of accretion from the cosmic web, with strong constraints anticipated with future all-sky catalogs such as by eROSITA. △ Less

Submitted 27 February, 2024; originally announced February 2024.

Comments: 4 pages, 5 figures + Appendices

Journal ref: A&A 686, L16 (2024)

Excess cataloged X-ray and radio sources at galaxy-cluster virial shocks

Authors: Gideon Ilani, Kuan-Chou Hou, Uri Keshet

Abstract: We detect a highly significant excess of X-ray (2RXS) and radio (NVSS, GMRT, VLSSr) catalog sources when stacked around MCXC galaxy clusters and groups, narrowly confined within $\lesssim100\mathrm{\,kpc}$ of the $\sim2.4 R_{500}$ virial shock radius (inferred from previous continuum stacking), with similar X-ray ($\sim4σ$ for $443$ clusters) and radio ($\sim4σ$ for $485$ clusters) characteristics… ▽ More We detect a highly significant excess of X-ray (2RXS) and radio (NVSS, GMRT, VLSSr) catalog sources when stacked around MCXC galaxy clusters and groups, narrowly confined within $\lesssim100\mathrm{\,kpc}$ of the $\sim2.4 R_{500}$ virial shock radius (inferred from previous continuum stacking), with similar X-ray ($\sim4σ$ for $443$ clusters) and radio ($\sim4σ$ for $485$ clusters) characteristics ($>5σ$ joint). The excess sources show $10-100$ kpc scales, $L_X(0.1-2.4\mbox{ keV})\simeq10^{42-43}\mathrm{\,erg\,s^{-1}}$ or $νL_ν(ν=1.4\mbox{ GHz}) \simeq 10^{40-41}\mathrm{\,erg\,s^{-1}}$ luminosities, and a preferentially radial radio-polarization. The narrow localization and properties of the excess identify these sources not as AGN, often invoked speculatively for excess X-ray sources at cluster outskirts, but rather as infalling gaseous clumps interacting with the virial shock, probably galactic halos and possibly outflow remnants. The local excess of such discrete, radio-to-$γ$-ray sources around an object can probe its virial shock also at high redshifts and sub-cluster scales. △ Less

Submitted 19 September, 2024; v1 submitted 26 February, 2024; originally announced February 2024.

Comments: Source tables added in Appendix H; to appear in JCAP

Detection of polarized Fermi-bubble synchrotron and dust emission

Authors: Uri Keshet

Abstract: The elusive polarized microwave signal from the Fermi bubbles is disentangled from the more extended polarized lobes, which similarly emanate from the Galactic plane but stretch farther west of the bubbles. The ~20% synchrotron polarization reveals magnetic fields preferentially parallel to the bubble edges, as expected downstream of a strong shock. The ~20% polarization of thermal dust emission i… ▽ More The elusive polarized microwave signal from the Fermi bubbles is disentangled from the more extended polarized lobes, which similarly emanate from the Galactic plane but stretch farther west of the bubbles. The ~20% synchrotron polarization reveals magnetic fields preferentially parallel to the bubble edges, as expected downstream of a strong shock. The ~20% polarization of thermal dust emission is similarly oriented, constraining grain alignment in an extreme environment. We argue that the larger lobes arise from an older Galactic-center, likely supermassive black-hole, outburst. △ Less

Submitted 1 January, 2024; originally announced January 2024.

Comments: 6 pages with 3 figures

Fermi-bubble bulk and edge analysis reveals dust, cooling breaks, and cosmic-ray diffusion, facilitating a self-consistent model

Authors: Uri Keshet, Ilya Gurwich, Assaf Lavi, Dina Avitan, Teodor Linnik

Abstract: The full, radio to $γ$-ray spectrum of the Fermi bubbles is shown to be consistent with standard strong-shock electron acceleration at the bubble edge, without the unnatural energy cutoffs and unrealistic electron cooling of previous studies, if the ambient interstellar radiation is strong; the $γ$-ray cooling break should then have a microwave counterpart, undetected until now. Indeed, a broadban… ▽ More The full, radio to $γ$-ray spectrum of the Fermi bubbles is shown to be consistent with standard strong-shock electron acceleration at the bubble edge, without the unnatural energy cutoffs and unrealistic electron cooling of previous studies, if the ambient interstellar radiation is strong; the $γ$-ray cooling break should then have a microwave counterpart, undetected until now. Indeed, a broadband bubble-edge analysis uncovers a pronounced downstream dust component, which masked the anticipated $\sim35$ GHz spectral break, and the same overall radio softening consistent with Kraichnan diffusion previously reported in $γ$-rays. △ Less

Submitted 26 March, 2024; v1 submitted 14 November, 2023; originally announced November 2023.

Comments: To appear in ApJ

Radio halos and relics from extended cosmic-ray ion distributions with strong diffusion in galaxy clusters

Authors: Uri Keshet

Abstract: A joint hadronic model is shown to quantitatively explain the observations of diffuse radio emission from galaxy clusters in the form of minihalos, giant halos, relics, and their hybrid, transitional stages. Cosmic-ray diffusion of order $D\sim 10^{31\text{--}32}\text{ cm}^2\text{ s}^{-1}$, inferred independently from relic energies, the spatial variability of giant-halo spectra, and the spectral… ▽ More A joint hadronic model is shown to quantitatively explain the observations of diffuse radio emission from galaxy clusters in the form of minihalos, giant halos, relics, and their hybrid, transitional stages. Cosmic-ray diffusion of order $D\sim 10^{31\text{--}32}\text{ cm}^2\text{ s}^{-1}$, inferred independently from relic energies, the spatial variability of giant-halo spectra, and the spectral evolution of relics, reproduces the observed spatio-spectral distributions, explains the recently discovered mega-halos as enhanced peripheral magnetization, and quenches electron (re)acceleration by weak shocks or turbulence. For instance, the hard-to-soft evolution along secondary-electron diffusion explains both the soft spectra in most halo peripheries and relic downstreams, and the hard spectra in most halo centres and relic edges, where the photon index can reach $α\simeq -0.5$ regardless of the Mach number $\mathcal{M}$ of the coincident shock. Such spatio-spectral modeling, recent $γ$-ray observations, and additional accumulated evidence are thus shown to support a previous claim (Keshet 2010) that the seamless transitions among minihalos, giant halos, and relics, their similar energetics, integrated spectra, and delineating discontinuities, the inconsistent $\mathcal{M}$ inferred from radio vs. X-rays in leptonic models, and additional observations, all indicate that these diffuse radio phenomena are manifestations of the same cosmic-ray ion population, with no need to invoke less natural alternatives. △ Less

Submitted 15 November, 2023; v1 submitted 14 March, 2023; originally announced March 2023.

Comments: To appear in MNRAS; Appendix E added

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 527, Issue 1, pp.1194-1215 (2024)

Synchrotron emission from virial shocks around stacked OVRO-LWA galaxy clusters

Authors: Kuan-Chou Hou, Gregg Hallinan, Uri Keshet

Abstract: Galaxy clusters accrete mass through large scale, strong, structure-formation shocks. Such a virial shock is thought to deposit fractions $ξ_e$ and $ξ_B$ of the thermal energy in cosmic-ray electrons (CREs) and magnetic fields, respectively, thus generating a leptonic virial ring. However, the expected synchrotron signal was not convincingly established until now. We stack low-frequency radio data… ▽ More Galaxy clusters accrete mass through large scale, strong, structure-formation shocks. Such a virial shock is thought to deposit fractions $ξ_e$ and $ξ_B$ of the thermal energy in cosmic-ray electrons (CREs) and magnetic fields, respectively, thus generating a leptonic virial ring. However, the expected synchrotron signal was not convincingly established until now. We stack low-frequency radio data from the OVRO-LWA around the 44 most massive, high latitude, extended MCXC clusters, enhancing the ring sensitivity by rescaling clusters to their characteristic, $R_{500}$ radii. Both high (73 MHz) and co-added low ($36\text{--}68\text{ MHz}$) frequency channels separately indicate a significant ($4\text{--}5σ$) excess peaked at $(2.4 \text{--} 2.6) R_{500}$, coincident with a previously stacked Fermi $γ$-ray signal interpreted as inverse-Compton emission from virial-shock CREs. The stacked radio signal is well fit (TS-test: $4$--$6σ$ at high frequency, $4$--$8σ$ at low frequencies, and $8$--$10σ$ joint) by virial-shock synchrotron emission from the more massive clusters, with $\dot{m}ξ_eξ_B\simeq (1\text{--}4)\times 10^{-4}$, where $\dot{m}\equiv \dot{M}/(MH)$ is the dimensionless accretion rate for a cluster of mass $M$ and a Hubble constant $H$. The inferred CRE spectral index is flat, $p \simeq 2.0 \pm 0.2$, consistent with acceleration in a strong shock. Assuming equipartition or using $\dot{m}ξ_e\sim0.6\%$ inferred from the Fermi signal yields $ξ_B\simeq (2\text{--}9)\%$, corresponding to $B \simeq (0.1\text{--}0.3)~μ\text{G}$ magnetic fields downstream of typical virial shocks. Preliminary evidence suggests non-spherical shocks, with factor $2$--$3$ elongations. △ Less

Submitted 13 March, 2023; v1 submitted 17 October, 2022; originally announced October 2022.

Comments: Revised version (added preliminary evidence for elongated shocks) to appear in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 521, Issue 4, pp.5786-5809 (2023)

Quasi-spiral solution to the mixed intracluster medium and the universal entropy profile of galaxy clusters

Authors: Uri Keshet, Itay Raveh, Arka Ghosh

Abstract: Well-resolved galaxy clusters often show a large-scale quasi-spiral structure in deprojected density $ρ$ and temperature $T$ fields, delineated by a tangential discontinuity known as a cold front, superimposed on a universal radial entropy profile with a linear $K(r)\propto Tρ^{-2/3}\propto r$ adiabat. We show that a spiral structure provides a natural quasi-stationary solution for the mixed intra… ▽ More Well-resolved galaxy clusters often show a large-scale quasi-spiral structure in deprojected density $ρ$ and temperature $T$ fields, delineated by a tangential discontinuity known as a cold front, superimposed on a universal radial entropy profile with a linear $K(r)\propto Tρ^{-2/3}\propto r$ adiabat. We show that a spiral structure provides a natural quasi-stationary solution for the mixed intracluster medium (ICM), introducing a modest pressure spiral that confines the locally buoyant or heavy plasma phases. The solution persists in the presence of uniform or differential rotation, and can accommodate both an inflow and an outflow. Hydrodynamic adiabatic simulations with perturbations that deposit angular momentum and mix the plasma thus asymptote to a self-similar spiral structure. We find similar spirals in Eulerian and Lagrangian simulations of 2D and 3D, merger and offset, clusters. The discontinuity surface is given in spherical coordinates $\{r,θ,φ\}$ by $φ(r,θ)\propto Φ(r)$, where $Φ$ is the gravitational potential, combining a trailing spiral in the equatorial ($θ=π/2$) plane and semicircles perpendicular to the plane, in resemblance of a snail shell. A local convective instability can develop between spiral windings, driving a modified global instability in sublinear $K(r)$ regions; evolved spirals thus imprint the observed $K\propto r$ onto the ICM even after they dissipate. The spiral structure brings hot and cold phases to close proximity, suggesting that the observed fast outflows could sustain the structure even in the presence of radiative cooling. △ Less

Submitted 22 April, 2023; v1 submitted 19 September, 2022; originally announced September 2022.

Comments: Revised version to appear in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 522, Issue 4, pp.4991-5012 (2023)

Nonlinear electromagnetic-wave interactions in pair plasma: (I) Non-relativistic regime

Authors: Arka Ghosh, Daniel Kagan, Uri Keshet, Yuri Lyubarsky

Abstract: This paper is the first in a series devoted to the numerical study of nonlinear interactions of electromagnetic waves with plasma. We start with non-magnetized pair plasmas, where the primary processes are induced (Compton) scattering and the filamentation instability. In this paper, we consider waves in which electron oscillations are non-relativistic. Here, the numerical results can be compared… ▽ More This paper is the first in a series devoted to the numerical study of nonlinear interactions of electromagnetic waves with plasma. We start with non-magnetized pair plasmas, where the primary processes are induced (Compton) scattering and the filamentation instability. In this paper, we consider waves in which electron oscillations are non-relativistic. Here, the numerical results can be compared to analytical theory, facilitating the development of appropriate numerical tools and framework. We distill the analytic theory, reconciling plasma and radiative transfer pictures of induced scattering and developing in detail the kinetic theory of modulation/filamentation instability. We carry out homogeneous numerical simulations using the particle-in-cell codes EPOCH and Tristan-MP, for both monochromatic waves and wave packets. We show that simulations of both processes are consistent with theoretical predictions, setting the stage for analyzing the highly nonlinear regime. △ Less

Submitted 31 October, 2021; originally announced November 2021.

Comments: 18 pages, 15 figures, Comments welcome

Journal ref: The Astrophysical Journal, Volume 930, Issue 2, id.106, 18 pp. (2022)

Fermi bubbles: the collimated outburst needed to explain forward-shock edges

Authors: Santanu Mondal, Uri Keshet, Kartick C. Sarkar, Ilya Gurwich

Abstract: The bipolar, nonthermal, high-latitude lobes known as the Fermi bubbles (FBs) are thought to originate from a massive energy release near the Galactic centre (GC). We constrain the FB engine and the circumgalactic medium (CGM) by analytically and numerically modeling the FB edges as strong forward shocks, as inferred from recent observations. A non-directed energy release produces shocks too spher… ▽ More The bipolar, nonthermal, high-latitude lobes known as the Fermi bubbles (FBs) are thought to originate from a massive energy release near the Galactic centre (GC). We constrain the FB engine and the circumgalactic medium (CGM) by analytically and numerically modeling the FB edges as strong forward shocks, as inferred from recent observations. A non-directed energy release produces shocks too spherical to account for observations even for a maximally massive Galactic disc, critical CGM rotation, or injection effectively offset from the GC. In contrast, collimated injection nearly perpendicular to the disc can account for observations in both ballistic (free expansion) and slowdown regimes, as we show using a simple stratified evolution model verified by hydrodynamic simulations. FBs still in their ballistic regime require injection (at $z\simeq100$ pc heights in our model) with a half-opening angle $θ\simeq4^\circ$, a normalized velocity $β_{-2}\equiv v/(0.01c)\gtrsim 0.4$, and an energy $E\gtrsim2β_{-2}^2\times 10^{55}$ erg, launched $\mathbb{T}\simeq 3.3β_{-2}^{-1}$ Myr ago, showing a distinctive low-pressure region behind the bubble head. Slowing-down (mass accumulated) FBs require a faster injection, a thinner jet, a smaller $E/(β_{-2}θ)^{2}$, and a comparable $\mathbb{T}$, and follow a ballistic stage that must reach a height $z_{s}\gtrsim 5$ kpc. △ Less

Submitted 8 September, 2021; originally announced September 2021.

Comments: 18 pages, 17 figures, 3 tables; comments welcome

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 514, Issue 2, pp.2581-2598 (2022)

Protruding bullet heads indicating dark matter pull

Authors: Uri Keshet, Itay Raveh, Yossi Naor

Abstract: A clump moving through the intracluster medium of a galaxy cluster can drive a bow shock trailed by a bullet-like core. In some cases, such as in the prototypical Bullet cluster, X-rays show a gas bullet with a protruding head and pronounced shoulders. We point out that these features, while difficult to explain without dark matter (DM), naturally arise as the head of the slowed-down gas is gravit… ▽ More A clump moving through the intracluster medium of a galaxy cluster can drive a bow shock trailed by a bullet-like core. In some cases, such as in the prototypical Bullet cluster, X-rays show a gas bullet with a protruding head and pronounced shoulders. We point out that these features, while difficult to explain without dark matter (DM), naturally arise as the head of the slowed-down gas is gravitationally pulled forward toward its unhindered DM counterpart. X-ray imaging thus provides a unique, robust probe of the offset, collisionless DM, even without gravitational lensing or other auxiliary data. Numerical simulations and a toy model suggest that the effect is common in major mergers, is often associated with a small bullet-head radius of curvature, and may lead to distinct bullet morphologies, consistent with observations. △ Less

Submitted 15 January, 2021; originally announced January 2021.

Comments: 9 pages, 3 figures, comments welcome

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 508, Issue 3, pp.3455-3462 (2021)

Maximally-hard spectra from diffusive shock-acceleration in pulsar-wind nebulae

Authors: Ofir Arad, Assaf Lavi, Uri Keshet

Abstract: The processes leading to the exceptionally hard radio spectra of pulsar-wind nebulae (PWNe) are not yet understood. Radio photon spectral indices among $29$ PWNe from the literature show an approximately normal, $α=0.2\pm0.2$ distribution. We present $\sim 3σ$ evidence for a distinct sub-population of PWNe, with a hard spectrum $α=0.01\pm0.06$ near the termination shock and significantly softer el… ▽ More The processes leading to the exceptionally hard radio spectra of pulsar-wind nebulae (PWNe) are not yet understood. Radio photon spectral indices among $29$ PWNe from the literature show an approximately normal, $α=0.2\pm0.2$ distribution. We present $\sim 3σ$ evidence for a distinct sub-population of PWNe, with a hard spectrum $α=0.01\pm0.06$ near the termination shock and significantly softer elsewhere, possibly due to a recent evacuation of the shock surroundings. Such spectra, especially in the hard sub-population, suggest a Fermi process, such as diffusive shock acceleration, at its extreme, $α=0$ limit. We show that this limit is approached for sufficiently anisotropic small-angle scattering, enhanced on either side of the shock for particles approaching the shock front. In the upstream, the spectral hardening is mostly associated with an enhanced energy gain, possibly driven by the same beamed particles crossing the shock. Downstream, the main effect is a diminished escape probability, but this lowers the acceleration efficiency to $\lesssim25\%$ for $α=0.3$ and $\lesssim1\%$ for $α=0.03$. △ Less

Submitted 1 December, 2020; originally announced December 2020.

Comments: 17 pages, 13 figures, 3 tables; comments welcome

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 504, Issue 4, pp.4952-4967 (2021)

Compressible potential flows around round bodies: Janzen-Rayleigh expansion inferences

Authors: Idan S. Wallerstein, Uri Keshet

Abstract: The subsonic, compressible, potential flow around a hypersphere can be derived using the Janzen-Rayleigh expansion (JRE) of the flow potential in even powers of the incident Mach number $\mathcal{M}_\infty$. JREs were carried out with terms polynomial in the inverse radius $r^{-1}$ to high orders in two dimensions (2D), but were limited to order $\mathcal{M}_\infty^4$ in three dimensions (3D). We… ▽ More The subsonic, compressible, potential flow around a hypersphere can be derived using the Janzen-Rayleigh expansion (JRE) of the flow potential in even powers of the incident Mach number $\mathcal{M}_\infty$. JREs were carried out with terms polynomial in the inverse radius $r^{-1}$ to high orders in two dimensions (2D), but were limited to order $\mathcal{M}_\infty^4$ in three dimensions (3D). We derive general JRE formulae to arbitrary order, adiabatic index, and dimension. We find that powers of $\ln(r)$ can creep into the expansion, and are essential in 3D beyond order $\mathcal{M}_\infty^4$. Such terms are apparently absent in the 2D disk, as we confirm up to order $\mathcal{M}_\infty^{100}$, although they do show in other dimensions (e.g. at order $\mathcal{M}_\infty^2$ in 4D) and in non-circular 2D bodies. This suggests that the disk, which was extensively used to study basic flow properties, has additional symmetry. Our results are used to improve the hodograph-based approximation for the flow in front of a sphere. The symmetry-axis velocity profiles of axisymmetric flows around different prolate spheroids are approximately related to each other by a simple, Mach-independent scaling. △ Less

Submitted 26 November, 2020; originally announced November 2020.

Comments: 24 pages, 6 figures, 5 tables, comments welcome

Journal ref: Journal of Fluid Mechanics, vol. 932, id. A6 (2022)

Diffusive shock acceleration in $N$ dimensions

Authors: Assaf Lavi, Ofir Arad, Yotam Nagar, Uri Keshet

Abstract: 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… ▽ More 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. △ Less

Submitted 25 February, 2020; originally announced February 2020.

Comments: 14 pages, 7 figures, comments welcome

Journal ref: The Astrophysical Journal, Volume 895, Issue 2, id.107 (2020)

Advection by large-scale spiral flows in galaxy clusters

Authors: Yossi Naor, Uri Keshet

Abstract: 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 t… ▽ More 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σ\text{-}4.2σ$) 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. △ Less

Submitted 5 February, 2020; originally announced February 2020.

Comments: 18 pages, 10 figures, 3 tables; comments welcome

Journal ref: The Astrophysical Journal, Volume 895, Issue 2, id.143 (2020)

Deprojecting galaxy-cluster cold fronts: evidence for bulk, magnetised spiral flows

Authors: Yossi Naor, Uri Keshet, Qian Wang, Ido Reiss

Abstract: 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 nontherma… ▽ More 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σ$) 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σ$) 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σ$ and $5.0σ$, 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σ$) with a $1.25^{+0.09}_{-0.08}$ metallicity drop across the CF. These findings suggest that GCs quite generally harbor extended spiral flows. △ Less

Submitted 3 February, 2020; originally announced February 2020.

Comments: 29 pages, 20 figures, 6 tables

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 495, Issue 4, pp.4392-4418 (2020)

Diffusive acceleration in relativistic shocks: particle feedback

Authors: Yotam Nagar, Uri Keshet

Abstract: 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 moti… ▽ More 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. △ Less

Submitted 22 October, 2019; originally announced October 2019.

Comments: 8 pages, 6 figures, comments welcome

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 501, Issue 1, pp.329-336 (2021)

Diffusive shock-acceleration: breakdown of spatial diffusion and isotropy

Authors: Uri Keshet, Ofir Arad, Yuri Lyubarski

Abstract: 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 c… ▽ More 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. △ Less

Submitted 17 October, 2019; originally announced October 2019.

Comments: 11 pages, 2 figures, comments welcome

Journal ref: The Astrophysical Journal, Volume 891, Issue 2, article id. 117 (2020)

Coincident SZ and $γ$-ray signals from cluster virial shocks

Authors: Uri Keshet, Ido Reiss, Guillaume Hurier

Abstract: 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… ▽ More 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 $γ$-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 $γ$-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. △ Less

Submitted 20 January, 2020; v1 submitted 4 January, 2018; originally announced January 2018.

Comments: Minor changes in presentation

Journal ref: The Astrophysical Journal, Volume 895, Issue 1, article id.72 (2020)

Detection of virial shocks in stacked Fermi-LAT clusters

Authors: Ido Reiss, Jonathan Mushkin, Uri Keshet

Abstract: 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… ▽ More 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 $γ$-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σ$ confidence level), probably dominated by active galactic nuclei, we identify ($5.9σ$) a bright, spectrally flat $γ$-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 $γ$-ray and radio backgrounds. △ Less

Submitted 18 December, 2017; originally announced December 2017.

Comments: Proceedings of the 7th International Fermi Symposium, Garmisch, Oct 15-20 2017

Journal ref: PoS(IFS2017)163

First detection of a virial shock with SZ data: implication for the mass accretion rate of Abell 2319

Authors: G. Hurier, R. Adam, U. Keshet

Abstract: 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σ$, using Planck thermal Sunyaev-Zel'dovich data. The sho… ▽ More 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σ$, 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$. △ Less

Submitted 15 December, 2017; originally announced December 2017.

Comments: 4 pages, 3 figures

X-ray to gamma-ray virial shock signal from the Coma cluster

Authors: Uri Keshet, Ido Reiss

Abstract: 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 $γ$-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σ$ LAT excess, and detect ($>5σ$) the expected signature in ROSAT bands R1 and… ▽ More 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 $γ$-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σ$ LAT excess, and detect ($>5σ$) 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. △ Less

Submitted 20 September, 2017; originally announced September 2017.

Comments: 14 pages, 13 figures, comments welcome

Journal ref: The Astrophysical Journal, Volume 869, Issue 1, article id. 53, 15 pp. (2018)

Analytic study of 1D diffusive relativistic shock acceleration

Authors: Uri Keshet

Abstract: 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[γ_d(1+β_d)\right]/\ln\left[γ_u(1+β_u)\right]$, where $β_u$ ($β_d)$ is the upstream (downstream) normalized velocity, and $γ$ is the respective Lor… ▽ More 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[γ_d(1+β_d)\right]/\ln\left[γ_u(1+β_u)\right]$, where $β_u$ ($β_d)$ is the upstream (downstream) normalized velocity, and $γ$ 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{ü}ttner-Synge equation of state) noticeably hardens with increasing $1<γ_u<57$, before logarithmically converging back to $p(γ_u\to\infty)=2$. The 1D spectrum is sensitive to drifts, but only in the downstream, and not in the ultra-relativistic limit. △ Less

Submitted 12 October, 2017; v1 submitted 19 July, 2017; originally announced July 2017.

Comments: Revised version to appear in JCAP

Journal ref: JCAP 10 (2017) 025

IceCube constraints on the Fermi Bubbles

Authors: Nimrod Sherf, Uri Keshet, Ilya Gurwich

Abstract: 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 $ξ_i<8\%$ upper limit on the cosmic-ray ion (CRI) acceleration efficiency, and an $η\equiv ξ_e/ξ_i \gtrsim0.006$ lower limit on the electron-to-ion ratio of accele… ▽ More 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 $ξ_i<8\%$ upper limit on the cosmic-ray ion (CRI) acceleration efficiency, and an $η\equiv ξ_e/ξ_i \gtrsim0.006$ lower limit on the electron-to-ion ratio of acceleration efficiencies (at the $2σ$ confidence level). For typical $ξ_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. △ Less

Submitted 22 September, 2017; v1 submitted 18 May, 2017; originally announced May 2017.

Comments: Revised version to appear in the Astrophysical Journal

Journal ref: The Astrophysical Journal, Volume 847, Issue 2, article id. 95, 4 pp. (2017)

Detection of virial shocks in stacked Fermi-LAT clusters

Authors: Ido Reiss, Uri Keshet

Abstract: 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 $γ$-ray virial ring. We stack Fermi-LAT data for the 112 most massive, high latitude, extended c… ▽ More 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 $γ$-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σ$ confidence level), probably dominated by AGN, we identify (at the $5.8σ$ confidence level) a bright, spectrally-flat $γ$-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 $γ$-ray and low-frequency radio backgrounds. △ Less

Submitted 7 October, 2018; v1 submitted 15 May, 2017; originally announced May 2017.

Comments: Revised version published in JCAP

Journal ref: JCAP10(2018)010

Fermi bubbles: high latitude X-ray supersonic shell

Authors: Uri Keshet, Ilya Gurwich

Abstract: The nature of the bipolar, $γ$-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.… ▽ More The nature of the bipolar, $γ$-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. △ Less

Submitted 17 April, 2017; originally announced April 2017.

Comments: 12 pages, comments welcome; high resolution version at http://physics.bgu.ac.il/~ukeshet/Online.html

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 480, Issue 1, p.223-235 (2018)

Compressible flow in front of an axisymmetric blunt object: analytic approximation and astrophysical implications

Authors: Uri Keshet, Yossi Naor

Abstract: 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 supersoni… ▽ More 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 $δ$ 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δ})$. △ Less

Submitted 18 November, 2016; originally announced November 2016.

Comments: 9 pages, 3 figures, ApJ accepted. arXiv admin note: text overlap with arXiv:1412.8406

Journal ref: The Astrophysical Journal, Volume 830, Issue 2, article id. 147, 9 pp. (2016)

Fermi bubble edges: spectrum and diffusion function

Authors: Uri Keshet, Ilya Gurwich

Abstract: Current measurements of the $γ$-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… ▽ More Current measurements of the $γ$-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 $η\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 $η$, indicate that the FB edge is a strong, Mach $\gtrsim5$, forward shock. △ Less

Submitted 15 April, 2017; v1 submitted 13 November, 2016; originally announced November 2016.

Comments: Revised version accepted by ApJ. High resolution version at http://physics.bgu.ac.il/~ukeshet/Online.html

Journal ref: The Astrophysical Journal, Volume 840, Issue 1, article id. 7, 12 pp. (2017)

Shock Revival in Core-Collapse Supernovae: A Phase-Diagram Analysis

Authors: Daniel Gabay, Shmuel Balberg, Uri Keshet

Abstract: 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 evolutio… ▽ More 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. △ Less

Submitted 14 August, 2015; originally announced August 2015.

Comments: 22 pages, 13 figures, submitted to ApJ. Comments welcome

Journal ref: The Astrophysical Journal, Volume 815, Issue 1, article id. 37, 20 pp. (2015)

Deep Chandra observation and numerical studies of the nearest cluster cold front in the sky

Authors: N. Werner, J. A. ZuHone, I. Zhuravleva, Y. Ichinohe, A. Simionescu, S. W. Allen, M. Markevitch, A. C. Fabian, U. Keshet, E. Roediger, M. Ruszkowski, J. S. Sanders

Abstract: 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… ▽ More 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. △ Less

Submitted 16 November, 2015; v1 submitted 21 June, 2015; originally announced June 2015.

Comments: Published in MNRAS 455, 846

Journal ref: 2016 MNRAS 455, 846

Self-similar ultra-relativistic jetted blast wave

Authors: Uri Keshet, Dani Kogan

Abstract: 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 $γ$ exceeds $\sim1/2$ of its maximal, nose value; a geometrically self-similar, expanding envelope slightly n… ▽ More 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 $γ$ 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 $γ$ (separating the core and the envelope), a sign reversal in $u$, and a minimal $γ$, 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 $γ$-ray burst. △ Less

Submitted 7 December, 2015; v1 submitted 15 June, 2015; originally announced June 2015.

Comments: Revised version (discussion added) accepted by ApJ

Journal ref: The Astrophysical Journal, Volume 815, Issue 2, article id. 100, 15 pp. (2015)

Relativistic Shocks: Particle Acceleration and Magnetization

Authors: Lorenzo Sironi, Uri Keshet, Martin Lemoine

Abstract: 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… ▽ More 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. △ Less

Submitted 5 June, 2015; originally announced June 2015.

Comments: 32 pages; 9 figures; invited topical review, comments welcome; submitted for publication in "The Strongest Magnetic Fields in the Universe" (Space Sciences Series of ISSI, Springer), Space Science Reviews

Journal ref: Space Science Reviews, Volume 191, Issue 1-4, pp. 519-544 (2015)

Magnetohydrodynamics using path or stream functions

Authors: Yossi Naor, Uri Keshet

Abstract: Magnetization in highly conductive plasmas is ubiquitous to astronomical systems. Flows in such media can be described by three path functions $Λ_α$, or, for a steady flow, by two stream functions $λ_κ$ and an additional field such as the mass density $ρ$, velocity $v$, or travel time $Δt$. While typical analyses of a frozen magnetic field $\boldsymbol{B}$ are problem-specific and involve nonlocal… ▽ More Magnetization in highly conductive plasmas is ubiquitous to astronomical systems. Flows in such media can be described by three path functions $Λ_α$, or, for a steady flow, by two stream functions $λ_κ$ and an additional field such as the mass density $ρ$, velocity $v$, or travel time $Δ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 $Λ$ or $λ$ space) solution $\boldsymbol{B}=(\partial\boldsymbol{x}/\partialΛ_α)_t \tilde{B}_αρ/\tildeρ$, 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ρ\boldsymbol{B}/ρ=(\partial\boldsymbol{x}/\partialλ_κ)_{Δt}\tilde{B}_κ+\boldsymbol{v}\tilde{B}_3/\tilde{v}$; here the electric field $\boldsymbol{E}\sim(\tilde{B}_2\boldsymbol{\nabla}λ_1-\tilde{B}_1\boldsymbol{\nabla}λ_2)/\tildeρ$ depends only on $λ_κ$ 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) △ Less

Submitted 16 August, 2015; v1 submitted 27 January, 2015; originally announced January 2015.

Comments: ApJ (accepted); high resolution version at https://sites.google.com/site/naoryossi7/publications

Journal ref: ApJ, 810, 152, 2015

A Dynamically Driven, Universal Thermal Profile of Galaxy Groups and Clusters

Authors: Ido Reiss, Uri Keshet

Abstract: 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 $τ\propto r^{0.72\pm0.01}$ ratio between cooling an… ▽ More 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 $τ\propto r^{0.72\pm0.01}$ ratio between cooling and free-fall times. Both profiles slightly flatten at small radii, as $τ$ 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. △ Less

Submitted 1 September, 2015; v1 submitted 29 December, 2014; originally announced December 2014.

Comments: ApJ (accepted)

Journal ref: ApJ 810 129 (2015)

Compressible flow in front of an axisymmetric blunt object

Authors: Uri Keshet, Yossi Naor

Abstract: 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 a… ▽ More 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. △ Less

Submitted 18 November, 2016; v1 submitted 29 December, 2014; originally announced December 2014.

Comments: This letter is superseded by a later, more detailed version with a parameter-free analysis

Preliminary evidence for a virial shock around the Coma galaxy cluster

Authors: Uri Keshet, Doron Kushnir, Abraham Loeb, Eli Waxman

Abstract: 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 $γ$-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'dovic… ▽ More 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 $γ$-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 $γ$-ray ring around the Coma cluster, elongated towards the large scale filament connecting Coma and Abell 1367, detected at the nominal $2.7σ$ confidence level ($5.1σ$ using control signal simulations). The $γ$-ray ring correlates both with a synchrotron signal and with the SZ cutoff, but not with Galactic tracers. The $γ$-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 $γ$-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. △ Less

Submitted 21 September, 2017; v1 submitted 4 October, 2012; originally announced October 2012.

Comments: Replaced with published version

Journal ref: The Astrophysical Journal, Volume 845, Issue 1, article id. 24, 14 pp. (2017)

Strong Magnetization Measured in the Cool Cores of Galaxy Clusters

Authors: Ido Reiss, Uri Keshet

Abstract: 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 magn… ▽ More 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. △ Less

Submitted 4 August, 2014; v1 submitted 31 July, 2012; originally announced August 2012.

Comments: PRL accepted, additional control tests added

Journal ref: Phys. Rev. Lett. 113, 071302 (2014)

The missing asymptotic sector of rotating black-hole spectroscopy

Authors: Uri Keshet, Arnon Ben-Meir

Abstract: 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 ω(n>>1) ~ mΩ- 2πiT(n-s), where Ωand T are the horizon angular velocity and temperature, s is the field spin, and m is an azimuthal eigenvalue. Some physical implicati… ▽ More 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 ω(n>>1) ~ mΩ- 2πiT(n-s), where Ω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. △ Less

Submitted 10 July, 2012; originally announced July 2012.

Comments: 4 pages, comments welcome

Journal ref: Physics Letters B, Volume 737, p. 379-382 (2014)

Interlayer Registry Determines the Sliding Potential of Layered Metal Dichalcogenides: The case of 2H-MoS2

Authors: Adi Blumberg, Uri Keshet, Inbal Zaltsman, Oded Hod

Abstract: We provide a simple and intuitive explanation for the interlayer sliding energy landscape of metal dichalcogenides. Based on the recently introduced registry index (RI) concept, we define a purely geometrical parameter which quantifies the degree of interlayer commensurability in the layered phase of molybdenum disulphide (2HMoS2). A direct relation between the sliding energy landscape and the cor… ▽ More We provide a simple and intuitive explanation for the interlayer sliding energy landscape of metal dichalcogenides. Based on the recently introduced registry index (RI) concept, we define a purely geometrical parameter which quantifies the degree of interlayer commensurability in the layered phase of molybdenum disulphide (2HMoS2). A direct relation between the sliding energy landscape and the corresponding interlayer registry surface of 2H-MoS2 is discovered thus marking the registry index as a computationally efficient means for studying the tribology of complex nanoscale material interfaces in the wearless friction regime. △ Less

Submitted 16 May, 2012; originally announced May 2012.

Comments: 13 pages, 7 figures

Journal ref: J. Phys. Chem. Lett. 3, 1936-1940 (2012)

Critical Conditions for Core-Collapse Supernovae

Authors: Uri Keshet, Shmuel Balberg

Abstract: 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 \… ▽ More 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. △ Less

Submitted 13 June, 2012; v1 submitted 10 February, 2012; originally announced February 2012.

Comments: PRL accepted

Journal ref: Phys. Rev. Lett. 108, 251101 (2012)

Spiral flows in cool-core galaxy clusters

Authors: Uri Keshet

Abstract: 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, imp… ▽ More 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, ρ\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. △ Less

Submitted 24 May, 2012; v1 submitted 9 November, 2011; originally announced November 2011.

Comments: Revised version to appear in the Astrophysical Journal

Journal ref: ApJ, 753, 120 (2012)

Common origin for radio relics and halos: galaxy cluster-wide, homogeneous cosmic-ray distribution, and evolving magnetic fields

Authors: Uri Keshet

Abstract: 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, r… ▽ More 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. △ Less

Submitted 2 November, 2010; originally announced November 2010.

Comments: 68 pages, 31 figures; comments welcome

Cold Fronts by Merging of Shocks

Authors: Yuval Birnboim, Uri Keshet, Lars Hernquist

Abstract: 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 expec… ▽ More 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. △ Less

Submitted 9 June, 2010; originally announced June 2010.

Comments: 16 pages, 13 figures, MNRAS accepted

Journal ref: MNRAS, 408, Issue 1, pp. 199-212 (2010)

Using Radio Halos and Minihalos to Measure the Distributions of Magnetic Fields and Cosmic-Rays in Galaxy Clusters

Authors: Uri Keshet, Abraham Loeb

Abstract: 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.… ▽ More 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. △ Less

Submitted 3 November, 2010; v1 submitted 5 March, 2010; originally announced March 2010.

Comments: Replaced with published version

Journal ref: The Astrophysical Journal, Volume 722, Issue 1, pp. 737-749 (2010)

Dynamics and Magnetization in Galaxy Cluster Cores Traced by X-ray Cold Fronts

Authors: Uri Keshet, Maxim Markevitch, Yuval Birnboim, Abraham Loeb

Abstract: 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… ▽ More 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. △ Less

Submitted 20 July, 2010; v1 submitted 18 December, 2009; originally announced December 2009.

Comments: Revised version to appear in Astrophys.J.Lett

Journal ref: Astrophys.J.719:L74-L78,2010

Cold Fronts from Shock Collisions

Authors: Yuval Birnboim, Uri Keshet, Lars Hernquist

Abstract: 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 eff… ▽ More 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. △ Less

Submitted 18 June, 2010; v1 submitted 10 December, 2009; originally announced December 2009.

Comments: 5 pages, 4 figures, this letter is superseded by a more detailed version in arXiv:1006.1892

Analytic study of mass segregation around a massive black hole

Authors: Uri Keshet, Clovis Hopman, Tal Alexander

Abstract: 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… ▽ More 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. △ Less

Submitted 21 December, 2009; v1 submitted 28 January, 2009; originally announced January 2009.

Comments: Revised version published in Astrophys. J. Lett

Journal ref: Astrophys.J.698:L64-L67,2009

Magnetic field evolution in relativistic unmagnetized collisionless shocks

Authors: Uri Keshet, Boaz Katz, Anatoly Spitkovsky, Eli Waxman

Abstract: 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 mag… ▽ More 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. △ Less

Submitted 9 March, 2009; v1 submitted 22 February, 2008; originally announced February 2008.

Comments: ApJ Letters, published (minor corrections, results unchanged)

Journal ref: Astrophys.J.693:L127-L130,2009; Astrophys.J.Lett.693:L127-L130,2009

Asymptotic Spectroscopy of Rotating Black Holes

Authors: Uri Keshet, Andrew Neitzke

Abstract: 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 semiclas… ▽ More 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. △ Less

Submitted 30 September, 2007; v1 submitted 11 September, 2007; originally announced September 2007.

Comments: References added

Journal ref: Phys.Rev.D78:044006,2008

Analytic Study of Rotating Black-Hole Quasinormal Modes

Authors: Uri Keshet, Shahar Hod

Abstract: 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.… ▽ More 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. △ Less

Submitted 24 August, 2007; v1 submitted 8 May, 2007; originally announced May 2007.

Comments: Shortened and typos fixed; to appear in PRD Rapid communications

Journal ref: Phys.Rev.D76:061501,2007