Reconnection rates and X line motion at the magnetopause : Global 2D-3V hybrid-Vlasov simulation results

We present results from a first study of the local reconnection rate and reconnection site motion in a 2D-3V global magnetospheric self-consistent hybrid-Vlasov simulation with due southward interplanetary magnetic field. We observe magnetic reconnection at multiple locations at the dayside magnetopause and the existence of magnetic islands, which are the 2-D representations of flux transfer events. The reconnection locations (the X lines) propagate over significant distances along the magnetopause, and reconnection does not reach a steady state. We calculate the reconnection rate at the location of the X lines and find a good correlation with an analytical model of local 2-D asymmetric reconnection. We find that despite the solar wind conditions being constant, the reconnection rate and location of the X lines are highly variable. These variations are caused by magnetosheath fluctuations, the effects of neighboring X lines, and the motion of passing magnetic islands.Peer reviewe

On the Importance of Spatial and Velocity Resolution in the Hybrid-Vlasov Modeling of Collisionless Shocks

In hybrid-Vlasov plasma modeling, the ion velocity distribution function is propagated using the Vlasov equation while electrons are considered a charge-neutralizing fluid. It is an alternative to particle-in-cell methods, one advantage being the absence of sampling noise in the moments of the distribution. However, the discretization requirements in up to six dimensions (3D position, 3V velocity) make the computational cost of hybrid-Vlasov models higher. This is why hybrid-Vlasov modeling has only recently become more popular and available to model large-scale systems. The hybrid-Vlasov model Vlasiator is the first to have been successfully applied to model the solar-terrestrial interaction. It includes in particular the bow shock and magnetosheath regions, albeit in 2D-3V configurations so far. The purpose of this study is to investigate how Vlasiator parameters affect the modeling of a plasma shock in a 1D-3V simulation. The setup is similar to the Earth's bow shock in previous simulations, so that the present results can be related to existing and future magnetospheric simulations. The parameters investigated are the spatial and velocity resolution, as well as the phase space density threshold, which is the key parameter of the so-called sparse velocity space. The role of the Hall term in Ohm's law is also studied. The evaluation metrics used are the convergence of the final state, the complexity of spatial profiles and ion distributions as well as the position of the shock front. In agreement with previous Vlasiator studies it is not necessary to resolve the ion inertial length and gyroradius in order to obtain kinetic phenomena. While the code remains numerically stable with all combinations of resolutions, it is shown that significantly increasing the resolution in one space but not the other leads to unphysical results. Past a certain level, decreasing the phase space density threshold bears a large computational weight without clear physical improvement in the setup used here. Finally, the inclusion of the Hall term shows only minor effects in this study, mostly because of the 1D configuration and the scales studied, at which the Hall term is not expected to play a major role.In hybrid-Vlasov plasma modeling, the ion velocity distribution function is propagated using the Vlasov equation while electrons are considered a charge-neutralizing fluid. It is an alternative to particle-in-cell methods, one advantage being the absence of sampling noise in the moments of the distribution. However, the discretization requirements in up to six dimensions (3D position, 3V velocity) make the computational cost of hybrid-Vlasov models higher. This is why hybrid-Vlasov modeling has only recently become more popular and available to model large-scale systems. The hybrid-Vlasov model Vlasiator is the first to have been successfully applied to model the solar-terrestrial interaction. It includes in particular the bow shock and magnetosheath regions, albeit in 2D-3V configurations so far. The purpose of this study is to investigate how Vlasiator parameters affect the modeling of a plasma shock in a 1D-3V simulation. The setup is similar to the Earth's bow shock in previous simulations, so that the present results can be related to existing and future magnetospheric simulations. The parameters investigated are the spatial and velocity resolution, as well as the phase space density threshold, which is the key parameter of the so-called sparse velocity space. The role of the Hall term in Ohm's law is also studied. The evaluation metrics used are the convergence of the final state, the complexity of spatial profiles and ion distributions as well as the position of the shock front. In agreement with previous Vlasiator studies it is not necessary to resolve the ion inertial length and gyroradius in order to obtain kinetic phenomena. While the code remains numerically stable with all combinations of resolutions, it is shown that significantly increasing the resolution in one space but not the other leads to unphysical results. Past a certain level, decreasing the phase space density threshold bears a large computational weight without clear physical improvement in the setup used here. Finally, the inclusion of the Hall term shows only minor effects in this study, mostly because of the 1D configuration and the scales studied, at which the Hall term is not expected to play a major role.In hybrid-Vlasov plasma modeling, the ion velocity distribution function is propagated using the Vlasov equation while electrons are considered a charge-neutralizing fluid. It is an alternative to particle-in-cell methods, one advantage being the absence of sampling noise in the moments of the distribution. However, the discretization requirements in up to six dimensions (3D position, 3V velocity) make the computational cost of hybrid-Vlasov models higher. This is why hybrid-Vlasov modeling has only recently become more popular and available to model large-scale systems. The hybrid-Vlasov model Vlasiator is the first to have been successfully applied to model the solar-terrestrial interaction. It includes in particular the bow shock and magnetosheath regions, albeit in 2D-3V configurations so far. The purpose of this study is to investigate how Vlasiator parameters affect the modeling of a plasma shock in a 1D-3V simulation. The setup is similar to the Earth's bow shock in previous simulations, so that the present results can be related to existing and future magnetospheric simulations. The parameters investigated are the spatial and velocity resolution, as well as the phase space density threshold, which is the key parameter of the so-called sparse velocity space. The role of the Hall term in Ohm's law is also studied. The evaluation metrics used are the convergence of the final state, the complexity of spatial profiles and ion distributions as well as the position of the shock front. In agreement with previous Vlasiator studies it is not necessary to resolve the ion inertial length and gyroradius in order to obtain kinetic phenomena. While the code remains numerically stable with all combinations of resolutions, it is shown that significantly increasing the resolution in one space but not the other leads to unphysical results. Past a certain level, decreasing the phase space density threshold bears a large computational weight without clear physical improvement in the setup used here. Finally, the inclusion of the Hall term shows only minor effects in this study, mostly because of the 1D configuration and the scales studied, at which the Hall term is not expected to play a major role.Peer reviewe

On the Importance of Spatial and Velocity Resolution in the Hybrid-Vlasov Modeling of Collisionless Shocks

In hybrid-Vlasov plasma modeling, the ion velocity distribution function is propagated using the Vlasov equation while electrons are considered a charge-neutralizing fluid. It is an alternative to particle-in-cell methods, one advantage being the absence of sampling noise in the moments of the distribution. However, the discretization requirements in up to six dimensions (3D position, 3V velocity) make the computational cost of hybrid-Vlasov models higher. This is why hybrid-Vlasov modeling has only recently become more popular and available to model large-scale systems. The hybrid-Vlasov model Vlasiator is the first to have been successfully applied to model the solar-terrestrial interaction. It includes in particular the bow shock and magnetosheath regions, albeit in 2D-3V configurations so far. The purpose of this study is to investigate how Vlasiator parameters affect the modeling of a plasma shock in a 1D-3V simulation. The setup is similar to the Earth's bow shock in previous simulations, so that the present results can be related to existing and future magnetospheric simulations. The parameters investigated are the spatial and velocity resolution, as well as the phase space density threshold, which is the key parameter of the so-called sparse velocity space. The role of the Hall term in Ohm's law is also studied. The evaluation metrics used are the convergence of the final state, the complexity of spatial profiles and ion distributions as well as the position of the shock front. In agreement with previous Vlasiator studies it is not necessary to resolve the ion inertial length and gyroradius in order to obtain kinetic phenomena. While the code remains numerically stable with all combinations of resolutions, it is shown that significantly increasing the resolution in one space but not the other leads to unphysical results. Past a certain level, decreasing the phase space density threshold bears a large computational weight without clear physical improvement in the setup used here. Finally, the inclusion of the Hall term shows only minor effects in this study, mostly because of the 1D configuration and the scales studied, at which the Hall term is not expected to play a major role

Evidence for transient, local ion foreshocks caused by dayside magnetopause reconnection

We present a scenario resulting in time-dependent behaviour of the bow shock and transient, local ion reflection under unchanging solar wind conditions. Dayside magnetopause reconnection produces flux transfer events driving fast-mode wave fronts in the magnetosheath. These fronts push out the bow shock surface due to their increased downstream pressure. The resulting bow shock deformations lead to a configuration favourable to localized ion reflection and thus the formation of transient, travelling foreshock-like field-aligned ion beams. This is identified in two-dimensional global magnetospheric hybrid-Vlasov simulations of the Earth's magnetosphere performed using the Vlasiator model (http://vlasiator.fmi.fi). We also present observational data showing the occurrence of dayside reconnection and flux transfer events at the same time as Geotail observations of transient foreshock-like field-aligned ion beams. The spacecraft is located well upstream of the fore-shock edge and the bow shock, during a steady southward interplanetary magnetic field and in the absence of any solar wind or interplanetary magnetic field perturbations. This indicates the formation of such localized ion foreshocks.Peer reviewe

Tail reconnection in the global magnetospheric context: Vlasiator first results

The key dynamics of the magnetotail have been researched for decades and have been associated with either three-dimensional (3-D) plasma instabilities and/or magnetic reconnection. We apply a global hybrid-Vlasov code, Vlasiator, to simulate reconnection self-consistently in the ion kinetic scales in the noon-midnight meridional plane, including both dayside and nightside reconnection regions within the same simulation box. Our simulation represents a numerical experiment, which turns off the 3-D instabilities but models ion-scale reconnection physically accurately in 2-D. We demonstrate that many known tail dynamics are present in the simulation without a full description of 3-D instabilities or without the detailed description of the electrons. While multiple reconnection sites can coexist in the plasma sheet, one reconnection point can start a global reconfiguration process, in which magnetic field lines become detached and a plasmoid is released. As the simulation run features temporally steady solar wind input, this global reconfiguration is not associated with sudden changes in the solar wind. Further, we show that lobe density variations originating from dayside reconnection may play an important role in stabilising tail reconnection

Vlasov methods in space physics and astrophysics

This paper reviews Vlasov-based numerical methods used to model plasma in space physics and astrophysics. Plasma consists of collectively behaving charged particles that form the major part of baryonic matter in the Universe. Many concepts ranging from our own planetary environment to the Solar system and beyond can be understood in terms of kinetic plasma physics, represented by the Vlasov equation. We introduce the physical basis for the Vlasov system, and then outline the associated numerical methods that are typically used. A particular application of the Vlasov system is Vlasia- tor, the world’s first global hybrid-Vlasov simulation for the Earth’s magnetic domain, the magnetosphere. We introduce the design strategies for Vlasiator and outline its numerical concepts ranging from solvers to coupling schemes. We review Vlasiator’s parallelisation methods and introduce the used high-performance computing (HPC) techniques. A short review of verification, validation and physical results is included. The purpose of the paper is to present the Vlasov system and introduce an example implementation, and to illustrate that even with massive computational challenges, an accurate description of physics can be rewarding in itself and significantly advance our understanding. Upcoming supercomputing resources are making similar efforts feasible in other fields as well, making our design options relevant for others facing similar challenges.This paper reviews Vlasov-based numerical methods used to model plasma in space physics and astrophysics. Plasma consists of collectively behaving charged particles that form the major part of baryonic matter in the Universe. Many concepts ranging from our own planetary environment to the Solar system and beyond can be understood in terms of kinetic plasma physics, represented by the Vlasov equation. We introduce the physical basis for the Vlasov system, and then outline the associated numerical methods that are typically used. Aparticular application of the Vlasov system is Vlasiator, the world’s first global hybrid-Vlasov simulation for the Earth’s magnetic domain, the magnetosphere. We introduce the design strategies for Vlasiator and outline its numerical concepts ranging from solvers to coupling schemes. We review Vlasiator’s parallelisation methods and introduce the used high-performance computing (HPC) techniques. A short review of verification, validation and physical results is included. The purpose of the paper is to present the Vlasov system and introduce an example implementation, and to illustrate that even with massive computational challenges, an accurate description of physics can be rewarding in itself and significantly advance our understanding. Upcoming supercomputing resources are making similar efforts feasible in other fields as well, making our design options relevant for others facing similar challenges.Peer reviewe