The Effect of Magnetic Turbulence Energy Spectra and Pickup Ions on the Heating of the Solar Wind C. S. Ng Geophysical Institute, University of Alaska Fairbanks A. Bhattacharjee, P. A. Isenberg, D. Munsi, and C. W. Smith Space Science Center, University of New Hampshire 9th Annual International Astrophysics Conference: Pickup Ions Throughout the Heliosphere and Beyond 15 March, 2010 Work partially supported by NASA, DOE, and NSF
MHD turbulence Energy cascade possible only if two Alfvén wave packets propagating in opposite directions collide Dimensionless parameter: Time scales: Eddy turn-over time τ N ~ 1/ kv k χ k 1/ 2 E k 1/ 2 V A 1 = v k / V A << 1 Alfvén time τ A ~ 1/ kv A = χτ N << τ N Energy cascade time τ E ~ τ A / χ 2 = τ N / χ >> τ N Kolmogorov cascade rate: ε K ~ v k 2 / τ N ~ k 5/ 2 E k 3/ 2 IK cascade rate: ε IK ~ v k 2 / τ E ~ ε K χ ~ k 3 E k 2 V A 1 << ε K [Iroshnikov 1963, Kraichnan 1965]
Solar wind turbulence model The steady state solar wind turbulence model developed by [Matthaeus et al. 1994, 1996] and later developments: Assumptions: Steady state Radially expanding solar wind with uniform speed V sw 1D (radial position r) Turbulence characterized by two fields: fluctuation velocity (Z) and correlation length (λ) Kolmogorov type cascade rate Solar wind (proton) temperature can be calculated passively
Solar wind turbulence model The steady state solar wind turbulence model developed by [Matthaeus et al. 1994, 1996] and later developments: dz 2 dr = AZ 2 r αz 3 λv SW + Q V SW dλ dr = Cλ r + βz V SW βλq V SW Z 2 dt dr = 4T 3r + mαz 3 3k B λv SW Z 2 : average turbulence energy with Z + 2 = Z 2 λ: turbulence correlation length T: solar wind proton temperature Q: turbulence generation rate due to pickup ions (more on the form of this operator later)
Solar wind turbulence model vs observations From [Smith et al. 2001]
Temperature comparison with pickup ions From [Isenberg et al. 2003]
Solar wind model with IK cascade The solar wind turbulence model changed to [Matthaeus et al. 1994, Hossain et al. 1995]: dz 2 dr = AZ 2 αz 4 + Q r λv SW V A V SW dλ dr = Cλ r + β Z 4 V SW V A Z 2 : average turbulence energy with Z 2 2 + = Z λ: turbulence correlation length T: solar wind temperature Q: turbulence generation rate due to pickup ions (also depends on the turbulence energy spectrum) 1/ 3 dt dr = 4T 3r + mαz 4 3k B λv SW V A [Ng et al., J. Geophys. Res., 115, A02101 (2010)] βλq V SW Z 2
Turbulence generation due to pickup ions Following Isenberg et al. [2003] and Isenberg [2005] Q = ζ V 2 SW n dn dt dn /dt = N 0 ν 0 (r E /r) 2 exp( L /r) V SW 2 / n : initial kinetic energy per pickup proton N 0 : neutral hydrogen density at the termination shock ν 0 : ionization rate L: scale of the ionization cavity (6 ~ 8 AU used in calculations) ζ: fraction of newly ionized pickup proton energy that generates waves (requires most theoretical considerations)
Wave energy generation fraction ζ(δ) =1 Δ + V 4 SW 2 Δ + V SW 1 Δ 1 Δ v 4 (µ)s(µ)dµ v 2 (µ)s(µ)dµ dv dµ = V j I j (k r ) 1 µv j j µv W j v j I j (k r ) 1 µv j µv W j v v: particle velocity µ: cosθ 2 1 Δ = Z / 3 1/ 2 V A v(µ = Δ) = V SW S(µ) = δv S(µ = Δ) = 1 k r = Ω V 3 j µvv 2 j + µvv 2 A = 0 W j = µν V j I(k) = A(r) k 5 / 3 ζ > 0 if v V SW <1, or I(k) = A(r) k 3 / 2 ζ 0 as dv dµ 0 2µvV A 2 2V j 2 2µvV j +V A 2 dependence on spectral index, but not A [Isenberg 2005]
Cyclotron resonance condition k r = Ω ω(k) = ±kv A 1+ ω Ω µν V j k 1 < k 2 dv dµ V ji j (k r ) 1 µv j µv W j v j < 0 usually in a turbulence cascade V 2 > 0 since I(k 1 ) > I(k 2 ) ζ is smaller for the IK spectrum than the Kolmogorov spectrum since dv/dµ is smaller due to I(k 2 ) getting closer to I(k 1 ) V 1 < 0 [Isenberg 2005]
ζ for Kolmogorov vs IK ζ Δ Δ Kolmogorov V A /V SW IK V A /V SW
ζ for Kolmogorov vs IK ζ [Isenberg 2005] ζ V A /V SW = 0.075 Kolmogorov IK Δ
Comparisons with observations IK with Q Adiabatic cooling Kolmogorov with Q IK without Q Kolmogorov without Q cf. [Isenberg et al. 2003]
Comparisons with observations IK with Q IK without Q Kolmogorov without Q Kolmogorov with Q cf. [Smith et al. 2001]
Comparisons with observations IK with Q IK without Q Kolmogorov without Q Kolmogorov with Q cf. [Smith et al. 2001]
Comparisons with observations Kolmogorov with Q IK with Q cf. [Smith et al. 2006] theory
Conclusion Although a solar wind model based on Kolmogorov cascade compares well with proton temperature observations when contributions from pickup ions are included, observations are consistent with IK cascade are less sensitive to the contribution of pickup ions. Pickup ions generate less turbulence under IK cascade. The prediction of the correlation length based on IK cascade appears to agree better with data than Kolmogorov cascade with pickup ions. Results obtained based on parameters/methods used in Smith et al. 2001, Isenberg et al. 2003, and Smith et al. 2006. Solar wind heating based on IK cascade is a possible alternative to Kolmogorov based theories and should be investigated more extensively. More precise description of how pickup ions generate turbulence might also change the prediction of solar wind temperature, especially in the outer heliosphere.