Context. Some of the most prominent sources for particle acceleration in our Solar System are large eruptions of magnetised plasma from the Sun called coronal mass ejections (CMEs). These accelerated particles can generate radio emission through various mechanisms. Aims. CMEs are often accompanied by a variety of solar radio bursts with different shapes and characteristics in dynamic spectra. Radio bursts directly associated with CMEs often show movement in the direction of CME expansion. Here, we aim to determine the emission mechanism of multiple moving radio bursts that accompanied a flare and CME that took place on 14 June 2012. Methods. We used radio imaging from the Nancay Radioheliograph, combined with observations from the Solar Dynamics Observatory and Solar Terrestrial Relations Observatory spacecraft, to analyse these moving radio bursts in order to determine their emission mechanism and three-dimensional (3D) location with respect to the expanding CME. Results. In using a 3D representation of the particle acceleration locations in relation to the overlying coronal magnetic field and the CME propagation, for the first time, we provide evidence that these moving radio bursts originate near the CME flanks and some that are possible signatures of shock-accelerated electrons following the fast CME expansion in the low corona. Conclusions. The moving radio bursts, as well as other stationary bursts observed during the eruption, occur simultaneously with a type IV continuum in dynamic spectra, which is not usually associated with emission at the CME flanks. Our results show that moving radio bursts that could traditionally be classified as moving type IVs can represent shock signatures associated with CME flanks or plasma emission inside the CME behind its flanks, which are closely related to the lateral expansion of the CME in the low corona.

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.

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.

This White Paper outlines the importance of addressing the fundamental science theme <<How are charged particles energized in space plasmas>> through a future ESA mission. The White Paper presents five compelling science questions related to particle energization by shocks, reconnection,waves and turbulence, jets and their combinations. Answering these questions requires resolving scale coupling, nonlinearity and nonstationarity, which cannot be done with existing multi-point observations. In situ measurements from a multi-point, multi-scale L-class plasma observatory consisting of at least 7 spacecraft covering fluid, ion and electron scales are needed. The plasma observatory will enable a paradigm shift in our comprehension of particle energization and space plasma physics in general, with very important impact on solar and astrophysical plasmas. It will be the next logical step following Cluster, THEMIS and MMS for the very large and active European space plasmas community. Being one of the cornerstone missions of the future ESA Voyage 2035-2050 science program, it would further strengthen the European scientific and technical leadership in this important field.

Using neutron monitor and space-borne data we performed an analysis of the second ground level enhancement of solar cycle 24, namely the event of 10 September 2017 (GLE 72) and derive the spectral and angular characteristics of GLE particles. We employ new neutron monitor yield function and a recently proposed model based on optimization procedure. The method consists of simulation of particle propagation in a model magnetosphere in order to derive the cut-off rigidity and neutron monitor asymptotic directions. Subsequently the rigidity spectrum and anisotropy of GLE particles are obtained in their dynamical evolution during the event on the basis of inverse problem solution. The derived angular distribution and spectra are briefly discussed.

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.

The Solar Accumulated and Peak Proton and Heavy Ion Radiation Environment (SAPPHIRE) model provides environment specification outputs for all aspects of the Solar Energetic Particle (SEP) environment. The model is based upon a thoroughly cleaned and carefully processed data set. Herein the evolution of the solar proton model is discussed with comparisons to other models and data. This paper discusses the construction of the underlying data set, the modelling methodology, optimisation of fitted flux distributions and extrapolation of model outputs to cover a range of proton energies from 0.1 MeV to 1 GeV. The model provides outputs in terms of mission cumulative fluence, maximum event fluence and peak flux for both solar maximum and solar minimum periods. A new method for describing maximum event fluence and peak flux outputs in terms of 1-in-x-year SPEs is also described. SAPPHIRE proton model outputs are compared with previous models including CREME96, ESP-PSYCHIC and the JPL model. Low energy outputs are compared to SEP data from ACE/EPAM whilst high energy outputs are compared to a new model based on GLEs detected by Neutron Monitors (NMs).

The aim of the study is to investigate the reason for the low productivity of high-energy SEPs in the present solar cycle. We employ scaling laws derived from diffusive shock acceleration theory and simulation studies including proton-generated upstream Alfv\'en waves to find out how the changes observed in the long-term average properties of the erupting and ambient coronal and/or solar wind plasma would affect the ability of shocks to accelerate particles to the highest energies. Provided that self-generated turbulence dominates particle transport around coronal shocks, it is found that the most crucial factors controlling the diffusive shock acceleration process are the number density of seed particles and the plasma density of the ambient medium. Assuming that suprathermal populations provide a fraction of the particles injected to shock acceleration in the corona, we show that the lack of most energetic particle events as well as the lack of low charge-to-mass ratio ion species in the present cycle can be understood as a result of the reduction of average coronal plasma and suprathermal densities in the present cycle over the previous one.

We study the link between an expanding coronal shock and the energetic particles measured near Earth during the Ground Level Enhancement (GLE) of 17 May 2012. We developed a new technique based on multipoint imaging to triangulate the 3-D expansion of the shock forming in the corona. It uses images from three vantage points by mapping the outermost extent of the coronal region perturbed by the pressure front. We derive for the first time the 3-D velocity vector and the distribution of Mach numbers, $M_{FM}$, of the entire front as a function of time. Our approach uses magnetic field reconstructions of the coronal field, full magneto-hydrodynamic simulations and imaging inversion techniques. We find that the highest $M_{FM}$ values appear near the coronal neutral line within a few minutes of the Coronal Mass Ejection (CME) onset; this neutral line is usually associated with the source of the heliospheric current and plasma sheet. We illustrate the variability of the shock speed, shock geometry and Mach number along different modeled magnetic field lines. Despite the level of uncertainty in deriving the shock Mach numbers, all employed reconstruction techniques show that the release time of GeV particles occurs when the coronal shock becomes super-critical ($M_{FM}>3$). Combining in-situ measurements with heliospheric imagery, we also demonstrate that magnetic connectivity between the accelerator (the coronal shock of 17 May 2012) and the near-Earth environment is established via a magnetic cloud that erupted from the same active region roughly five days earlier.

Solar Orbiter is a joint ESA-NASA mission planed for launch in October 2018. The science payload includes remote-sensing and in-situ instrumentation designed with the primary goal of understanding how the Sun creates and controls the heliosphere. The spacecraft will follow an elliptical orbit around the Sun, with perihelion as close as 0.28 AU. During the late orbit phase the orbital plane will reach inclinations above 30 degrees, allowing direct observations of the solar polar regions. The Energetic Particle Detector (EPD) is an instrument suite consisting of several sensors measuring electrons, protons and ions over a broad energy interval (2 keV to 15 MeV for electrons, 3 keV to 100 MeV for protons and few tens of keV/nuc to 450 MeV/nuc for ions), providing composition, spectra, timing and anisotropy information. We present an overview of Solar Orbiter from the energetic particle perspective, summarizing the capabilities of EPD and the opportunities that these new observations will provide for understanding how energetic particles are accelerated during solar eruptions and how they propagate through the Heliosphere.

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.

Transport and acceleration of charged particles in turbulent media is a topic of great interest in space physics and interstellar astrophysics. These processes are dominated by the scattering of particles off magnetic irregularities. The scattering process itself is usually described by small-angle scattering with the pitch-angle coefficient $D_{\mu\mu}$ playing a major role. Since the diffusion coefficient $D_{\mu\mu}$ can be determined analytically only for the approximation of quasi-linear theory, the determination of this coefficient from numerical simulations has, therefore, become more important. So far these simulations yield particle tracks for small-scale scattering, which can then be interpreted using the running diffusion coefficients. This method has a limited range of validity. This paper presents two new methods that allow for the calculation of the pitch-angle diffusion coefficient from numerical simulations. These methods no longer analyse particle trajectories, but the change of particle distribution functions. It is shown that they provide better resolved results and allow for the analysis of strong turbulence. The application of these methods to Monte Carlo simulations of particle scattering and hybrid MHD-particle simulations is presented. Both analysis methods are able to recover the diffusion coefficients used as input for the Monte Carlo simulations and provide better results in MHD simulations especially for stronger turbulence.

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.

We study the differences in the heavy ion composition of solar energetic particle (SEP) events between solar cycles 23 and 24. We have surveyed the SOHO/ERNE heavy ion data from the beginning of solar cycle 23 until the end of June 2015, that is, well into the declining phase of cycle 24. We used this long observation period to study the properties of heavy ions (from C to Fe) and to compare the two solar cycles in this respect. We surveyed the data for SEP events with enhancements in the Fe/C and Fe/O intensity ratios in the energy range 5-15 MeV per nucleon, and associated the events with solar flare and coronal mass ejections (CME) when possible. We studied the properties of heavy ions in these events and compared the average relative abundances of heavy ions between the two solar cycles. We found that fewer days had C and O intensities higher than ~10$^{-3}$ cm$^{-2}$sr$^{-1}$s$^{-1}$(MeVn$^{-1}$)$^{-1}$ during solar cycle 24 than during cycle 23. For Fe this difference was clear even at lower intensities. We also found that fewer days had Fe/(C+O) > 0.183 during cycle 24. We identified 86 SEP events with at least one Fe-rich day, 65 of which occurred during cycle 23 and only 21 during cycle 24. We found that impulsive events have been almost completely absent during cycle 24. Mean abundances of heavy ions in the events were found to be significantly lower during cycle 24 than in cycle 23. Our results reflect the reduced solar activity in cycle 24 and indicate lower efficiency of particle acceleration processes for both gradual and impulsive SEP events in cycle 24.

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.

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.

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.

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.

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.

Three-wave interactions of plasma waves propagating parallel to the mean magnetic field at frequencies below the electron cyclotron frequency are considered. We consider Alfv\'en--ion-cyclotron waves, fast-magnetosonic--whistler waves, and ion-sound waves. Especially the weakly turbulent low-beta plasmas like the solar corona are studied, using the cold-plasma dispersion relation for the transverse waves and the fluid-description of the warm plasma for the longitudinal waves. We analyse the resonance conditions for the wave frequencies $\omega$ and wavenumbers $k$, and the interaction rates of the waves for all possible combinations of the three wave modes, and list those reactions that are not forbidden.

Aims: Solar radio type II bursts are rarely seen at frequencies higher than a few hundred MHz. Since metric type II bursts are thought to be signatures of propagating shock waves, it is of interest to know how these shocks, and the type II bursts, are formed. In particular, how are high-frequency, fragmented type II bursts created? Are there differences in shock acceleration or in the surrounding medium that could explain the differences to the "typical" metric type IIs? Methods: We analyse one unusual metric type II event in detail, with comparison to white-light, EUV, and X-ray observations. As the radio event was associated with a flare and a coronal mass ejection (CME), we investigate their connection. We then utilize numerical MHD simulations to study the shock structure induced by an erupting CME in a model corona including dense loops. Results: Our simulations show that the fragmented part of the type II burst can be formed when a coronal shock driven by a mass ejection passes through a system of dense loops overlying the active region.To produce fragmented emission, the conditions for plasma emission have to be more favourable inside the loop than in the interloop area. The obvious hypothesis, consistent with our simulation model, is that the shock strength decreases significantly in the space between the denser loops. The later, more typical type II burst appears when the shock exits the dense loop system and finally, outside the active region, the type II burst dies out when the changing geometry no longer favours the electron shock-acceleration.

A high-speed halo-type coronal mass ejection (CME), associated with a GOES M4.6 soft X-ray flare in NOAA AR 0180 at S12W29 and an EIT wave and dimming, occurred on 9 November 2002. A complex radio event was observed during the same period. It included narrow-band fluctuations and frequency-drifting features in the metric wavelength range, type III burst groups at metric--hectometric wavelengths, and an interplanetary type II radio burst, which was visible in the dynamic radio spectrum below 14 MHz. To study the association of the recorded solar energetic particle (SEP) populations with the propagating CME and flaring, we perform a multi-wavelength analysis using radio spectral and imaging observations combined with white-light, EUV, hard X-ray, and magnetogram data. Velocity dispersion analysis of the particle distributions (SOHO and Wind in situ observations) provides estimates for the release times of electrons and protons. Our analysis indicates that proton acceleration was delayed compared to the electrons. The dynamics of the interplanetary type II burst identify the burst source as a bow shock created by the fast CME. The type III burst groups, with start times close to the estimated electron release times, trace electron beams travelling along open field lines into the interplanetary space. The type III bursts seem to encounter a steep density gradient as they overtake the type II shock front, resulting in an abrupt change in the frequency drift rate of the type III burst emission. Our study presents evidence in support of a scenario in which electrons are accelerated low in the corona behind the CME shock front, while protons are accelerated later, possibly at the CME bow shock high in the corona.

We apply a semi-classical approach of handling waves as quasiparticle gas in a slowly varying flow -- analogous to ray tracing -- to calculate the Alfven wave transmission parameters, the resulting cross-helicity of the waves and the scattering-centre compression ratio, for cases where the shock thickness is large enough for the turbulent waves in the plasma to see the transition of the background flow parameters as smooth and slowly varying. For nonrelativistic shocks the wave transmission produces similar effects on the downstream turbulence and the scattering-centre compression ratio as does the transmission through a step shock: the downstream Alfven waves propagate predominantly towards the shock in the local plasma frame and, thus, the scattering-centre compression ratio is larger than the gas compression ratio. For thick relativistic shocks, however, we find qualitative differences with respect to the step-shock case: for low-Alfvenic-Mach-number shocks the downstream waves propagate predominantly away from the shock, and the scattering-centre compression ratio is lower than that of the gas. Thus, when taken into account, the Alfven wave transmission can decrease the efficiency of the first-order Fermi acceleration in a thick relativistic shock.

We report studies on first-order Fermi acceleration in parallel modified shock waves with a large scattering center compression ratio expected from turbulence transmission models. Using a Monte Carlo technique we have modeled particle acceleration in shocks with a velocity ranging from nonrelativistic to ultrarelativistic and a thickness extending from nearly steplike to very wide structures exceeding the particle diffusion length by orders of magnitude. The nonrelativistic diffusion approximation is found to be surprisingly accurate in predicting the spectral index of a thick shock with large compression ratio even in the cases involving relativistic shock speeds.

We present results of test-particle simulations on both the first- and the second-order Fermi acceleration for relativistic parallel shock waves. Our studies suggest that the role of the second-order mechanism in the turbulent downstream of a relativistic shock may have been underestimated in the past, and that the stochastic mechanism may have significant effects on the form of the particle spectra and its time evolution.

(abridged) We present results of test-particle simulations on both the first and the second order Fermi acceleration at relativistic parallel shock waves. We consider two scenarios for particle injection: (i) particles injected at the shock front, then accelerated at the shock by the first order mechanism and subsequently by the stochastic process in the downstream region; and (ii) particles injected uniformly throughout the downstream region to the stochastic process. We show that regardless of the injection scenario, depending on the magnetic field strength, plasma composition, and the employed turbulence model, the stochastic mechanism can have considerable effects on the particle spectrum on temporal and spatial scales too short to be resolved in extragalactic jets. Stochastic acceleration is shown to be able to produce spectra that are significantly flatter than the limiting case of particle energy spectral index -1 of the first order mechanism. Our study also reveals a possibility of re-acceleration of the stochastically accelerated spectrum at the shock, as particles at high energies become more and more mobile as their mean free path increases with energy. Our findings suggest that the role of the second order mechanism in the turbulent downstream of a relativistic shock with respect to the first order mechanism at the shock front has been underestimated in the past, and that the second order mechanism may have significant effects on the form of the particle spectra and its evolution.

It is generally accepted that relativistic jet outflows power the nonthermal emission from active galactic nuclei (AGN). The composition of these jets -- leptonic versus hadronic -- is still under debate. We investigate the microphysical details of the conversion process of the kinetic energy in collimated relativistic pair outflows into radiation through interactions with the ambient interstellar medium. Viewed from the coordinate system comoving with the pair outflow, the interstellar protons and electrons represent a proton-electron beam propagating with relativistic speed in the pair plasma. We demonstrate that the beam excites both electrostatic and low-frequency magnetohydrodynamic Alfven-type waves via a two-stream instability in the pair background plasma, and we calculate the time evolution of the distribution functions of the beam particles and the generated plasma wave turbulence power spectra. For standard AGN jet outflow and environment parameters we show that the initial beam distributions of interstellar protons and electrons quickly relax to plateau-distributions in parallel momentum, transferring thereby one-half of the initial energy density of the beam particles to electric field fluctuations of the generated electrostatic turbulence. On considerably longer time scales, the plateaued interstellar electrons and protons will isotropise by their self-generated transverse turbulence and thus be picked-up in the outflow pair plasma. These longer time scales are also characteristic for the development of transverse hydromagnetic turbulence from the plateaued electrons and protons. This hydromagnetic turbulence upstream and downstream is crucial for diffusive shock acceleration to operate at external or internal shocks associated with pair outflows.