SOLAR WIND ION AND ELECTRON DISTRIBUTION FUNCTIONS AND THE TRANSITION FROM FLUID TO KINETIC BEHAVIOR JUSTIN C. KASPER HARVARD-SMITHSONIAN CENTER FOR ASTROPHYSICS GYPW01, Isaac Newton Institute, July 2010
Overview The solar wind as a laboratory to understand plasma dynamics As a function of beta As a function of collision/dynamic timescales When can we describe these plasmas with Fluid MHD (equal temperatures, velocities, ) Kinetic equations (non-maxwellian VDFs, wave-particle coupling, ) Example of current research Instabilities Ion Heating Ion Electron Interactions
Plasma physics in the heliosphere The corona is not in hydrostatic equilibrium and a supersonic solar wind is generated. The solar wind is highly structured, with streams and shocks In interplanetary space, our spacecraft are much smaller than the Debye length, so we can measure the VDFs of ions and e- s. Credit: SOHO (ESA/NASA) We cannot prepare shots, but we can treat the solar wind as an ergodic system that explores a wide range of states. Our starting point is a massive database of accessible plasma conditions
A fleet of spacecraft explore the heliosphere
The Wind spacecraft Launched in fall 1994 and still operating Understand solar wind upstream of Earth Lots of extra fuel new locations (L1, L2, distant prograde orbits, double lunar swing-by, ) Spins out of ecliptic plane at 3 seconds Excellent complement of stable, accurate, well-understood in situ instruments Ions and Electrons (>5.3M 3D VDFs) Electromagnetic fluctuations
Kinetic physics of the solar wind Earliest solar wind measurements showed that ion and electron velocity distribution functions (phase space density) are almost, but not quite, Maxwellian What keeps them from being Maxwellian? Very low densities and high temperatures at 1 AU a typical proton experiences a weak Coulomb interaction once per AU Interactions between electromagnetic fields and particles of a particular speed, not as a fluid OK, then why are they even close to Maxwellian? Kinetic physics of the solar wind When is MHD valid and when do we need to include details of the velocity distribution? What are VDF signatures of wave-particle coupling, heating, dissipation, acceleration? Testing and pushing theory with large and accurate samples of the solar wind (millions of individual spectral measurements) Schwen and Marsch Kasper et al., 2008
Collisional relaxation in the solar wind
0.3 AU 0.5 AU 1 AU Kinetic properties of the solar wind Slow Fast Kinetic or non-thermal features Field-aligned anisotropies Beams Different temperatures Evolution Stronger closer to Sun Stronger in fast wind Paradigm Slow wind is Maxwellian and behaves as single fluid Fast wind has strong nonthermal aspects Slow/fluid, Fast/kinetic Schwen and Marsch
Kinetic properties ordered by speed Consider three non-thermal effects: Alpha/proton temperature ratio Proton temperature anisotropy Differential flow 2D histograms of these parameters as a function of solar wind speed Clearly some dependence on speed But why sort by speed? Consider a thermalization timescale instead.
Coulomb collisional age Coulomb collisional age A c is the number of Coulomb collisions experienced by typical particle over propagation from Sun R A c U n 3/ 2 T Fast wind is hotter and less dense Previous work has shown that large A c leads to smaller differential flow and T p ~ T α (Mckenzie 1987, Marsch 1983, Neugebauer 1976) > 4 orders of magnitude A c
Wind speed or collisional age?
Proton temperature anisotropy
What regulates temperature anisotropy? Double adiabatic expansion (Chew, Goldberger & Low, 1965) T B, T 2 n / B 2 R j T T j j B n 3 2 Measured anisotropy Predicted variation Clack et al., 2003
What regulates temperature anisotropy? T Double adiabatic expansion B, T 2 n / B Coulomb relaxation A c 2 n 3/ 2 T R j T T j j B n 3 2
What regulates temperature anisotropy? T Double adiabatic expansion B, T 2 n / Coulomb relaxation Instabilities (Gary, Vinas, Hellinger, Marsch) Firehose R B Cyclotron, Mirror 2 n A c T p R p 1 1 3/ 2 R S p p p S p p p j T T j j B n 3 2 If R p exceeds threshold, instability begins to grow, generating waves that make the velocity distribution function more isotropic The larger the anisotropy, the stronger the wave growth and scattering Instability Limit is curve for a predefined growth rate of the instability
What regulates temperature anisotropy? T Double adiabatic expansion B, T 2 n / Coulomb relaxation Instabilities (Gary, Vinas, Hellinger, Marsch) Firehose R B Mirror, Cyclotron 2 n A c T p R p 1 1 3/ 2 R S p p p S p p p T T j j Cyclotron-resonant heating R p 1 j B n 3 2 R p 1 Cyclotron resonant heating Cyclotron Instability Expansion Mirror Instability Firehose Instability p
Distribution of observed anisotropy at 1 AU Cyclotron resonant heating R p Cyclotron Instability Mirror Instability 1 Firehose Instabilit y p Expansion Kasper et al., 2003
Range of anisotropy limited by instabilities R p 1 S p p p Kasper et al., 2003 R p 1 S p p p
Identify the dominant instabilities using the statistical distribution of the solar wind Cyclotron Mirror Hellinger et al., 2006 Parallel Firehose Oblique Firehose
Power in magnetic fluctuations Instabilities should generate EM waves Look at average power of magnetic fluctuations across R p - p plane Bale et al. (2009)
Multiple sources of power at small scales Bale et al. (2005) Bale et al. (2009)
Anisotropic heating of protons
r L 2 /BT /B In situ signatures of heating process Large perpendicular temperature anisotropies Proton magnetic moment not conserved as the solar wind expands Marsch et al., 1982 Marsch, 1991
Radial evolution and heating T B, T 2 n / B 2 R j B n 3 2 Matteini et al., 2007
Maruca et al. (2009) Proton temperature near instability thresholds 2.0 R p 1.0 0.0 Liu et al. (2006)
Proton temperature near instability thresholds
Components of proton temperature near instability thresholds Kasper et al. (in prep)
Relative heating of ions
Helium/hydrogen temperature ratio T T p T 4T p T / / Tp m m p T 4T p T / / Tp m m p Kasper et al. (2008)
The Alfven-Cyclotron Resonant Interaction V B /C A -0.2 1.0 0.5 0.0 0.3 0.8 After Hollweg and Isenberg (2002) Pitchangle diffusion V B /C A Turbulence in the corona and solar wind transports power to shorter scales, creating a spectrum of Alfven waves Several models of dissipation have focused on the absorption of Alfvencyclotron fluctuations in a plasma with multiple ion species (Hu 1999, Cranmer 2003, Hellinger 2005, Isenberg 2007) In a frame moving with the wave, particles will stochastically diffuse along the surface v B2 + V B2 = const Heavier ions resonate with slightly faster waves so if diffusion is sufficiently fast a heavy ion can achieve T i >(m i /m p )T p
Test for Alfvén-cyclotron resonant heating in a two species plasma V protons C A alphas V Relative heating may also be regulated by differential flow (Gary 2004, Gary 2005, Hellinger 2005) Significant because heavier ions often flow faster than protons. In general the differential flow V p = V -V p reaches but rarely exceeds the Alfven speed, C A (Neugebauer 1976, Marsch 1982, Reisenfeld 2001) Prediction For V p /C A ~0, helium will have a stronger cyclotron resonance than hydrogen, and thus will be heated faster and develop larger temperature anisotropies. On the other hand, if V p /C A is increased, the helium ions come out of resonance with the waves and the protons are heated instead.
Test for Alfvén-cyclotron resonant heating in a two species plasma V V protons C A C A alphas V V Relative heating may also be regulated by differential flow (Gary 2004, Gary 2005, Hellinger 2005) Significant because heavier ions often flow faster than protons. In general the differential flow V p = V -V p reaches but rarely exceeds the Alfven speed, C A (Neugebauer 1976, Marsch 1982, Reisenfeld 2001) Prediction For V p /C A < 0.15, helium will have a stronger cyclotron resonance than hydrogen, and thus will be heated faster and develop larger temperature anisotropies. On the other hand, if V p /C A is increased, the helium ions come out of resonance with the waves and the protons are heated instead.
T α /T p as a function of A c and V αp For fixed A c, T α /T p increases as V αp drops At small A c and V αp, the average value of T α /T p ~ 6 In fact, higher T α /T p is seen even at high A c
Anisotropy of heating and dependence on Vαp/CA
Ion-Electron Interactions
Relative heating of ions and electrons Cranmer et al. (2008)
Prescription Start with empirical dependence of temperature and heat flux on distance: Insert into equations for internal energy: Extract Q p and Q e
Relative heating with distance Red curve is best-fit fraction of power dissipated into protons as a function of distance from the Sun Change assumed mean speed Raise Coulomb collision frequency R [A.U.] R [A.U.]
What explains this partition?
Next Steps Models Data Spacecraft
Theoretical models and data analysis Models Is the mirror instability really more important than the cyclotron instability at regulating temperature anisotropy? Predictions of wave power Data analysis Ion-electron interactions Plasma at higher cadence Combine electromagnetic fields with plasma
Solar Probe Plus (SP+)
Advantages of SPP orbit
Rich kinetic dataset
Conclusions We can use large sets of precise solar wind measurements to conduct sensitive and quantitative kinetic physics experiments Microphysics is important at all solar wind speeds Instabilities Experienced by protons, electrons (alphas, minor ions?) Limit anisotropy, generate waves, heat plasma Heating mechanisms Alfven-ion cyclotron heating is occurring in solar wind Upcoming missions to the inner heliosphere such as Solar Orbiter and Solar Probe can conduct groundbreaking kinetic physics experiments Kinetic physics will be important at all solar wind speeds once we get close in Need to be able to make fast and accurate measurements of bulk properties of solar wind ions and electrons, electric and magnetic field Radial gradients important operate these missions beyond encounters