Source author record

T. Neukirch

T. Neukirch appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2020arXiv

A family of Vlasov-Maxwell equilibrium distribution functions describing a transition from the Harris sheet to the force-free Harris sheet

We discuss a family of Vlasov-Maxwell equilibrium distribution functions for current sheet equilibria that are intermediate cases between the Harris sheet and the force-free (or modified) Harris sheet. These equilibrium distribution functions have potential applications to space and astrophysical plasmas. The existence of these distribution function had been briefly discussed in by Harrison and Neukirch 2009, but here it is shown that their approach runs into problems in the limit where the guide field goes to zero. The nature of this problem will be discussed and an alternative approach will be suggested that avoids the problem. This is achieved by considering a slight variation of the magnetic field profile, which allows a smooth transition between the Harris and force-free Harris sheet cases.

preprint2020arXiv

Kinetic models of tangential discontinuities in the solar wind

Kinetic-scale current sheets observed in the solar wind are frequently approximately force-free despite the fact that their plasma $β$ is of the order of one. In-situ measurements have recently shown that plasma density and temperature often vary across the current sheets, while the plasma pressure is approximately uniform. In many cases these density and temperature variations are asymmetric with respect to the center of the current sheet. To model these observations theoretically we develop in this paper equilibria of kinetic-scale force-free current sheets that have plasma density and temperature gradients. The models can also be useful for analysis of stability and dissipation of the current sheets in the solar wind.

preprint2016arXiv

From one-dimensional fields to Vlasov equilibria: Theory and application of Hermite polynomials

We consider the theory and application of a solution method for the inverse problem in collisionless equilibria, namely that of calculating a Vlasov-Maxwell equilibrium for a given macroscopic (fluid) equilibrium. Using Jeans' Theorem, the equilibrium distribution functions are expressed as functions of the constants of motion, in the form of a Maxwellian multiplied by an unknown function of the canonical momenta. In this case it is possible to reduce the inverse problem to inverting Weierstrass transforms, which we achieve by using expansions over Hermite polynomials. A sufficient condition on the pressure tensor is found which guarantees the convergence and the boundedness of the candidate solution, when satisfied. This condition is obtained by elementary means, and it is clear how to put it into practice. We also argue that for a given pressure tensor for which our method applies, there always exists a positive distribution function solution for a sufficiently magnetised plasma. Illustrative examples of the use of this method with both force-free and non-force-free macroscopic equilibria are presented, including the full verification of a recently derived distribution function for the force-free Harris Sheet (Allanson $\textit{et al., Phys. Plasmas}$, vol. 22 (10), 2015, 102116). In the effort to model equilibria with lower values of the plasma beta, solutions for the same macroscopic equilibrium in a new gauge are calculated, with numerical results presented for $β_{pl}=0.05$.

preprint2016arXiv

Neutral and non-neutral collisionless plasma equilibria for twisted flux tubes: The Gold-Hoyle model in a background field

We calculate exact one-dimensional collisionless plasma equilibria for a continuum of flux tube models, for which the total magnetic field is made up of the `force-free' Gold-Hoyle magnetic flux tube embedded in a uniform and anti-parallel background magnetic field. For a sufficiently weak background magnetic field, the axial component of the total magnetic field reverses at some finite radius. The presence of the background magnetic field means that the total system is not exactly force-free, but by reducing its magnitude the departure from force-free can be made as small as desired. The distribution function for each species is a function of the three constants of motion; namely the Hamiltonian and the canonical momenta in the axial and azimuthal directions. Poisson's Equation and Ampère's Law are solved exactly, and the solution allows either electrically neutral or non-neutral configurations, depending on the values of the bulk ion and electron flows. These equilibria have possible applications in various solar, space and astrophysical contexts, as well as in the laboratory.

preprint2016arXiv

Particle dynamics in a non-flaring solar active region model

The aim of this work is to investigate and characterise particle behaviour in an (observationally-driven) 3D MHD model of the solar atmosphere above a slowly evolving, non-flaring active region. We use a relativistic guiding-centre particle code to investigate the behaviour of selected particle orbits, distributed throughout a single snapshot of the 3D MHD simulation. Two distinct particle acceleration behaviours are recovered, which affect both electrons and protons: (i) direct acceleration along field lines and (ii) tangential drifting of guiding centres with respect to local magnetic field. However, up to 40\% of all particles actually experience a form of (high energy) particle trap, because of changes in the direction of the electric field and unrelated to the strength of the magnetic field; such particles are included in the first category. Additionally, category (i) electron and proton orbits undergo surprisingly strong acceleration to non-thermal energies ($\lesssim42$MeV), because of the strength and extent of super-Dreicer electric fields created by the MHD simulation. Reducing the electric field strength of the MHD model does not significantly affect the efficiency of the (electric field-based) trapping mechanism, but does reduce the peak energies gained by orbits. We discuss the implications for future experiments, which aim to simulate non-flaring active region heating and reconnection.

preprint2016arXiv

Particular solutions of the inverse problem for 1D Vlasov-Maxwell equilibria using Hermite polynomials

We present the solution to an inverse problem arising in the context of finding a distribution function for a specific collisionless plasma equilibrium. The inverse problem involves the solution of two integral equations, each having the form of a Weierstrass transform. We prove that inverting the Weierstrass transform using Hermite polynomials leads to convergent infinite series. We also comment on the non-negativity of the distribution function, with more detail on this in Allanson $\textit{et al., Journal of Plasma Physics}$, vol. 82 (03), 2016. Whilst applied to a specific magnetic field, the inversion techniques used in this paper (as well as the derived convergence criteria and discussion of non-negativity) are of a general nature, and are applicable to other smooth pressure functions.

preprint2015arXiv

An exact collisionless equilibrium for the Force-Free Harris Sheet with low plasma beta

We present a first discussion and analysis of the physical properties of a new exact collisionless equilibrium for a one-dimensional nonlinear force-free magnetic field, namely the Force-Free Harris Sheet. The solution allows any value of the plasma beta, and crucially below unity, which previous nonlinear force-free collisionless equilibria could not. The distribution function involves infinite series of Hermite Polynomials in the canonical momenta, of which the important mathematical properties of convergence and non-negativity have recently been proven. Plots of the distribution function are presented for the plasma beta modestly below unity, and we compare the shape of the distribution function in two of the velocity directions to a Maxwellian distribution.

preprint2015arXiv

Magneto-static modelling of the mixed plasma Beta solar atmosphere based on SUNRISE/IMaX data

Our aim is to model the 3D magnetic field structure of the upper solar atmosphere, including regions of non-negligible plasma beta. We use high-resolution photospheric magnetic field measurements from SUNRISE/IMaX as boundary condition for a magneto-static magnetic field model. The high resolution of IMaX allows us to resolve the interface region between photosphere and corona, but modelling this region is challenging for the following reasons. While the coronal magnetic field is thought to be force-free (the Lorentz-force vanishes), this is not the case in the mixed plasma $β$ environment in the photosphere and lower chromosphere. In our model, pressure gradients and gravity forces are taken self-consistently into account and compensate the non-vanishing Lorentz-force. Above a certain height (about 2 Mm) the non-magnetic forces become very weak and consequently the magnetic field becomes almost force-free. Here we apply a linear approach, where the electric current density consists of a superposition of a field-line parallel current and a current perpendicular to the Sun's gravity field. We illustrate the prospects and limitations of this approach and give an outlook for an extension towards a non-linear model.

preprint2014arXiv

Particle acceleration at a reconnecting magnetic separator

While the exact acceleration mechanism of energetic particles during solar flares is (as yet) unknown, magnetic reconnection plays a key role both in the release of stored magnetic energy of the solar corona and the magnetic restructuring during a flare. Recent work has shown that special field lines, called separators, are common sites of reconnection in 3D numerical experiments. To date, 3D separator reconnection sites have received little attention as particle accelerators. We investigate the effectiveness of separator reconnection as a particle acceleration mechanism for electrons and protons. We study the particle acceleration using a relativistic guiding-centre particle code in a time-dependent kinematic model of magnetic reconnection at a separator. The effect upon particle behaviour of initial position, pitch angle and initial kinetic energy are examined in detail, both for specific (single) particle examples and for large distributions of initial conditions. The separator reconnection model contains several free parameters and we study the effect of changing these parameters upon particle acceleration, in particular in view of the final particle energy ranges which agree with observed energy spectra.

preprint2012arXiv

A Systematic Examination of Particle Motion in a Collapsing Magnetic Trap Model for Solar Flares

Context. It has been suggested that collapsing magnetic traps may contribute to accelerating particles to high energies during solar flares. Aims. We present a detailed investigation of the energization processes of particles in collapsing magnetic traps, using a specific model. We also compare for the first time the energization processes in a symmetric and an asymmetric trap model. Methods. Particle orbits are calculated using guiding centre theory. We systematically investigate the dependence of the energization process on initial position, initial energy and initial pitch angle. Results. We find that in our symmetric trap model particles can gain up to about 50 times their initial energy, but that for most initial conditions the energy gain is more moderate. Particles with an initial position in the weak field region of the collapsing trap and with pitch angles around 90 degrees achieve the highest energy gain, with betatron acceleration of the perpendicular energy the dominant energization mechanism. For particles with smaller initial pitch angle, but still outside the loss cone, we find the possibility of a significant increase in parallel energy. This increase in parallel energy can be attributed to the curvature term in the parallel equation of motion and the associated energy gain happens in the center of the trap where the field line curvature has its maximum. We find qualitatively similar results for the asymmetric trap model, but with smaller energy gains and a larger number of particles escaping from the trap.

preprint2008arXiv

Computing nonlinear force free coronal magnetic fields

Knowledge 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.

preprint2008arXiv

Helmet Streamers with Triple Structures: Simulations of resistive dynamics

Recent observations of the solar corona with the LASCO coronagraph on board of the SOHO spacecraft have revealed the occurrence of triple helmet streamers even during solar minimum, which occasionally go unstable and give rise to large coronal mass ejections. There are also indications that the slow solar wind is either a combination of a quasi-stationary flow and a highly fluctuating component or may even be caused completely by many small eruptions or instabilities. As a first step we recently presented an analytical method to calculate simple two-dimensional stationary models of triple helmet streamer configurations. In the present contribution we use the equations of time- dependent resistive magnetohydrodynamics to investigate the stability and the dynamical behaviour of these configurations. We particularly focus on the possible differences between the dynamics of single isolated streamers and triple streamers and on the way in which magnetic reconnection initiates both small scale and large scale dynamical behaviour of the streamers. Our results indicate that small eruptions at the helmet streamer cusp may incessantly accelerate small amounts of plasma without significant changes of the equilibrium configuration and might thus contribute to the non-stationary slow solar wind. On larger time and length scales, large coronal eruptions can occur as a consequence of large scale magnetic reconnection events inside the streamer configuration. Our results also show that triple streamers are usually more stable than a single streamer.

preprint2008arXiv

Including stereoscopic information in the reconstruction of coronal magnetic fields

We present a method to include stereoscopic information about the three dimensional structure of flux tubes into the reconstruction of the coronal magnetic field. Due to the low plasma beta in the corona we can assume a force free magnetic field, with the current density parallel to the magnetic field lines. Here we use linear force free fields for simplicity. The method uses the line of sight magnetic field on the photosphere as observational input. The value of $α$ is determined iteratively by comparing the reconstructed magnetic field with the observed structures. The final configuration is the optimal linear force solution constrained by both the photospheric magnetogram and the observed plasma structures. As an example we apply our method to SOHO MDI/EIT data of an active region. In the future it is planned to apply the method to analyse data from the SECCHI instrument aboard the STEREO mission.

preprint1997arXiv

Helmet streamers with triple structure: Weakly two-dimensional stationary states

Recent observations of the solar corona with the LASCO coronagraph on board of the SOHO spacecraft have revealed the occurrence of triple helmet streamers even during solar minimum, which occasionally go unstable and give rise to particularly huge coronal mass ejections. We present a method to calculate (semi-)analytically self-consistent stationary configurations of triple helmet streamers which can serve as input for stability considerations and dynamical calculations. The method is based on an asymptotic expansion procedure using the elongated structure of the streamers. The method is very flexible and can be used in both Cartesian and spherical geometry. We discuss the effects of magnetic shear, gravity and field-aligned flow on open field lines. Example solutions illustrating the influence of each of these features on the solution are presented.