Magnetic Reconnection: Some Answers and Open Questions

1. Magnetic Reconnection: Some Answers and Open Questions Amitava Bhattacharjee Space Science Center and Center for Magnetic Self-Organization University of New Hampshire NESSC Meeting, Durham, New Hampshire, January 9, 2007

2. Collaborators Naoki Bessho, John Dorelli, Terry Forbes, Kai Germaschewski, Chung-Sang Ng, David Pontin, Joachim Raeder, Hong-Ang Yang, UNH John Greene, GA Shuanghui Hu, UC-Irvine Zhiwei Ma, Unversity of Iowa Xiaogang Wang, Beijing University CMSO : John Finn, Chris Hegna, Hantao Ji, Yi-Min Hwang, Vladimir Mirnov, Stewart Prager, Dalton Schnack, Carl Sovinec, Masaaki Yamada, Ping Zhu, Ellen Zweibel

3. Outline • Classical steady-state models of Sweet-Parker and Petschek. • Impulsive, time-dependent reconnection dynamics: comparative studies of magnetospheric substorms, and solar flares. • Scaling of the reconnection rate. • 3D reconnection: Separator reconnection Current sheet formation and reconnection in geometries without nulls or closed field lines. I will not discuss the important subjects of particle acceleration, and turbulent reconnection.

4. Classical (2D) Steady-State Models of Reconnection Sweet-Parker [Sweet 1958, Parker 1957] Geometry of reconnection layer : Y-points [Syrovatskii 1971] Length of the reconnection layer is of the order of the system size >> width  Reconnection time scale Solar flares: Too long to account for solar flares!

5. Q. Why is Sweet-Parker reconnection so slow? A. Geometry Conservation relations of mass, energy, and flux Petschek [1964] Geometry of reconnection layer: X-point Length (<< L) is of the order of the width Solar flares:

6. Numerical Simulations of the Petschek Model [Sato and Hayashi 1979, Ugai 1984, Biskamp 1986, Forbes and Priest 1987, Scholer 1989, Yan, Lee and Priest 1993, Ma et al. 1995, Uzdensky and Kulsrud 2000, Breslau and Jardin 2003, Malyshkin, Linde and Kulsrud 2005] • Conclusions • Petschek model is not realizable in high-S plasmas, unless the resistivity is locally strongly enhanced at the X-point. • In the absence of such anomalous enhancement, the reconnection layer evolves dynamically to form Y-points and realize a Sweet-Parker regime.

7. 2D coronal loop : high-Lundquist number resistive MHD simulation T = 0 T = 30 [Ma, Ng, Wang, and Bhattacharjee 1995]

8. Biskamp’s Critique of the Petschek Model Magnetic energy accumulates upstream of the current sheet (flux pileup). increases with increasing Lundquist number; this directly contradicts the Petschek model, which requires  Flux pile-up: Litvinenko, Forbes and Priest [1996], Craig, Henton, and Rickard[1993], Dorelli and Birn [2001], Simakov, Chacon, and Knoll [2006]

9. Impulsive Reconnection: The Trigger Problem Dynamics exhibits an impulsiveness, that is, a sudden change in the time-derivative of the reconnection rate. The magnetic configuration evolves slowly for a long period of time, only to undergo a sudden dynamical change over a much shorter period of time. Dynamics is characterized by the formation of near-singular current sheets in finite time. Examples Sawtooth oscillations in tokamaks Magnetospheric substorms Impulsive solar flares

10. Sawtooth events in MST [Almagri and the MST group 2003]

12. Current Disruption in the Near-Earth Magnetotail [Ohtani, Kokubun, and Russell 1992]

13. Time Profile of Reconnection Rates for X3 Flareobserved by TRACE Linear Plot Log Plot [Saba, Tarbell, and Gaeng 2003]

14. Hall MHD (or Two-Fluid) Model and the Generalized Ohm’s Law In high-S plasmas, when the width of the thin current sheet ( ) satisfies (or if there is a guide field) “collisionless” terms in the generalized Ohm’s law cannot be ignored. Generalized Ohm’s law (dimensionless form) Electron skin depth Ion skin depth Electron beta

15. Forced Magnetic Reconnection Due to Inward Boundary Flows Magnetic field Inward flows at the boundaries , Two simulations:Resistive MHD versus Hall MHD [Ma and Bhattacharjee 1996] For another perspective, with similar conclusions, see Horiuchi and Sato [1994] and more recently, Cassak, Shay, and Drake [2005]

18. ….. Hall ___ Resistive

19. Transition from Collisional to Collisionless Regimes in MRX

20. “Standard Model”

22. Scaling of Reconnection • Observations: Yokoyama et al. (2001), Isobe et al. (2002) This Workshop: Lin et al. (2004), Noglik et al, (2004) Reconnection rates : inflow velocity/Alfven velocity~0.001-0.1. Hall MHD reconnection scaling depends on the ion skin depth, the Lundquist number (weakly), and system size. The range of reconnection rates in Hall MHD and PIC simulations ~ 0.01-0.1.

23. Observations (Ohtani et al. 1992) Hall MHD Simulation [Ma and Bhattacharjee 1998]

24. Scaling of the collisionless reconnection rate How does it scale with ion/electron skin depth, resistivity, plasma beta, and system size? y Vin - + i e - e + x i out-of-plane magnetic field, Bz

25. Scaling is a controversial subject • Is the reconnection rate insensitive to the details of the electron layer (current sheet layer), and controlled by ions? • The GEM Challenge Perspective • [Birn et al. 2001; Mandt et al. 1994, Shay and Drake 1998, Hesse et al. 1999, Rogers et al. 2001, Pritchett 2001, Ricci et al. 2002] • Reconnection is insensitive to the mechanism that breaks field lines (electron inertia or resistivity). • In the presence of Hall currents, whistler waves mediate reconnection. The characteristic outflow speed is the whistler phase speed (based on the upstream magnetic field). • The inflow velocity where (= system size). This rate is independent of .

26. Scaling is a controversial subject • The whistler waves generate an out-of plane quadrupolar magnetic field (seen in MRX/PPPL, SSX/Swarthmore as well as in situ satellite observations in space). • The ratio of the horizontal electron outflow to the horizontal magnetic perturbation scales as k for the dispersive whistler (or kinetic Alfven) wave. • How does the reconnection rate scale with the system size? The length of the reconnection layer . Reconnection rate is a “universal constant”, [Shay et al. 2001, 2004, Huba and Rudakov 2004].

27. Scaling is a controversial subject The GEM perspective is not universally accepted. An alternate point of view provides evidence that: • Reconnection is not a universal constant, and depends on system parameters (such as ion/electron skin depth, plasma beta, boundary conditions) • Reconnection rate is not independent of the system size, and in fact, often decreases as the system size increases. [Wang et al. 2001, 2006, Fitzpatrick 2004, Bhattacharjee et al. 2005, Daughton et al. 2006] Three examples : (1) Forced reconnection without guide field (2) Undriven reconnection with guide field (3) Undriven reconnection with open boundaries

28. Scaling of steady-state reconnection rate in adriven Harris sheet (without guide field) Two-fluid simulations as di is varied [Wang et al., 2006] Confirms analytical scaling law [Wang, Bhattacharjee, and Ma 2001]

29. Effect of Open Boundary Conditions and System Size [Daughton, Scudder and Karimabadi 2006]

30. Hall currents cancel exactly in electron-positron plasmas: no whistlers Generalized Ohm’s law does not contain the Hall term. [Bessho and Bhattacharjee 2006]

31. Open questions on scaling • While collisionless reconnection models do provide a clear mechanism for the onset of impulsive and fast reconnection, issues of scaling remain wide open. • Is reconnection controlled entirely by ions? Or by electrons? Or is it a hybrid of ion and electron parameters, as well as boundary conditions? • What is the role of periodic and open boundary conditions? • Is the reconnection rate independent of the system size? If it isn’t, the theory faces great challenges in working for solar coronal plasmas while it may be adequate for magnetotail and fusion plasmas. • Is the Hall current essential for fast reconnection? • Role of turbulence--the Eyink-Aluie theorem (2006): “Frozen-in-flux condition in a turbulent quasi-ideal plasma will occur if the sites of current and vortex sheets intersect one another.”

32. 3D Reconnection: Geometry and Dynamics

33. Resistive Tearing Modes in 2D/2.5D Geometry x Neutral line at y=0 y Time Scales Lundquist Number Tearing modes [Furth, Killeen and Rosenbluth 1963, Coppi, Galvao, Pellat, Rosenbluth, and Rutherford 1976]

34. Magnetic nulls in 3D play the role of X-points in 2D Spine Fan [Greene 1988, Lau and Finn 1990]

35. Dungey’s Model for Southward and Northward IMF “Magnetopause phenomena are more complicated as a result of merging. This is why I no longer work on the magnetopause.” -- J. W. Dungey [Dungey 1961, 1963]

36. Observation of magnetic null in space[Xiao, Pu, Wang et al., Nature, 2006]

37. Equilbrium Current Density, J [Greene and Miller 1995]

38. Equilbrium B-field

39. B-lines in a spherically tearing state [Hu, Bhattacharjee, Dorelli, and Greene 2004]

40. Equilibrium Perturbed Field lines penetrating the spherical tearing surface

41. Magnetopause reconnection topology in global MHD Greene, J., Locating three-dimensional roots by a bisection method, JCP, 1992. [Dorelli, Bhattacharjee, and Raeder 2006]

42. Magnetopause reconnection geometry [Dorelli, Bhattacharjee, and Raeder, 2006] Current density is concentrated at the intersection of the two separatrix surfaces; the X line has collapsed into a double Y line topology. Such ribbons were discussed by Longcope and Cowley [1996], but for force-free fields.

43. Solar corona astron.berkeley.edu/~jrg/ ay202/img1731.gif www.geophys.washington.edu/ Space/gifs/yokohflscl.gif

44. Solar corona: heating problem photosphere corona Temperature Density Time scale Magnetic fields (~100G) --- role in heating? ~ ~  Alfvén wave  current sheets

45. Parker's Model (1972) Straighten a curved magnetic loop Photosphere

46. Reduced MHD equations low  limit of MHD

47. Magnetostatic equilibrium Key Question: What is the nature of the solutions, given a sufficiently complicated footpoint mapping? Objections to Parker’s claim: van Ballegooijen [1985], Longcope and Strauss [1994]. Cowley et al. [1997] and others.

48. A theorem onParker's model For any given footpoint mapping connected with the identity mapping, there is at most one smooth equilibrium. Caveat: A proof based on reduced MHD equations, periodic boundary condition in x [Related results by Aly 2005 and Low 2006]

49. Implication An unstable but smooth equilibrium cannot relax to a second smooth equilibrium, hence must have current sheets.

50. Possible current sheet topology in line-tied geometry at Quasi-Separatrix Layers (QSL)