# Effects of self-generated turbulence on Galactic Cosmic Ray propagation and associated diffuse -ray emission

###### Abstract

Recent results obtained by analyzing diffuse -ray emission detected by Fermi-LAT show a substantial variation of the CR spectrum as a function of the distance from the Galactic Center. For energies up to tens of GeV, the CR proton density in the outer Galaxy appears to be weakly dependent upon the galactocentric distance while the density in the central region of the Galaxy was found to exceed the value measured in the outer Galaxy. At the same time, Fermi-LAT data suggest a gradual spectral softening while moving outward from the center of the Galaxy to its outskirts. These findings represent a challenge for standard calculations of CR propagation based on assuming a uniform diffusion coefficient within the Galactic volume. Here we present a model of non-linear CR propagation in which transport is due to particle scattering and advection off self-generated turbulence. We will show that for a realistic distribution of CR sources following the spatial distribution of supernova remnants and the space dependence of the magnetic field on galactocentric distance, both the spatial profile of CR density and the spectral softening can easily be accounted for.

###### keywords:

cosmic rays: general, gamma rays: diffuse background, ISM: general^{†}

^{†}volume:

00 \journalnameNuclear and Particle Physics Proceedings \runauthG. Morlino \jidnppp \jnltitlelogoNuclear and Particle Physics Proceedings

## 1 Introduction

One of the key aspect of the Cosmic Ray (CR) physics concerns their diffusive behavior. Traditionally the average diffusion coefficient in the Galaxy is determined from local measurements of the ratio between secondary and primary CRs (mainly Boron/Carbon). The inferred diffusion coefficient is then assumed to hold in the whole Galactic volume and is often used to interpret the diffuse -ray emission resulting from interactions of CRs with the interstellar medium (ISM).

The diffusion properties of CRs is the Galaxy could be tightly connected with the “radial gradient problem”, concerning the dependence of CR intensity on Galactocentric distance: the CR density as measured from -ray emission in the Galactic disc is much more weakly dependent upon galactocentric distance than the spatial distribution of the alleged CR sources, modeled following pulsar and supernova remnant (SNR) catalogues. This result was obtained for the first time from the analysis of the -ray emissivity in the Galactic disk derived from the COS-B data (Bhat86; Bloemen86) and then confirmed by later work based on EGRET data (Hunter97; Strong96). In more recent years the existence of a “gradient problem” in the external region of our Galaxy has also been confirmed by data collected by Fermi-LAT (Ackermann11; Ackermann12). The results of a more detailed analysis of Fermi-LAT data, including also the inner part of the Galaxy, was published in two independent papers (Acero16; Yang16) based on data accumulated over seven years. This study highlighted a more complex situation: while in the outer Galaxy the density gradient problem has been confirmed once more, the density of CR protons in the inner Galaxy turns out to be appreciably higher than in the outer regions of the disk. Moreover this analysis suggests a softening of the CR spectrum with Galactocentric distance, with a slope ranging from 2.6, at a distance of kpc, to 2.9 in the external regions.

The emerging scenario is difficult to reconcile with the standard approach to CR propagation, which is based upon the assumption that the diffusive properties are the same in the whole propagation volume (see, e.g., Berezinskii90). Within the context of this approach, several proposals have been put forward to explain the radial gradient problem. Among them: a) assuming a larger halo size or b) a flatter distribution of sources in the outer Galaxy (Ackermann11); c) accounting for advection effects due to the presence of a Galactic wind (Bloemen93); d) assuming a sharp rise of the CO-to-H ratio in the external Galaxy (Strong04); e) speculating on a possible radial dependence of the injected spectrum (Erlykin16). None of these ideas, taken individually, can simultaneously account for both the spatial gradient and the spectral behavior of CR protons. Moreover, many of them have issues in accounting for other observables (see, e.g., the discussion in Evoli12).

A different class of solutions invoke the breakdown of the hypothesis of a spatially constant diffusion coefficient. For instance, Evoli12 proposed a correlation between the diffusion coefficient parallel to the Galactic plane and the source density in order to account for both the CR density gradient and the small observed anisotropy of CR arrival directions. Gaggero15a; Gaggero15b followed the same lines of thought and showed that a phenomenological scenario where the transport properties (both diffusion and convection) are position-dependent can account for the observed gradient in the CR density. It is however unsatisfactory that these approaches do not provide a convincing physical motivation for the assumed space properties of the transport parameters.

In the present paper, following the results presented in (Recchia2016b), we discuss the possibility that diffusion and advection in self-generated waves produced by CR-streaming could play a major role in determining the CR radial density and spectrum. The effects of self-generated diffusion has been shown to provide a viable explanation to the hardening of CR proton and helium spectra observed by PAMELA (Adriani2011) and AMS-02 (Aguilar15), supporting the idea that below GeV, particle transport at the Sun’s location may be dominated by self-generated turbulence (Aloisio15; 2012PhRvL.109f1101B). We suggest that this could be the case everywhere in the Galaxy and explore the implications of this scenario. In the assumption that the sources of Galactic CRs trace the spatial distribution of SNRs and that the magnetic field drops at large galactocentric distances, the density of CRs and their spectrum are well described if CRs are allowed to diffuse and advect in self-produced waves.

The paper is structured as follows. In § 2 we discuss the importance of self-generation while in § 3 we develop a solution of the CR transport in a 1D slab model with constant advection and with purely self-generated diffusion. In § 4 we discuss the distribution of sources and the behavior of the magnetic field in the Galactic plane and we compare our results for the CR proton spectrum with the data obtained by Fermi-LAT . Finally we summarize in § 5.

## 2 The role of self-generated turbulence

Since the first works on the shock acceleration theory (Skilling1975; Bell1978) it became clear that the Alfvén waves self generated by streaming of CRs should dominate the CR scattering in the acceleration region around the shock. In fact, in absence of self generation, the maximum energy reached by particles would be ridiculously low (Lagage-Cesarsky1983) and the shock acceleration could not be invoked to explain the CR spectrum. When magnetic perturbations are weak (), one can use quasi-linear theory to determine the diffusion coefficient in the direction parallel to the large scale magnetic field , which is usually written as

(1) |

where is the Bohm diffusion coefficient, with the Larmor radius and the particle’s speed. is the normalized energy density per unit logarithmic wavenumber , calculated at the resonant wavenumber . In the limit , the diffusion perpendicular to the field lines can be neglected, being suppressed by a factor of with respect to the one parallel to (see, e.g. Drury1983), hence the problem usually reduces to one spatial dimention.

Always in quasi-linear theory the growth rate due to the CR-streaming instability is (Skilling1975):

(2) |

where the Alfvén speed is . It is worth stressing two relevant properties of equation (2): the first is that , implying that the amplification is more effective in regions where the is small; the second is that the growth rate is proportional to the CR spatial gradient along the direction of . This implies that wherever the CR distribution is not uniform the amplification is turned on. Hence the importance of self-amplification is not limited to the case of shocks and, in fact, in recent studies it has been proven to be relevant in several contexts. Among them: a) CR escaping from Galactic sources Malkov2013; Nava2016; Dangelo2016; b) propagation of CRs close to molecular clouds Skilling-Strong1976; Cesarsky-Volk1978; Morlino-Gabici2015; c) diffusion in the Galactic halo Blasi-Amato-Serpico2012; Recchia2016a; Recchia2016b; d) CR escaping from galaxies Blasi-Amato-Dangelo2015.

In all those situations, the final amplitude of amplified waves can be reduced by their damping which can take place through different mechanisms whose relative importance depends on the properties of the plasma we are considering. A strong damping mechanism is due to ion-neutral collisions which transfer energy from ions (which oscillates with the waves) to neutrals (which do not oscillate) (Kulsrud-Cesarsky1971; Drury-Duffy-Kirk1996). Here we are mainly interested in describing the diffusion in the Galactic halo where the density of neutral atoms is small enough that the ion-neutral damping can be safely neglected.

The dominant damping process in a region where the background gas is totally ionized, is the non-linear Landau Damping (NLLD) due to wave-wave interactions. It occurs at a rate (Zhou90):

(3) |

with . In the rest of the paper we consider only the NLLD, even if another relevant mechanism has been proposed by Farmer and Goldreich (farmer2004). They suggested that the waves generated by streaming instability interact with oppositely directed turbulent wave packets characterizing a pre-existing MHD turbulence and, as a consequence, the wave energy dissipates through a cascade towards smaller scales. The analytical expression for the FG damping rate is

(4) |

where is the scale at which the turbulence is injected. Here we neglect such damping because we are not accounting for any external turbulence. Nevertheless, we warn the reader that in realistic cases the FG mechanism could dominate the damping rate when the level of self-generated turbulence drops below some threshold. In fact, assuming that is of the order of the coherence scale of the background magnetic field, i.e. pc (beck2016), we can estimate the relative importance of with respect to comparing Eq.(4) with Eq.(3). One gets that if (where the resonant condition has been used). It is clear that for large enough energies the FG damping becomes the dominant one but in this work we are mainly interested in particles with energy up to few tens of GeV where is at most comparable with .

## 3 CR transport in the Galactic halo with self-generated turbulence

The simplified cylindrical model for the Galaxy that we use here is sketched in Figure 1. We assume a purely poloidal magnetic field orthogonal to the Galactic plane, , where the strength is constant along but can be a function of the Galactocentric distance, . Hence, for any given , the transport can be described as one-dimensional, so that particles diffuse and advect only along the direction. The transport equation for CR protons can then be written as follows:

(5) |

where the advection is only due to the motion of Alfvén waves directed away from the disk, i.e. . We further assume that the ion density is constant everywhere in the halo, which extends out to .

The injection of particles occurs only in the Galactic disk (at ) and is a power law in momentum:

(6) |

Here erg is the total kinetic energy released by a single supernova explosion, is the fraction of such energy channeled into CRs, is the SN explosion rate per unit area and is a function of the galactocentric distance. Finally, the normalization constant is .

A standard technique to solve the transport equation (5) is to integrate between and and between and with the boundary conditions and (see, e.g., 2012PhRvL.109f1101B; Recchia2016a). One gets the following result for :

(7) |

where and the distribution function in the disk is:

(8) |

This solution is only implicit because the diffusion coefficient is a function of through the self generation mechanism, Eq.(2). The local value of is determined by its transport equation where the growth of the waves is balanced by their damping and advection along ,

(9) |

As discussed in § 2 we only consider the NLLD. In general the damping and the growth are much faster than the wave advection at the Alfvén speed, as we will show in a while. Hence the rhs in equation (10) can be neglected and a good approximation for the wave distribution is obtained by equating with , which returns the following implicit form for the wave spectrum:

(10) |

Inserting the diffusive flux, (from equation (7)) into equation (10), we can derive the diffusion coefficient using equation (1):

(11) |

where is the diffusion coefficient at and reads:

(12) |

It follows that the diffusion coefficient is maximum at and decreases linearly with .

Now, we can verify a posteriori that the advection term in Eq.(10) is always negligible. Using Eq.(11) and recalling that , we have , hence the ratio between the advection and the damping terms can be written as

(13) |

which, for any realistic value of the Alfvén speed in the Galactic halo, is always for .

The exact solutions, equations (7)-(3) and (11)-(12), can be written in an explicit form in the two opposite limits of diffusion-dominated and advection-dominated transport. In particular, it is straightforward to verify that in the diffusion dominated case (i.e. when ) the standard solution is recovered. In fact in this limit and becomes constant in , namely

(14) |

In the same limit equation (7) reduces to the well known result

(15) |

Replacing this last expression for into equation (14) we obtain explicit expressions for both and , which read

(16) |

and

(17) |

respectively. Notice that Eq.(16) is telling us that the self-generated turbulence is given by , where is the energy flux injected into CRs and is the energy density of the background magnetic field.

In the opposite limit, when , we have that , hence and from equation (3) we recover:

(18) |

On the other hand, we see from equation (11) that the diffusion coefficient in the disk behaves like a constant in momentum, namely . This happen because for small , , hence . Clearly this dependence is restricted to the momenta for which diffusion becomes comparable to advection, typically below GeV/c (see below). We refer to this regime as advection dominated regime, although particles never reach a fully advection dominated transport because diffusion and advection time are of the same order.

We are interested in describing the dependence of the CR spectrum on the Galactocentric distance, which enters the calculation only through the injection term and the magnetic field strength. Comparing equations (17) and (18) we see that the CR density has the following scalings with quantities depending on :

(19) |

In a more general case, equations (3), (11) and (12) can be solved numerically using an iterative technique. We start by choosing a guess function for (for instance the expression given by equation (16), obtained without advection) and then we iterate until convergence is reached, a procedure which usually requires only few iterations. Notice that the general case of a transport equation (5) where the advection speed may depend on the -coordinate has been recently discussed by (Recchia2016a) and used to describe CR-induced Galactic winds. For the sake of simplicity, and to retain the least number of parameters, here we assume that the advection velocity is simply the Alfvén speed of self-generated waves and we assume that is independent of .

## 4 Results

### 4.1 Fitting the local CR spectrum

The rate of injection of CRs per unit surface can be calibrated to reproduce the energy density and spectrum of CRs as observed at the Earth. In all our calculations, following most of current literature, we choose a size of the halo kpc, while the ion density in the halo is fixed as cm, a value consistent with the density of the warm ionized gas component (see, e.g, Ferriere01). The magnetic field at the location of the Sun is assumed to be G (since we are only describing the propagation in the direction, this should be considered as the component of the field perpendicular to the disc). Notice that, given and , the value of the Alfvén speed is also fixed, and this is very important in that it also fixes the momentum where the transition from advection propagation to diffusion dominated propagation takes place for a given injection spectrum and product of injection efficiency times the local SN explosion rate, . Following 2012PhRvL.109f1101B; 2013JCAP...07..001A; Aloisio15 we adopted a slope at injection . Then, by requiring that the local CR density at GeV is equal to the observed one, we get .

It is worth stressing that the CR spectrum in the energy region GeV, may be heavily affected by either pre-existing turbulence 2012PhRvL.109f1101B; 2013JCAP...07..001A; Aloisio15 or a -dependent diffusion coefficient (Tomassetti). Both possibilities have been proposed to explain the spectral hardening observed in both the protons’ and helium spectrum at rigidities above GV. For this reason, a model including only self-generated diffusion can be considered as reliable only below GeV. In the following, we limit our attention to CRs that are responsible for the production of -rays of energy GeV, as observed by Fermi-LAT , namely protons with energy of order GeV. This threshold is sufficiently low that the slope derived by ignoring the high energy spectral hardening can be considered reliable.

The injection parameters (efficiency and spectrum) found by fitting the CR density and spectrum at the Sun’s location are assumed to be the same for the whole Galaxy. As discussed in §4.2, the rate of injection of CRs per unit surface is then proportional to the density of SNRs which is, in turn, inferred from observations.

### 4.2 CR spectrum in the Galactic disk

The SNR distribution is usually inferred based on two possible tracers: radio SNRs and pulsars. Here we adopt the distribution of SNRs recently obtained by Green15 from the analysis of bright radio SNRs. He adopted a cylindrical model for the Galactic surface density of SNRs as a function of the Galactocentric radius, in the form:

(20) |

where the position of the Sun is assumed to be at kpc. For the best fit Green15 obtained and , so that the distribution is peaked at kpc. However, as noted by Green15, it is worth keeping in mind that the best-fitting model is not very well defined, as there is some level of degeneracy between the parameters and .

In a previous work CB98 largely cited in the literature, the authors also adopted a fitting function as in equation (20) but obtained their best fit for and , resulting in a distribution peaked at kpc and broader for larger values of with respect to the one of Green15. In CB98 the source distances is estimated using the so called ‘–D’ relation, that is well known to be affected by large uncertainties. Moreover in Green15 it is argued that the –D used by CB98 appears to have been derived incorrectly.

An important caveat worth keeping in mind is that the SNR distribution derived in the literature is poorly constrained for large galactocentric radii. For instance Green15 used a sample of 69 bright SNRs but only two of them are located at galactic latitude . Similarly, CB98 used a larger sample with 198 SNRs, but only 7 of them are located at kpc and there are no sources beyond 16 kpc.

The distribution of pulsars is also expected to trace that of SNRs after taking into account the effect of birth kick velocity, that can reach km/s. These corrections are all but trivial, (see, e.g. Faucher06), hence in what follows we adopt the spatial distribution as inferred by Green15.

One last ingredient needed for our calculation is the magnetic field strength, , as a function of galactocentric distance . While there is a general consensus that the magnetic field in the Galactic disk is roughly constant in the inner region, in particular in the so-called “molecular ring”, between 3 and 5 kpc (Jansson12; stanev97), much less is known about what the trend is in the very inner region around the Galactic center, and in the outer region, at kpc. Following the prescription of Jansson12 (see also stanev97), we assume the following radial dependence:

(21) |

where the normalization is fixed at the Sun’s position, that is G. Using this prescription we calculate the CR spectrum as a function of the Galactocentric distance, as discussed in §3. In the left panel of Figure 2 we plot the density of CRs with energy GeV (dashed green line) and compare it with the same quantity as derived from Fermi-LAT data. Our results are in remarkably good agreement with data, at least out to a distance of kpc. At larger distances, our predicted CR density drops faster than the one inferred from data, thereby flagging again the well known CR gradient problem. In fact, the non-linear theory of CR propagation, in its most basic form (dashed line) makes the problem even more severe: where there are more sources, the diffusion coefficient is reduced and CRs are trapped more easily, but where the density of sources is smaller the corresponding diffusion coefficient is larger and the CR density drops. A similar situation can be seen in the trend of the spectral slope as a function of , plotted in the right panel of the same Figure. The dashed line reproduces well the slope inferred from Fermi-LAT data out to a distance of kpc, but not in the outer regions where the predicted spectrum is steeper than observed. It is important to understand the physical motivation for such a trend: at intermediate values of , where there is a peak in the source density, the diffusion coefficient is smaller and the momenta for which advection dominates upon diffusion is higher. This implies that the equilibrium CR spectrum is closer to the injection spectrum, (harder spectrum). On the other hand, for very small and for large values of , the smaller source density implies a larger diffusion coefficient and a correspondingly lower momentum where advection dominates upon diffusion. As a consequence the spectrum is steeper, namely closer to . In fact, at distances kpc, the spectrum reaches the full diffusive regime, hence , meaning that the slope in Figure (2), right panel, is 3.6. As pointed out in §3, the non-linear propagation is quite sensitive to the dependence of the magnetic field on .

Both the distribution of sources and the magnetic field strength in the outer regions of the Galaxy are poorly known. Hence, we decided to explore the possibility that the strength of the magnetic field may drop faster than at large galactocentric distances. As a working hypothesis we assumed the following form for the dependence of on R, at kpc:

(22) |

where the scale length, , is left as a free parameter. We found that using kpc, both the resulting CR density and spectral slope describe very well the Fermi-LAT data in the outer Galaxy. The results of our calculations for this case are shown in Fig. 2 with solid red lines.

The diffusion coefficient in the Galactic plane resulting from the non-linear CR transport in the Galaxy, calculated as in §3, is illustrated in Fig. 3, for different galactocentric distances. It is interesting to notice that at all values of (and especially at the Sun’s position) is almost momentum independent at GeV/c. This reflects the fact that at those energies the transport is equally contributed by both advection and diffusion, as discussed above. This trend, that comes out as a natural consequence of the calculations, is remarkably similar to the one that in numerical approaches to CR transport is imposed by hand in order to fit observations.

Contrary to a naive expectation, in the case in which drops exponentially, the diffusion coefficient becomes smaller in the external Galaxy than in the inner part, in spite of the smaller number of sources in the outer Galaxy. This counterintuitive result is due to the fact that the growth rate due to streaming instability is while the damping rate . As a consequence the amplification is more effective in region where the background magnetic field has a smaller strength. Such trends immediately reflect in the result of equation (16), i.e. , hence if drops faster than , the diffusion coefficient decreases at large , as occurs in or model. Clearly, this result loses validity when approaches unity and the amplification enters the non linear regime. Using equation (11), such condition in the disk can be written as which, for 1 GeV particles occurs for kpc (red-dashed lines in Figure (2)). In any case, the density of CRs at large galactocentric distances drops down.

## 5 Conclusions

Understanding the transport properties of CRs in the Galactic halo as well as in the Galactic disc remain an open issue in the CR physics and is of paramount importance to correctly interpret the diffuse -ray emission. In the standard picture, largely used in the literature, a too simplistic model, where the diffusion coefficient is constant everywhere in the Galaxy, is used. The CR density recently inferred from Fermi-LAT observations of the diffuse Galactic -ray emission appears to be all but constant with galactocentric distance Acero16; Yang16. In the inner kpc from the Galactic center, such density shows a pronounced peak around kpc, while it drops with for kpc, but much slower than what one would expect based on the distribution of SNRs, as possible sources of Galactic CRs. Moreover, the inferred slope of the CR spectrum shows a gradual steepening in the outer regions of the Galaxy. This puzzling CR gradient is hard to accommodate if the diffusion were the same in the whole Galaxy.

Here we showed that both the gradient and the spectral shape of CR spectrum can be explained in a simple model of non-linear CR transport: CRs excite waves through streaming instability in the ionized Galactic halo and are advected with such Alfvén waves. In this model, the diffusion coefficient is smaller where the source density is larger and this phenomenon enhances the CR density in the inner Galaxy. In the outer Galaxy, the data can be well explained only by assuming that the background magnetic field drops exponentially at kpc, with a suppression scale of kpc. This scenario also fits well the spectral slope of the CR spectrum as a function of , as a result of the fact that at different the spectrum at a given energy may be dominated by advection (harder spectrum) or diffusion (softer spectrum). Our analysis is limited to CR energy GeV because data on the -ray diffuse emission are available mainly for photon energies few GeV. The knowledge of the spectral shape at larger energies would be especially interesting. In fact a simple prediction of our calculations is that the spectral hardening should disappear at GeV, where transport is diffusion dominated at all galactocentric distances.

## References

- A LaTeX macro package for message sequence charts—maintenance document—describing \mscpack version \mscversion. Note: Included in MSC macro package distribution
- A LaTeX macro package for message sequence charts—reference manual—describing \mscpack version \mscversion. Note: Included in MSC macro package distribution
- A LaTeX macro package for message sequence charts—user manual—describing \mscpack version \mscversion. Note: Included in MSC macro package distribution
- The LaTeX Graphics Companion. Addison-Wesley.
- ITU-TS Recommendation Z.120: Message Sequence Chart (MSC). Geneva.
- LaTeX—a document preparation system—user’s guide and reference manual. 2nd edition, Adsison-Wesley. Note: Updated for LaTeX2e
- Tutorial on message sequence charts (MSC’96). In FORTE,