COSMIC RAY ACCELERATION and TRANSPORT LECTURE 2

1. COSMIC RAY ACCELERATION andTRANSPORTLECTURE 2 Pasquale Blasi INAF/Arcetri Astrophysical Observatory 4th School on Cosmic Rays and Astrophysics UFABC - Santo André - São Paulo – Brazil

2. Acceleration of charged particles The presence of non-thermal particles is deduced in a myriad of situations in Nature (from the solar wind to the AGNs, from SNRs to GRBs, from pulsars to mQSO) PARTICLE ACCELERATION BUT usually (through not always) in the same regions there is evidence for Thermalized plasmas, therefore the questions arises WHICH PROCESSES DETERMINE WHETHER A PARTICLE IS GOING TO BE ACCELERATED OR RATHER BE THERMALIZED ?

3. Acceleration of charged particles All acceleration processes we are aware of are of electro-magnetic nature – but magnetic fields DO NOT MAKE WORK on charged particles WHAT IS THE ORIGIN OF THE ELECTRIC FIELDS THAT PRODUCE ACCELERATION? ACCELERATION MECHANISMS ARE CLASSIFIED ACCORDING WITH THE ORIGIN OF THE ELECTRIC FIELDS REGULAR ACCELERATION STOCHASTIC ACCELERATION

4. REGULAR ACCELERATION Large mean scale electric fields are produced on some spatial scale Lreg DIFFICULT TO CREATE NET ELECTRIC FIELD IN ASTROPHYSICS BECAUSE OF HIGH CONDUCTIVITY, BUT SOME EXCEPTIONS: Unipolar Inductor Magnetic reconnection

5. STOCHASTIC ACCELERATION Most astrophysical acceleration processes belong to this class The stochastic electric field may result from random fluctuations on a typical scale Lst but with random orientations so that on average the field vanishes. If both regular and stochastic acceleration occur:

6. 2nd order Fermi Acceleration (Fermi, 1949) LOSSES AND GAINS ARE PRESENT BUT DO NOT COMPENSATE EXACTLY PROBABILITY OFENCOUNTER

7. WHY WOULD MAGNETIC CLOUDS ACCELERATE PARTICLES? WHERE ARE THE ELECTRIC FIELDS? In the Fermi example the electric fields are induced by the motion of the magnetized moving clouds In reality we need to go back to our example of motion of a charged particles in a group of Alfven waves…what if we do not sit in the reference frame of the waves? As usual:

8. Where you should recall that: Therefore: The time for diffusion in momentum space is then: DIFFUSION IN SPACE IMPLIES THAT A (2nd ORDER) DIFFUSION IN MOMENTUM TAKES PLACE (ACCELERATION)

9. A PRIMER ON SHOCK WAVES Forσ~10-25cm2 and density n~1 cm-3 the typical interaction length is ~3 Mpc >> than the typical size of astrophysical objects and even Larger than the Galaxy!!! COLLISIONLESS SHOCKS UPSTREAM DOWNSTREAM -∞ 0 +∞ U1 U2

10. STATIONARY SHOCKS 4 M→∞ M→∞ M→∞ SHOCK WAVES ARE MAINLY HEATING MACHINES!

11. BOUNCING BETWEEN APPROACHING MIRRORS UPSTREAM DOWNSTREAM V=U1-U2>0 Relative velocity INITIAL ENERGY DOWNS: E -1< μ <0 -∞ 0 +∞ U1 U2 0< μ’ <1 TOTAL FLUX

12. FIRST ORDER A FEW IMPORTANT POINTS: • There are no configurations that lead to losses • The mean energy gain is now first order in V • The energy gain is basically independent of any detail on how particles scatter back and forth!

13. RETURN PROBABILITIES AND SPECTRUM OF ACCELERATED PARTICLES UPSTREAM DOWNSTREAM -∞ 0 +∞ Return Probability from Downstream U1 U2 HIGH PROBABILITY OF RETURN FROM DOWNSTREAM BUTTENDS TO ZERO FOR HIGH U2

14. ENERGY GAIN: E0→ E1 → E2 → --- →EK=[1+(4/3)V]K E0 N0 → N1=N0*Pret → --- → NK=N0*PretK

15. Putting these two expressions together we get: Therefore: THE SLOPE OF THE DIFFERENTIAL SPECTRUM WILL BE γ+1=(r+2)/(r-1)→ 2 FOR r→4 (STRONG SHOCK)

16. THE TRANSPORT EQUATION APPROACH Integrating around the shock: Integrating from upstr. infinity to 0-: and requiring homogeneity downstream: UP DOWN -∞ 0- 0+ +∞ U1 U2

17. THE TRANSPORT EQUATION APPROACH INTEGRATION OF THIS SIMPLE EQUATION GIVES: NOTE THAT THIS IS IN P SPACE NAMELY N(p)dp=4π p2 f(p)dp Therefore the slope is 3r/(r-1) • THE SPECTRUM OF ACCELERATED PARTICLES IS A POWER LAW EXTENDING TO INFINITE MOMENTA • THE SLOPE DEPENDS UNIQUELY ON THE COMPRESSION FACTOR AND IS INDEPENDENT OF THE DIFFUSION PROPERTIES • INJECTION IS TREATED AS A FREE PARAMETER WHICH DETERMINES THE NORMALIZATION

18. TEST PARTICLE SPECTRUM

19. SOME IMPORTANT COMMENTS • THE STATIONARY PROBLEM DOES NOT ALLOW TO HAVE A MAX MOMENTUM! • THE NORMALIZATION IS ARBITRARY THEREFORE THERE IS NO CONTROL ON THE AMOUNT OF ENERGY IN CR • AND YET IT HAS BEEN OBTAINED IN THE TEST PARTICLE APPROXIMATION • THE SOLUTION DOES NOT DEPEND ON WHAT IS THE MECHANISM THAT CAUSES PARTICLES TO BOUNCE BACK AND FORTH • FOR STRONG SHOCKS THE SPECTRUM IS UNIVERSAL AND CLOSE TO E-2 • IT HAS BEEN IMPLICITELY ASSUMED THAT WHATEVER SCATTERS THE PARTICLES IS AT REST (OR SLOW) IN THE FLUID FRAME

20. A FREE ESCAPE BOUNDARY CONDITION THE ESCAPE OF PARTICLES AT X=X0 CAN BE SIMULATED BY TAKING UP DOWN THIS REFLECTS IN AN EXP CUTOFF AT SOME MAX MOMENTUM x0 ESCAPE FLUX TOWARDS UPSTREAM INFINITY!!!

21. ESCAPE FLUX IN TEST PARTICLE THEORY FOR D(E) PROPORTIONAL TO E (BOHM DIFFUSION):

22. SOME FOOD FOR THOUGHT • WHAT DETERMINES THE MAX MOMENTUM IN REALITY? • IF THE RETURN PROBABILITY FROM UPSTREAM IS UNITY, WHAT ARE COSMIC RAYS MADE OF? • ARE WE SURE THAT THE 10-20% EFFICIENCY WE NEED FOR SNR TO BE THE SOURCES OF GALACTIC CR ARE STILL COMPATIBLE WITH THE TEST PARTICLE REGIME?

23. MAXIMUM MOMENTUM OF ACCELERATED PARTICLES THE ACCELERATION TIME IS GIVEN BY: AND SHOULD BE COMPARED WITH THE AGE OF THE ACCELERATOR, FOR INSTANCE A SUPERNOVA REMNANT AS AN ESTIMATE: Emax

24. IF THE SHOCK IS PROPAGATING IN THE ISM ONE WOULD BE TEMPTED TO ASSUME D(E)=Dgal(E) WHERE TYPICALLY: A=(1-10) 1027 cm2/s α=0.3-0.5 FOR ALL CHOICES OF PARAMETERS THE MAX ENERGY OBTAINED IN THIS WAY IS FRACTIONS OF GeV, THEREFORE IRRELEVANT !!! …BUT IT WOULD BE HIGHER IF D(E) WERE MUCH SMALLER…CAN IT HAPPEN?

25. DIFFERENT PHASES OF A SNR THERE IS AN INITIAL PERIOD DURING WHICH THE SHELL OF THE SN EXPANDS FREELY (FREE EXPANSION PHASE -BALLISTIC MOTION): MASS OF THE EJECTA: Mej TOTAL KINETIC ENERGY: E51 FREE EXPANSION VELOCITY:

26. BUT THE SHOCK SWEEPS THE MATERIAL IN FRONT OF IT AND AT SOME POINT IT ACCUMULATES ENOUGH MATERIAL TO SLOW DOWN THE EXPANDING SHELL: SEDOV PHASE: The sound speed in the ISM is about 106 cm/s STRONG SHOCK

27. Simple implications During free expansion the shock fron moves with constant speed Therefore its position scales with t The diffusion front moves proportional to t1/2 During the free expansion phase the particles are not allowed to Leave the acceleration box, which is the reason why the maximum Energy increases

28. During the Sedov-Taylor expansion the radius of the blast waves Grows as t2/5, slower than the diffusion front The Maximum energy of accelerated particles decreases with Time This is the phase during which the particles can possibly Become Cosmic rays

29. OVERLAP OF ESCAPE FLUXES: A SIMPLE ESTIMATE BE VERY CAREFUL…THIS IS JUST A WAY TO SHOW HOW YOU GET ROUGHLY A POWER LAW BUT SUMMING NON-POWER LAWS. MORE DETAILED CALC’S SHOW DEPARTURES FROM THIS SIMPLE TREND