Three-dimensional analytical magnetohydrostatic equilibria of rigidly rotating magnetospheres in cylindrical geometry

We present three-dimensional solutions of the magnetohydrostatic equations in the co-rotating frame of reference outside a magnetized rigidly rotating cylinder. We make no symmetry assumption for the magnetic field, but to be able to make analytical progress we neglect outflows and specify a particular form for the current density. The magnetohydrostatic equations can then be reduced to a single linear partial differential equation for a pseudo-potential $U$, from which the magnetic field can be calculated by differentiation. The equation for $U$ can be solved by standard methods. The solutions can also be used to determine the plasma pressure, density and temperature as functions of all three spatial coordinates. Despite the obvious limitations of this approach, it can for example be used as a simple tool to create three-dimensional models for the closed field line regions of rotating magnetospheres without rotational symmetry.Comment: 13 pages, 2 figures, accepted for publication by Geophysical and Astrophysical Fluid Dynamic

Collisionless distribution function for the relativistic force-free Harris sheet

A self-consistent collisionless distribution function for the relativistic analogue of the force-free Harris sheet is presented. This distribution function is the relativistic generalization of the distribution function for the non-relativistic collisionless force-free Harris sheet recently found by Harrison and Neukirch [Phys. Rev. Lett. 102, 135003 (2009)], as it has the same dependence on the particle energy and canonical momenta. We present a detailed calculation which shows that the proposed distribution function generates the required current density profile (and thus magnetic field profile) in a frame of reference in which the electric potential vanishes identically. The connection between the parameters of the distribution function and the macroscopic parameters such as the current sheet thickness is discussed. (C) 2012 American Institute of Physics. [doi: 10.1063/1.3677268]PostprintPeer reviewe

Particle energisation in a collapsing magnetic trap model : the relativistic regime

The authors acknowledge financial support by the UK’s Science and Technology Facilities Council through a Doctoral Training Grant (SEO) and Consolidated Grant ST/K000950/1 (SEO and TN).Context. In solar flares, a large number of charged particles is accelerated to high energies. By which physical processes this is achieved is one of the main open problems in solar physics. It has been suggested that during a flare, regions of the rapidly relaxing magnetic field can form a collapsing magnetic trap (CMT) and that this trap may contribute to particle energisation. Aims. In this Research Note we focus on a particular analytical CMT model based on kinematic magnetohydrodynamics. Previous investigations of particle acceleration for this CMT model focused on the non-relativistic energy regime. It is the specific aim of this Research Note to extend the previous work to relativistic particle energies. Methods. Particle orbits were calculated numerically using the relativistic guiding centre equations. We also calculated particle orbits using the non-relativistic guiding centre equations for comparison. Results. For mildly relativistic energies the relativistic and non-relativistic particle orbits mainly agree well, but clear deviations are seen for higher energies. In particular, the final particle energies obtained from the relativistic calculations are systematically lower than the energies reached from the corresponding non-relativistic calculations, and the mirror points of the relativistic orbits are systematically higher than for the corresponding non-relativistic orbits. Conclusions. While the overall behaviour of particle orbits in CMTs does not differ qualitatively when using the relativistic guiding centre equations, there are a few systematic quantitative differences between relativistic and non-relativistic particle dynamics.Publisher PDFPeer reviewe

Computing nonlinear force free coronal magnetic fields

International audienceKnowledge of the structure of the coronal magnetic field is important for our understanding of many solar activity phenomena, e.g. flares and CMEs. However, the direct measurement of coronal magnetic fields is not possible with present methods, and therefore the coronal field has to be extrapolated from photospheric measurements. Due to the low plasma beta the coronal magnetic field can usually be assumed to be approximately force free, with electric currents flowing along the magnetic field lines. There are both observational and theoretical reasons which suggest that at least prior to an eruption the coronal magnetic field is in a nonlinear force free state. Unfortunately the computation of nonlinear force free fields is way more difficult than potential or linear force free fields and analytic solutions are not generally available. We discuss several methods which have been proposed to compute nonlinear force free fields and focus particularly on an optimization method which has been suggested recently. We compare the numerical performance of a newly developed numerical code based on the optimization method with the performance of another code based on an MHD relaxation method if both codes are applied to the reconstruction of a semi-analytic nonlinear force-free solution. The optimization method has also been tested for cases where we add random noise to the perfect boundary conditions of the analytic solution, in this way mimicking the more realistic case where the boundary conditions are given by vector magnetogram data. We find that the convergence properties of the optimization method are affected by adding noise to the boundary data and we discuss possibilities to overcome this difficulty

Negative Specific Heat of a Magnetically Self-Confined Plasma Torus

It is shown that the thermodynamic maximum entropy principle predicts negative specific heat for a stationary magnetically self-confined current-carrying plasma torus. Implications for the magnetic self-confinement of fusion plasma are considered.Comment: 10p., LaTeX, 2 eps figure file

Particle-in-cell simulations of collisionless magnetic reconnection with a non-uniform guide field

Results are presented of a first study of collisionless magnetic reconnection starting from a recently found exact nonlinear force-free Vlasov–Maxwell equilibrium. The initial state has a Harris sheet magnetic field profile in one direction and a non-uniform guide field in a second direction, resulting in a spatially constant magnetic field strength as well as a constant initial plasma density and plasma pressure. It is found that the reconnection process initially resembles guide field reconnection, but that a gradual transition to anti-parallel reconnection happens as the system evolves. The time evolution of a number of plasma parameters is investigated, and the results are compared with simulations starting from a Harris sheet equilibrium and a Harris sheet plus constant guide field equilibrium

Force-free collisionless current sheet models with non-uniform temperature and density profiles

We present a class of one-dimensional, strictly neutral, Vlasov-Maxwell equilibrium distribution functions for force-free current sheets, with magnetic fields defined in terms of Jacobian elliptic functions, extending the results of Abraham-Shrauner [Phys. Plasmas 20, 102117 (2013)] to allow for non-uniform density and temperature profiles. To achieve this, we use an approach previously applied to the force-free Harris sheet by Kolotkov et al. [Phys. Plasmas 22, 112902 (2015)]. In one limit of the parameters, we recover the model of Kolotkov et al. [Phys. Plasmas 22, 112902 (2015)], while another limit gives a linear force-free field. We discuss conditions on the parameters such that the distribution functions are always positive and give expressions for the pressure, density, temperature, and bulk-flow velocities of the equilibrium, discussing the differences from previous models. We also present some illustrative plots of the distribution function in velocity space

Collisionless distribution functions for force-free current sheets: using a pressure transformation to lower the plasma beta

So far, only one distribution function giving rise to a collisionless nonlinear force-free current sheet equilibrium allowing for a plasma beta less than one is known (Allanson et al., Phys. Plasmas, vol. 22 (10), 2015, 102116; Allanson et al., J. Plasma Phys., vol. 82 (3), 2016a, 905820306). This distribution function can only be expressed as an infinite series of Hermite functions with very slow convergence and this makes its practical use cumbersome. It is the purpose of this paper to present a general method that allows us to find distribution functions consisting of a finite number of terms (therefore easier to use in practice), but which still allow for current sheet equilibria that can, in principle, have an arbitrarily low plasma beta. The method involves using known solutions and transforming them into new solutions using transformations based on taking integer powers (N) of one component of the pressure tensor. The plasma beta of the current sheet corresponding to the transformed distribution functions can then, in principle, have values as low as 1/N. We present the general form of the distribution functions for arbitrary N and then, as a specific example, discuss the case for N=2 in detail