cosmic rays - university of iowaastro.physics.uiowa.edu/~kaaret/heastro10s/l15_cosmicray.pdftotal...

Cosmic Rays

• Discovery of cosmic rays

• Local measurements

• Gamma-ray sky (and radio sky)

• Origin of cosmic rays

Discovery of Cosmic Rays

Problem that electroscopes would always lose their charge.

In 1912 Hess flew electroscopes in balloons (up to 17,500 feet) and showed that the rate of loss increased with altitude, thus showing that the particles causing the loss of charge were produced external to the Earth. He called them cosmic radiation.

γ = 2.7

γ ≈ 3.1

γ ≈ 2.7

Cosmic raysIsotropic

CR: 2% electrons, 98% hadrons. Hadrons:89% H, 10% He, 1% heavier elements.

Energy density CR ~ 1 eV/cm3

Starlight ~ 0.3 eV/cm3

B-field ~ 0.2 eV/cm3

CMB ~ 0.3 eV/cm3

“Knee” at 3× 1015 eV

Cosmic ray spectrum

Confinement

Milky Way B ~ 3 µG.

Lamour radius rg = p/qB, for protons, rg = 1012 cm (E/GeV).

Scale height of Galactic disk is ~5× 1020 cm, thus, protons with energies up to about 1017 eV can be confined.

At low energies, heliosphere affects trajectories.

Greisen-Zatsepin-Kuzmin Cutoff

• Cosmic rays will interact with cosmic microwave background: p+γ → p+π

• Only occurs when proton has enough energy to produce pion, E ~ 5× 1019 eV

• Detection of particles above this energy requires “local” sources (or new physics).

Cosmic ray abundances

Rare elements and radioactive isotopes are over abundant due to spallation.

Spallation

p 816O 3

7 L i X

Cosmic rays collide with nuclei in ISM, changing the composition. For example:

Observed CR composition depends on:● Initial CR composition● CR path length/life time● CR energy

Measurement of CR composition, including isotopes, allows one to constrain these quantities.Find path length traversed by CR by comparing abundant elements to those produced by spallation (i.e. B vs C, Cr vs Fe). Path length ~ 50 kg m-2 (with some energy dependence).

Radioactive Clocks

● If there is a cosmic ray nucleus with a radioactive isotope and we can determine the ratios of the different isotopes produced via spallation, then we can figure out the age of cosmic rays.● 10Be is about 10% of all Be produced by spallation, has a half-life of 1.6× 106 years.● The measured fraction of 10Be in Be cosmic rays is about 3%.● Conclude that cosmic rays are confined in the Galaxy for about 107 years before escaping.

Leaky Box“Leaky box” model – CRs diffuse inside a containment volume (volume of Milky Way disk) and are reflected at boundaries with some probability of escape. Typical path length ~ 50 kg m-2 and typical lifetime ~ 107 years.

Confinement is assumed to be done by Galactic magnetic field.

Note that highest energy CR are not confined by magnetic fields of Milky Way.

Implies roughly uniform distribution of cosmic rays through out the Galaxy.

Photon production by cosmic rays

• Pion production: p + N → π 0 + X π 0 has total charge = 0, baryon number = 0 so π 0 can decay via π 0 → γ γ

Neutron pion mass = 135 MeV, Decay produces two photons of ~70 MeV

• Electron bremsstrahlung – cosmic ray electrons on ISM• Inverse Compton – cosmic ray electrons on star light

( )dduu −=2

10π

Gamma-Ray Sky

Gamma-ray spectrum

Modeling of spatial distribution and spectrum requires 3-d models of cosmic ray, matter, and star light distributions in Milky Way.

Need to multiply by cross sections and convolve along lines of sight.

Suggestion that CR spectrum is harder towards Galactic center.

Photon production by cosmic rays

• Electron cosmic rays produce radio emission via synchrotron radiation in the Galactic magnetic field• Radiated spectrum peaks at a frequency

≈

G

BE

µν

3GeV 10GHz) 2.3(

2

As for gamma-rays, need to convolve the electron CR spectrum with the Galactic magnetic field distribution along each line of sight.

Milky Way at 408 MHz

Total Power in Cosmic Rays

Volume of Galactic disk V ~ π R2d.

For R ~ 15 kpc, d ~ 200 pc, find V ~ 4× 1066 cm3.

Power in cosmic rays L ~ Vρ/τ .ρ = energy density ~ 1 eV/cm3τ = lifetime of CR ~ 107 years.

Find L ~ 1041 erg/s in high energy particles.

Power Source: Supernovae?

• Supernovae –– E = Mechanical energy ~ 1051 erg

– R = rate 1/100 years– η = efficiency for conversion of mechanical energy

into relativistic particles ~ 10% (?)

• LSN ~ ηER ~ 2× 1041 erg/s

• Need mechanism for acceleration, need to know if acceleration is really 10% efficient.

Power Source: Massive Star Winds?

• O and B star winds –– Winds have speeds of 2000-4000 km/s

– Expect multiple stars within OB associations

– Cyg OB2 has a wind luminosity ~ 1039 erg/s, would need 1000 such associations in the Milky Way to power cosmic rays

• OB associations are bright in gamma-rays

Cosmic Ray Map

Power Source: X-Ray Binaries?

• Jets from X-ray binaries known to contain relativistic particles –– Only SS 433 is know to accelerate hadrons, and that

jet is not ultrarelativistic– Integrated output of X-ray binaries appears to be too

low to power full CR population, but may contribute a few percent