Cosmic (Particle) Accelerators II
- Sources & Mechanisms -
Max Planck Institut für Kernphysik Heidelberg, Germany
Frank M. Rieger ISAPP School Heidelberg, May 28, 2019
ITA Univ. Heidelberg
Cosmic (Particle) Accelerators II - Sources & Mechanisms - Frank - - PowerPoint PPT Presentation
Cosmic (Particle) Accelerators II - Sources & Mechanisms - Frank M. Rieger ISAPP School Heidelberg, May 28, 2019 Max Planck Institut ITA Univ. Heidelberg fr Kernphysik Heidelberg, Germany Outline radio (VLA) Particle
Max Planck Institut für Kernphysik Heidelberg, Germany
Frank M. Rieger ISAPP School Heidelberg, May 28, 2019
ITA Univ. Heidelberg
2
radio (VLA)
3
radio (VLA)
(1st order)
(2nd order)
wave-particle interactions
charge density? topology? transparency? limited in size? efficiency? spectral shape? efficiency (Γs, σ)? localized? efficiency? spectral shape? efficiency? energetic seeds? AGN & Pulsars… AGN, PWN… AGN, SNRs, PWN.. AGN… AGN…
4
Goldreich & Julian 1969
⃗ B ∝ (2 cos θ ⃗ e r + sin θ ⃗ eθ) / r3 ⇒ ρGJ(r) ∝ (sin2 θ − 2 cos2 θ) / r3
RLC=c/Ω
5
Possible sites of particle acceleration
(Credits: A. Harding)
(e.g., Ruderman & Sutherland 1975; Cheng et al. 1985; Muslimov & Harding 2003)
6
Levinson & Segev 2017 (e.g., Blandford & Znajek 1977; Beskin et al. 1992, Hirotani & Okamoto 1998)
7
Koide+
8
9
9
EII H h
9
dE|| dh ' 4π ρGJ ) E||(h) = 4π ρGJh + const E||(h = H) = 0 ) const = 4πρGJH Thus : E||(h) = E0 (H h) H , where E0 = 4πρGJH
EII H h
9
dE|| dh ' 4π ρGJ ) E||(h) = 4π ρGJh + const E||(h = H) = 0 ) const = 4πρGJH Thus : E||(h) = E0 (H h) H , where E0 = 4πρGJH
EII H h
9
dE|| dh ' 4π ρGJ ) E||(h) = 4π ρGJh + const E||(h = H) = 0 ) const = 4πρGJH Thus : E||(h) = E0 (H h) H , where E0 = 4πρGJH dE|| dh ' 4π d(ρ ρGJ) dh |h=H/2 (h H/2) ) E||(h) = EA h(H h) H2 with EA =2π d(ρ ρGJ) dh H2
EII H h
10
radio (VLA)
e.g., Blandford & Znajek 1982, Levinson 2000 Levinson & FR 2011 e.g., Hirotani & Pu 2016 Katsoulakos & FR 2018
dE|| dh = 4π (ρe − ρGJ) ”Gauss0 law”
10
radio (VLA)
Jet power constraints can become relevant
e.g., Blandford & Znajek 1982, Levinson 2000 Levinson & FR 2011 e.g., Hirotani & Pu 2016 Katsoulakos & FR 2018
dE|| dh = 4π (ρe − ρGJ) ”Gauss0 law”
11
5 6 7 8 9 10 11 12 13 −6 −4 −2 2 4 6 8 10
Lorentz Factor, log10(γe) Characteristic time scale, log10(τ)
τcur τic τacc (η=1.0, ν=1.0) τacc (η=1.0, ν=2.0) τacc (η=1/6, ν=3.0)
Katsoulakos & FR 2018
Parameters: M9=5, ṁ=10-4 (ADAF), h/rg=0.5
11
5 6 7 8 9 10 11 12 13 −6 −4 −2 2 4 6 8 10
Lorentz Factor, log10(γe) Characteristic time scale, log10(τ)
τcur τic τacc (η=1.0, ν=1.0) τacc (η=1.0, ν=2.0) τacc (η=1/6, ν=3.0)
Katsoulakos & FR 2018
can reach Lorentz factors ɣ~1010
Parameters: M9=5, ṁ=10-4 (ADAF), h/rg=0.5
12
1 2 3 4 14 13 12 11 10 logEGeV logE2dNdEerg cm2s1
(cf. FR & Levinson 2018 for review and references)
Cen A M87 IC 310
VHE flare in Nov 2012
13
−7 −6 −5 −4 −3 −2 −1 −10 −9 −8 −7 −6 −5 −4 −3 −2 −1 1
Accretion rate, log10( ˙ m) log10[10−48LgapM −1
9 ( h rg )−β]
η = 1 η = 1/6
Critical accretion rate
M87 (β=1) IC310 (β=1)
IC310 (β=4)
Katsoulakos & FR 2018
Gap Power Limit
accretion: B ∝ ṁ 1/2
variability: H ~ c Δt
13
−7 −6 −5 −4 −3 −2 −1 −10 −9 −8 −7 −6 −5 −4 −3 −2 −1 1
Accretion rate, log10( ˙ m) log10[10−48LgapM −1
9 ( h rg )−β]
η = 1 η = 1/6
Critical accretion rate
M87 (β=1) IC310 (β=1)
IC310 (β=4)
Katsoulakos & FR 2018
Gap ɣ-ray origin: M87 possible IC310 not expected
Gap Power Limit
accretion: B ∝ ṁ 1/2
variability: H ~ c Δt
14
super-exponential cut-off in gamma-ray emission
in gamma-ray emission
Vela) : even show sub-exponential cut-off
Fermi LAT Collab. 2009
Crab Pulsar
Phase-averaged spectrum, Ecut~ 5 GeV
F(E)*exp(-E/Ecut)b , b ≃1
15
16
radio (VLA)
(1st order)
(2nd order)
wave-particle interactions
charge density? topology? transparency? limited in size? efficiency? spectral shape? efficiency (Γs, σ)? localized? efficiency? spectral shape? efficiency? energetic seeds? AGN & Pulsars… AGN, PWN… AGN, SNRs, PWN.. AGN… AGN…
17
Fermi, Phys. Rev. 75, 578 [1949]
Magnetic cloud
Vc E < E
f i
E i E > E
f i
E i
18
19
c /c2 −
19
c /c2 −
can we do better?
20
radio (VLA)
, l , l , l , l , l , l
21
radio (VLA)
22
radio (VLA)
(residence time tc = diffusion length / c , with diffusion length l ~ √𝛌 t, and t ~ l / Vs )
tacc = E (dE/dt) ≃ E ΔE × tc ∼ ( c Vs) κ Vs c ∼ κ V2
s
∝ λ V2
s
tacc = E (dE/dt) ≃ E ΔE × τ ∼ ( c VA)
2 λ
c ∝ λ V2
A
22
radio (VLA)
(residence time tc = diffusion length / c , with diffusion length l ~ √𝛌 t, and t ~ l / Vs )
tacc = E (dE/dt) ≃ E ΔE × tc ∼ ( c Vs) κ Vs c ∼ κ V2
s
∝ λ V2
s
tacc = E (dE/dt) ≃ E ΔE × τ ∼ ( c VA)
2 λ
c ∝ λ V2
A
also second
speed !
23
2
2
2
24
Radio (VLA) Optical (HST) X-ray (Chandra)
1 arcsec ~ 0.1 kpc (0.081 kpc)
Marshall+ 2002
M87
Relativistic particles throughout whole jet
HST Chandra
SED can be fitted by broken power-law
(B = 3x10-4 G, ɣb~106, ɣmax~108, Pjet~ 1043 erg/s, Δ𝛃 ~ 2)
Synchrotron origin: spectral index radio ≠ X-ray (should be similar in IC-CMB)
VLA (14 GHz) ~1400 ksec (2000-2016)
25
Parameters: B = 3μG, vj,max ~ 0.4c, rj ~ 30 pc, βA~ 0.007, Δr ~ rj/10, mean free path λ = ξ-1 rg (rg/𝝡max)1-q ∝ ɣ2-q, q=5/3 (Kolmogorov), ξ=0.1 synchrotron flux particle distribution
(Liu, FR, Aharonian 2017)
4/3 2/3
26
(e.g., Drury 1983; Kirk 1994; Duffy & Blundell 2005; FR+ 2007)
tacc ∼ ( c VA )
2 λ
c ∝ λ V2
A
tacc ∼ ( c Vs) κ Vs c ∼ κ V2
s
∝ λ V2
s
2
27
radio (VLA)
Chandra X-ray emission (Credits: NASA+)
Mixture of line radiation (hot plasma) & synchrotron continuum (relativistic electrons). For electron synchrotron in (amplified) magnetic field of ~ 0.1-1 mG:
(but: degeneracy in B & ɣ)
28
radio (VLA)
Vs = c / 50
R ~ 1019 cm
tacc ≃ 8 κ V2
s
= 8 3 λ c V2
s
SNR radius
(Lagage & Cesarsky 1983) }
28
radio (VLA)
Vs = c / 50
R ~ 1019 cm
tacc ≃ 8 κ V2
s
= 8 3 λ c V2
s
SNR radius
(Lagage & Cesarsky 1983) } need amplified magnetic field (e.g., Lucek & Bell 2000)
28
radio (VLA)
Vs = c / 50
R ~ 1019 cm
tacc ≃ 8 κ V2
s
= 8 3 λ c V2
s
SNR radius
(Lagage & Cesarsky 1983) } need amplified magnetic field (e.g., Lucek & Bell 2000)
but limitations due to self-regulated CR escape implying Emax < 1 PeV for Tycho, Cas A, Kepler (e.g., Bell+ 2013)
29
radio (VLA)
tacc/tescape ; if synchrotron-loss limited, relativistic Maxwellian distributions may occur…
(e.g., isotropization upstream not guaranteed; relativistic shocks are generically quasi- perpendicular as B⟂ = 3 Γs B⟂’…)
(e.g., Sironi+ 2013, Lemoine & Pelletier 2017; Bell+ 2018, Webb+ 2018, FR 2019)
30
radio (VLA)