5. SYNCHROTRON RADIATION 1

Transcription

1 5. SYNCHROTRON RADIATION Charge motions in a static magnetic field Charged particles moving inside a static magnetic field continuously accelerate due to the Lorentz force and continuously emit radiation. In the case of non-relativistic motions, the emitted radiation is named cyclotron radiation and is due to the spiraling motion of particles and consequent acceleration. The radiation frequency is just that of the gyration qb/mc, where q and m are the electric charge and mass of the particle. This is a line emission, if we neglect broadening factors. The nature of the emission changes radically if the particle becomes relativistic, and the emission is then named synchrotron emission. The motion of a generic charge is ruled by the Lorentz dynamical equation of motion [5.0] and the energy conservation equation 1 This Section particularly benefited, in addition to the Rybicki & Lightman book, by lecture notes on High Energy Astrophysics by A. Cavaliere. 5.1

2 (the second term has been neglected because of the scalar product and because an assumption about the lack of significant large scale electric fields - that would be immediately erased by charge motions - was made). Then we use the second eq. into the former and get (see the figure below for an illustration of the relevant quantities) [5.1] The solution is then an helicoidal motion (see graph below), with the velocity along the direction of the B field keeping constant, and a uniform circular motion around that direction with rotational pulsation given by 5.

3 ω B qb γ mc =. [5.] 5. Total synchrotron emitted power The total emitted power by the particle is easily studied in the instantaneous reference frame of the electron (K'). Here the particle will not remain at rest for long time, but just for a short time interval and the electron motion in that frame will in any case be non-relativistic during this time. So we can use here the Larmor formula for the emitted power. Alternatively, one might use the full relativistic eq. (1.1), but this would be a quite more complicated derivation (see Sect. 5.7 below). Because energy and time transform in the same way from K' to the system K of the observer (and of the B field), we have dw = γdw '; dt = γt ' dw dw ' = P = P' = dt dt ' q d P = P' = a' 3 3c We need to consider how the vector acceleration transforms from one system (K) to the other (K'). It can be seen (e.g. Rybicki & Lightman) that the following relativistic rules apply if the particle is at rest in K', where parallel and orthogonal refer to the direction of the instantaneous velocity vector From this and [5.] we have a a' = aγ ; a' = a γ 3. qb = v = ω r = ω v γ mc B B q qb P= = r c B 3c γ mc 3 4 γ v 3 0 β γ v β β a c e 0 = ; sin = = r mc [5.3] 5.3

4 r 0 is the classical electron radius, q= e the electron charge. α is the pitch angle, that is the angle between the velocity and B field vectors, as shown in the figure: For an isotropic distribution of charge velocities, we average over α as: β β β = sin ad, such that 4π Ω= 3 P= r cβγb 3 0 or, if we use the Thomson cross section the magnetic field U B σ T = 8πr 3 and the energy density of 4 P = σcβγu, with U = B / 8π [5.4] 0 3 T B B For the cyclotron emission we have instead 4 v P = σ U 3 T c B. 5.3 Aberration, beaming, angular distribution of radiation Recall that in the lab reference system, in which the electron is seen moving with relativistic velocity, the emission is strongly enhanced (beamed) in the direction of the velocity. According to the relativistic aberration law, we have 5.4

5 where α is in the lab frame (K) and α' in the electron frame (K'). This implies that, except when. We can see this by a power series development If this condition does not verify, then (see below figure for an illustration of effect) Case of [5.5] Case of 5.5

6 5.4 The transition from the cyclotron to the synchrotron spectrum Let us consider here the evolution of a synchrotron spectrum at increasing electron velocity and energy. This makes a transition from a cyclotron line emission to the synchrotron continuum one that is illustrated by the following figures. 5.6

7 For another graphical representation we refer to Ghisellini (High Energy Astrophysics): From cyclo to synchro: if the emitting particle has a very small velocity, the observer sees a sinusoidal (in time) electric field E(t). Increasing the velocity the pattern becomes asymmetric, and the second harmonic appears. For 0 < β 1 the power in the second harmonic is a factor β less than the power in the first, see eq. (5.4). For relativistic particles, the pattern becomes strongly beamed, the emission is concentrated in the time ta. As a consequence the Fourier transformation of E(t) must contain many harmonics, and the power is concentrated in the harmonics of frequencies ν 1/ ta. Broadening of the harmonics due to several effects ensures that the spectrum in this case becomes continuous. Note that the fundamental harmonic becomes smaller at increasing γ (since νb 1/γ). The spectrum of emission of the first 0 harmonics of mildly relativistic cyclotron radiation for an electron with v = 0.4c (Bekefi, 1966). 5.7

8 Hard X-ray spectrum of the neutron star X observed with Beppo-SAX, showing the first four cyclotron harmonics. (From Santangelo et al. 1999) Cyclotron lines are observed more frequently in absorption in the high-energy spectra of neutron stars, like in the above figure, where the first 4 harmonics are observed. From this, the value of the (strong) stellar magnetic field can be derived from : where the last factor accounts for the gravitational redshift. 5.5 Qualitative behavior of the synchrotron spectrum of a single pulse We have seen that the spectral emission in the case of a cyclotron is concentrated essentially in a line of frequency. All the electron orbit is seen by the observer. Instead, in the synchrotron case, because of the relativistic beaming, the radiation field seen by the observer is concentrated inside a narrow cone of width around the electron velocity vector. As a consequence, the observer will see photon emission only from a small fraction of the orbit: therefore, the emitted spectrum will 5.8

9 include a much wider range of photon frequencies. All this is easily understood in terms of the Fourier transformation. Let us consider the above figure (from Ribicki & Lightman) representing the portion of the orbit seen by the observer. a is the radius of the orbit. The points 1 and correspond to the extremes of the sector seen: in 1 the emission cone enters the lineof-sight, in it exits it. The emission cone is shown in the two points. The angle Δθ is that subtending the sector, and is given by Δs between positions 1 and is then. The distance [5.6] Now the radius of the orbit can be obtained from [5.7] where is the gyration frequency of the relativistic electron from [5.] and where α is the pitch angle (obviously, the radius a increases when α decreases, that is the velocity approaches the B direction). So, from [5.6] and [5.7] [5.8] 5.9

10 s = v( t t ) In the observer frame K, the transit times t 1 and t are such that 1 t = t t1 = γω sinα This time interval however does not correspond to the times of arrival to the observer of the emitted photons, because these will be shorter by an amount equal to s/ c, the time that radiation emitted from 1 takes to travel from 1 to. So the arrival times will be arr v 1 t = 1 sin c sin 3 γωb a γ ωb a because (1-v/c) 1/γ. In conclusion, the observed pulse will have a temporal duration 3 1 γ times the period of the electron gyration. B From consideration of the Fourier transform, we then expect the presence of characteristic photon frequencies in the observed spectrum proportionally larger by a 3 factor γ, that is γ qb n ω γω a a γn π π πγmc 3 3 = = sin = sin = char B cycl We then define a critical frequency for the synchrotron emission as: 3 nc γωcycl sinα = [5.9] 5.10

11 So we expect significant emission to happen at this frequency. We will see later that indeed the peak synchrotron emission happen at this critical frequency. Let us see now some further scaling relations concerning the spectrum of a single radiation pulse by a single orbiting electron. Due to the relativistic beaming and the enormously increased emitted power in the direction of the electron velocity vector, the dependence of the electric field intensity E on the angle θ (the polar angle between the instantaneous velocity vector and the observer line-of-sight) goes under the product γθ. So also concerning the temporal dependence we have Et () g[ gθ () t], where t is the time in the K observer frame. This means that a given value of E will be obtained for large values of θ if γ is proportionally small, or vice-versa small θ if γ is large. Let us set the origin of times t and paths s when θ = 0. We will have: s θ s/ a; t s s (1 v / c) v c v γ s γθ γ s/ a ωb sina. Since v s tγ, so that γθ γ ( γ ωbsin a) t ωct v (arrival time of the signal) [5.10] So for the time dependences of the E field we can write where [ ] 1 v/c 1/γ, we have from [5.10] [ ω ], Et () f ct f x is an a-dimensional function, so that every dependence on t goes under the product ω. The Fourier transform of E will be c t + ( ω) ( ω )exp( ω ) x = ω E f t i t dt and by transforming t t c 1 + E ω ( ω) f( x)exp i x dx ω c ωc For the power spectrum of the pulse from the single charge we will have (after defining T as the electron rotation period): dw 1 dw = = P ( ω) = CF ω dt dω T dω ω c c [5.11] 5.11

12 that is, every dependence of the spectrum on frequency will go under the ratio ω. ω c F is an a-dimensional function and C is a quantity that does not depend on the frequency ω. We can easily determine the value of this constant and verify that it does not depend on γ either. Indeed if we write the total power P = P( ω) dω = C F ω dω ωcc F( x) dx; x ω ω = c ωc [5.1] if we now compare this relation with [5.3] and [5.4], in the ultra-relativistic limit we get qbγ β sin P = 3 3mc 3γ qbsinα wc = mc 4 From [5.1], ω ( ) c 0 α with the definition C = P F x dx and we finally obtain 3 3 qbsinα P( ω) = F ω π mc ω c [5.13] We need later to find an expression for the function F(x), which requires a rather complex calculation (Sect. 5.7 below). However, what we have found already allows us to obtain the key results of the synchrotron emission without going into complex calculations, just using the general expression [5.13]. 5.6 Emission by a non-thermal electron energy distribution We have seen in Sect. 5.4 that the emission spectrum of a single charge can be expressed as a function of an a-dimensional function F of the ratio ωω, that is the c whole dependence of the individual spectrum on electron energy is completely contained in that ratio. We need now to define the distribution of energies for a set of emitting electrons. This choice will be such that later comparison of predicted spectra will bring them into agreement with observations. We anticipate that this requires to introduce a power-law energy distribution of non-thermal electrons: 5.1

13 [5.14] Of course, the normalization constant C 1 may or may not depend on direction. Then We change the integration variable to and, remembering that ; ;, the Jacobian gets, so where the limits will depend on γ 1 and γ but if these limits are sufficiently wide, we have x 1 =0 and x =, such that [5.15] log S(n) -(p-1)/ log n Figure a) Emission spectrum for a power-law energy distribution for emitting particles. b) Observed radio spectrum of a radio-source, showing indeed the power-law behaviour, in addition to a low-frequency cutoff due to synchrotron self-absorption. The spectral index α of the emitted radiation is simply related to that of the energy distribution of electrons, and the synchrotron spectrum will be a simple power-law as in figure. 5.13

14 Typical values of the energy index p are typically between and 3. Then the radio spectral index will assume values around 0.7, as we will see in some detail below. How does this compare with data on the flux spectrum of cosmic rays in our Galaxy will be discussed in Sect. 7. As we see from the figures, cosmic rays in the Galaxy show indeed energy spectral slopes between and 3 over an enormous interval of particle energies. This is completely consistent with the observed synchrotron spectra of radio-galaxies (similar acceleration mechanisms appear to operate on such different scales). 5.14

15 More information on the cosmic rays that we directly collected near the Earth, as well as those inferred in Galactic and extragalactic radio-sources are reported in the following Section 6, including their presumed origin and composition. 5.7 The complete analysis So far we have developed an heuristic analysis of the synchrotron emission that allowed us to obtain in a simple way the main results. We need to mention here a more complete and rigorous treatment, that we take from Rybicki & Lightman. The geometry of the process that we adopt is illustrated in the figure below. The time considered is the retarded time as usual. It is assumed to be t'=0 at the moment the particle transits the origin of the system. The instantaneous velocity v is taken to be along the x axis. is the line-of-sight versor, is a versor orthogonal to v at t'=0 within the plane (x,y) containing the electron orbit (with curvature radius a) at that time. We define. Then the calculation starts from the general equation of electrodynamics in the "wave zone", eq. [1.1]. We can show that the total emitted power by the particle can be expressed as the sum of two contributions, one with the electric field parallel to and the other parallel to. We then start from the famous general eq. of electrodynamics, eq. (1.1) of Sect.1 (considering the acceleration term only, of course): 5.15

16 5.16

17 5.17

18 5.18

19 5.8 Synchrotron self-absorption and spectral cutoffs If there aren't cutoffs in the electron energy distribution that might affect the spectrum, it could happen in some circumstances - e.g. for high particle densities - that the photons emitted by synchrotron are re-absorbed by the synchrotron process itself before exiting the source. This implies a synchrotron optical depth τ>1. If τ>>1, we can make the argument extreme, by stating that the intensity of the photon field produced is I(ν)=flux of radiant energy (number of independent states of the photon field) (average energy of each state ε). From the black-body theory 5.19

20 3 we borrow the first factor being given by ns = n / c, while for calculating ε we can consider it to correspond to an equilibrium situation between the radiation field and the emitting particles. The appropriate value for ε is that of electrons emitting at frequency ν: ν ν = γω ε c B such that the characteristic energy of the emitting photons is the self-absorption regime, τ>>1, the spectrum will be: S( ν) ( ν) ν ν ν ε ν 1/. Overall, in 0.5 5/ I. [5.0] Figure 6. Overall spectrum of the synchrotron emission, including synchrotron selfabsorption at low frequencies. We immediately verify that the absorption coefficient a ν for synchrotron [defined by the transfer equation di(ν) = a ν I(ν) ds], can be obtained from the relation of the source function S(ν)=j ν /a ν with S(ν) ν.5, j ν ν -(p-1)/, hence α ν + ( p 4)/ ν [5.1] that, considering the typical values of p, shows a strong inverse dependence on frequency a ν ν 3. The meaning of this is that high values of optical depth, or the self-absorption, happen at low frequencies (Figs. 6). 5.9 Synchrotron polarization Let us first consider the simplest case of a single non-relativistic particle radiating cyclotron emission. The geometrical situation is illustrated in the figure, assuming for simplicity the particle simply rotating in a plane orthogonal to B. 5.0

21 The emitted radiation is described as usual in terms of the Larmor equation and of the time variations of the electric dipole. In our case, the dipole d d makes a rotation in the x-y plane, a motion that can be represented as the sum of two orthogonal oscillating dipoles 1 d d and d d. At the point A of the figure only the oscillating dipole d 1contributes radiation: here d d and d d dipoles contribute, hence the polarization is just linear. At point B both 1 the polarization is circular. At intermediate directions the polarization is intermediate, i.e. it is elliptical. Elliptical polarizations as observed above and below the orbit are opposite. All this is useful to understand the more complex synchrotron situation. In the case of the synchrotron emission, the flux is linearly polarized looking in the direction of maximal emission, that is of v, and is elliptically polarized when looking outside this direction. If the line-of-sight is outside the emission cone of the particle, the emission is essentially null. The situation is well illustrated by the figure below. The direction of the B field and of the particle motion are indicated, as well as the path made on the sky by the velocity vector and the trace of the emission cone in shaded grey. In analogy with what happens for the cyclotron, the polarization inside the emission cone in the inner and outer shaded region, with respect to the trace of the velocity vector, is elliptical, in one sense and in the opposite, respectively. So, for a roughly uniform distribution of electrons, these two elliptical polarization components cancel out, and only the linear polarization survives. If we want to quantitatively estimate the degree of polarization, we need to refer to the analysis in Sect. 5.7, where we have seen that the emitted wave in the observer 5.1

22 frame can be split into two components, one with the electric field ε parallel to the projection of the B field on the celestial sphere, and the other component ε orthogonal to it. The corresponding radiation powers are P ( ω) = dw dω and P ( ω) = dw dω. The degree of polarization can then be immediately obtained from P ( ω) P ( ω) Gx ( ) P ( ω) = = P ( ω) + P( ω) F( x) [5.4] where G(x) and F(x) are the two functions defined in Sect. 5.7, and F(x) corresponds to the total emitted power. Emission cone of the single particle, where the line-of-sight and direction of the magnetic field are indicated. Because of the geometry of the emission, the synchrotron light is polarized linearly along the direction of the velocity vector (i.e. along the cone). Polarization is almost frequency independent. Once integrated over all frequencies, we obtain an average (linear) polarization fraction of about P 75%, a rather high value. For an entire population of particles with power-law distribution function as in [5.14], we get Π= p + 1 p + 7/3 [5.5] 5.

23 which makes still a very large fraction assuming a typical value for the spectral index p=.5, P 70%. It has a slow dependence on the energy spectral index p, which is illustrated in the following graph. P p orientation of the B field on the sky. Note however that all this assumes a completely ordered B field, which is not a realistic assumption; the field may vary even chaotically inside the source or on its surface, which diminishes the polarization. Values of about a few tens of percent are often observed in radio-galaxies. The direction of polarization is orthogonal to the direction of the projection of the B field on the celestial sphere: polarization studies then allow us to estimate the The observational evidence for polarized light offers important indications about the physical origin of radiations in cases in which this might not be obvious. One such case is the X-ray emission by the Crab nebula (that is a very important astrophysical source and often used as a cosmic calibrator at high energies). This source will be further discussed in Sect Observations of the degree of linear polarization of these X-rays have demonstrated that the origin of this flux is synchrotron emission by a population of extremely energetic electrons. Note in any case that the presence of polarized light is not always evidence for nonthermal processes in operation: other phenomena can originate it, e.g. scattering and reflections, or even anisotropic absorption by polarized orientation of dust particles Validity limits of our treatment Two basic assumptions have been made in treating synchrotron. From the quantum mechanical limit condition that hν<<w, since,, we have that The calculation is simply:, after change of variable. 5.3

24 13 γ B < 5 10 Gauss [5.6] as the condition such that quantum corrections are negligible. Another limitation to our treatment comes from considering that the radiation reaction can become important and modify the helical orbit, hence the spectrum, if the particle in an orbit loses a too large fraction of its energy: or Since B πγ E = P t(1 orbit) = σtcbγ << γmc 4π eb mc =, γ B πmcν e γ B << r0 0 Gauss 3 πmcν e mhc << hν << 70 MeV! [5.7] e r e (See implications of this when discussing about the Crab Nebula emission below.) 5.11 Energy losses and synchrotron spectral evolution Generalities We consider here how the synchrotron spectrum evolves as an effect of the energy losses of the emitting particles. Every electron produces photons from transformation of its kinetic energy according to (note that the B component considered in the following is just that orthogonal to the line-of-sight): de B dγ σ B = σtc γ = γ dt 4π dt mc 4π γ dγ σ t T B = dt γ 0 γ mc 4π 0 T ; whose solution details how the electron energy changes with time: 5.4

25 [5.8] for sufficiently long times t. Vice versa, the characteristic time for the energy loss through radiation emission is:. [5.9] From all this we see that particles with higher energy loose energy faster than the others. These energy losses produce a progressive modification of the emitted spectrum: the electrons with the highest energies producing the photons of highest frequencies loose energy fast and move to lower energies in an energy-differential way, as shown in the figure. Consequently, we expect a bending of the synchrotron spectrum at the highest frequencies, progressing towards the lower energies with time. From [5.18] we get the synchrotron frequency where the energy loss effect starts to manifest itself: [5.30] Synchrotron spectral evolution The exact solution of the problem of estimating the spectral evolution by synchrotron emission requires to solve the energy continuity equation for the electron population, where N is the phase-space density of electrons: 5.5

26 [5.31] The detailed process then depends on the source function S, which details how the energetic electrons are injected into the volume, if this happens in a continuous way or impulsively at the start time. It can be shown (see details in Blumenthal & Gould 1974) that if the injection is continuous we have that the energy distribution spectrum bends by -1, hence the synchrotron spectrum bends by (see figure below, showing the evolution in the case of a continuous injection). 5.6

27 Instead, in the case of an impulsive injection of electrons at t=0, we have a qualitatively similar situation, in which the bending in the photon spectrum gets larger: instead of -0.5, it is now α = ( p + 5) This difference may be intuitively understood, considering that in this latter case there is no compensation to the energy loss by the continuous electron injection as in the previous case The case of the Crab Nebula A very interesting illustration of the energy loss effect by synchrotron emission is offered by a particular galactic object, which may be taken as a prototypical synchrotron source: the Crab Nebula. We have already mentioned that the evidence for linear polarization demonstrates the synchrotron nature of this emission up to the X-rays at least. The source (in particular the diffuse component within the network of filaments that constitute the nebula) has been observed in any detail at all frequencies. The Crab total broad-band spectrum is shown in the figure. 70 MeV Note among other things that the strong spectral break at about 10 Hz corresponds to the limit of validity of the synchrotron treatment of 70 MeV in [5.7]. For photons at higher energies the electrons lose their energy almost immediately during one single orbital round. This spectrum is composed of various sections, all with a power-law behavior. Let us concentrate on the region between radio and the optical-ir: above a power-law 5.7

28 spectrum with a rigorously constant index α. 0.3 extending from 10 7 to 10 1 Hz, we see a first bending at about Hz that can be interpreted as due to energy losses in the electrons. On the other hand, we know exactly the time of the birth of the nebula, A.D.1054 (from Chinese records), so from this we can calculate the magnetic field intensity permeating the nebula, from eq. [5.30]: B 5 10 / 10 ( ) /3 4 Gauss (note that this value depends weakly on the cutoff frequency and its uncertainty, * 1/3 B ν ). The bending of the spectrum is approximately α. 0.5, which is consistent with a continuous injection 3. There is in any case a stronger argument in favor of a continuous injection since the formation of the remnant 960 years ago: if the B field is B Gauss, energetic electrons emitting in X-rays loose energy very quickly: to emit at ν= Hz, they need to have but then they lose energy in a time 1/ ν mc Hz γ eb B 10 Gauss t (5 10 ) sec 1 yrs! Such energetic electrons need to be continuously refurbished to the nebula. This argument demonstrates that there is an engine in the source doing the job, and this is the pulsar (indeed, this argument was known even before the discovery of the pulsar 7 in 1966). Electrons with γ 5 10 have an energy of MeV=5 TeV, or about 100 erg, the kinetic energy of a flying mosquito. Instead the energies of the most energetic cosmic rays appearing in Fig. pag.14 correspond to that of a fast tennis ball! None of these can be reached in terrestrial accelerators (not even the lastgeneration LHC!, in particular 5 TeV is times more than the highest energy achieved in LHC in proton collisions note it is much more difficult to accelerate electrons, see Sect. 6), the main reason being the energy losses of the accelerated electrons due to their synchrotron emission: to reach higher energies we would need much lower B fields to constrain them, which would mean however accelerating rings of much larger radii. 0.3 A synchrotron spectrum with P( ν) ν would mean a flow of electrons injected 1.6 with an energy distribution N( ε). Cε. 1 1/ 3 There is also a second change in the spectral index in the far-uv that is possibly related to the acceleration mechanism. 5.8