Kinetic Vlasov Simulations of collisionless magnetic Reconnection

A fully kinetic Vlasov simulation of the Geospace Environment Modeling (GEM) Magnetic Reconnection Challenge is presented. Good agreement is found with previous kinetic simulations using particle in cell (PIC) codes, confirming both the PIC and the Vlasov code. In the latter the complete distribution functions $f_k$ ($k=i,e$) are discretised on a numerical grid in phase space. In contrast to PIC simulations, the Vlasov code does not suffer from numerical noise and allows a more detailed investigation of the distribution functions. The role of the different contributions of Ohm's law are compared by calculating each of the terms from the moments of the $f_k$. The important role of the off--diagonal elements of the electron pressure tensor could be confirmed. The inductive electric field at the X--Line is found to be dominated by the non--gyrotropic electron pressure, while the bulk electron inertia is of minor importance. Detailed analysis of the electron distribution function within the diffusion region reveals the kinetic origin of the non--gyrotropic terms

The interaction between transpolar arcs and cusp spots

Transpolar arcs and cusp spots are both auroral phenomena which occur when the interplanetary magnetic field is northward. Transpolar arcs are associated with magnetic reconnection in the magnetotail, which closes magnetic flux and results in a "wedge" of closed flux which remains trapped, embedded in the magnetotail lobe. The cusp spot is an indicator of lobe reconnection at the high-latitude magnetopause; in its simplest case, lobe reconnection redistributes open flux without resulting in any net change in the open flux content of the magnetosphere. We present observations of the two phenomena interacting--i.e., a transpolar arc intersecting a cusp spot during part of its lifetime. The significance of this observation is that lobe reconnection can have the effect of opening closed magnetotail flux. We argue that such events should not be rare

Toward detailed prominence seismology - I. Computing accurate 2.5D magnetohydrodynamic equilibria

Context. Prominence seismology exploits our knowledge of the linear eigenoscillations for representative magnetohydro- dynamic models of filaments. To date, highly idealized models for prominences have been used, especially with respect to the overall magnetic configurations. Aims. We initiate a more systematic survey of filament wave modes, where we consider full multi-dimensional models with twisted magnetic fields representative of the surrounding magnetic flux rope. This requires the ability to compute accurate 2.5 dimensional magnetohydrodynamic equilibria that balance Lorentz forces, gravity, and pressure gradients, while containing density enhancements (static or in motion). Methods. The governing extended Grad-Shafranov equation is discussed, along with an analytic prediction for circular flux ropes for the Shafranov shift of the central magnetic axis due to gravity. Numerical equilibria are computed with a finite element-based code, demonstrating fourth order accuracy on an explicitly known, non-trivial test case. Results. The code is then used to construct more realistic prominence equilibria, for all three possible choices of a free flux-function. We quantify the influence of gravity, and generate cool condensations in hot cavities, as well as multi- layered prominences. Conclusions. The internal flux rope equilibria computed here have the prerequisite numerical accuracy to allow a yet more advanced analysis of the complete spectrum of linear magnetohydrodynamic perturbations, as will be demonstrated in the companion paper.Comment: Accepted by Astronomy & Astrophysics, 15 pages, 15 figure

Magnetohydrodynamics dynamical relaxation of coronal magnetic fields. II. 2D magnetic X-points

We provide a valid magnetohydrostatic equilibrium from the collapse of a 2D X-point in the presence of a finite plasma pressure, in which the current density is not simply concentrated in an infinitesimally thin, one-dimensional current sheet, as found in force-free solutions. In particular, we wish to determine if a finite pressure current sheet will still involve a singular current, and if so, what is the nature of the singularity. We use a full MHD code, with the resistivity set to zero, so that reconnection is not allowed, to run a series of experiments in which an X-point is perturbed and then is allowed to relax towards an equilibrium, via real, viscous damping forces. Changes to the magnitude of the perturbation and the initial plasma pressure are investigated systematically. The final state found in our experiments is a "quasi-static" equilibrium where the viscous relaxation has completely ended, but the peak current density at the null increases very slowly following an asymptotic regime towards an infinite time singularity. Using a high grid resolution allows us to resolve the current structures in this state both in width and length. In comparison with the well known pressureless studies, the system does not evolve towards a thin current sheet, but concentrates the current at the null and the separatrices. The growth rate of the singularity is found to be tD, with 0 < D < 1. This rate depends directly on the initial plasma pressure, and decreases as the pressure is increased. At the end of our study, we present an analytical description of the system in a quasi-static non-singular equilibrium at a given time, in which a finite thick current layer has formed at the null

The Structure and Evolution of Magnetized Cloud Cores in a Zero--Density Background

Molecular-line observations of star-forming cloud cores indicate that they are not the flattened structures traditionally considered by theory. Rather, they are elongated, perhaps in the direction of their internal magnetic field. We are thus motivated to consider the structure and evolution of axisymmetric, magnetized clouds that start from a variety of initial states, both flattened (oblate) and elongated (prolate). We devise a new technique, dubbed the $q$-method, that allows us to construct magnetostatic equilibria of any specified shape. We find, in agreement with previous authors, that the field lines in oblate clouds bend inward. However, those in prolate clouds bow outward, confining the structures through magnetic tension. We next follow the quasi-static evolution of these clouds via ambipolar diffusion, under the assumption of constant core mass. An oblate cloud either relaxes to a magnetically force-free sphere or, if sufficiently massive, flattens along its polar axis as its central density runs away. A prolate cloud always relaxes to a sphere of modest central density. We finally consider the evolution of an initially spherical cloud subject to the tidal gravity of neighboring bodies. Although the structure constricts equatorially, it also shortens along the pole, so that it ultimately flattens on the way to collapse. In summary, none of our initial states can evolve to the point of collapse while maintaining an elongated shape. We speculate that this situation will change once we allow the cloud to gain mass from its environment.Comment: 19 pages, plus 20 postscript figures. Accepted by Ap

In situ evidence for the structure of the magnetic null in a 3D reconnection event in the Earth's magnetotail

Magnetic reconnection is one of the most important processes in astrophysical, space and laboratory plasmas. Identifying the structure around the point at which the magnetic field lines break and subsequently reform, known as the magnetic null point, is crucial to improving our understanding reconnection. But owing to the inherently three-dimensional nature of this process, magnetic nulls are only detectable through measurements obtained simultaneously from at least four points in space. Using data collected by the four spacecraft of the Cluster constellation as they traversed a diffusion region in the Earth's magnetotail on 15 September, 2001, we report here the first in situ evidence for the structure of an isolated magnetic null. The results indicate that it has a positive-spiral structure whose spatial extent is of the same order as the local ion inertial length scale, suggesting that the Hall effect could play an important role in 3D reconnection dynamics.Comment: 14 pages, 4 figure

Observations of whistler mode waves with nonlinear parallel electric fields near the dayside magnetic reconnection separatrix by the Magnetospheric Multiscale mission

We show observations from the Magnetospheric Multiscale (MMS) mission of whistler mode waves in the Earth's low-latitude boundary layer (LLBL) during a magnetic reconnection event. The waves propagated obliquely to the magnetic field toward the X line and were confined to the edge of a southward jet in the LLBL. Bipolar parallel electric fields interpreted as electrostatic solitary waves (ESW) are observed intermittently and appear to be in phase with the parallel component of the whistler oscillations. The polarity of the ESWs suggests that if they propagate with the waves, they are electron enhancements as opposed to electron holes. The reduced electron distribution shows a shoulder in the distribution for parallel velocities between 17,000 and 22,000 km/s, which persisted during the interval when ESWs were observed, and is near the phase velocity of the whistlers. This shoulder can drive Langmuir waves, which were observed in the high-frequency parallel electric field data