Source author record

Rami Vainio

Rami Vainio 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

17works
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

17 published item(s)

preprint2026arXiv

Turbulence-Driven Corrugation of Collisionless Fast-Magnetosonic Shocks

Collisionless fast-magnetosonic shocks are often treated as smooth, planar boundaries, yet observations point to organized corrugation of the shock surface. A plausible driver is upstream turbulence. Broadband fluctuations arriving at the front can continually wrinkle it, changing the local shock geometry and, in turn, conditions for particle injection and radiation. We develop a linear-MHD formulation that treats the shock as a moving interface rather than a fixed boundary. In this approach the shock response can be summarized by an effective impedance determined by the Rankine-Hugoniot base state and the shock geometry, while the upstream turbulence enters only through its statistics. This provides a practical mapping from an assumed incident spectrum to the corrugation amplitude, its drift along the surface, and a coherence scale set by weak damping or leakage. The response is largest when the transmitted downstream fast mode propagates nearly parallel to the shock in the shock frame, which produces a Lorentzian-type enhancement controlled by the downstream normal group speed. We examine how compression, plasma $β$, and obliquity affect these corrugation properties and discuss implications for fine structure in heliospheric and supernova-remnant shock emission.

preprint2022arXiv

CMEs and SEPs During November-December 2020: A Challenge for Real-Time Space Weather Forecasting

Predictions of coronal mass ejections (CMEs) and solar energetic particles (SEPs) are a central issue in space weather forecasting. In recent years, interest in space weather predictions has expanded to include impacts at other planets beyond Earth as well as spacecraft scattered throughout the heliosphere. In this sense, the scope of space weather science now encompasses the whole heliospheric system, and multi-point measurements of solar transients can provide useful insights and validations for prediction models. In this work, we aim to analyse the whole inner heliospheric context between two eruptive flares that took place in late 2020, i.e. the M4.4 flare of November 29 and the C7.4 flare of December 7. This period is especially interesting because the STEREO-A spacecraft was located ~60° east of the Sun-Earth line, giving us the opportunity to test the capabilities of "predictions at 360°" using remote-sensing observations from the Lagrange L1 and L5 points as input. We simulate the CMEs that were ejected during our period of interest and the SEPs accelerated by their shocks using the WSA-Enlil-SEPMOD modelling chain and four sets of input parameters, forming a "mini-ensemble". We validate our results using in-situ observations at six locations, including Earth and Mars. We find that, despite some limitations arising from the models' architecture and assumptions, CMEs and shock-accelerated SEPs can be reasonably studied and forecast in real time at least out to several tens of degrees away from the eruption site using the prediction tools employed here.

preprint2022arXiv

Observation-based modelling of the energetic storm particle event of 14 July 2012

We model the energetic storm particle (ESP) event of 14 July 2012 using the energetic particle acceleration and transport model named PARADISE, together with the solar wind and coronal mass ejection (CME) model named EUHFORIA. The simulation results illustrate both the capabilities and limitations of the utilised models. We show that the models capture some essential structural features of the ESP event; however, for some aspects the simulations and observations diverge. We describe and, to some extent, assess the sources of errors in the modelling chain of EUHFORIA and PARADISE and discuss how they may be mitigated in the future. The PARADISE model evolves energetic particle distributions in a background solar wind generated by the ideal MHD module of EUHFORIA. The CME generating the ESP event is simulated by using the spheromak model of EUHFORIA, which approximates the CME's flux rope as a linear force-free spheroidal magnetic field. In addition, a tool was developed to trace CME-driven shock waves in the EUHFORIA simulation domain. This tool is used in PARADISE to (i) inject 50 keV protons continuously at the CME-driven shock and (ii) include a foreshock and a sheath region, in which the energetic particle parallel mean free path, $λ_\parallel$, decreases towards the shock wave. The value of $λ_\parallel$ at the shock wave is estimated from in situ observations of the ESP event. For energies below 1 MeV, the simulation results agree well with both the upstream and downstream components of the ESP event observed by the Advanced Composition Explorer (ACE). This suggests that these low-energy protons are mainly the result of interplanetary particle acceleration. In the downstream region, the sharp drop in the energetic particle intensities is reproduced at the entry into the following magnetic cloud, illustrating the importance of a magnetised CME model.

preprint2022arXiv

The Probabilistic Solar Particle Event foRecasting (PROSPER) Model

The Probabilistic Solar Particle Event foRecasting (PROSPER) model predicts the probability of occurrence and the expected peak flux of Solar Energetic Particle (SEP) events. Predictions are derived for a set of integral proton energies (i.e. E$>$10, $>$30 and $>$100 MeV) from characteristics of solar flares (longitude, magnitude), coronal mass ejections (width, speed) and combinations of both. Herein the PROSPER model methodology for deriving the SEP event forecasts is described and the validation of the model, based on archived data, is presented for a set of case studies. The PROSPER model has been incorporated into the new operational Advanced Solar Particle Event Casting System (ASPECS) tool to provide nowcasting (short term forecasting) of SEP events as part of ESA's future SEP Advanced Warning System (SAWS). ASPECS also provides the capability to interrogate PROSPER for historical cases via a run on demand functionality.

preprint2020arXiv

Electron acceleration by magnetosheath jet-driven bow waves

Magnetosheath jets are localized fast flows with enhanced dynamic pressure. When they supermagnetosonically compress the ambient magnetosheath plasma, a bow wave or shock can form ahead of them. Such a bow wave was recently observed to accelerate ions and possibly electrons. The ion acceleration process was previously analyzed, but the electron acceleration process remains largely unexplored. Here we use multi-point observations by Time History of Events and Macroscale during Substorms from three events to determine whether and how magnetosheath jet-driven bow waves can accelerate electrons. We show that when suprathermal electrons in the ambient magnetosheath convect towards a bow wave, some electrons are shock-drift accelerated and reflected towards the ambient magnetosheath and others continue moving downstream of the bow wave resulting in bi-directional motion. Our study indicates that magnetosheath jet-driven bow waves can result in additional energization of suprathermal electrons in the magnetosheath. It implies that magnetosheath jets can increase the efficiency of electron acceleration at planetary bow shocks or other similar astrophysical environments.

preprint2020arXiv

Statistical study of magnetosheath jet-driven bow waves

When a magnetosheath jet (localized dynamic pressure enhancements) compresses ambient magnetosheath at a (relative) speed faster than the local magnetosonic speed, a bow wave or shock can form ahead of the jet. Such bow waves or shocks were recently observed to accelerate particles, thus contributing to magnetosheath heating and particle acceleration in the extended environment of Earth bow shock. To further understand the characteristics of jet-driven bow waves, we perform a statistical study to examine which solar wind conditions favor their formation and whether it is common for them to accelerate particles. We identified 364 out of 2859 (13%) magnetosheath jets to have a bow wave or shock ahead of them with Mach number typically larger than 1.1. We show that large solar wind plasma beta, weak interplanetary magnetic field (IMF) strength, large solar wind Alfven Mach number, and strong solar wind dynamic pressure present favorable conditions for their formation. We also show that magnetosheath jets with bow waves or shocks are more frequently associated with higher maximum ion and electron energies than those without them, confirming that it is common for these structures to accelerate particles. In particular, magnetosheath jets with bow waves have electron energy flux enhanced on average by a factor of 2 compared to both those without bow waves and the ambient magnetosheath. Our study implies that magnetosheath jets can contribute to shock acceleration of particles especially for high Mach number shocks. Therefore, shock models should be generalized to include magnetosheath jets and concomitant particle acceleration.

preprint2019arXiv

Calibration of RADMON Radiation Monitor Onboard Aalto-1 CubeSat

RADMON is a small radiation monitor designed and assembled by students of the University of Turku and the University of Helsinki. It is flown on-board Aalto-1, a 3-unit CubeSat in low Earth orbit at about 500 km altitude. The detector unit of the instrument consists of two detectors, a Si solid-state detector and a CsI(Tl) scintillator, and utilizes the \textDelta{E}-E technique to determine the total energy and species of each particle hitting the detector. We present the results of the on-ground and in-flight calibration campaigns of the instrument, as well as the characterization of its response through extensive simulations within the Geant4 framework. The overall energy calibration margin achieved is about 5\%. The full instrument response to protons and electrons is presented and the issue of proton contamination of the electron channels is quantified and discussed.

preprint2019arXiv

Particle Telescope aboard FORESAIL-1: simulated performance

The Particle Telescope (PATE) of FORESAIL-1 mission is described. FORESAIL-1 is a CubeSat mission to polar Low Earth Orbit. Its scientific objectives are to characterize electron precipitation from the radiation belts and to observe energetic neutral atoms (ENAs) originating from the Sun during the strongest solar flares. For that purpose, the 3-unit CubeSat carries a particle telescope that measures energetic electrons in the nominal energy range of 80--800 keV in seven energy channels and energetic protons at 0.3--10 MeV in ten channels. In addition, particles penetrating the whole telescope at higher energies will be measured in three channels: one $>$800 keV electron channel, two integral proton channels at $>$10 MeV energies. The instrument contains two telescopes at right angles to each other, one measuring along the spin axis of the spacecraft and one perpendicular to it. During a spin period (nominally 15 s), the rotating telescope will, thus, deliver angular distributions of protons and electrons, at 11.25-degree clock-angle resolution, which enables one to accurately determine the pitch-angle distribution and separate the trapped and precipitating particles. During the last part of the mission, the rotation axis will be accurately pointed toward the Sun, enabling the measurement of the energetic hydrogen from that direction. Using the geomagnetic field as a filter and comparing the rates observed by the two telescopes, the instrument can observe the solar ENA flux for events similar to the only one so far observed in December 2006. We present the Geant4-simulated energy and angular response functions of the telescope and assess its sensitivity showing that they are adequate to address the scientific objectives of the mission.

preprint2016arXiv

Assessment of Source and Transport Parameters of Relativistic SEPs Based on Neutron Monitor Data

As part of the HESPERIA Horizon 2020 project, we developed a software package for the direct inversion of Ground Level Enhancements (GLEs) based on data of the worldwide network of Neutron Monitors (NMs). The new methodology to study the release processes of relativistic solar energetic particles (SEPs) makes use of several models, including: the propagation of relativistic SEPs from the Sun to the Earth, their transport in the Earth's magnetosphere and atmosphere, as well as the detection of the nucleon component of the secondary cosmic rays by ground based NMs. The combination of these models allows to compute the expected ground-level NM counting rates caused by a series of instantaneous particle releases from the Sun. The proton release-time profile at the Sun and the interplanetary transport conditions are then inferred by fitting NM observations with modeled NM counting rates. In the paper the used models for the different processes, the software and first findings with the new software are presented.

preprint2016arXiv

Self-consistent Monte Carlo simulations of proton acceleration in coronal shocks: Effect of anisotropic pitch-angle scattering of particles

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

preprint2014arXiv

Fundamental and harmonic plasma emission in different plasma environments

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.

preprint2013arXiv

Acceleration of Energetic Particles Through Self-Generated Waves in a Decelerating Coronal Shock

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én waves at gyro-resonant wavenumbers. The Alfvén 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én 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én 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.

preprint2013arXiv

Heavy-ion Acceleration and Self-generated Waves in Coronal Shocks

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

preprint2013arXiv

Injection of thermal and suprathermal seed particles into coronal shocks of varying obliquity

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.

preprint2013arXiv

Particle scattering in turbulent plasmas with amplified wave modes

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 $μ=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.

preprint2012arXiv

Emission of Type II Radio Bursts - Single-Beam versus Two-Beam Scenario

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.

preprint2012arXiv

Nonlinear Wave Interactions as Emission Process of Type II Radio Bursts

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.