Cosmological Shocks Research Articles

Cosmological shocks around galaxy clusters: a coherent investigation with DES, SPT, and ACT

ABSTRACT We search for signatures of cosmological shocks in gas pressure profiles of galaxy clusters using the cluster catalogues from three surveys: the Dark Energy Survey (DES) Year 3, the South Pole Telescope (SPT) SZ survey, and the Atacama Cosmology Telescope (ACT) data releases 4, 5, and 6, and using thermal Sunyaev–Zeldovich (SZ) maps from SPT and ACT. The combined cluster sample contains around 105 clusters with mass and redshift ranges $10^{13.7} \lt M_{\rm 200m}/\, {\rm M}_\odot \lt 10^{15.5}$ and 0.1 < z < 2, and the total sky coverage of the maps is $\approx 15\, 000 \deg ^2$. We find a clear pressure deficit at R/R200m ≈ 1.1 in SZ profiles around both ACT and SPT clusters, estimated at 6σ significance, which is qualitatively consistent with a shock-induced thermal non-equilibrium between electrons and ions. The feature is not as clearly determined in profiles around DES clusters. We verify that measurements using SPT or ACT maps are consistent across all scales, including in the deficit feature. The SZ profiles of optically selected and SZ-selected clusters are also consistent for higher mass clusters. Those of less massive, optically selected clusters are suppressed on small scales by factors of 2–5 compared to predictions, and we discuss possible interpretations of this behaviour. An oriented stacking of clusters – where the orientation is inferred from the SZ image, the brightest cluster galaxy, or the surrounding large-scale structure measured using galaxy catalogues – shows the normalization of the one-halo and two-halo terms vary with orientation. Finally, the location of the pressure deficit feature is statistically consistent with existing estimates of the splashback radius.

Cosmic-Ray Acceleration and Nonthermal Radiation at Accretion Shocks in the Outer Regions of Galaxy Clusters

Cosmology models predict that external accretion shocks form in the outer region of galaxy clusters owing to supersonic gas infall from filaments and voids in the cosmic web. They are characterized by high sonic and Alfvénic Mach numbers, M s ∼ 10–102 and M A ∼ 102–103, and propagate into weakly magnetized plasmas of β ≡ P g /P B ≳ 102. Although strong accretion shocks are expected to be efficient accelerators of cosmic rays (CRs), nonthermal signatures of shock-accelerated CRs around clusters have not been confirmed, and detailed acceleration physics at such shocks has yet to be understood. In this study, we first establish through two-dimensional particle-in-cell simulations that at strong high-β shocks electrons can be pre-energized via stochastic Fermi acceleration owing to the ion Weibel instability in the shock transition region, possibly followed by injection into diffusive shock acceleration. Hence, we propose that the models derived from conventional thermal leakage injection may be employed for the acceleration of electrons and ions at accretion shocks as well. Applying these analytic models to numerical shock zones identified in structure formation simulations, we estimate nonthermal radiation, such as synchrotron and inverse Compton (IC) emission due to CR electrons and π 0-decay γ-rays due to CR protons, around simulated clusters. Our models with injection parameter Q ≈ 3.5–3.8 predict synthetic synchrotron maps, which seem consistent with recent radio observations of the Coma Cluster. However, the detection of nonthermal IC X-rays and γ-rays from accretion shocks would be quite challenging. We suggest that the proposed analytic models may be adopted as generic recipes for CR production at cosmological shocks.

Shocks in the stacked Sunyaev-Zel’dovich profiles of clusters II: Measurements from SPT-SZ + Planck Compton-y map

ABSTRACT We search for the signature of cosmological shocks in stacked gas pressure profiles of galaxy clusters using data from the South Pole Telescope (SPT). Specifically, we stack the latest Compton-y maps from the 2500 deg2 SPT-SZ survey on the locations of clusters identified in that same data set. The sample contains 516 clusters with mean mass $\langle M_{\rm 200m}\rangle = 10^{14.9} \, {\rm M}_\odot$ and redshift 〈z〉 = 0.55. We analyse in parallel a set of zoom-in hydrodynamical simulations from the three hundred project. The SPT-SZ data show two features: (i) a pressure deficit at R/R200m = 1.08 ± 0.09, measured at 3.1σ significance and not observed in the simulations, and; (ii) a sharp decrease in pressure at R/R200m = 4.58 ± 1.24 at 2.0σ significance. The pressure deficit is qualitatively consistent with a shock-induced thermal non-equilibrium between electrons and ions, and the second feature is consistent with accretion shocks seen in previous studies. We split the cluster sample by redshift and mass, and find both features exist in all cases. There are also no significant differences in features along and across the cluster major axis, whose orientation roughly points towards filamentary structure. As a consistency test, we also analyse clusters from the Planck and Atacama Cosmology Telescope Polarimeter surveys and find quantitatively similar features in the pressure profiles. Finally, we compare the accretion shock radius ($R_{\rm sh,\, acc}$) with existing measurements of the splashback radius (Rsp) for SPT-SZ and constrain the lower limit of the ratio, $R_{\rm sh,\, acc}/R_{\rm sp}\gt 2.16 \pm 0.59$.

Exploring the role of cosmological shock waves in the Dianoga simulations of galaxy clusters

ABSTRACT Cosmological shock waves are ubiquitous to cosmic structure formation and evolution. As a consequence, they play a major role in the energy distribution and thermalization of the intergalactic medium (IGM). We analyze the Mach number distribution in the Dianoga simulations of galaxy clusters performed with the SPH code gadget-3. The simulations include the effects of radiative cooling, star formation, metal enrichment, supernova, and active galactic nuclei feedback. A grid-based shock-finding algorithm is applied in post-processing to the outputs of the simulations. This procedure allows us to explore in detail the distribution of shocked cells and their strengths as a function of cluster mass, redshift, and baryonic physics. We also pay special attention to the connection between shock waves and the cool-core/non-cool-core (CC/NCC) state and the global dynamical status of the simulated clusters. In terms of general shock statistics, we obtain a broad agreement with previous works, with weak (low-Mach number) shocks filling most of the volume and processing most of the total thermal energy flux. As a function of cluster mass, we find that massive clusters seem more efficient in thermalizing the IGM and tend to show larger external accretion shocks than less massive systems. We do not find any relevant difference between CC and NCC clusters. However, we find a mild dependence of the radial distribution of the shock Mach number on the cluster dynamical state, with disturbed systems showing stronger shocks than regular ones throughout the cluster volume.

Collisionless shock heating of heavy ions in SN 1987A

Astrophysical shocks at all scales, from those in the heliosphere up to the cosmological shock waves, are typically "collisionless", because the thickness of their jump region is much shorter than the collisional mean free path. Across these jumps, electrons, protons, and ions are expected to be heated at different temperatures. Supernova remnants (SNRs) are ideal targets to study collisionless processes because of their bright post-shock emission and fast shocks. Although optical observations of Balmer-dominated shocks in young SNRs showed that the post-shock proton temperature is higher than the electron temperature, the actual dependence of the post-shock temperature on the particle mass is still widely debated. We tackle this longstanding issue through the analysis of deep multi-epoch and high-resolution observations of the youngest nearby supernova remnant, SN 1987A, made with the Chandra X-ray telescope. We introduce a novel data analysis method by studying the observed spectra in close comparison with a dedicated full 3-D hydrodynamic simulation. The simulation is able to reproduce self-consistently the whole broadening of the spectral lines of many ions altogether. We can therefore measure the post shock temperature of protons and selected ions through comparison of the model with observations. We have obtained information about the heating processes in collisional shocks by finding that the ion to proton temperature ratio is always significantly higher than one and increases linearly with the ion mass for a wide range of masses and shock parameters.

On the interplay between cosmological shock waves and their environment

Cosmological shock waves are tracers of the thermal history of the structures in the Universe. They play a crucial role in redistributing the energy within the cosmic structures and are also amongst the main ingredients of galaxy and galaxy cluster formation. Understanding this important function requires a proper description of the interplay between shocks and the different environments where they can be found. In this paper, an Adaptive Mesh Refinement (AMR) Eulerian cosmological simulation is analysed by means of a shock-finding algorithm that allows to generate shock wave maps. Based on the population of dark matter halos and on the distribution of density contrast in the simulation, we classify the shocks in five different environments. These range from galaxy clusters to voids. The shock distribution function and the shocks power spectrum are studied for these environments dynamics. We find that shock waves on different environments undergo different formation and evolution processes, showing as well different characteristics. We identify three different phases of formation, evolution and dissipation of these shock waves, and an intricate migration between distinct environments and scales. Shock waves initially form at external, low density regions and are merged and amplified through the collapse of structures. Shock waves and cosmic structures follow a parallel evolution. Later on, shocks start to detach from them and dissipate. We also find that most of the power that shock waves dissipate is found at scales of k sim0.5~mbox{Mpc}^{-1}, with a secondary peak at k sim8~mbox{Mpc}^{-1}. The evolution of the shocks power spectrum confirms that shock waves evolution is coupled and conditioned by their environment.

RADIO AND X-RAY SHOCKS IN CLUSTERS OF GALAXIES

Radio relics detected in the outskirts of galaxy clusters are thought to trace radio-emitting relativistic electrons accelerated at cosmological shocks. In this study, using the cosmological hydrodynamic simulation data for the large-scale structure formation and adopting a diffusive shock acceleration (DSA) model for the production of cosmic-ray (CR) electrons, we construct mock radio and X-ray maps of simulated galaxy clusters that are projected in the sky plane. Various properties of shocks and radio relics, including the shock Mach number, radio spectral index and luminosity are extracted from the synthetic maps and compared with observations. A substantial fraction of radio and X-ray shocks identified in these maps involve multiple shock surfaces along line of sights (LoSs), and the morphology of shock distributions in the maps depends on the projection direction. Among multiple shocks in a given LoS, radio observations tend to pick up stronger shocks with flatter radio spectra, while X-ray observations preferentially select weaker shocks with larger kinetic energy flux. As a result, the shock Mach numbers and locations derived from radio and X-ray observations could differ from each other in some cases. We also find that the distributions of the spectral index and radio power of the synthetic radio relics are somewhat inconsistent with those of observed real relics; a bit more radio relics have been observed closer to the cluster core and with steeper spectral indices. We suggest the inconsistency could be explained, if very weak shocks with $M_s \la 2$ accelerate CR electrons more efficiently, compared to the DSA model adopted here.

Shock finding on a moving mesh – I. Shock statistics in non-radiative cosmological simulations

Cosmological shock waves play an important role in hierarchical structure formation by dissipating and thermalizing kinetic energy of gas flows, thereby heating the universe. Furthermore, identifying shocks in hydrodynamical simulations and measuring their Mach number accurately is critical for calculating the production of non-thermal particle components through diffusive shock acceleration. However, shocks are often significantly broadened in numerical simulations, making it challenging to implement an accurate shock finder. We here introduce a refined methodology for detecting shocks in the moving-mesh code AREPO, and show that results for shock statistics can be sensitive to implementation details. We put special emphasis on filtering against spurious shock detections due to tangential discontinuities and contacts. Both of them are omnipresent in cosmological simulations, for example in the form of shear-induced Kelvin-Helmholtz instabilities and cold fronts. As an initial application of our new implementation, we analyse shock statistics in non-radiative cosmological simulations of dark matter and baryons. We find that the bulk of energy dissipation at redshift zero occurs in shocks with Mach numbers around ${\cal M}\approx2.7$. Furthermore, almost $40\%$ of the thermalization is contributed by shocks in the warm hot intergalactic medium (WHIM), whereas $\approx60\%$ occurs in clusters, groups and smaller halos. Compared to previous studies, these findings revise the characterization of the most important shocks towards higher Mach numbers and lower density structures. Our results also suggest that regions with densities above and below $\delta_b=100$ should be roughly equally important for the energetics of cosmic ray acceleration through large-scale structure shocks.

Probing the cosmic-ray content of galaxy clusters by stackingFermi-LAT count maps

Observations in radio have shown that galaxy clusters are giant reservoirs of cosmic rays (CR). Although a gamma- ray signal from the cluster volume is expected to arise through interactions of CR protons with the ambient plasma, a confirming observation is still missing. We search for a cumulative gamma-ray emission in direction of galaxy clusters by analysing a collection of stacked Fermi-LAT count maps. Additionally, we investigate possible systematic differences in the emission between cool-core and non-cool core cluster populations. Making use of a sample of 53 clusters selected from the HIFLUGCS catalog, we do not detect a significant signal from the stacked sample. The upper limit on the average flux per cluster derived for the total stacked sample is at the level of a few 1e-11 ph cm-2 s-1 at 95% confidence level in the 1-300 GeV band, assuming power-law spectra with photon indices 2.0, 2.4, 2.8 and 3.2. Separate stacking of the cool core and non-cool core clusters in the sample lead to similar values of around 5e-11 ph cm-2 s-1 and 2e-11 ph cm-2 s-1, respectively. Under the assumption that decaying pi0, produced in collisions between CRs and the ambient thermal gas, are responsible for the gamma-ray emission, we set upper limits on the average CR content in galaxy clusters. For the entire cluster population, our upper limit on the gamma-ray flux translates into an upper limit on the average CR-to-thermal energy ratio of 4.6% for a photon index of 2.4, although it is possible for individual systems to exceed this limit. Our 95% upper limits are at the level expected from numerical simulations, which likely suggests that the injection of CR at cosmological shocks is less efficient than previously assumed.

DIFFUSIVE SHOCK ACCELERATION AT COSMOLOGICAL SHOCK WAVES

We reexamine nonlinear diffusive shock acceleration (DSA) at cosmological shocks in the large scale structure of the Universe, incorporating wave-particle interactions that are expected to operate in collisionless shocks. Adopting simple phenomenological models for magnetic field amplification (MFA) by cosmic-ray (CR) streaming instabilities and Alfv'enic drift, we perform kinetic DSA simulations for a wide range of sonic and Alfv'enic Mach numbers and evaluate the CR injection fraction and acceleration efficiency. In our DSA model the CR acceleration efficiency is determined mainly by the sonic Mach number Ms, while the MFA factor depends on the Alfv'enic Mach number and the degree of shock modification by CRs. We show that at strong CR modified shocks, if scattering centers drift with an effective Alfv'en speed in the amplified magnetic field, the CR energy spectrum is steepened and the acceleration efficiency is reduced significantly, compared to the cases without such effects. As a result, the postshock CR pressure saturates roughly at ~ 20 % of the shock ram pressure for strong shocks with Ms>~ 10. In the test-particle regime (Ms<~ 3), it is expected that the magnetic field is not amplified and the Alfv'enic drift effects are insignificant, although relevant plasma physical processes at low Mach number shocks remain largely uncertain.

COSMOLOGICAL SHOCK WAVES IN THE LARGE SCALE STRUCTURE OF THE UNIVERSE

Based on the cosmological hydrodynamic simulation, we study the properties of shock waves formed during the formation of the large scale structure (LSS) of the universe, and investigate their contribution to the cosmic ray (CR) fraction in the intergalactic medium (IGM). It is found that while strong accretion shocks prevail at high redshift, weak internal shocks become dominant in the intracluster medium (ICM) as galaxy clusters form and virialize at low redshift, z &lt; 1. The accumulated CR proton energy is likely to be less than 10 % of the thermal energy in the ICM, since weak shocks of M ≲ 3 are most abundant. This is consistent with the upper limit constrained by radio and gamma-ray observations of galaxy clusters. In the warm-hot medium (WHIM) inside filaments, CRs and gas could be almost in energy equipartition, since relatively stronger shocks of 5 ≲ M ≲ 10 are dominant there. We suggest that the non-thermal emissions from the CR electrons and protons accelerated by cosmological shock waves could provide a new way to detect the WHIM of the universe.

Cosmological shock waves: clues to the formation history of haloes

Shock waves developed during the formation and evolution of cosmic structures encode crucial information on the hierarchical formation of the Universe. We analyze an Eulerian AMR hydro + N-body simulation in a $\Lambda$CDM cosmology focused on the study of cosmological shock waves. The combination of a shock-capturing algorithm together with the use of a halo finder allows us to study the morphological structures of the shock patterns, the statistical properties of shocked cells, and the correlations between the cosmological shock waves appearing at different scales and the properties of the haloes harbouring them. The shocks in the simulation can be split into two broad classes: internal weak shocks related with evolutionary events within haloes, and external strong shocks associated with large-scale events. The shock distribution function contains information on the abundances and strength of the different shocks, and it can be fitted by a double power law with a break in the slope around a Mach number of 20. We introduce a generalised scaling relation that correlates the average Mach numbers within the virial radius of haloes and their virial masses. In this plane, haloes occupy different areas according to their early evolutionary histories: those with a quiet evolution have an almost constant Mach number independently of their masses, whereas haloes undergoing significant merger events very early in their evolution show a linear dependence with their masses. At high redshift, the halo distribution in this scaling relation forms a L-like pattern that changes due to the evolution of the haloes. The analysis of the propagation speed and size of the shock waves around haloes could give some hints on the formation time and main features of the haloes. (Abridged)

Gamma-ray signatures of cosmic ray acceleration, propagation, and confinement in the era of CTA

Galactic cosmic rays are commonly believed to be accelerated at supernova remnants via diffusive shock acceleration. Despite the popularity of this idea, a conclusive proof for its validity is still missing. Gamma-ray astronomy provides us with a powerful tool to tackle this problem, because gamma rays are produced during cosmic ray interactions with the ambient gas. The detection of gamma rays from several supernova remnants is encouraging, but still does not constitute a proof of the scenario, the main problem being the difficulty in disentangling the hadronic and leptonic contributions to the emission. Once released by their sources, cosmic rays diffuse in the interstellar medium, and finally escape from the Galaxy. The diffuse gamma-ray emission from the Galactic disk, as well as the gamma-ray emission detected from a few galaxies is largely due to the interactions of cosmic rays in the interstellar medium. On much larger scales, cosmic rays are also expected to permeate the intracluster medium, since they can be confined and accumulated within clusters of galaxies for cosmological times. Thus, the detection of gamma rays from clusters of galaxies, or even upper limits on their emission, will allow us to constrain the cosmic ray output of the sources they contain, such as normal galaxies, AGNs, and cosmological shocks. In this paper, we describe the impact that the Cherenkov Telescope Array, a future ground-based facility for very-high energy gamma-ray astronomy, is expected to have in this field of research.

Experience and the Ultimacy of God: In Memoriam John Edwin Smith

Experience and the Ultimacy of God: In Memoriam John Edwin Smith

MULTI-FREQUENCY STUDIES OF RADIO RELICS IN THE GALAXY CLUSTERS A4038, A1664, AND A786

We present a multi-frequency study of radio relics associated with the galaxy clusters A4038, A1664 and A786. Radio images, integrated spectra, spectral index maps and fits to the integrated spectra in the framework of the adiabatic compression model are presented. Images of the relic in A4038 at 150, 240 and 606 MHz with the Giant Metrewave Radio Telescope have revealed extended ultra-steep spectrum (\alpha ~ -1.8 to -2.7) emission of an extent 210 x 80 kpc^2. The model of passively evolving radio lobes compressed by a shock best fits the integrated spectrum. The relic with a circular morphology at the outskirts of the cluster A1664 has an integrated spectral index of -1.10+/-0.06 and is best fit by the model of radio lobes lurking for ~ 4x10^7 yr. The relic near A786 has a curved spectrum and is best fit by a model of radio lobes lurking for ~ 3x10^7 yr. At 4.7 GHz, a compact radio source, possibly the progenitor of the A786 relic, is detected near the center of the radio relic. The A786 radio relic is thus likely a lurking radio galaxy rather than a site of cosmological shock as has been considered in earlier studies.

GALAXY CLUSTER RADIO RELICS IN ADAPTIVE MESH REFINEMENT COSMOLOGICAL SIMULATIONS: RELIC PROPERTIES AND SCALING RELATIONSHIPS

Cosmological shocks are a critical part of large-scale structure formation, and are responsible for heating the intracluster medium in galaxy clusters. In addition, they are also capable of accelerating non-thermal electrons and protons. In this work, we focus on the acceleration of electrons at shock fronts, which is thought to be responsible for radio relics - extended radio features in the vicinity of merging galaxy clusters. By combining high resolution AMR/N-body cosmological simulations with an accurate shock finding algorithm and a model for electron acceleration, we calculate the expected synchrotron emission resulting from cosmological structure formation. We produce synthetic radio maps of a large sample of galaxy clusters and present luminosity functions and scaling relationships. With upcoming long wavelength radio telescopes, we expect to see an abundance of radio emission associated with merger shocks in the intracluster medium. By producing observationally motivated statistics, we provide predictions that can be compared with observations to further improve our understanding of magnetic fields and electron shock acceleration.

Cosmic Rays Accelerated at Cosmological Shock Waves

Based on hydrodynamic numerical simulations and diffusive shock acceleration model, we calculated the ratio of cosmic ray (CR) to thermal energy. We found that the CR fraction can be less than ~0.1 in the intracluster medium, while it would be of order unity in the warm-hot intergalactic medium.

RE-ACCELERATION OF NON-THERMAL PARTICLES AT WEAK COSMOLOGICAL SHOCK WAVES

We examine diffusive shock acceleration (DSA) of the pre-exisiting as well as freshly injected populations of nonthermal, cosmic-ray (CR) particles at weak cosmological shocks. Assuming simple models for thermal leakage injection and Alfv\'enic drift, we derive analytic, time-dependent solutions for the two populations of CRs accelerated in the test-particle regime. We then compare them with the results from kinetic DSA simulations for shock waves that are expected to form in intracluster media and cluster outskirts in the course of large-scale structure formation. We show that the test-particle solutions provide a good approximation for the pressure and spectrum of CRs accelerated at these weak shocks. Since the injection is extremely inefficient at weak shocks, the pre-existing CR population dominates over the injected population. If the pressure due to pre-existing CR protons is about 5 % of the gas thermal pressure in the upstream flow, the downstream CR pressure can absorb typically a few to 10 % of the shock ram pressure at shocks with the Mach number $M \la 3$. Yet, the re-acceleration of CR electrons can result in a substantial synchrotron emission behind the shock. The enhancement in synchrotron radiation across the shock is estimated to be about a few to several for $M \sim 1.5$ and $10^2-10^3$ for $M \sim 3$, depending on the detail model parameters. The implication of our findings for observed bright radio relics is discussed.

Cosmological Shocks in Adaptive Mesh Refinement Simulations and the Acceleration of Cosmic Rays

We present new results characterizing cosmological shocks within adaptive mesh refinement N-body/hydrodynamic simulations that are used to predict nonthermal components of large-scale structure. This represents the first study of shocks using adaptive mesh refinement. We propose a modified algorithm for finding shocks from those used on unigrid simulations that reduces the shock frequency of low Mach number shocks by a factor of ~3. We then apply our new technique to a large, (512 h−1 Mpc)3, cosmological volume and study the shock Mach number () distribution as a function of preshock temperature, density, and redshift. Because of the large volume of the simulation, we have superb statistics that result from having thousands of galaxy clusters. We find that the Mach number evolution can be interpreted as a method to visualize large-scale structure formation. Shocks with 20 generally follow accretion onto filaments and galaxy clusters, respectively. By applying results from nonlinear diffusive shock acceleration models using the first-order Fermi process, we calculate the amount of kinetic energy that is converted into cosmic-ray protons. The acceleration of cosmic-ray protons is large enough that in order to use galaxy clusters as cosmological probes, the dynamic response of the gas to the cosmic rays must be included in future numerical simulations.

Shock waves in the large-scale structure of the Universe

Cosmological shock waves are induced during hierarchical formation of large-scale structure in the universe. Like most astrophysical shocks, they are collisionless, since they form in the tenuous intergalactic medium through electromagnetic viscosities. The gravitational energy released during structure formation is transferred by these shocks to the intergalactic gas as heat, cosmic-rays, turbulence, and magnetic fields. Here we briefly described the properties and consequences of the shock waves in the context of the large-scale structure of the universe.