Cosmic Rays and Thermal Instability

1. Cosmic Rays and Thermal Instability T. W. Hartquist, A. Y. Wagner, S. A. E. G. Falle, J. M. Pittard, S. Van Loo

2. Outline • The thermal instability of a non-magnetized, uniform static fluid, in the absence of cosmic rays • The thermal instability and cloud formation • The effect of cosmic rays on the thermal instability of a uniform medium • Instability of radiative, non-magnetic shocks • The effect of cosmic rays on radiative shocks

3. Uniform, Static, No CRs (Field 1965) The difference between the heating rate per unit volume and the cooling rate per unit volume Thermal equilibrium

4. Thermal Equilibrium P – n Relationship

5. Isobaric Perturbations Stable if

6. Isentropic Perturbations Stable If

7. Isobaric and Isentropic If L = Λ(T)n, then the above criterion gives stability if d ln(Λ)/d ln(T) > 1 There is also a criterion for the growth of isentropic perturbations (i.e. sound waves); for the same assumption about L, it gives stability if d ln(Λ)/d ln (T) > -3/2 for γ = 5/3

8. Shock-induced formation of • Giant Molecular Clouds Sven Van Loo Collaborators: Sam Falle and Tom Hartquist

9. Overview • Introduction • Thermal properties of ISM • Formation of molecular clouds • Conclusions

10. Introduction • Hierarchical density structure in molecular clouds • Emission line maps of the Rosette Molecular Cloud (Blitz 1987) • MCs that do not harbour any young stars are rare • Old stellar associations (few Myr) are devoid of molecular gas • ⇒ Cloud and core formation are entangled • Not homogeneous,but highly structured • Stars embedded in dense cores

11. Cloud formation • Compression + • Thermal processes in diffuse atomic gas: • Heating: photoelectric heating, • cosmic rays, soft X-rays, … • Cooling: fine-structure lines, • electron recombination, • resonance lines, … • 2 stable phases in which • heating balances cooling: • Rarefied, warm gas (w; T > 6102 K) • Dense, cold gas (c; T < 313 K) Net cooling w c Net heating (Wolfire et al. 1995; Sanchez-Salcedo et al. 2002)

12. Cloud formation: flow-driven Flow-driven formation or colliding streams e.g. expanding and colliding supershells Heitsch, Stone & Hartmann (2009) Hennebelle et al. (2008) • Collision region prone to instabilities, i.e. KH, RT, NTSI • Turbulent shocked layer • Fragmentation into cold clumps • Structure depends strongly on magnetic field (both orientation and magnitude) //v ⊥v B//v

13. Cloud formation: shock-induced • Shock-induced formation • e.g. shocks and winds sweeping up material • Similar processes as flow-driven • Can explain different cloud morphologies • e.g. filamentary, head-tail,… • ⇒ Shock-cloud interaction Lim,Falle & Hartquist (2005) Inutsuka & Koyama (2006) Van Loo et al. (2007) W3 GMC (Bretherton 2003) Previous work: 2D: adiabatic: MacLow et al. (1994), Nakamura et al. (2006) radiative: Fragile et al. (2005), Van Loo et al. (2007) 3D: adiabatic: Stone & Norman (1992), Shin, Stone & Snyder (2008) radiative: Leão et al. (2009) (nearly isothermal), Van Loo et al. (2010)

14. Numerical simulation • Interaction of shock with initially warm, thermally stable cloud (n = 0.45 cm-3, T = 6788K, R = 200pc) which is in pressure equilibrium with hot ionised gas (n = 0.01 cm-3, T = 282500K) and β = 1. • Numerics: • Ideal MHD code with AMR (Falle 1991): • 2nd order Godunov scheme with linear Riemann solver • + divergence cleaning algorithm (Dedner et al. 2002) • Include cooling as source function • Resolution: 640/120 cells (2D/3D) across initial cloud radius • (120 cells = resolution for adiabatic convergence in 2 and 3D)

15. Slow-mode shock Fast-mode shock Dynamical evolution: 2D Mach 2.5 (but similar for other values)

16. Dynamical evolution: 2D From 2D simulation 12CO • Typical GMC values: n ≈ 20 cm-3 & R ≈ 50 pc • High-mass clumps in boundary and low-mass • clumps inside cloud  precursors of stars • Similar to observations of e.g. W3 GMC (Bretherton 2003)

17. Results: 2D Dependency on Mach number • Weak shocks (M ≤ 2): • NOT magnetically dominated • Strong shocks (M > 4): • formation time too short, because time-scale for formation of H2 is a few Myr Volume fraction of cloud for which β=Pg/Pm < 0.1 M = 2.5 M = 5 M = 1.5 • only moderate-strength shocks can produce clouds similar to GMCs(obs: β ~ 0.03-0.6)

18. Small (constant ram pressure) Large (significant decrease in ram pressure ) Results: small vs. large • Pressure decrease behind shock, e.g. blast waves

19. Stronger Shocks Possible

20. Phase diagram log(p/k) log(n) ⇒ Rapid condensation at boundary Dynamical evolution: 3D parallel Parallel shock Geometry

21. Dynamical evolution:3D oblique Oblique shock ~45o Geometry Phase diagram log(p/k) log(n) ⇒ Condensation along equilibrium curve

22. Results: 3D • Cloud properties: • large differences between parallel and oblique/perpendicular • Oblique/perpendicular → HI clouds; Parallel → molecular clouds • Ideal conditions (β < 1) for MHD waves to produce large density contrasts

23. Results: 3D • Column density • Large column density >1021 cm-2 • Some filaments, but not much substructure Parallel Oblique

24. Substructure • Effect of increasing resolution: overall the same but more detail 640 cpr 120 cpr

25. Future work • Shock interacting with multiple clouds • Low resolution simulation (60 cpr) of 2 identical clouds overrun by an oblique shock (~45o) • Qualitative differences: • Shape • Density structure • Still need further study…

26. Conclusions • Magnetically-dominated clouds form due to thermal instability and compression by weak or moderately-strong shocks • The time-lag between cloud and core formation is short

27. Uniform Media – Incoporating Cosmic Rays • A. Y. Wagner, S. A. E. G. Falle, T. W. Hartquist, J. M. Pittard (2005) • CR Pressure Gradient Term in Momentum Density Eqn. and Corresponding Term in Energy Density Eqn.

28. Additional Parameters

29. First of Three Conditions • Similar to condition for isobaric perturbations

30. Other Two Conditions Analogous to conditions for isentropic perturbations Obviously Complicated

31. Limit of Large Φ and Small Diffusion Coefficient ϕ big compared to 1 and absolute value of any other wavenumber divided by cosmic ray diffusion wavenumber (a/χ). Stable if Roughly satisfied for all values of k with small enough magnitudes

32. Limit of Small Φ and Diffusion Coefficient Compared to 1; magnitudes of ratios of all other wavenumbers to the cosmic ray wavenumber are small compared to 1 (corresponds to big diffusion coefficient). Stable if

33. Thermal Instability • Falle (1975); Langer, Chanmugum, and Shaviv (1981); Imamura, Wolfe, and Durisen (1984) showed that single fluid, non-magnetic, radiative shocks are unstable if the logarithmic temperature derivative, α, of the energy radiated per unit time per unit volume is less than a critical value • Pittard, Dobson, Durisen, Dyson, Hartquist, and O’Brien (2005) investigated the dependence of thermal stability on Mach number and boundary conditions

34. Alpha = -1.5, M = 1.4, 2, 3, and 5

35. Do Magnetic Fields Affect the Themal Instability? • Interstellar magnetic pressure is comparable to interstellar thermal pressure (about 1 eV/cc) • Immediately behind a strong shock propagating perpendicular to the magnetic field, the magnetic pressure increases by a factor of 16 • Immediately behind a strong shock the thermal pressure increases by roughly the Mach number squared

36. Magnetic pressure limits the ultimate compression behind a strong radiative shock, but it does not affect the thermal instability

37. How About Cosmic Rays? • In interstellar medium the pressure due to roughly GeV protons is comparable to the thermal pressure. • Krymskii (1977); Axford et al. (1977); Blandford and Ostriker (1978); Bell (1978) showed that shocks are the sites of first order Fermi acceleration of cosmic rays. • Studies were restricted to adiabatic shocks but indicated that cosmic ray pressure is great enough to modify the thermal fluid flow.

38. Two Fluid Model of Cosmic Ray Modified Adiabatic Shocks • Völk, Drury, and McKenzie (1984) used such a model to study the possible cosmic ray acceleration efficiency • Thermal fluid momentum equation includes the gradient of the cosmic ray pressure • Thermal fluid equation for its entire energy includes a corresponding term containing cosmic ray pressure

39. Equation governing cosmic ray pressure derived from appropriate momentum moment of cosmic ray transport equation including diffusion – diffusion coefficient is a weighted mean • Concluded that for a large range of parameter space most ram pressure is converted into cosmic ray pressure and that the compression factor is 7 rather than 4 behind a strong shock

40. Two Fluid Model of Cosmic Ray Modified Radiative Shocks • Developed by Wagner, Falle, Hartquist, and Pittard (2006)

41. Cosmic Ray Pressure Held Constant Over Whole Grid Until t = 0

42. Problems • Compression is much less than observed • Too high of a fraction of ram pressure goes into cosmic ray pressure which is inconsistent with comparable interstellar themal and cosmic ray pressures

43. Possible Solution • Drury and Falle (1986) showed that if the length scale over which the cosmic ray pressure changes is too small compared to the diffusion length an acoustic instability occurs • Wagner, Falle, and Hartquist (2007, 2009) assumed that energy transfer from cosmic rays to thermal fluid then occurs

44. Including Acoustic Instability Induced Energy Transfer

45. Tycho Optical FeaturesWagner, Lee, Falle, Hartquist, Raymond (2009)