Context. Solar Energetic Particles (SEPs) with energy in the GeV range can propagate to Earth from their acceleration region near the Sun and produce Ground Level Enhancements (GLEs). The traditional approach to interpreting and modelling GLE observations assumes particle propagation only parallel to the magnetic field lines of interplanetary space, i.e. it is spatially 1D. Recent measurements by PAMELA have characterised SEP properties at 1 AU for the ~100 MeV-1 GeV range at high spectral resolution. Aims. We model the transport of GLE-energy solar protons through the Interplanetary Magnetic Field (IMF) using a 3D approach, to assess the effect of the Heliospheric Current Sheet (HCS) and drifts associated to the gradient and curvature of the Parker spiral. The latter are influenced by the IMF polarity. We derive 1 AU observables and compare the simulation results with data from PAMELA. Methods. We use a 3D test particle model including a HCS. Monoenergetic populations are studied first to obtain a qualitative picture of propagation patterns and numbers of crossings of the 1 AU sphere. Simulations for power law injection are used to derive intensity profiles and fluence spectra at 1 AU. A simulation for a specific event, GLE 71, is used to compare with PAMELA data. Results. Spatial patterns of 1 AU crossings and the average number of crossings are strongly influenced by 3D effects, with significant differences between periods of A+ and A- polarities. The decay time constant of 1 AU intensity profiles varies depending on the polarity and position of the observer, and it is not a simple function of the mean free path as in 1D models. Energy dependent leakage from the injection flux tube is particularly important for GLE energy particles, in many cases resulting in a roll-over in the fluence spectrum.

Solar energetic particles (SEPs) are an important product of solar activity. They are connected to solar active regions and flares, coronal mass ejections (CMEs), EUV waves, shocks, Type II and III radio emissions, and X-ray bursts. These phenomena are major probes of the partition of energy in solar eruptions, as well as for the organization, dynamics, and relaxation of coronal and interplanetary magnetic fields. Many of these phenomena cause terrestrial space weather, posing multiple hazards for humans and their technology from space to the ground. Since particular flares, shocks, CMEs, and EUV waves produce SEP events but others do not, since propagation effects from the low corona to 1 AU appear important for some events but not others, and since Type II and III radio emissions and X-ray bursts are sometimes produced by energetic particles leaving these acceleration sites, it is necessary to study the whole system with a multi-frequency and multi-instrument perspective that combines both in-situ and remote observations with detailed modelling of phenomena. This article demonstrates this comprehensive approach, and shows its necessity, by analysing a trio of unusual and striking solar eruptions, radio and X-ray bursts, and SEP events that occurred on 4 November 2015. These events show both strong similarities and differences from standard events and each other, despite having very similar interplanetary conditions and only two are sites and CME genesis regions. They are therefore major targets for further in-depth observational studies, and for testing both existing and new theories and models. Based on the very limited modelling available we identify the aspects that are and are not understood, and we discuss ideas that may lead to improved understanding of the SEP, radio, and space-weather events.

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.

The source of high-energy protons (>500 MeV) responsible for the so-called ground level enhancements (GLEs) remains an open question in solar physics. One of the candidates is a shock wave driven by a coronal mass ejection, which is thought to accelerate particles via diffusive-shock acceleration. We perform physics-based simulations of proton acceleration using information on the shock and ambient plasma parameters derived from the observation of a real GLE event. We analyse the simulation results with the aim to find out which of the parameters are significant in controlling the acceleration efficiency and to get a better understanding of the conditions under which the shock can produce relativistic protons. We use results of the recently developed technique to determine the shock and ambient plasma parameters, applied to the 17 May 2012 GLE event, and carry out proton acceleration simulations with the Coronal Shock Acceleration model. We have performed proton acceleration simulations for nine individual magnetic field lines characterised by various plasma conditions. Analysis of the simulation results shows that the acceleration efficiency of the shock, i.e., its ability to accelerate particles to high energies, tends to be higher for those shock portions that are characterised by larger values of the scattering-centre compression ratio and/or the fast-mode Mach number. At the same time, the acceleration efficiency can be strengthened due to enhanced plasma density in the flux tube. Analysis of the delays between the flare onset and the production times of protons of 1 GV rigidity for different field lines in our simulations, and a subsequent comparison of those with the observed values indicate a possibility that quasi-perpendicular portions of the shock play the main role in producing relativistic protons.

Forecasting Solar Energetic Particle (SEP) fluence, as integrated over an SEP event, is an important element when estimating the effect of solar eruptions on humans and technology in space. Current real-time estimates are based on SEP measurements at a single location in space. However, the interplanetary magnetic field corotates with the Sun approximately 13$^\circ$ each day with respect to Earth, thus in 4 days a near-Earth spacecraft will have changed their connection about 60$^\circ$ from the original SEP source. We estimate the effect of the corotation on particle fluence using a simple particle transport model, and show that ignoring corotation can cause up to an order of magnitude error in fluence estimations, depending on the interplanetary particle transport conditions. We compare the model predictions with STEREO observations of SEP events.

The injection, propagation and arrival of solar energetic particles (SEPs) during eruptive solar events is an important and current research topic of heliospheric physics. During the largest solar events, particles may have energies up to a few GeVs and sometimes even trigger ground-level enhancements (GLEs) at Earth. We study the first GLE-event of solar cycle 24, from 17th May 2012, using data from multiple spacecraft (SOHO, GOES, MSL, STEREO-A, STEREO-B and MESSENGER). These spacecraft are located throughout the inner heliosphere, at heliocentric distances between 0.34 and 1.5 astronomical units (au), covering nearly the whole range of heliospheric longitudes. We present and investigate sub-GeV proton time profiles for the event at several energy channels, obtained via different instruments aboard the above spacecraft. We investigate issues due to magnetic connectivity, and present results of three-dimensional SEP propagation simulations. We gather virtual time profiles and perform qualitative and quantitative comparisons with observations, assessing longitudinal injection and transport effects as well as peak intensities. We distinguish different time profile shapes for well-connected and weakly connected observers, and find our onset time analysis to agree with this distinction. At select observers, we identify an additional low-energy component of Energetic Storm Particles (ESPs). Using well-connected observers for normalisation, our simulations are able to accurately recreate both time profile shapes and peak intensities at multiple observer locations. This synergetic approach combining numerical modeling with multi-spacecraft observations is crucial for understanding the propagation of SEPs within the interplanetary magnetic field. Our novel analysis provides valuable proof of the ability to simulate SEP propagation throughout the inner heliosphere, at a wide range of longitudes.

Understanding the transport of Solar Energetic Particles (SEPs) from acceleration sites at the Sun into interplanetary space and to the Earth is an important question for forecasting space weather. The Interplanetary Magnetic Field (IMF), with two distinct polarities and a complex structure, governs energetic particle transport and drifts. We analyse for the first time the effect of a wavy Heliospheric Current Sheet (HCS) on the propagation of SEPs. We inject protons close to the Sun and propagate them by integrating fully 3D trajectories within the inner heliosphere in the presence of weak scattering. We model the HCS position using fits based on neutral lines of magnetic field source surface maps (SSMs). We map 1 au proton crossings, which show efficient transport in longitude via HCS, depending on the location of the injection region with respect to the HCS. For HCS tilt angles around $30^\circ-40^\circ$, we find significant qualitative differences between A+ and A$-$ configurations of the IMF, with stronger fluences along the HCS in the former case but with a distribution of particles across a wider range of longitudes and latitudes in the latter. We show how a wavy current sheet leads to longitudinally periodic enhancements in particle fluence. We show that for an A+ IMF configuration, a wavy HCS allows for more proton deceleration than a flat HCS. We find that A$-$ IMF configurations result in larger average fluences than A+ IMF configurations, due to a radial drift component at the current sheet.

Solar Energetic Particles (SEPs), a major component of space weather, propagate through the interplanetary medium strongly guided by the Interplanetary Magnetic Field (IMF). In this work, we analyse the implications a flat Heliospheric Current Sheet (HCS) has on proton propagation from SEP release sites to the Earth. We simulate proton propagation by integrating fully 3-D trajectories near an analytically defined flat current sheet, collecting comprehensive statistics into histograms, fluence maps and virtual observer time profiles within an energy range of 1--800 MeV. We show that protons experience significant current sheet drift to distant longitudes, causing time profiles to exhibit multiple components, which are a potential source of confusing interpretation of observations. We find that variation of current sheet thickness within a realistic parameter range has little effect on particle propagation. We show that IMF configuration strongly affects deceleration of protons. We show that in our model, the presence of a flat equatorial HCS in the inner heliosphere limits the crossing of protons into the opposite hemisphere.

The event-averaged charge state of heavy ion Solar Energetic Particles (SEPs), measured at 1 AU from the Sun, typically increases with the ions' kinetic energy. The origin of this behaviour has been ascribed to processes taking place within the acceleration region. In this paper we study the propagation through interplanetary space of SEP Fe ions, injected near the Sun with a variety of charge states that are uniformly distributed in energy, by means of a 3D test particle model. In our simulations, due to gradient and curvature drifts associated with the Parker spiral magnetic field, ions of different charge propagate with very different efficiencies to an observer that is not magnetically well connected to the source region. As a result we find that, for many observer locations, the 1 AU event-averaged charge state <Q>, as obtained from our model, displays an increase with particle energy E, in qualitative agreement with spacecraft observations. We conclude that drift-associated propagation is a possible explanation for the observed distribution of <Q> versus E in SEP events, and that the distribution measured in interplanetary space cannot be taken to represent that at injection.

Context. Solar energetic particles observed in association with coronal mass ejections (CMEs) are produced by the CME-driven shock waves. The acceleration of particles is considered to be due to diffusive shock acceleration (DSA). Aims. We aim at a better understanding of DSA in the case of quasi-parallel shocks, in which self-generated turbulence in the shock vicinity plays a key role. Methods. We have developed and applied a new Monte Carlo simulation code for acceleration of protons in parallel coronal shocks. The code performs a self-consistent calculation of resonant interactions of particles with Alfv\'en waves based on the quasi-linear theory. In contrast to the existing Monte Carlo codes of DSA, the new code features the full quasi-linear resonance condition of particle pitch-angle scattering. This allows us to take anisotropy of particle pitch-angle scattering into account, while the older codes implement an approximate resonance condition leading to isotropic scattering.We performed simulations with the new code and with an old code, applying the same initial and boundary conditions, and have compared the results provided by both codes with each other, and with the predictions of the steady-state theory. Results. We have found that anisotropic pitch-angle scattering leads to less efficient acceleration of particles than isotropic. However, extrapolations to particle injection rates higher than those we were able to use suggest the capability of DSA to produce relativistic particles. The particle and wave distributions in the foreshock as well as their time evolution, provided by our new simulation code, are significantly different from the previous results and from the steady-state theory. Specifically, the mean free path in the simulations with the new code is increasing with energy, in contrast to the theoretical result.

Context: Diffusive shock acceleration in the solar corona can accelerate solar energetic particles to very high energies. Acceleration efficiency is increased by entrapment through self-generated waves, which is highly dependent on the amount of accelerated particles. This, in turn, is determined by the efficiency of particle injection into the acceleration process. Aims: We present an analysis of the injection efficiency at coronal shocks of varying obliquity. We assessed injection through reflection and downstream scattering, including the effect of a cross-shock potential. Both quasi-thermal and suprathermal seed populations were analysed. We present results on the effect of cross-field diffusion downstream of the shock on the injection efficiency. Methods: Using analytical methods, we present applicable injection speed thresholds that were compared with both semi-analytical flux integration and Monte Carlo simulations, which do not resort to binary thresholds. Shock-normal angle and shock-normal velocity were varied to assess the injection efficiency with respect to these parameters. Results: We present evidence of a significant bias of thermal seed particle injection at small shock-normal angles. We show that downstream isotropisation methods affect the theta-dependence of this result. We show a non-negligible effect caused by the cross-shock potential, and that the effect of downstream cross-field diffusion is highly dependent on boundary definitions. Conclusions: Our results show that for Monte Carlo simulations of coronal shock acceleration a full distribution function assessment with downstream isotropisation through scatterings is necessary to realistically model particle injection. Based on our results, seed particle injection at quasi-parallel coronal shocks can result in significant acceleration efficiency, especially when combined with varying field-line geometry.

High-energy particles stream during coronal mass ejections or flares through the plasma of the solar wind. This causes instabilities, which lead to wave growth at specific resonant wave numbers, especially within shock regions. These amplified wave modes influence the turbulent scattering process significantly. In this paper, results of particle transport and scattering in turbulent plasmas with excited wave modes are presented. The method used is a hybrid simulation code, which treats the heliospheric turbulence by an incompressible magnetohydrodynamic approach separately from a kinetic particle description. Furthermore, a semi-analytical model using quasilinear theory (QLT) is compared to the numerical results. This paper aims at a more fundamental understanding and interpretation of the pitch-angle scattering coefficients. Our calculations show a good agreement of particle simulations and the QLT for broad-band turbulent spectra; for higher turbulence levels and particle beam driven plasmas, the QLT approximation gets worse. Especially the resonance gap at $\mu=0$ poses a well-known problem for QLT for steep turbulence spectra, whereas test-particle computations show no problems for the particles to scatter across this region. The reason is that the sharp resonant wave--particle interactions in QLT are an oversimplification of the broader resonances in test-particle calculations, which result from nonlinear effects not included in the QLT. We emphasise the importance of these results for both numerical simulations and analytical particle transport approaches, especially the validity of the QLT.

Context: Acceleration in coronal mass ejection driven shocks is currently considered the primary source of large solar energetic particle events. Aims: The solar wind, which feeds shock-accelerated particles, includes numerous ion populations, which offer much insight into acceleration processes. We present first simulations of shock-accelerated minor ions, in order to explore trapping dynamics and acceleration timescales in detail. Methods: We have simulated diffusive shock acceleration of minor ions (3He2+, 4He2+, 16O6+ and 56Fe14+) and protons using a Monte Carlo method, where self-generated Alfv\'enic turbulence allows for repeated shock crossings and acceleration to high energies. Results: We present the effect of minor ions on wave generation, especially at low wavenumbers, and show that it is significant. We find that maximum ion energy is determined by the competing effects of particle escape due to focusing in an expanding flux tube and trapping due to the amplified turbulence. We show the dependence of cut-off energy on the particle charge to mass ratio to be approximately (Q/A)^1.5. Conclusions: We suggest that understanding the acceleration of minor ions at coronal shocks requires simulations which allow us to explore trapping dynamics and acceleration timescales in detail, including evolution of the turbulent trapping boundary. We conclude that steady-state models do not adequately describe the acceleration of heavy ions in coronal shocks.

We have developed a simulation model of particle acceleration in coronal shock waves. The model is based on a Monte Carlo method, where particles are traced in prescribed large-scale electromagnetic fields utilizing the guiding center approximation. The particles are scattered in the turbulence according to quasilinear theory, with the scattering amplitude directly proportional to the intensity of Alfv\'en waves at gyro-resonant wavenumbers. The Alfv\'en waves are traced simultaneously with the particles, so that the wave field is propagated outwards from the Sun using WKB propagation supplemented with a phenomenological wavenumber diffusion term and a growth rate computed from the net flux of the accelerated particles. We consider initial wave amplitudes small enough to allow rapid escape of particles from the shock to the ambient medium. Thus, in our model the Alfv\'en waves responsible for the diffusive acceleration of particles are generated by the accelerated particles themselves. In this work, we study the effects of non-constant shock velocity and non-monotonic Alfv\'en velocity on particle acceleration scenarios. We report in particular how the deceleration of a shock affects particle intensity and turbulence power evolution in the vicinity of the shock.