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

Aims: Emission of radio waves from plasmas through plasma emission with fundamental and harmonic frequencies is a familiar process known from solar type II radio bursts. Current models assume the existence of counterstreaming electron beam populations excited at shocks as sources for these emission features, which limits the plasma parameters to reasonable heliospheric shock conditions. However, situations in which counterstreaming electron beams are present can also occur with different plasma parameters, such as higher magnetisation, including but not limited to our Sun. Similar radio emissions might also occur from these situations. Methods: We used particle-in-cell simulations, to compare plasma microphysics of radio emission processes from counterstreaming beams in different plasma environments that differed in density and magnetization. Results: Although large differences in wave populations are evident, the emission process of type II bursts appears to be qualitatively unaffected and shows the same behaviour in all environments.

Generation of radio bursts in CME foreshock regions and turbulent cascades in the solar wind are assumed to be results of three-wave interaction processes of dispersive plasma modes. Using our Particle in Cell code ACRONYM, we have studied the behaviour of kinetic wavemodes in the presence of beamed electron populations, with a focus on type II radio burst emission processes. We discuss the numerical challenges in generating and analyzing self-consistently evolving wave coupling processes with a PiC-Code and present preliminary results of said project.

The foreshock region of a CME shock front, where shock accelerated electrons form a beam population in the otherwise quiescent plasma is generally assumed to be the source region of type II radio bursts. Nonlinear wave interaction of electrostatic waves excited by the beamed electrons are the prime candidates for the radio waves' emission. To address the question whether a single, or two counterpropagating beam populations are a requirement for this process, we have conducted 2.5D particle in cell simulations using the fully relativistic ACRONYM code. Results show indications of three wave interaction leading to electromagnetic emission at the fundamental and harmonic frequency for the two-beam case. For the single-beam case, no such signatures were detectable.

The emission of fundamental and harmonic frequency radio waves of type II radio bursts are assumed to be products of three-wave interaction processes of beam-excited Langmuir waves. Using a particle-in-cell code, we have performed simulations of the assumed emission region, a CME foreshock with two counterstreaming electron beams. Analysis of wavemodes within the simulation shows self-consistent excitation of beam driven modes, which yield interaction products at both fundamental and harmonic emission frequencies. Through variation of the beam strength, we have investigated the dependence of energy transfer into electrostatic and electromagnetic modes, confirming the quadratic dependence of electromagnetic emission on electron beam strength.

Almost all sources of high energy particles and photons are associated with jet phenomena. Prominent sources of such highly relativistic outflows are pulsar winds and Active Galactic Nuclei. The current understanding of these jets assumes diluted plasmas which are best described as kinetic phenomena. In this kinetic description particle acceleration to ultra-relativistic speeds can occur in completely unmagnetized and neutral plasmas through insetting effects of instabilities. Even though the morphology and nature of particle spectra are understood to a certain extent, the composition of the jets is not known yet. While Poynting-flux dominated jets are certainly composed of electron-positron plasmas, the understanding of the governing physics in AGN jets is mostly unclear. In this article we investigate how the constituting elements of an electron-positron-proton plasma behave differently under the variation of the fundamental mass-ratio m_p/m_e. We studied initially unmagnetized counterstreaming plasmas using fully relativistic three-dimensional particle-in-cell simulations to investigate the influence of the mass-ratio on particle acceleration and magnetic field generation in electron-positron-proton plasmas. We covered a range of mass-ratios m_p/m_e between 1 and 100 with a particle number composition of n_{p^+}/n_{e^+} of 1 in one stream, only protons are injected in the other, whereas electrons are present in both to guarantee charge neutrality in the simulation box. We find that with increasing proton mass the instability takes longer to develop and for mass-ratios > 20 the particles seem to be accelerated in two phases which can be accounted to the individual instabilities of the different species. This means that for high mass ratios the coupling between electrons/positrons and the heavier protons, which occurs in low mass-ratios, disappears.