Plasma, a gas of charged particles exhibiting collective behaviour, can be found everywhere in our vast Universe. The characteristics of plasma in very distant parts of the Universe can be similar to characteristics in our solar system and near-Earth space. We can therefore gain an understanding of what happens in astrophysical plasmas by studying processes occurring in near Earth space, an environment much easier to reach. Large volumes in space are filled with plasma and when different plasmas interact distinct boundaries are often created. Many important physical processes, for example particle acceleration, occur at these boundaries. Thus, it is very important to study and understand such boundaries. In Paper I we study magnetic nulls, regions of vanishing magnetic fields, that form inside boundaries separating plasmas with different magnetic field orientations. For the first time, a statistical study of magnetic nulls in the Earth’s nightside magnetosphere has been done by using simultaneous measurements from all four Cluster spacecraft. We find that magnetic nulls occur both in the magnetopause and the magnetotail. In addition, we introduce a method to determine the reliability of the type identification of the observed nulls. In the manuscript of Paper II we study a different boundary, the shocked solar wind plasma in the magnetosheath, using the new Magnetospheric Multiscale mission. We show that a region of strong current in the form of a current sheet is forming inside the turbulent magnetosheath behind a quasi-parallel shock. The strong current sheet can be related to the jets with extreme dynamic pressure, several times that of the undisturbed solar wind dynamic pressure. The current sheet is also associated with electron acceleration parallel to the background magnetic field. In addition, the current sheet satisfies the Walén relation suggesting that plasmas on both sides of the current region are magnetically connected. We speculate on the formation mechanisms of the current sheet and the physical processes inside and around the current sheet.

The work in this thesis focuses primarily on equilibrium and stability properties of collisionless current sheet models, in particular of the force-free Harris sheet model. A detailed investigation is carried out into the properties of the distribution function found by Harrison and Neukirch (Physical Review Letters 102, 135003, 2009) for the force-free Harris sheet, which is so far the only known nonlinear force-free Vlasov-Maxwell equilibrium. Exact conditions on the parameters of the distribution function are found, which show when it can be single or multi-peaked in two of the velocity space directions. This is important because it may have implications for the stability of the equilibrium. One major aim of this thesis is to find new force-free equilibrium distribution functions. By using a new method which is different from that of Harrison and Neukirch, it is possible to find a complete family of distribution functions for the force-free Harris sheet, which includes the Harrison and Neukirch distribution function (Physical Review Letters 102, 135003, 2009). Each member of this family has a different dependence on the particle energy, although the dependence on the canonical momenta remains the same. Three detailed analytical examples are presented. Other possibilities for finding further collisionless force-free equilibrium distribution functions have been explored, but were unsuccessful. The first linear stability analysis of the Harrison and Neukirch equilibrium distribution function is then carried out, concentrating on macroscopic instabilities, and considering two-dimensional perturbations only. The analysis is based on the technique of integration over unperturbed orbits. Similarly to the Harris sheet case (Nuovo Cimento, 23:115, 1962), this is only possible by using approximations to the exact orbits, which are unknown. Furthermore, the approximations for the Harris sheet case cannot be used for the force-free Harris sheet, and so new techniques have to be developed in order to make analytical progress. Full analytical expressions for the perturbed current density are derived but, for the sake of simplicity, only the long wavelength limit is investigated. The dependence of the stability on various equilibrium parameters is investigated.

In this thesis we investigate current sheets in the solar corona. The well known 1D model for the tearing mode instability is presented, before progressing to 2D where we introduce a non-uniform resistivity. The effect this has on growth rates is investigated and we find that the inclusion of the non-uniform term in η cause a decrease in the growth rate of the dominant mode. Analytical approximations and numerical simulations are then used to model current sheet formation by considering two distinct experiments. First, a magnetic field is sheared in two directions, perpendicular to each other. A twisted current layer is formed and we find that as we increase grid resolution, the maximum current increases, the width of the current layer decreases and the total current in the layer is approximately constant. This, together with the residual Lorentz force calculated, suggests that a current sheet is trying to form. The current layer then starts to fragment. By considering the parallel electric field and calculating the perpendicular vorticity, we find evidence of reconnection. The resulting temperatures easily reach the required coronal values. The second set of simulations carried out model an initially straight magnetic field which is stressed by elliptical boundary motions. A highly twisted current layer is formed and analysis of the energetics, current structures, magnetic field and the resulting temperatures is carried out. Results are similar in nature to that of the shearing experiment.

In this thesis examples of translationally invariant one-dimensional (1D) Vlasov-Maxwell (VM) equilibria are investigated. The 1D VM equilibrium equations are equivalent to the motion of a pseudoparticle in a conservative pseudopotential, with the pseudopotential being proportional to one of the diagonal components of the plasma pressure tensor. A necessary condition on the pseudopotential (plasma pressure) to allow for force-free 1D VM equilibria is formulated. It is shown that linear force-free 1D VM solutions correspond to the case where the pseudopotential is an attractive central potential. The pseudopotential for the force-free Harris sheet is found and a Fourier transform method is used to find the corresponding distribution function. The solution is extended to include a family of equilibria that describe the transition between the Harris sheet and the force-free Harris sheet. These equilibria are used in 2.5D particle-in-cell simulations of magnetic reconnection. The structure of the diffusion region is compared for simulations starting from anti-parallel magnetic field configurations with different strengths of guide field and self-consistent linear and non-linear force-free magnetic fields. It is shown that gradients of off-diagonal components of the electron pressure tensor are the dominant terms that give rise to the reconnection electric field. The typical scale length of the electron pressure tensor components in the weak guide field case is of the order of the electron bounce widths in a field reversal. In the strong guide field case the scale length reduces to the electron Larmor radius in the guide magnetic field.

No description available.