4 Transport of cosmic rays in the Galaxy

Transcription

1 4 Transport of cosmic rays in the Galaxy Note, the following discussion follows quite closely the Wefel overview article in Cosmic rays, supernovae and the interstellar medium, eds. M.M. Shapiro, R. Silberberg, and J.P. Wefel, p44, Dordrecht: Kluver Academic Publishers. A summary of data is given in Simpson s article from Many of the calculations are also discussed in Sections 9 and 20 of the Longair book. More recent data in the review by Strong A.W., Moskalenko, I.V., and Ptuskin, V.S. Ann.Rev.Nucl.Part.Sci. 57, , 2007, as preprint arxiv:astro-ph/ In order to set the stage, let us review what we know from measurements of the local cosmic-ray population: 1. The spectral energy distribution of the all-particle spectrum in dn = N(0) (E/E 0 ) Γ with Γ 2.7. The broad-band power-law stretches from the energy of the heliospheric cut-off up to the so-called knee at an energy of a few PeV (10 15 ev). At higher energies, the energy spectrum softens to Γ 3.1 (see Fig. 1). 2. The arrival direction is (almost) isotropic. Up to knee energies, the isotropy is good within 0.1 %, while at higher energies, the anisotropy increases, reaching a value of about 1 % at energies of ev (see Fig. 2). 3. The composition (chemical abundance) of the cosmic rays arriving at Earth is at low energies similar (but not equal) to the solar metallicity. The cosmic-ray abundance of rare elements like Lithium (Li), Berrylium (Be), and Boron (B) is significantly higher than the solar abundance (see Fig. 3). Furthermore, isotopic composition deviates significantly in the neutron rich isotopes. 4. The cosmic-ray composition changes with energy and the position of the knee seems to be correlated with mass or charge (see Fig. 4). These are the most important observational findings. We will re-visit these issues in the following and derive an ansatz for an interpretation. We are mainly interested in the following questions: What is the confinement volume of cosmic-rays (as a function of energy)? 1

2 Energies and rates of the cosmic-ray particles E 2 dn/de (GeV cm -2 sr -1 s -1 ) 10 0 protons only electrons positrons antiprotons all-particle CAPRICE AMS BESS98 Ryan et al. Grigorov JACEE Akeno Tien Shan MSU KASCADE CASA-BLANCA DICE HEGRA CasaMia Tibet Fly Eye Haverah Yakutsk AGASA HiRes E kin (GeV / particle) Figure 1: Energy spectrum of cosmic rays from a recent compilation. Note, the energy dependent scaling on the flux. What is the confinement time of cosmic-rays (as a function of energy)? 2

3 Figure 2: Observed relative anisotropy of cosmic-rays (Ambrosio et al. PRD 2003), compared with predictions. What is the total power required to sustain the cosmic-ray population? Note, these questions are mainly related to the nucleonic component of cosmic-rays. The origin of primary cosmic-ray electrons (and possibly positrons) requires a different treatment 1 We can identify different approaches to address the questions given above: Analysis of cosmic-ray spallation leading to the observed secondary nuclei Be, Li, B Analysis of radio-isotopes (and other rare isotopes) and the relative abundance of cosmic-ray clocks 1 The main difference is the importance of radiative losses for electrons which can be neglected for cosmic-ray nuclei. 3

4 Figure 3: Comparison of solar abundances (open symbols) and cosmic-ray abundances (filled symbols). Note, the logarithmic scale. Measurement of Photons and neutrinos from cosmic-ray interaction. Again, cosmic-ray electrons require a separate treatment as a consequence of radiative losses. The overall scheme on how to proceed is shown schematically in the cartoon (Fig. 5). 4

5 Figure 4: Compilation of broad band energy spectra for H, He, Fe (top to bottom). The compilation combines direct measurements with indirect (air shower based) measurements. Taken from Hörandel (2005). 5

6 Figure 5: Cartoon of the cosmic ray transport in our Galaxy- taken from Moskalenko et al Origin of light elements (Li, Be, B) in cosmic-rays Obviously, the cosmic-ray interaction with the interstellar medium requires some consideration of nuclear physics. The interaction of a cosmic-ray nucleus with a hydrogren atom (and to some extent and much rarer an interaction with a He or a metall in the interstellar medium) leads to a spallation reaction in which the original nucleus breaks up in potentially excited nuclear fragments. The cross section for these processes are in principle accessible to lab experiments, however, the data-base of spallation cross sections is not complete. As it turns out, some of the analysis of cosmic-ray data is limited by the lack of available cross section data. For now, we assume, that cross sections for processes like σ(c + p He + X) are known or phenomenologically calculated. The general equation for a particle distribution of type i with uniform spatial diffusion (neglecting convection and any form of momentum diffusion) is N i t = D 2 N i + E (b in i ) + Q i N i + P ji N j. (1) τ i j>i τ j 6

7 In the following, we will argue to use a simplified version of this equation. The energy losses are negligible for typical gas densities of n = 1 cm 3 and for the following discussion, we neglect the source term. This corresponds to an instantaneous injection at t = 0. Furthermore, we neglect the diffusion which would correspond to a homogeneous distribution in the considered medium. Generally, this is not valid, but sufficient for the purpose of demonstrating the basic properties of the effect of spallation (and the limitation of this approximation). Furthermore, we change the variable from time to the traversed column density ξ = ρ x = ρ v t. In this case, after substituting, we arrive at the following equation: N i ξ = N i ξ i + j>i P ji ξ j N j. (2) In this case, we assume, that all nuclei traverse the same column density (slab), this approximation is therefore known as the single slab model (we will see shortly, that this model is too simple). For now, let us consider the group of light (L) nuclei (Li, Be, B) and the most abundant group of medium (M) nuclei (C, N, O). The boundary condition can be simply set to N L (ξ = 0) = 0, because the light nuclei are not enriched by stellar fusion processes (Li is only produced in the primordial nucleosynthesis). For the M-group, we ignore the additional contribution from spallation of heavier elements (we will see shortly, that this assumed decoupling is a reasonable approximation): dn M dξ dn L dξ = N M ξ M (3) = N L ξ L + P ML ξ M N M. (4) The solution to Eqn.3 with the boundary condition that N M (ξ = 0) = N M (0) is N M (ξ) = N M (0) exp( ξ/ξ M ). (5) We can find a solution to Eqn.4 after substituting Eqn.5, expanding with exp(ξ/ξ L ), and integrating (exercise). The ratio of the groups is then simply: N L (ξ) N M (ξ) = P [ ( ML ξ L ξ exp ξ ) ] 1. (6) ξ L ξ M ξ M ξ L 7

8 Before comparing with the measurements, let us consider briefly which values to use for ξ L, ξ M, P ML. The column density ξ is related to the integral cross section σ and the molar mass m A : ξ = m A N A σ, (7) with N A the Loschmidt/Avogadro number ( mol 1 ). As an example, consider a typical cross section of 30 mbarn for inelastic pp-scattering. In this case, ξ pp = 55 g/cm 2. The total cross section for the spallation of an element of the M-group is simply the sum of all partial cross sections (e.g. σ(c + p B + X) etc.). The values of the partial cross sections are listed in Table 1 in units of mbarn. In order to calculate the values for ξ M and P ML, we calculate the abundance weighted cross section. The abundances are the actual values measured for the individual species C, N, and O (using the values from Table 2). The averaged cross section is simply σ = i w i σ i / i w i with the weights taken from the relative abundance ( MeV/n) and σ i the total cross section for species i {C, N, O}. The resulting values are σ M = 280 mbarn and for σ L = 200 mbarn 2. The probability for production of an element from the L-group is estimated by taking the ratio of P ML = i {Li,Be,B} σ Mi / σ M = The measurement (see also Table 2) indicates for the ratio N L N M = 0.25 ξ = 4.8 g cm 2. (8) This result is obviously consistent with the naive expectation that the column density is of similar magnitude as ξ M. When looking into the relative abundance of the individual elements in the L-group ([Li] = 136, [Be] = 67, [B] = ), the ratio of elements is consistent with the respective weighted production cross sections (σ(m Li) = 24 mbarn, σ(m Be) = 16.4 mbarn, σ(m B) = 35 mbarn). This is reassuring and shows that the simple, one slab model is sufficient to explain the observed abundance of light elements. The same exercise can be done to estimate the slab thickness for the production of 3 He which is again roughly 5 g/cm 2. 2 note, that the values for L are not included in Table 1 3 The abundances are always considered relative to the abundance of Si 8

9 The simple one slab model works surprisingly well for the spallation of the M-group elements into the light nuclei. However, when looking into the details of the spallation of iron nuclei, a simple one slab model does not work. The spallation cross section for iron nuclei is very high (see Table 1): σ F e = mbarn, this corresponds to ξ F e = 2.2 g/cm 2. If the slab thickness were constant at 5 g/cm 2, the iron nuclei abundance observed in cosmic rays should be strongly depleted with respect to the average abundance. Independent of the absolute value of the abundance, the ratio of spallation products to the observed abundance of iron in cosmic rays should be strongly dominated by the spallation products: [products] [primaries] = 1 exp( ξ/ξ F e) exp( ξ/ξ F e ) = 8.7 (9) for ξ = 5 g/cm 2. This is clearly not consistent with the observed abundances: The spallation of iron results mainly in the production of nuclei in the range from Cl...V (see also the partial cross sections listed in Table 1). Using the values given in the appendix, the observations are [Cl,..., V ] [F e] = 1.5 (10) clearly in contradiction with the higher value expected. The discrepancy of the simple one slab model can be resolved, when treating the transport of cosmic rays more realistically. There are in principle two ways of treating the cosmic ray transport. The most simple approach is to consider a distribution of path lengths and try to find the distribution that matches the observed ratio of secondaries to primaries (phenomenological approach). A different approach is based upon a more elaborate, physics oriented model of cosmic ray transport in the Galaxy which then in return results in a prediction of the path length distribution (theoretical approach). In principle, both approaches have strengths and short comings. Let us for the sake of clarity consider two extreme cases for the distribution of path lengths: 1. Leaky box: In this approach, the confinement volume of free streaming cosmic rays is defined by a boundary region with a finite escape probability. In this case the entire population of cosmic rays can be described by N t = N τ e (E). (11) 9

10 In the most simple case, the escape time would be energy independent (we will see later, that the data suggest an energy dependence) which would simplify the solution to be N exp( t/τ e ) exp( ξ/ξ e ). Exponential path length distribution. 2. Diffusion in infinite volume: The other extreme would be an infinite escape time τ e with a diffusive transport: N t = D 2 N (12) with a solution which corresponds to a Gaussian distribution of path lengths. Before comparing the expectations from various models (variants of the theme suggested above), let us summarize the observational results for secondary to primary ratios. The most important measurements of the secondary to primary ratio are the B/C (Boron to Carbon) and (Sc+Ti+V)/Fe-measurements. The main result of both measurements is an energy dependence of the secondary/primary ratio which indicates that the path length distribution changes with energy (see Fig. 6 for a recent compilation of measurements). The energy dependence is such that the secondary/primary ratio increases with increasing energy until it reaches a peak and subsequently for increasing energy it drops. The behaviour for the two measurements is very similar, indicating that the path length distribution is similar for heavy and for medium group. In order to appreciate the relevance of the observations, we take a step back and consider the path length distribution required to match the observational data. Again, we simplify the underlying equation by assuming a steady-state case (i.e. N = 0). In this case, we can simply solve the following equation: t N L ξ e (E) + P ML N M (ξ) N L = 0 (13) ξ L ξ M N L = P ML N M (ξ)/ξ M ξ e (E) 1 ξ 1 L which we simplify by assuming (realistically) ξ e ξ L :, (14) N L (ξ) N M (ξ) = P ML ξe(e) ξ M. (15) 10

11 This effectively means that the secondary-to-primary ratio represents the effective path-length as a function of energy (ξ M is fairly constant over energy). Returning to the measured secondary-to-primary ratio we clearly see that for large energies an exponential drop-off in the path-length with energy is a fairly good description of the model. However, at the low energy end as well as at the high energy end, there are markable differences visible. Let us summarize the main results: The increase of the secondary production with increasing energy at the low energy end is not naively expected. It requires an ad-hoc assumption in the leaky box model. The exponential path length distribution is however a good approximation at high energies. More realistic models have been considered in the literature (see e.g. Jones, Lukasiak, and Ptuskin 2001). The increase in path length is a natural consequence in a model where the cosmic rays are injected in an infinitesimally thin disk with a halo in which cosmic rays move with a convective or turbulent wind. The wind speed can be tuned such that the cosmic rays are removed from the disk before they can interact with the gas in the disk. This naturally leads to the peak in the secondary/primary ratio. The observed drop off in secondary-to-primary ratio indicates that cosmic rays with increasing energy leave the galactic disk and traverse on average a smaller column density on the way to the observer on Earth. The general result of the stable secondaries indicates that the mean path length is of the order of 5 g/cm 2, decreasing with increasing energy roughly as ξ e (E) E When considering a value of 5g/cm 2, what geometrical path length does this amount to? As a rough estimate, consider a medium with average density ρ (in units of g/cm 3 ), the traversed column density relates to the spatial distance travelled x: ξ = x ρ. With a typical gas density of 1 proton/cm 3, the corresponding mass density ρ = g/cm 3. For ξ = 5 g/cm 2, the corresponding path length is x = cm 1 Mpc. When comparing this with the typical radius of the gas disk in our Galaxy of 10 kpc, it is obvious that the propagation of cosmic rays is not rectilinear/free-streaming but is clearly related to a diffusive transport of cosmic-rays in the entangled magnetic field 11

12 Figure 6: Compilation of secondary to primary ratios (from Strong, Moskalenko, and Ptuskin 2008). The left panel shows the B/C ratio representative for the spallation induced light nuclei abundance, the right panel shows the production of heavy nuclei in spallation from iron. of the interstellar medium. Similar numbers can be estimated by considering the propagation time. The particles move roughly with the speed of light. Therefore, the time is given by τ e = n/(1 cm 3 ) yrs. We will compare this number with the estimate derived from the cosmic clocks. 4.2 Cosmic-ray clocks Now, what is the confinement volume? In order to investigate this point, the most elegant approach is the use of cosmic-ray clocks. The best-measured cosmic-clock is the radio active isotope 10 Be. The lifetime of 10 Be is τ( 10 Be) = yrs for the decay of 10 Be 10 B +e + ν e. 10 Be is mainly produced in the spallation of C and O. By measuring the ratio of stable 7 Be to 10 Be, it is possible to measure the average time for the propagation of cosmic-rays from the source to the observer. Quantitatively, we consider again a leaky box model (steady-state). In order to shorten the notation, we define C i := j>i N i τ e (i) + C i N i τ spal. (i) N i τ r (i) 12 P ij τ j N j (16) = 0. (17)

13 The most important change is the introduction of the last term which takes into account the decay of the radio isotopes. The equation for the stable isotopes: Combining both equations: N k τ e (k) + C k N k τ spal. (k) = 0. (18) N( 10 Be) N( 7 Be) = C(10 Be) C( 7 Be) τ e ( 7 Be) 1 + τ spal ( 7 Be) 1, (19) τ e ( 10 Be) 1 + τ spal ( 10 Be) 1 + τ r ( 10 Be) 1 assuming that τ spal τ e, this simplifies to: N( 10 Be) N( 7 Be) = C(10 Be) C( 7 Be) τ e ( 7 Be) 1. (20) τ e ( 10 Be) 1 + τ r ( 10 Be) 1 Inserting the measurements, which show a slight variation with energy (see Fig. 7), the typical value for τ e derived is τ e 10 7 years for [ 10 Be]/[ 7 Be + 9 Be+ 10 Be] When comparing this with the column density traversed, cosmic-rays obviously are on average propagating in a medium which has a density of n 0.3 cm 3, considerably smaller than the average density in the Galactic disk. This indicates that cosmic rays very likely spend a considerable time propagating outside the disk in an extended halo. It is (at this point) not possible to conclude the actual extent of the halo but it certainly extends beyond the scale height of cold molecular gas. 4.3 Comments on particle diffusion Diffusive transport of particles is a well-known phenomena that has been studied intensively in the general context of transportation processes (e.g. heat conductivity). The diffusive transport of cosmic-rays is however a more complicated process as it requires the treatment of charged particles moving through a (partially) ionised medium with magnetic field. This is a highly non-linear problem which we do not introduce here in all its depth. Let us highlight a number of important issus. The general principle expressed in Fick s second law applies to the general problem of diffusive transport: N t = D xx N (21) 13

14 Measured values for the ratio of radioactive Be to stable Be- Figure 7: isotopes. which describes the temporal change of a particle (concentration, density) that is subject to a diffusive (random) transport. The spatial diffusion coefficent D xx, given in units of length 2 time 1, is related to the mean free path length λ: D xx = v λ 3. (22) The mean free path length is commonly estimated to be the distance travelled at which the particle has changed the direction considerably and the general distribution of pitch angles is isotropic. The most common value used is the gyro-radius of a particle moving in a homogeneous magnetic field. Note however, that the motion of a particle is usually not as simple as a gyrating motion. The gyro-radius r g = p/(zeb) can be readily calculated with the following equation: ( ) 1 r G = 0.4 pc Z 1 E B. (23) ev 3 µg The net motion of diffusing particles is towards regions with smaller values of density (see also Eqn. 21). The distance R travelled in a time τ is given 14

15 by R 2 = D τ, (24) which is characteristic for diffusive motion. Before we estimate the diffusion coefficient of cosmic-ray propagation, we discuss the microscopic properties of diffusive transport. Scattering of cosmic-rays on self-generated Alfvén waves: The interstellar medium can be simplified by assuming that a static, background magnetic field B 0 is present throughout the medium perturbed by an additional fluctuating field component B 1. The fluctuation B 1 can be described by a spectrum of turbulences, where the power (B1) 2 is a function of the length scale (λ). For cosmic-rays, the background field B 0 leads to a regular motion of charged particles which singles out a characteristic length scale similar to r G. For λ r G or λ r G, the effect of B 1 is averaged either over many cycles (λ r G ) or there are not sufficient cycles (λ r G ) to affect the motion of the particle. If however, particles and irregularities are close in phase and wave-length, a resonant scattering takes place, where energy is transferred from plasma waves to cosmic-rays or vice versa. The scattering of the particle will lead after one cycle to a change of the angle δϕ B 1 B 0. (25) Since the change of the angle is random, it will take N = (δφ) 2 scattering events to change the angle by one radian. During the N cycles, the particle will have moved ( ) 2 λ sc Nλ r G (δφ) 2 B0 r G. (26) B 1 In the following, we want to show that these irregularities are produced by the cosmic-rays themselves (this is why frequently, cosmic-rays are called self-confining particles). The irregularities turn out to be Alfvén (plasma) waves. Let us turn for a moment to plasma waves. A plasma is defined as a medium which is (partially) ionised. On average (or large scales), the plasma is electrically neutral. In astrophysics, the plasma is very dilute 15

16 such that actual particle-particle collisions are negligible. The long-distance Coulomb interaction however still leads to scattering of particles in electrical and magnetic fields. This type of plasma is often referred to as a collisionless plasma. There is a large number of plasma waves as well as instabilities known. For our purposes, we consider only one particular (and relevant) type of plasma waves which are called Alfvén-waves 4 This particular wave propagates in the plasma with density ρ = N p m p with a velocity (note the use of the Gaussian system of units where B cgs = B SI 4π/µ 0 ) 5 : v A = B 0 4πρ0. (27) For the typical interstellar medium with N p = 1 cm 3 and a magnetic field B 0 = 3 µg, we can calculate v A = 8 km/s. The growth rate of Alfvénic waves depend on the streaming velocity of cosmic-rays. Without deriving the result (see Longair, Section 20.4), the energy density in Alfvén waves is U = U 0 exp(γ t) with a growth factor Γ(λ) = Ω 0 N(> E(λ)) N p ( 1 + v ) v A. (28) The particles N(> E) which resonantly scatter with the wave length λ lead to a growth in the perturbance spectrum for streaming velocities larger than v A. For cosmic-rays N(> E) E 1.7 implies that the time-scale τ = Γ 1 for the growth of Alfvén waves increases with E 1.7. This implies that particles with increasing energy are less efficiently confined. We have so far neglected processed leading to wave-damping, which effectively requires energy to be removed from the waves. Among these processes, interactions with the neutral phase of the interstellar medium lead to a damping on time-scales short in comparison with the growth time scale (Kulsrud and Pierce, 1969). The interstellar medium in the Galactic plane consists of a mixture of cool, 4 Biographical note: Hannes Alfvén (*1908, 1995) was a Swedish Physicist who was awarded the Nobel price in 1970, recognizing his fundamental work in Magnetohydrodynamics and plasma physics. 5 For further details on conversion and also on plasma physics, the naval research lab (NRL) publishes regularly a review, the NRL Plasma Formulary, available at their web site 16

17 mostly molecular but also neutral atomic (hydrogen) gas with embedded regions with a high level of ionisation. Outside the disk, the gas is mostly ionised. This leads to the following, simple picture for cosmic-ray transport in our Galaxy: Inside the (mostly neutral) gas in the Galactic disk, cosmic-rays are mostly free-streaming, possibly drifting along field lines. When encountering a region of highly ionized gas (as e.g. outside the Galactic plane), the cosmic-rays start moving diffusively, with a net streaming velocity approaching the Alfvén velocity. At the interface between the two regions, particle conservation should hold, and therefore: N int c = N ext v A. (29) This readily translates into the time it takes for cosmic rays to leave the disk with thickness L (this can only take place at the velocity v A ): τ L = yrs L ( N p v A kpc 0.1 cm 3 ) 1/2 B 3 µg. (30) Again, this value is fairly close to the value we have derived using the cosmicray clocks. The remaining difference could in principle be an indication that the magnetic field may be on average smaller along the propagation path. Estimate of the Diffusion coefficient for Galactic cosmic-rays: The height of the halo that is occupied by cosmic-rays can be estimated from the already discussed average density of the traversed interstellar medium. The total column density of (neutral) gas in the direction of the Galactic poles is N H cm 2. Strictly speaking this is a lower limit as the ionised gas is not accounted for. The height of the halo is then H = N H n = cm = 1.6 kpc. (31) Using this value in combination with the estimated time travelled of 10 7 yrs: D = H2 τ 28 cm s. (32) 17

18 The energy dependence of the B/C-secondary-to-primary-ratio implies, that this value obviously changes with energy. A common parameterization (independent of any underlying model) is: ( ) E α D(E) = D 0 (33) 10 GeV with α = (derived from the energy dependence of the B/C-ratio). 4.4 Anisotropy of cosmic-rays The anisotropy of cosmic-rays is an important ingredient on the way to trace the sources of cosmic-rays. The measurements indicate the presence of a relative anisotropy δ < 10 3, with the following definition of δ: δ := I max I min I max + I min. (34) Generally, the measurement of such a small anisotropy is experimentally very challenging. Especially indirect (air-shower) measurements show variations of the detection rate orders of magnitude larger than the level of anisotropy searched for. Furthermore, a large scale (e.g. dipole) anisotropy requires a reasonably homogeneous exposure of the entire sky, which again is not trivial for an individual detector located on Earth. Generally, the observed level of anisotropy should be (cautiously) considered as an upper limit. Nevertheless, we can at least qualitatively discuss the expectations. The degree of anisotropy depends on the location of the observer with respect to the sources. We would start with the assumption that the accelerators of cosmic rays are probably following the general mass distribution of the Galaxy. Most of the accelerators would therefore be located within the disk and fairly close to the inner Galaxy. The solar system is located at the exterior part of the Galaxy and therefore, we would expect an anisotropy with I max located grossly in the direction of the Galactic center and I min probably in the opposite hemisphere or even directed towards the Galactic poles. The fact, that the observed relative anisotropy is well below unity, implies already that cosmic-rays do not propagate rectilinear but obviously diffusively (no surprise after the discussion in the previous sections). For the diffusive transport, we can readily estimate the relative anisotropy. Based upon the 18

19 estimate of the diffusion coefficient in the Galaxy (D cm 2 s 1 ), the net streaming velocity of the cosmic-ray gas is therefore V D R = 107 cm s = 100 km s c. (35) The cosmic-rays are streaming in the direction of the gradient of the cosmicray density. In the rest frame of the streaming motion, the cosmic-ray distribution is isotropic. Given that the nuclei travel roughly with the speed of light, the relative anisotropy for an observer δ V/c A more detailed treatment would have to take into account that the relative motion leads to a shift in the observed energy such that δ = v/c (p + 2) with p the power-law index of the differential energy spectrum. For the cosmic-rays, p = 2.7 which leads to δ = , which in turn is fairly close to the observed anisotropy. The observed degree of anisotropy is in agreement with a diffusive transport of cosmic-rays in our Galaxy. However, there are a number of buts which (in the light of the good agreement of the simple estimate with the measurement) should not be forgotten: The sun is not at the edge of the cosmic-ray source distribution. We can certainly expect, that cosmic ray sources are also active in the direction of the Galactic anti-center. Even more important is the fact that the sun is located close to the center of a local active star forming region called Gould belt. Here, we have e.g. a high number of massive stars as well as supernova remnants which in principle inject cosmic rays in the local medium. The sun is located in the so-called local bubble, ie. the interstellar neigbhorhood of the sun 6 is characterized by a hot, ionised, and tenuous gas, possibly heated up by multiple supernova explosions. The presence of spiral arms. In the spiral arm structure the magnetic field lines are more ordered (as we can derive from the polarization direction of radio-synchrotron emission seen predominantly in spiral galaxies). This ordered magnetic field component leads to a drift of the cosmic-rays parallel to the field lines (the particles gyrate along the 6 As the name implies, the interstellar medium in the solar environment is not typical for the interstellar medium. 19

20 field line and move parallel to the field). This intrinsically should lead to an anisotropy larger than the rough estimate given above. The motion of the sun around the Galactic center as well as the motion of Earth around the sun: The net streaming velocity of the cosmic rays is smaller than the rotational velocity of the sun (> 200 km/s). This leads to an additional dipole-moment that is annualy modulated by the motion of the Earth around the sun ( 30 km/s) that is either parallel or anti-parallel to the orbital motion of the solar system. Both effects need to be removed from the data to infer the streaming velocity. This is commonly done by considering the anisotropy perpendicular to the direction of motion. 4.5 Anti-matter from cosmic-rays in the Galaxy The search for anti-matter is mainly of interest for two reasons: (i)accelerated anti-matter could be the result of acceleration of anti-nuclei in a anti-matter galaxy and (ii) anti-protons (as well as positrons) are produced in cosmicray interactions with gas and provide insights into the cosmic-ray transport. Deviations from the p-flux expected from secondary production may lead us to conclude that additional, primary anti-proton-production is necessary as e.g. suggested in models of evaporating black holes as well as from selfannihilating dark matter (e.g. WIMP-dark matter). 4.6 On the total power required to sustain the cosmicray population Pulling all pieces together, we are ready to estimate the total power required to sustain the cosmic-ray population. Starting with a cylindrical volume of V = πr 2 H with R = 10 kpc, H 1.5 kpc: V cm 3. The energy density of cosmic rays is simply taken to be the locally measured energy density (with a correction of the effect of solar modulation). The canonical value is usually taken to be u cr 1 ev/cm 3. The total energy stored in the form of cosmic rays in the Galaxy is therefore E = u cr V = ev ergs. The power required to balance the escape losses is P = E/τ esc = ergs/s= 10 7 L (in words 10 million solar luminosities) in the form of cosmic-rays. 20

21 Generally, the most attractive class of objects suggested to explain the cosmicray population are supernova remnants. Each supernova remnant releases roughly ergs in the form of kinetic energy in the interstellar medium. Only a fraction η < 1 is released in the form of cosmic-rays. Together with the super-nova remnant rate ṅ snr = 1/100 yrs, we can estimate the efficiency η required to match L cr : L cr = η E SNR ṅ snr (36) ( ) ( ) 1 ( ) L cr ESNR 1/ṅ η = 9 % (37) ergs/s ergs 100 yrs 21

22 5 Appendix Table 1: Partial cross section for spallation reactions. The table is taken from Longair, Chapter 5. 22

23 Table 2: Measured abundances of cosmic rays, solar system, and local Galactic. Table is taken from Longair, chapter 9. 23