HIGH ENERGY ASTROPHYSICS - Lecture 5 PD Frank Rieger ITA & MPIK Heidelberg Wed 1 Synchrotron Emission I 1 Overview • Radiation of charged particles gyrating in magnetic fields. • Total emitted power for single particle (Larmor formula), spectral distribu- tion. • Characteristic Emission at ν / γ2B, nearly linearly polarised. • For a non-thermal, power law electron distribution, n(γ) / γ−p, emitted −(p−1)=2 spectrum a power law jν / ν too. • Believed to be the fundamental radiation process at the lower end of the HE spectrum. 2 2 Synchrotron Radiation "Magnetobremsstrahlung" = Radiation emitted by relativistic charged particles (mostly electrons) due to acceleration ("gyro-motion") in a static magnetic field. Note: In the sub-relativistic regime this process is called "Cyclotron emission" Approach: Semi-quantitative analysis of spectral features as detailed calcula- tions are lengthy (see, e.g., Rybicki & Lightman, Chap. 6). Structure: • Motion of electrons in magnetic fields • Look at emission from a single electron • Consider electron distribution and opacity effects to obtain final spectrum 3 3 Motion of a Charged Particle in Static Magnetic Field Remember: Particle gyrates around magnetic field with angular frequency (Lar- mor frequency) qB Ω = L γmc 2 and associated radius of gyration orbit (Larmor radius), γmv?=rL = (q=c) v×B, 2 v? γmcv sin θ γmc rL = = ' ΩL qB qB Numerically: B m 1 Ω ' 1:8 × 107 e Hz L 1 G m γ c 3 1 G m rL ' ' 1:7 × 10 γ cm ΩL B me 10−6 G m ' 1:7 × 109 γ cm B me −6 ) rL is small on cosmic scales (B ∼ 10 G in interstellar medium, B ∼ 1 G in center of AGN [but away from BH]). 4 4 Total Emitted Power for a Single Particle • Remember: Larmor's Formula for relativistically moving particle (lecture 4): 2 4 dE 2 q γ 2 2 2 = 3 (a? + γ ak) dt K 3 c with a?; ak = acceleration components perpendicular and parallel to direction of motion (not to magnetic field!). 2 • For gyro-motion, ak = 0 and a? = v?=rL = ΩLv? (centrifugal), with v? = v sin θ, so 2 e2γ4 v2 e2B2 2 β2e4B2 P = ? = γ2 sin2 θ syn 3 c3 γ2m2c2 3 m2c3 • For isotropic distribution of particles, average pitch angle 1 Z 4π 1 Z 2π Z π < sin2 θ >= sin2 θdΩ = dφ sin2 θ sin θdθ 4π 0 4π 0 0 1 2 = [cos(3θ) − 9 cos(θ)]π=12 = 2 0 3 5 • Average single particle synchrotron power [erg/s] for β ! 1 4 e4B2 4 m 2 hP i = γ2 = cσ e γ2 u syn 9 m2c3 3 T m B with magnetic energy density B2 u := ; B 8π and Thomson cross-section 2 2 8π e −25 2 σT := 2 = 6:65 × 10 cm 3 mec • Characteristic electron cooling timescale: 2 E γmec 8 1 tcool(γ) := = = 7:8 × 10 2 sec jdE=dtj < Psyn > B γ Note: 2 1. Psyn / 1=m ) Synchrotron radiation from charged particles with larger mass (e.g. protons) is usually negligible. 2 2 2. E = γmc ) Psyn / E uB, so more energetic particles radiate more. 6 5 Energy Evolution for Electrons due to Synchrotron Losses Start at t0 with initial energy of electron = E(t0): 2 • From Psyn = −dE=dt and energy E = γmec 2 dE 2 4 1 B −3 2 = −c1E with c1 := cσT 2 2 = 1:6 × 10 B [1=erg=s] dt 3 (mec ) 8π 2 • Integrate (dE=dt)=E = −c1 from t0 to t (t > t0) to give 1 t 1 1 − = −c (t − t ) viz: − = c (t − t ) E 1 0 E(t) E(t ) 1 0 t0 0 E(t0) − E(t) = c1(t − t0) E(t)E(t0) so finally E(t ) 1 E(t) = 0 / for large times 1 + c1(t − t0)E(t0) t • Note: Timescale tcool for cooling down to 1/2 of initial energy E0 at t0 = 0: E0 E0 1 E(tcool) = = to give tcool(E0) = . 2 1+c1tcoolE0 c1E0 7 6 Single Electron Spectrum Remember: Radiated power as a function of time, P (t) / j~a(t)j2, reflects time- dependence of acceleration a(t). Distribution of this radiation over frequency, 2 Pν is called spectrum and reflects power (Fourier) spectrum of a (t). For synchrotron emission, observed radiation varies with time due to relativistic beaming (in addition to fundamental gyro-motion), so expect observed spectrum to contain power at frequencies much higher than ΩL=2π. � � � Figure 1: A charged particle in slow gyro-motion around a magnetic field will emit an approximate dipole pattern with maximum in the direction of motion v. For higher speeds, aberration (beaming) changes pattern to an asymmetric shape, with most emission contained in forward cone of half opening angle ∼ 1/γ. 8 • Aberration (Lecture 4): Radiation of relativistic particle is strongly beamed in forward direction, focused into cone of opening angle∆ α ' 2/γ [rad]. �� r ��� 1 2 Observer • Observer only sees radiation once per obit when particle velocity is within 2/γ of line of sight. ∆tobs:=duration for which beam stays within line of sight ) Observer sees narrow pulse of duration ∆tobs once per gyro-period ) Most of radiation power must appear at frequency !c ∼ 1=∆tobs. (Basic motion is periodic, t / 1=ΩL (discrete spectrum), but frequency spacing ΩL between successive harmonics becomes so narrow as to be negligible, so can treat as continuous spectrum) 9 Estimate Characteristic Frequency νc or !c = 2πνc: • Path length ∆s for which we see radiation: ∆s = r ∆α = r (2/γ). • Radius of curvature r of path given by equation of motion: ∆~v q γm = ~v × B~ ∆t c With j∆~vj ' v ∆α and ∆s ' v∆t, i.e. ∆t = ∆s=v, this gives: ∆α qB sin θ ' ∆s γmcv so that ∆s γmcv r = ' ∆α qB sin θ (Note: 3-dim r is different from projected/perpendicular Larmor radius) 2vmc • Path length ∆s ' r ∆α = r (2/γ) = qB sin θ. 10 • Times t1 and t2 at which particles passes point 1 and 2 related by ∆s = v(t2 − t1), thus 2mc t − t ' 2 1 qB sin θ • Travel time effects: Arrival time difference is less than (t2 − t1) by ∆s=c: v ∆s h vi 2mc h vi 1 + c mc ∆tobs = (t2−t1)− = (t2−t1) 1 − = 1 − v ' 2 c c qB sin θ c 1 + c γ qB sin θ ) Observed pulse has duration shortened by1 /γ2 Note: Converting from rad to Hz, observed duration would be = 2π∆tobs • Spectrum (in terms of frequencies) will be broad with characteristic angular frequency (!ch = 2πνch) at which significant part of power is radiated: 1 2 qB 2 3 !ch ' = γ sin θ = γ Ω0 sin θ = γ ΩL sin θ tobs mc with cyclotron frequency Ω0 := qB=(mc). 11 7 Resulting Observed Electric Field y t i s n e t n i width: � ������� time ����� � The observed time-dependent E-Field, E(t), from an electron is a sequence of 1 1 pulses of width∆ tobs ' 2 = 3 , separated in time by ∆t ∼ 2π=ΩL. Note: γ Ω0 γ ΩL Larmor frequency ΩL = eB/γmc =: Ω0/γ. 12 8 Single Particle Synchrotron Spectrum Pν(γ) Approach: Notation alert: Note E here refers to electric field, not particle energy! • Want spectrum=power per unit frequency, Pν(γ) [erg/s/Hz], R 1 2 β2e4B2 2 2 satisfying 0 Pν(γ)dν = Psyn = 3 m2c3 γ sin θ. ) Find dimensionless frequency distribution function F~(ν/νc) with 1 Pν(γ) = Psyn F~(ν/νc) νc satisfying Z 1 Z 1 1 νc ! F~(ν/νc)dν = F~(x)dx = 1: νc 0 νc 0 • Start by Fourier-transforming E(t) / g(!ct) (get from full L-W potentials), 1 Z 1 E^(!) = E(t)ei!tdt ; 2π −∞ R 1 ^ −i!t where inverse transform E(t) = −∞ E(!)e d!. 13 • Remember Parseval' theorem: Z 1 Z 1 E2(t)dt = 2π jE^(!)j2d! −∞ −∞ • Total energy W per unit area in pulse (Poynting theorem): dW c Z 1 Z 1 = E2(t)dt = c jE^(!)j2d! dA 4π −∞ 0 so total energy per area per unit time per unit frequency: dW 1 dW c ≡ = jE^(!)j2 dAdtd! T dAd! T • Spectrum Pν(γ) = 2πP!(γ) from 1 dW c Z P (γ) ≡ = r2jE^(!)j2dΩ ! T d! T 2 where r distance, T = 2π=ΩL period and solid angle element dΩ = dA=r , using ∆θ ∼ 1/γ. 14 ) Result: Total synchrotron power per unit frequency [erg/s/Hz] for a single electron with Lorentz factor γ (relativistic case β ' 1) p e3B sin θ ν Pν(γ) = 3 2 F mc νc R 1 0 0 0:3 −x with F (x) := x x K5=3(x )dx ' 1:8 x e , K5=3 modified Bessel function of order 5=3, and γ entering via 3 eB 3 ν = γ2 sin θ = γ2 Ω sin θ c 4π mc 4π 0 1 Note: By convention, F is normalized somewhat differently compared to F~, i.e. R F (x)dx = 8pπ , introducing the 0 9 3 slightly different numerical value in Pν. 1 1 0.8 0.6 ) ) c c ! ! / / ! ! F( F( 0.1 0.4 0.2 0.01 0 0.001 0.01 0.1 1 0 0.5 1 1.5 2 2.5 3 3.5 4 !/!c !/!c Figure 2: Single Particle Spectrum F (ν/νc) as function of normalized frequency ν/νc. Left: Log-log repre- sentation. Right: Normal representation. The spectrum is continuous and has a maximum at νmax = 0:29νc. 15 9 Example: Synchrotron Emission from the Crab Nebula Figure 3: Crab Nebula (d ∼ 2 kpc) caused by SN explosion in 1054 A.D.