Turbulence and Particle Acceleration in Giant Radio Halos: the Origin of Seed Electrons
About of X-ray-luminous clusters show smooth, unpolarized radio emission on Mpc scales, known as giant radio halos. One promising model for radio halos is Fermi II acceleration of seed relativistic electrons by turbulence of the intracluster medium (ICM); Coulomb losses prohibit acceleration from the thermal pool. However, the origin of seed electrons has never been fully explored. Here, we integrate the Fokker-Planck equation of the cosmic ray (CR) electron and proton distributions in a cosmological simulations of cluster formation. For standard assumptions, structure formation shocks lead to a seed electron population which produces too centrally concentrated radio emission. Instead, we present three realistic scenarios that each can reproduce the spatially flat radio emission observed in the Coma cluster: (1) the ratio of injected turbulent energy density to thermal energy density increase significantly with radius, as seen in cosmological simulations. This generates a flat radio profile even if the seed population of CRs is steep with radius. (2) Self-confinement of energetic CR protons can be inefficient, and CR protons may stream at the Alfvén speed to the cluster outskirts when the ICM is relatively quiescent. A spatially flat CR proton distribution develops and produces the required population of secondary seed electrons. (3) The CR proton to electron acceleration efficiency is assumed to be larger than in our Galaxy (), due to the magnetic geometry at the shock. The resulting primary electron population dominates. Due to their weaker density dependence compared to secondary electrons, these primaries can also reproduce radio observations. These competing non-trivial solutions provide incisive probes of non thermal processes in the high- ICM.
About one third of X-ray-luminous clusters show smooth, unpolarized radio emission on Mpc scales, known as giant radio halos (RHs) (Brunetti and Jones, 2014). They appear only in disturbed, merging clusters and the RH luminosity correlates with the X-ray luminosity (Govoni et al., 2001; Feretti et al., 2012) and the Compton -parameter (Planck Collaboration et al., 2013). The RHs show that CRs and magnetic fields permeate a large volume fraction of the intra-cluster medium (ICM). The dominant CR source, given the smoothness and enormous extent of RHs, is thought to be structure formation shocks (Miniati et al., 2001; Pfrommer, 2008). At the same time, plasma processes, the origin of magnetic fields and particle acceleration in a turbulent, high- plasma like the ICM are not well understood. Radio halos thus provide an incisive probe of non-thermal processes in the ICM.
One promising model for RHs is re-energetization of seed suprathermal electrons by Fermi II acceleration when ICM turbulence becomes transonic during mergers (Brunetti et al., 2001; Brunetti and Lazarian, 2007). Due to the short radiative cooling time of high-energy relativistic electrons, the cluster synchrotron emission quickly fades away after a merger, which naturally explains the observed bimodality of RHs.
However, there is a salient piece missing in the turbulent reacceleration model. It relies heavily on the assumption of an abundant, volume-filling population of seed suprathermal electrons; direct Fermi II acceleration from the thermal pool is precluded by strong Coulomb losses (Petrosian and East, 2008; Chernyshov et al., 2012). These seeds are presumed to be either fossil CR electrons (CRes) accelerated by diffusive shock acceleration (DSA) during structure formation (Sarazin, 1999), or secondaries injected by hadronic interaction of CR protons (CRps) with thermal protons (Brunetti and Lazarian, 2011). While analytic estimates have been made, there has been no ab initio demonstration that structure formation can lead to the required abundance of seed electrons with the correct spatial and spectral characteristics. This is a non-trivial requirement: Coulomb cooling in dense cluster cores is severe, and DSA fossil electrons may not survive. On the other hand, for secondaries to constitute the seed population, the CRp population required in the best-studied case of the Coma cluster must have a very broad and flat (or even slightly inverted) spatial profile (Brunetti et al., 2012), in contrast with the thermal plasma whose energy density declines steeply with radius. In this Letter we show that such a distribution is not predicted by cosmological simulations (see lower right panel of Fig. 1) (see also Pinzke and Pfrommer, 2010; Vazza et al., 2014).
Indeed, arriving at a seed population with the required characteristics is highly constraining, and has the potential to teach us much about the origin of CRps/CRes in clusters. We consider 3 new possibilities: (i) Our model M-turbulence: a significantly flatter turbulent profile than what was adopted in Brunetti et al. (2012), which allows seed CRps to follow the steep profile that is suggested by structure formation simulations. (ii) Our model M-streaming: streaming CRps that produce flat distributions of CRs in the ICM (Enßlin et al., 2011; Wiener et al., 2013), which also flattens the secondary electron distribution. (iii) Our model M-primaries: if the acceleration efficiency of CRps is below about % in weak (perpendicular) shocks and the ratio of injected electrons-to-protons , this yields a dominant primary population with a flat spatial distribution, since primaries have a weaker density dependence than secondaries. In this work we pursue these three possibilities further. We employ cosmological simulations of CRs in clusters, in tandem with new insights from our recent work on DSA generated fossil electrons (Pinzke et al., 2013), to generate the first quantitative calculation of primary and secondary seed electrons.
Method. The transport of relativistic electrons and protons in the ICM is a complex process that depends both on the details of the thermal component (gas density, temperature, and pressure) as well as non-thermal component (turbulence, magnetic fields, fossil CRs). We use high resolution galaxy cluster simulations to derive the thermal and fossil CR properties (shock accelerated primary CRes and CRps, as well as secondary CRes produced in p-p collisions) Pfrommer et al. (2007); Pfrommer (2008); Pinzke and Pfrommer (2010); Pinzke et al. (2013). In this Letter we focus on our simulated cluster, g72a, which is a massive cluster that experienced a merger about 1-2 Gyrs ago. Since the cluster mass, density and temperature profiles are all similar to the well studied Coma cluster Pfrommer et al. (2007); Pinzke and Pfrommer (2010), we will compare our calculations to radio and gamma-ray observations of Coma.
In our Galaxy, the CRe-to-CRp ratio at a few GeV is . Hence, we adopt this as a fiducial value for the CRe-to-CRp acceleration efficiency (see Pinzke et al. (2013) for more discussion). However, as recent PIC simulations have shown, this is likely very different at weak shocks, with electrons efficiently accelerated at perpendicular shocks (Guo et al., 2014a, b), and ions efficiently accelerated at parallel shocks (Caprioli and Spitkovsky, 2014). Thus, depending on magnetic geometry, could be either larger or smaller. In this work we use a simple test-particle model for the CRp acceleration Kang and Ryu (2011); Pinzke et al. (2013). The ratio of accelerated proton-to-dissipated energy in the downstream of strong shocks varies from 1-10%, depending on the adopted model (for more details, see the Results section), and is a factor 10-100 lower for weak shocks. However, some observations of radio relics suggest higher values of , due to the absence of gamma-ray emission, which probes the CRp population Vazza and Brüggen (2014). This suggests primary CRes as a viable alternative scenario to secondary CRes as seeds for the giant RHs. In our M-primaries scenario, we adopt (viable for primarily perpendicular shocks) to test this possibility.
As previously noted, secondaries produced by shock accelerated CRp have the wrong spatial profile to explain RH observations; because they arise from a two body process, they are too centrally concentrated. They also produce -ray emission in excess of Fermi-LAT upper limits (Arlen et al., 2012; Brunetti et al., 2012; Ackermann et al., 2014). However, if CRps stream in the ICM, then their spatial profile could potentially flatten sufficiently (Enßlin et al., 2011; Wiener et al., 2013). This scenario is very attractive: it generates seed electrons with the right spatial footprint, and by removing CRps from the core, obeys gamma-ray constraints. Turbulence plays two opposing roles: Alfvénic turbulence damps waves generated by the CR streaming instability (Yan and Lazarian, 2002; Farmer and Goldreich, 2004), thus reducing self-confinement; but compressible fast modes scatter CRs directly. Turbulent damping is still efficient for highly subsonic conditions (Wiener et al., 2013), while we assume compressible fast modes only provide effective spatial confinement during the periods of transonic, highly super-Alfvénic () turbulence associated with mergers. Thus, CRs can stream out when the cluster is kinematically quiescent. Furthermore, even Alfvénic streaming timescales are relatively short ( Gyr; Wiener et al. (2013)) compared to the timescale on which the CRp population is built up. Based on these findings, we adopt a toy model for our M-streaming scenario in which CR streaming quickly produces flat CRp profiles. We assume that CRs cannot stream significantly past perpendicular -fields at the accretion shock, so that the total number of CRs is conserved.
Given a seed population of CRs, we adopt essentially the same set of plasma physics assumptions as the reacceleration model for RHs Brunetti and Lazarian (2007, 2011), leaving exploration of parameter space to future work. We solve the isotropic, gyro-phase averaged Fokker-Planck equation (via a Crank-Nicholson scheme) for the time evolution of the CRe distribution in the Lagrangian frame Brunetti and Lazarian (2007, 2011):
Here is the one-dimensional distribution in position (suppressed for clarity), momentum and time (which is normalized such that the number density is given by ), is the Lagrangian derivative, is the gas velocity, represents radiative (r) and Coulomb (C) losses, is the momentum space diffusion coefficient, and denotes the injection rate of primary and secondary electrons in the ICM. The terms represent adiabatic gains and losses. During post-processing of our Coma-like cluster simulation, we solve the Fokker-Planck equation over a redshift interval from to 0. The simulated cluster undergoes a major merger over the last 1-2 Gyrs that injects large turbulent eddies. After about 1 Gyr those have decayed down to the scale needed to reaccelerate particles. In all our calculations we assume that turbulent reacceleration is efficiently accelerating particles for 650 Myrs and that during this turbulent phase CR streaming and spatial diffusion can be neglected. In our M-streaming model, CR streaming and diffusion are incorporated separately during kinematically quiescent times that precede the merger. As a result, flat CRp profiles are produced on relatively short timescales ( Gyr).
The time evolution of the spectral energy distribution of CRps, , is similarly given by:
where denotes the injection rate of shock accelerated CRps as a function of momentum and time , and is the timescale of hadronic losses that produce pions via CRp collisions with thermal protons of the ICM (e.g. Brunetti and Lazarian, 2011). We incorporate momentum diffusion for electrons and protons from transit-time-damping (TTD) resonance with compressible magneto-hydrodynamic (MHD) turbulence, to model Fermi-II reacceleration Brunetti and Lazarian (2007, 2011). The TTD resonance requires the wave frequency , where and are the parallel (projected along the magnetic field) wavenumber and particle velocity, respectively. This implies that the particle transit time across the confining wave region matches the wave period, . The resonance changes the component of particle momentum parallel to seed magnetic fields, which over time leads to increasing anisotropy in the particle distribution that decreases the efficiency of reacceleration with time. As in ref. Brunetti and Lazarian (2011), we assume that there exists a mechanism—such as the firehose instability—that isotropizes the CR distribution function at the gyroscale and on the reacceleration time scale, which ensures sustained efficient reacceleration with time. The particle pitch-angle averaged momentum-diffusion coefficient of isotropic particles that couple to fast magnetosonic modes via TTD resonance is Brunetti and Lazarian (2007) (Eqn. 47):
where is the thermal-to-magnetic pressure ratio, and is the speed of light. The energy density of a mode in a magnetized plasma stems from both electromagnetic fields and resonant particles. For a high- plasma, the pitch angle averaged ratio of beta-weighted magnetic-to-total energy density saturates to (see figure 2 in Brunetti and Lazarian (2007)). The pitch angle of the CR momentum with the magnetic field orientation is given by , and . Here is the phase velocity of the fast magnetosonic waves given approximately by the sound speed, . For a sound speed typical for the ICM of 1000 km/s, . As in Brunetti and Lazarian (2007), we initially assume that the velocity of turbulent eddies is throughout the cluster. This gives a turbulent acceleration time scale, , that is typically few 100 Myrs in the ICM.
We adopt a simplified isotropic Kraichnan MHD turbulent spectrum for the fast modes per elemental range of the form
for , where we assume an injection scale for the turbulence, . The volumetric injection rate of turbulent energy, , is fixed by requiring that the total turbulent energy density on the largest scales , where is the thermal energy density Brunetti and Lazarian (2007, 2011). In this work we investigate different spatial models for injected turbulence. We assume that , where for M-turbulence, for M-streaming, and for M-primaries (note that in previous work, was adopted Brunetti et al. (2012)). Our flatter turbulent profiles are motivated by fits to cosmological simulations Lau et al. (2009); Shaw et al. (2010); Battaglia et al. (2012) and the range indicates uncertainties of the turbulent profile in Coma. Future observations (by Astro-H) and simulations will help to clarify this issue. Provided dissipation of turbulence in the ICM is collisionless, turbulent cascades of compressible modes become suppressed when thermal and relativistic particles resonantly interact with magnetosonic waves via TTD on a timescale that is approaching the cascading timescale given by . Here the wave-wave diffusion coefficient of magnetosonic modes is given by
Thus, the cascade is suppressed for wave numbers above:
where kpc in the ICM. This constitutes an effective mean free path for CRs, unless plasma instabilities can mediate interactions between turbulence and particles on smaller scales Brunetti and Lazarian (2011). In this work we only consider damping via TTD due to thermal electrons, and neglect subdominant damping with thermal protons and relativistic particles. The latter will be subdominant in the ICM for a CR to thermal energy density ratio Brunetti and Lazarian (2007), which is always satisfied. The azimuthally averaged turbulent damping rate from thermal electrons Brunetti and Lazarian (2007) in a high- plasma is , where . The magnetic field in the ICM is typically . To compute the synchrotron surface brightness profiles, we use the profile of the magnetic field strength derived from Faraday rotation observations of Coma Bonafede et al. (2010) in combination with the density profile derived from X-ray measurements Briel et al. (1992).
Results and Discussion. Let us first consider the two models which rely on secondary electrons. After turbulent reacceleration, the volume-weighted, relative CRp energy density and CRp number density inside the RH for M-turbulence (M-streaming), are found to be 3 (2) % and (), respectively. As we will see later, these densities are just of the right order of magnitude to reproduce radio observations in the Coma cluster. In addition we predict the gamma-ray flux within the virial radius of the Coma cluster from CRps that produce decaying neutral pions for M-turbulence (M-streaming) with . Both fluxes are well below current limits set by Fermi-LAT Ackermann et al. (2014), and will be challenging to probe in the near future. The spectral index of the CRp distribution is relatively steep for the CRp energies GeV that are relevant for the injection of radio-emitting secondary CRes. The steep spectrum is ultimately a consequence of our test particle model for Fermi-I acceleration Pinzke et al. (2013), where we steepen the spectral index to avoid acceleration efficiencies above .
In Fig. 1, we find that all three scenarios in which the seeds undergo Fermi-II reacceleration can reproduce the Coma RH profile at 352 MHz. In the panel Brunetti et al. (2012) we show that without CR streaming or a flat turbulent profile, our simulations of reaccelerated CRs produce radio profiles that are too steep. Indeed, even using the assumptions of previous work – where complete freedom in the seed population was allowed – it is not possible to reproduce observations in both frequencies. This signals that this problem is generic and requires either additional modifications to the plasma physics of Fermi-II acceleration or a better understanding of potential observational systematics.
In principle, reacceleration via TTD leads to spectral steepening with
particle energy due to the inefficiency of the acceleration process to
counter the stronger cooling losses with increasing energy. Since
synchrotron emission peaks at frequency , this translates into a
spectral steepening of the total radio spectrum (see
Fig. 2). A given radio window samples higher
energy electrons for a decreasing field strength in the cluster
outskirts. Hence, the spectral steepening with energy should translate
into a radial spectral steepening Brunetti et al. (2012). However, because
of the weak dependence of the electron Lorentz factor on emission
frequency (), this effect is only
visible in our simulations for GHz. Most
importantly, our simulated fluid elements at a given radius sample a
broad distribution of shock history, density and temperature, which
implies very similar synchrotron brightness profiles at
MHz and 1.4 GHz. The discrepancy of the observed
and simulated 1.4 GHz profiles could instead be due to systematic flux
calibration error in single dish observations. Interestingly, we can
match the 1.4 GHz data if we reduce the zero point by adding 10% of
the central flux to every data point; this flattens the outer
In Fig. 2 we show that our three models that include Fermi-II reacceleration can individually reproduce the convexly curved total radio spectrum found in the Coma cluster. Seed CRs that do not experience turbulent reacceleration have a power-law spectrum in disagreement with observations. In order to match both the spatial and spectral profiles in Coma, we can constrain the acceleration efficiency for the strongest shocks in our three models M-streaming, M-turbulence, and M-primaries to , , and , respectively. Following the Mach number ()-dependence of the acceleration efficiency suggested in Pinzke et al. (2013), the efficiency in weak shocks () that dominates the CR distribution function, has an acceleration efficiency for protons , and for electrons .
Conclusions. The standard reacceleration model for RHs requires a population of seed electrons to undergo turbulent reacceleration. These seeds are generally thought to be secondary electrons from hadronic CRp interactions. In this work we use cosmological simulations to derive a population of seed CRps originating from structure formation shocks and merger shocks during the cluster build up. The resulting secondary population is inconsistent with RH observations. We propose 3 possible solutions that all produce gamma-ray emission below current upper limits and that reproduce both the spectrum and the surface brightness profiles of the Coma RH: (i) injected turbulence that is flatter than in previously adopted models, (ii) streaming CRs, (iii) shock accelerated CR electrons injected with . We will pursue further implications and distinguishing characteristics of these competing models in future work.
Acknowledgments. We thank Josh Wiener for discussions on CR streaming. We are also grateful to Lawrence Rudnick for discussion on uncertainties in the 1.4 GHz radio data. A.P. is grateful to the Swedish research council for financial support. S.P.O. thanks NASA grant NNX12AG73G for support. C.P. gratefully acknowledges financial support of the Klaus Tschira Foundation.
- Lawrence Rudnick, private communication.
- G. Brunetti and T. W. Jones, International Journal of Modern Physics D 23, 1430007 (2014).
- F. Govoni, T. A. Enßlin, L. Feretti, and G. Giovannini, A&A 369, 441 (2001).
- L. Feretti, G. Giovannini, F. Govoni, and M. Murgia, AAPR 20, 54 (2012).
- Planck Collaboration, P. A. R. Ade, N. Aghanim, M. Arnaud, M. Ashdown, F. Atrio-Barandela, J. Aumont, C. Baccigalupi, A. Balbi, A. J. Banday, et al., A&A 554, A140 (2013).
- F. Miniati, D. Ryu, H. Kang, and T. W. Jones, Astrophys. J. 559, 59 (2001).
- C. Pfrommer, MNRAS 385, 1242 (2008).
- G. Brunetti, G. Setti, L. Feretti, and G. Giovannini, MNRAS 320, 365 (2001).
- G. Brunetti and A. Lazarian, MNRAS 378, 245 (2007).
- V. Petrosian and W. E. East, Astrophys. J. 682, 175 (2008).
- D. O. Chernyshov, V. A. Dogiel, and C. M. Ko, Astrophys. J. 759, 113 (2012).
- C. L. Sarazin, Astrophys. J. 520, 529 (1999).
- G. Brunetti and A. Lazarian, MNRAS 412, 817 (2011).
- G. Brunetti, P. Blasi, O. Reimer, L. Rudnick, A. Bonafede, and S. Brown, MNRAS 426, 956 (2012).
- A. Pinzke and C. Pfrommer, MNRAS 409, 449 (2010).
- F. Vazza, C. Gheller, and M. Brüggen, MNRAS 439, 2662 (2014).
- T. Enßlin, C. Pfrommer, F. Miniati, and K. Subramanian, A&A 527, A99+ (2011).
- J. Wiener, S. P. Oh, and F. Guo, MNRAS 434, 2209 (2013).
- A. Pinzke, S. P. Oh, and C. Pfrommer, MNRAS 435, 1061 (2013).
- C. Pfrommer, T. A. Enßlin, V. Springel, M. Jubelgas, and K. Dolag, MNRAS 378, 385 (2007).
- X. Guo, L. Sironi, and R. Narayan, Astrophys. J. 794, 153 (2014a).
- X. Guo, L. Sironi, and R. Narayan, Astrophys. J. 797, 47 (2014b).
- D. Caprioli and A. Spitkovsky, Astrophys. J. 783, 91 (2014).
- H. Kang and D. Ryu, Astrophys. J. 734, 18 (2011).
- F. Vazza and M. Brüggen, MNRAS 437, 2291 (2014).
- T. Arlen, T. Aune, M. Beilicke, W. Benbow, A. Bouvier, J. H. Buckley, V. Bugaev, K. Byrum, A. Cannon, A. Cesarini, et al., Astrophys. J. 757, 123 (2012).
- M. Ackermann, M. Ajello, A. Albert, A. Allafort, W. B. Atwood, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, K. Bechtol, et al., Astrophys. J. 787, 18 (2014).
- H. Yan and A. Lazarian, Physical Review Letters 89, B1102+ (2002).
- A. J. Farmer and P. Goldreich, Astrophys. J. 604, 671 (2004).
- E. T. Lau, A. V. Kravtsov, and D. Nagai, Astrophys. J. 705, 1129 (2009).
- L. D. Shaw, D. Nagai, S. Bhattacharya, and E. T. Lau, Astrophys. J. 725, 1452 (2010).
- N. Battaglia, J. R. Bond, C. Pfrommer, and J. L. Sievers, Astrophys. J. 758, 74 (2012).
- A. Bonafede, L. Feretti, M. Murgia, F. Govoni, G. Giovannini, D. Dallacasa, K. Dolag, and G. B. Taylor, A&A 513, A30 (2010).
- U. G. Briel, J. P. Henry, and H. Boehringer, A&A 259, L31 (1992).
- S. Brown and L. Rudnick, MNRAS 412, 2 (2011).
- B. M. Deiss, W. Reich, H. Lesch, and R. Wielebinski, A&A 321, 55 (1997).
- R. F. Pizzo, Ph.D. thesis, University of Groningen (2010).