Cosmic (Particle) Accelerators II - Sources & Mechanisms - Frank - - PowerPoint PPT Presentation

cosmic particle accelerators ii
SMART_READER_LITE
LIVE PREVIEW

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


slide-1
SLIDE 1

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

slide-2
SLIDE 2

Outline

2

radio (VLA)

  • ptical (HST)
  • Particle Acceleration Mechanisms
  • Gap-type particle acceleration (pulsars, black holes)
  • concept & relevance
  • Fermi-type particle acceleration
  • stochastic 2nd order Fermi
  • shock acceleration - 1st order Fermi (SNR)
  • shear acceleration (AGN)
  • Conclusions
slide-3
SLIDE 3

Possible Acceleration Processes & Sites (not exhaustive)

3

radio (VLA)

  • ptical (HST)

“vacuum” gap shock

(1st order)

reconnection stochastic

(2nd order)

shear “Fermi-type”

wave-particle interactions

“one-shot”

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…

slide-4
SLIDE 4

The Occurrence of Gaps in Pulsar Magnetospheres I

4

  • ptical (HST)
  • in vacuum: e EII >> Fgrav at surface
  • vacuum conditions cannot exist
  • if enough charges, force-free

conditions possible:

  • Goldreich-Julian charge density:
  • co-rotating dipole magnetic field

defines null charge surface

  • no particle acceleration (E|| =0)

⃗ E = − ( ⃗ v × ⃗ B) / c = − ([ ⃗ Ω × ⃗ r] × ⃗ B) /c

Goldreich & Julian 1969

ρGJ = ⃗ ∇ ⋅ ⃗ E 4π ≃ − ⃗ Ω ⋅ ⃗ B 2πc

⃗ B ∝ (2 cos θ ⃗ e r + sin θ ⃗ eθ) / r3 ⇒ ρGJ(r) ∝ (sin2 θ − 2 cos2 θ) / r3

RLC=c/Ω

slide-5
SLIDE 5

The Occurrence of Gaps in Pulsar Magnetospheres II

5

  • ptical (HST)
  • ideal MHD in most of magneto-

sphere:

  • deficient charge supply:

⇒ particle acceleration

  • Solve Gauss’ law:

⃗ ∇ ⋅ ⃗ E = 4π (ρ − ρGJ)

⃗ E ⋅ ⃗ B = 0

⃗ E ⋅ ⃗ B ≠ 0

Possible sites of particle acceleration

(Credits: A. Harding)

(e.g., Ruderman & Sutherland 1975; Cheng et al. 1985; Muslimov & Harding 2003)

slide-6
SLIDE 6

The Occurrence of Gaps in BH Magnetospheres

6

  • ptical (HST)
  • Null surface in Kerr Geometry (r ~ rg≣GM/c2)

for force-free magnetosphere, vanishing of poloidal electric field Ep ∝ (ΩF-ω) ∇Ψ = 0, ω =Lense-Thirring ⇒ 𝜍GJ changes sign, “gap” may easily develop

  • Stagnation surface (r ~ few rg)

Inward flow of plasma below due to gravitation field,

  • utward motion above

⇒ charges need to be continuously replenished

Levinson & Segev 2017 (e.g., Blandford & Znajek 1977; Beskin et al. 1992, Hirotani & Okamoto 1998)

slide-7
SLIDE 7

The Conceptual Relevance of BH Gaps

7

  • ptical (HST)
  • BH-driven jets (Blandford-Znajek)
  • Self-consistency: Plasma source needed to

ensure force-free MHD

  • Non-thermal Particle Acceleration
  • Implication: efficient (direct) acceleration of

electrons & positrons

  • Radiation & Pair Cascade…..
  • Features: expect ɣ-ray production,
  • ɣɣ-absorption triggers pair cascade
  • generating charge multiplicity
  • ensuring electric field screening (closure)

Koide+

slide-8
SLIDE 8

Gamma-Ray Emission from AGN Magnetospheres

8

  • ptical (HST)
  • Direct electric field acceleration:

Rate of energy gain for electron:

dɣ/dt ∝ e Δϕgap·(c/h)

  • Curvature & Inverse Compton:

HE ɣ-rays via curvature: 𝓦~(0.2c) (ɣ3/Rc) VHE ɣ-rays via IC: h𝓦 ≲ ɣ me c2

  • Accretion environment (RIAF):

Radiatively inefficient needed to facilitate escape of VHE photons

  • Maximum Gap luminosity:

Lgap ∝ nGJ (Volume) (dɣ/dt)

slide-9
SLIDE 9

9

  • ptical (HST)

Characterizing the Magnetospheric Potential

Possible boundary conditions in the pulsar case :

  • “non-free escape” (Ruderman): EII(h=0)≠0, E||(h=H)=0, ρe << ρGJ :
  • “free escape” (Arons): EII(h=0)=0, E||(h=H)=0, ρe ~ ρGJ (ρe ≠ ρGJ ≡ΩB cosθb ) :

dE|| dh = 4π (ρe − ρGJ) ”Gauss0 law”

slide-10
SLIDE 10

9

  • ptical (HST)

Characterizing the Magnetospheric Potential

Possible boundary conditions in the pulsar case :

  • “non-free escape” (Ruderman): EII(h=0)≠0, E||(h=H)=0, ρe << ρGJ :
  • “free escape” (Arons): EII(h=0)=0, E||(h=H)=0, ρe ~ ρGJ (ρe ≠ ρGJ ≡ΩB cosθb ) :

dE|| dh = 4π (ρe − ρGJ) ”Gauss0 law”

EII H h

slide-11
SLIDE 11

9

  • ptical (HST)

Characterizing the Magnetospheric Potential

Possible boundary conditions in the pulsar case :

  • “non-free escape” (Ruderman): EII(h=0)≠0, E||(h=H)=0, ρe << ρGJ :
  • “free escape” (Arons): EII(h=0)=0, E||(h=H)=0, ρe ~ ρGJ (ρe ≠ ρGJ ≡ΩB cosθb ) :

dE|| dh = 4π (ρe − ρGJ) ”Gauss0 law”

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

slide-12
SLIDE 12

9

  • ptical (HST)

Characterizing the Magnetospheric Potential

Possible boundary conditions in the pulsar case :

  • “non-free escape” (Ruderman): EII(h=0)≠0, E||(h=H)=0, ρe << ρGJ :
  • “free escape” (Arons): EII(h=0)=0, E||(h=H)=0, ρe ~ ρGJ (ρe ≠ ρGJ ≡ΩB cosθb ) :

dE|| dh = 4π (ρe − ρGJ) ”Gauss0 law”

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

slide-13
SLIDE 13

9

  • ptical (HST)

Characterizing the Magnetospheric Potential

Possible boundary conditions in the pulsar case :

  • “non-free escape” (Ruderman): EII(h=0)≠0, E||(h=H)=0, ρe << ρGJ :
  • “free escape” (Arons): EII(h=0)=0, E||(h=H)=0, ρe ~ ρGJ (ρe ≠ ρGJ ≡ΩB cosθb ) :

dE|| dh = 4π (ρe − ρGJ) ”Gauss0 law”

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

slide-14
SLIDE 14

10

radio (VLA)

  • ptical (HST)

Magnetospheric Potential & Jet Power in AGN - Differences

  • Gap potential:
  • Δϕgap ~ aspin rg B (H/rg)2
  • Constraining losses:
  • Curvature, IC…
  • Jet power:
  • LVHE ~ Ljet x (H/rg)2 …

highly under-dense: 𝞻e<<𝞻GJ weakly under-dense: 𝞻e~𝞻GJ

  • Gap potential:
  • Δϕgap ~ aspin rg B (H/rg)3
  • Constraining losses:
  • IC, curvature…
  • Jet power:
  • LVHE ~ Ljet x (H/rg)4 …

Solving Gauss’ laws depending on different boundaries

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”

slide-15
SLIDE 15

10

radio (VLA)

  • ptical (HST)

Magnetospheric Potential & Jet Power in AGN - Differences

  • Gap potential:
  • Δϕgap ~ aspin rg B (H/rg)2
  • Constraining losses:
  • Curvature, IC…
  • Jet power:
  • LVHE ~ Ljet x (H/rg)2 …

highly under-dense: 𝞻e<<𝞻GJ weakly under-dense: 𝞻e~𝞻GJ

  • Gap potential:
  • Δϕgap ~ aspin rg B (H/rg)3
  • Constraining losses:
  • IC, curvature…
  • Jet power:
  • LVHE ~ Ljet x (H/rg)4 …

Jet power constraints can become relevant

Solving Gauss’ laws depending on different boundaries

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”

slide-16
SLIDE 16

Timescales (example)

11

  • ptical (HST)

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)

Loss time scales:

𝜐cur ∝1/ɣ3, 𝜐IC ∝1/ɣ𝛃 (𝛃=1 Thompson, 𝛃<0 KN)

Energy gain

Katsoulakos & FR 2018

Parameters: M9=5, ṁ=10-4 (ADAF), h/rg=0.5

slide-17
SLIDE 17

Timescales (example)

11

  • ptical (HST)

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)

Loss time scales:

𝜐cur ∝1/ɣ3, 𝜐IC ∝1/ɣ𝛃 (𝛃=1 Thompson, 𝛃<0 KN)

Energy gain

Katsoulakos & FR 2018

can reach Lorentz factors ɣ~1010

Parameters: M9=5, ṁ=10-4 (ADAF), h/rg=0.5

slide-18
SLIDE 18

Example: Phenomenological Relevance of Gaps in AGN

12

  • ptical (HST)
  • Gamma-Ray Emission from Radio Galaxies:

misaligned jets: moderate Doppler boosting of jet emission only ⇒ gap IC & curvature emission may show up at hard HE-VHE gamma-rays

  • Possibly related to observable AGN features in:
  • M87 (d ~17 Mpc): day-scale VHE variability, radio-VHE outburst correlation…
  • Cen A (d ~ 4 Mpc): spectral hardening of core emission above ~5 GeV…
  • IC 310 (d ~ 80 Mpc): rapid (5 min) VHE variability, huge power (Lɣ ~ 1044 erg/sec)
  • 1

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

slide-19
SLIDE 19

Example: Maximum Gap Power Constraints

13

  • ptical (HST)

−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

  • Magnetic field limited by

accretion: B ∝ ṁ 1/2

  • Gap size limited by

variability: H ~ c Δt

slide-20
SLIDE 20

Example: Maximum Gap Power Constraints

13

  • ptical (HST)

−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

  • Magnetic field limited by

accretion: B ∝ ṁ 1/2

  • Gap size limited by

variability: H ~ c Δt

slide-21
SLIDE 21

Preference for Outer Gap Acceleration in Pulsars ?

14

  • ptical (HST)
  • Polar Cap Acceleration:
  • absorption via magnetic pair creation,

super-exponential cut-off in gamma-ray emission

  • Outer Gap Acceleration:
  • curvature radiation, exponential cut-off

in gamma-ray emission

  • Fermi-LAT HE observations:
  • super-exponential cutoff excluded
  • brightest pulsars (Crab,

Vela) : even show sub-exponential cut-off

➡ superposition (states & sites) ?

  • cut-offs in narrow band Ecut ~1-5 GeV
  • compatible with curvature radiation
  • origin of TeV (IC)?

Fermi LAT Collab. 2009

Crab Pulsar

Phase-averaged spectrum, Ecut~ 5 GeV

F(E)*exp(-E/Ecut)b , b ≃1

slide-22
SLIDE 22

Gap-type Particle Acceleration - Summary

15

  • ptical (HST)
  • gaps (“unscreened parallel electric fields”) are to be expected in the mag-

netospheres of pulsars, and may occur around supermassive black holes

  • most efficient (“direct - one-shot”) particle acceleration mechanism
  • energy gain dE/dt ≃ e ɸ (c/H)
  • acceleration timescales can be as short as tacc ~ ɣ m c / (eB)
  • unavoidable max. cutoff due to curvature radiation
  • pulsars : ɣmax ~ 107-8 (e+e-)
  • AGN : ɣmax ~ 1010 (e, p)
  • Development of pair cascade may limit size of gap & lead to closure
slide-23
SLIDE 23

Possible Acceleration Processes & Sites (not exhaustive)

16

radio (VLA)

  • ptical (HST)

“vacuum” gap shock

(1st order)

reconnection stochastic

(2nd order)

shear “Fermi-type”

wave-particle interactions

“one-shot”

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…

slide-24
SLIDE 24

17

  • ptical (HST)

Fermi-type Particle Acceleration

Kinematic effect resulting from scattering off magnetic inhomogeneities

Fermi, Phys. Rev. 75, 578 [1949]

⇒ energy gain as results of multiple scatterings (stochastic process) _Ingredients: in frame of scattering centre

  • momentum magnitude conserved
  • particle direction randomised

Magnetic cloud

Vc E < E

f i

E i E > E

f i

E i

slide-25
SLIDE 25

18

  • ptical (HST)

Fermi Acceleration - energy change in elastic scattering event I

slide-26
SLIDE 26

19

  • ptical (HST)

Fermi Acceleration - energy change in elastic scattering event II

_Characteristic energy change per scattering (non-relativistic Vc):

energy gain for head-on (p Vc < 0), loss for following collision (p Vc > 0)

  • stochastic: average energy gain 2nd order: < ΔE > ~ (Vc / c) 2 E

ΔE = Ef − Ei = 2 (Ei V2

c /c2 −

⃗ p i ⋅ ⃗ V c)

slide-27
SLIDE 27

19

  • ptical (HST)

Fermi Acceleration - energy change in elastic scattering event II

_Characteristic energy change per scattering (non-relativistic Vc):

energy gain for head-on (p Vc < 0), loss for following collision (p Vc > 0)

  • stochastic: average energy gain 2nd order: < ΔE > ~ (Vc / c) 2 E

ΔE = Ef − Ei = 2 (Ei V2

c /c2 −

⃗ p i ⋅ ⃗ V c)

can we do better?

slide-28
SLIDE 28

20

radio (VLA)

  • ptical (HST)

Fermi Acceleration @ shocks

, l , l , l , l , l , l

slide-29
SLIDE 29

21

radio (VLA)

  • ptical (HST)

Fermi Acceleration @ shocks

_For particles crossing the shock, scattering is always head-on:

  • shock: spatial diffusion, gain on crossing is 1st order: < ΔE > ~ (Δv / c) E
slide-30
SLIDE 30

22

radio (VLA)

  • ptical (HST)

Fermi Acceleration Timescales

_Acceleration timescale ~ particle energy / (rate of energy change):

  • stochastic: 𝜐 = λ / c “mean scattering time” (λ = mean free path):
  • shock: spatial diffusion process 𝜐 = tc ~ 𝛌 / (Vs c) “crossing time”

(residence time tc = diffusion length / c , with diffusion length l ~ √𝛌 t, and t ~ l / Vs )

tacc = E (dE/dt) ≃ E ΔE × τ

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

slide-31
SLIDE 31

22

radio (VLA)

  • ptical (HST)

Fermi Acceleration Timescales

_Acceleration timescale ~ particle energy / (rate of energy change):

  • stochastic: 𝜐 = λ / c “mean scattering time” (λ = mean free path):
  • shock: spatial diffusion process 𝜐 = tc ~ 𝛌 / (Vs c) “crossing time”

(residence time tc = diffusion length / c , with diffusion length l ~ √𝛌 t, and t ~ l / Vs )

tacc = E (dE/dt) ≃ E ΔE × τ

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

  • rder in shock

speed !

slide-32
SLIDE 32

_Gradual shear flow with frozen-in scattering centres:

  • like 2nd Fermi, stochastic process with average gain:

using characteristic effective velocity:

, where λ = particle mean free path

  • leads to
  • seed from acceleration @ shock or stochastic….
  • easier for protons….

23

Fermi Acceleration @ shear flows

⌅ u = uz(x) ⌅ ez

non-relativistic

< ΔE > E ∝ ( V c )

2

= 1 c2 ( ∂ux ∂x )

2

λ2

V = Δu = ( ∂uz ∂x ) λ tacc = E (dE/dt) ≃ E ΔE × τ ∼ ( c [∂uz/∂x] λ )

2

λ c ∝ 1 λ

slide-33
SLIDE 33

24

  • ptical (HST)

_Emission from large-scale jets

  • extended X-ray electron synchrotron emission
  • electron Lorentz factors ɣ ~108
  • short cooling timescale tcool ∝ 1/ɣ ; cooling length c tcool << kpc
  • distributed acceleration mechanism required (Sun,Yang, FR+ 2018 for M87)

ν ∝ ɣ2 B

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)

Example: Stochastic & shear acceleration in large-scale AGN jets I

slide-34
SLIDE 34

25

Ansatz: Fokker-Planck equation for f(t,p) including stochastic, shear and synchrotron for cylindrical jet.

  • from 2nd Fermi (tacc∝ λ) to shear

(tacc ∝ 1 / λ)…

  • electron acceleration up to ɣ~109

possible

  • formation of multi-component

particle distribution

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

Radiative-loss-limited acceleration in mildly relativistic flows

(Liu, FR, Aharonian 2017)

4/3 2/3

Example: Stochastic & shear acceleration in large-scale AGN jets II

slide-35
SLIDE 35

_”1st order” Fermi - standard shock (non-relativistic): with shock crossing time tc ~ κ /(us c), where κ ~ λ c _”2nd order” Fermi (stochastic): with scattering time τ ~ λ/c _Shear - gradual (non-relativistic):

26

(e.g., Drury 1983; Kirk 1994; Duffy & Blundell 2005; FR+ 2007)

Fermi Acceleration Timescales - Summary

tacc ∼ ( c VA )

2 λ

c ∝ λ V2

A

tacc ∼ ( c Vs) κ Vs c ∼ κ V2

s

∝ λ V2

s

tacc ∼ ( c [∂uz/∂x] λ)

2

λ c ∝ 1 λ

slide-36
SLIDE 36

27

radio (VLA)

  • ptical (HST)

Example: Shocks in SNRs - historical shell SNR

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:

  • radio (GHz): ɣe ~ 103-4
  • X-rays (keV): ɣe ~ 107-8

(but: degeneracy in B & ɣ)

slide-37
SLIDE 37

28

radio (VLA)

  • ptical (HST)

Example: Efficient Cosmic Ray (PeV) Acceleration @ SNR shocks ?

  • Acceleration timescale:
  • with spatial diffusion coefficient:
  • smallest possible mean free path: λ ≈ rgyro = E / (e B)

⇒ Limit on maximum CR energy:

  • typical for young SNR:

ISM mag field: few μG

Vs = c / 50

R ~ 1019 cm

tacc ≃ 8 κ V2

s

= 8 3 λ c V2

s

κ = λ c / 3 }

λ ≤ 3 8 Vs c R

tacc ≤ tage ≃ R / Vs implies:

SNR radius

Emax ≤ 3 8 Vs c R e B Emax ≲ 1014 (B / 5μG) eV

(Lagage & Cesarsky 1983) }

slide-38
SLIDE 38

28

radio (VLA)

  • ptical (HST)

Example: Efficient Cosmic Ray (PeV) Acceleration @ SNR shocks ?

  • Acceleration timescale:
  • with spatial diffusion coefficient:
  • smallest possible mean free path: λ ≈ rgyro = E / (e B)

⇒ Limit on maximum CR energy:

  • typical for young SNR:

ISM mag field: few μG

Vs = c / 50

R ~ 1019 cm

tacc ≃ 8 κ V2

s

= 8 3 λ c V2

s

κ = λ c / 3 }

λ ≤ 3 8 Vs c R

tacc ≤ tage ≃ R / Vs implies:

SNR radius

Emax ≤ 3 8 Vs c R e B Emax ≲ 1014 (B / 5μG) eV

(Lagage & Cesarsky 1983) } need amplified magnetic field (e.g., Lucek & Bell 2000)

slide-39
SLIDE 39

28

radio (VLA)

  • ptical (HST)

Example: Efficient Cosmic Ray (PeV) Acceleration @ SNR shocks ?

  • Acceleration timescale:
  • with spatial diffusion coefficient:
  • smallest possible mean free path: λ ≈ rgyro = E / (e B)

⇒ Limit on maximum CR energy:

  • typical for young SNR:

ISM mag field: few μG

Vs = c / 50

R ~ 1019 cm

tacc ≃ 8 κ V2

s

= 8 3 λ c V2

s

κ = λ c / 3 }

λ ≤ 3 8 Vs c R

tacc ≤ tage ≃ R / Vs implies:

SNR radius

Emax ≤ 3 8 Vs c R e B Emax ≲ 1014 (B / 5μG) eV

(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)

slide-40
SLIDE 40

29

radio (VLA)

  • ptical (HST)

Some Issues Concerning Fermi-type Particle Acceleration

  • Stochastic particle acceleration:
  • generates no unique power-law particle distribution, e.g., index depends on ratio of

tacc/tescape ; if synchrotron-loss limited, relativistic Maxwellian distributions may occur…

  • slow process unless the scattering center speed is high (Alfven speed; AGN jets)…
  • Shock acceleration:
  • highly relativistic shocks (PWN, GRBs etc) are not expected to be efficient accelerators

(e.g., isotropization upstream not guaranteed; relativistic shocks are generically quasi- perpendicular as B⟂ = 3 Γs B⟂’…)

  • no longer a unique power law…
  • Shear acceleration:
  • only efficient in relativistic shear flows
  • particle transport across flow still to be understood
  • destruction of flow (KH/shear instabilities)?

(e.g., Sironi+ 2013, Lemoine & Pelletier 2017; Bell+ 2018, Webb+ 2018, FR 2019)

slide-41
SLIDE 41

30

radio (VLA)

  • ptical (HST)

The END

Thank you!