Particle acceleration at multiple internal relativistic shocks

Relativistic shocks provide an efficient method for high-energy particle acceleration in many astrophysical sources. Multiple shock systems are even more effective and of importance, for example, in the internal shock model of gamma-ray bursts. We investigate the reacceleration of pre-existing energetic particles at such relativistic internal shocks by the first order Fermi process of pitch angle scattering. We use a well established eigenfunction method to calculate the resulting spectra for infinitely thin shocks. Implications for GRBs and relativistic jets are discussed. Paul Dempsey would like to thank IRCSET for their financial support.     

Langmuir waves associated with collisionless shocks; a review

Progress in understanding the Langmuir waves which accompany collisionless shocks everywhere in the solar system is briefly reviewed, with some emphasis on the discovery papers, and with discussion and illustrative examples of the most recent progress.     

Particle spectra from acceleration at forward and reverse shocks of young Type Ia Supernova Remnants

We study cosmic-ray acceleration in young Type Ia Supernova Remnants (SNRs) by means of test-particle diffusive shock acceleration theory and 1-D hydrodynamical simulations of their evolution. In addition to acceleration at the forward shock, we explore the particle acceleration at the reverse shock in the presence of a possible substantial magnetic field, and consequently the impact of this acceleration on the particle spectra in the remnant. We investigate the time evolution of the spectra for various time-dependent profiles of the magnetic field in the shocked region of the remnant. We test a possible influence on particle spectra of the Alfvénic drift of scattering centers in the precursor regions of the shocks. In addition, we study the radiation spectra and morphology in a broad band from radio to gamma-rays. It is demonstrated that the reverse shock contribution to the cosmic-ray particle population of young Type Ia SNRs may be significant, modifying the spatial distribution of particles and noticeably affecting the volume-integrated particle spectra in young SNRs. In particular spectral structures may arise in test-particle calculations that are often discussed as signatures of non-linear cosmic-ray modification of shocks. Therefore, the spectrum and morphology of emission, and their time evolution, differ from pure forward-shock solutions.     

Trajectories of ions specularly reflected from non-planar collisionless shocks

An important contribution to the thermalization of the solar wind ions at the Earth's bow shock for high Mach numbers comes from the reflection of a fraction of these ions from the shock. Previous studies have examined the trajectories of the reflected ions assuming the shock to be an infinite plane. In this paper a model is developed to describe the trajectories of particles after reflection for a variety of shock geometries. Of particular interest are the initial conditions which allow the particle to return to the shock with a greater normal velocity than at first encounter, or to return to the shock at all. The effects of the magnetic field direction and the curvature of the shock on particle trajectories are discussed for cylindrical and spherical shock geometries and compared to those for a planar shock.     

Particle acceleration in shocks

The mechanism of particle acceleration by shocks is reviewed, and the resulting modifications of the shock structure contrasted with those which occur in molecular shocks with magnetic precursors. The danger associated with the use of steady models is emphasised and the applicability of the theory to cosmic ray acceleration in supernova remnants briefly discussed.     

Fast magnetic reconnection and energetic particle acceleration

Our numerical simulations show that the reconnection of magnetic field becomes fast in the presence of weak turbulence in the way consistent with the Lazarian and Vishniac (1999) model of fast reconnection. We trace particles within our numerical simulations and show that the particles can be efficiently accelerated via the first order Fermi acceleration. We discuss the acceleration arising from reconnection as a possible origin of the anomalous cosmic rays measured by Voyagers.     

Particle acceleration by ultrarelativistic shocks: theory and simulations

We consider the acceleration of charged particles near ultrarelativistic shocks, with Lorentz factor     . We present simulations of the acceleration process and compare these with results from semi-analytical calculations. We show that the spectrum that results from acceleration near ultrarelativistic shocks is a power law,     , with a nearly universal value     for the slope of this power law.
We confirm that the ultrarelativistic equivalent of the Fermi acceleration at a shock differs from its non-relativistic counterpart by the occurrence of large anisotropies in the distribution of the accelerated particles near the shock. In the rest frame of the upstream fluid, particles can only outrun the shock when their direction of motion lies within a small loss cone of opening angle     around the shock normal.
We also show that all physically plausible deflection or scattering mechanisms can change the upstream flight direction of relativistic particles originating from downstream by only a small amount:     . This limits the energy change per shock crossing cycle to     , except for the first cycle where particles originate upstream. In that case the upstream energy is boosted by a factor     for those particles that are scattered back across the shock into the upstream region.     

Diffusion in momentum space as a picture of second-order Fermi acceleration

Energetic particles in a turbulent medium can be subject to second-order Fermi acceleration due to scattering on moving plasma waves. This mechanism leads to growing particle momentum dispersion and, at the same time, increases the mean particle energy. In the most frequently met situations both processes can be represented by a single momentum diffusion term in the particle kinetic equation. In the present paper we discuss the conditions allowing the additional term for regular acceleration to arise. For forward-backward asymmetric scattering centres, besides the diffusive term one should explicitly consider the regular acceleration term in momentum space, which can consist of the first-order (∝ V), as well as the second-order (∝ V²) part in the wave velocity V. We derive the condition for the scattering probability in the wave rest frame requied for vanishing the regular acceleration term and provide a simple mechanical example illustrating the theoretical concepts. Finally, we address its possible role in cosmic ray acceleration processes.     

Stochastic Fermi acceleration in turbulent fields with non-vanishing wave helicities

In the previous paper [Astropart. Phys. 10 (1999) 121] we showed that the opposite helicity circularly polarized Alfvén waves of finite amplitudes provide conditions to forward–backward asymmetry of particle scattering. In this letter we present an analytic solution of kinematic equation proving the enhancement of stochastic acceleration efficiency due to regular (asymmetry) term. The process is controlled by the ratio of the regular and the ordinary diffusion term.     

The role of suprathermal particle measurements in CrossScale studies of collisionless plasma processes

The CrossScale mission will advance our understanding of fundamental plasma processes in collisionless plasmas. It will exploit the excellent natural plasma laboratory provided by the Earth’s magnetosphere and the near-Earth solar wind and, in particular, carry out multi-scale studies that will strongly complement plasma studies in ground-based laboratories. Previous studies of collisionless plasmas in space environments across the solar system have shown the ubiquitous nature of suprathermal particles and that these particles exhibit a power-law energy spectrum. In this paper we discuss the great significance of these suprathermal particles for CrossScale studies. We show that the presence of these particles is a natural consequence of the collisionless regime as they can propagate across the heliosphere with little spectral change and are not thermalised by collisions. They are a key indicator of the non-equilibrium nature of collisionless plasmas and an important source of free energy that can drive plasma processes. We discuss how these suprathermal particles influence the overall properties of the plasma. In particular, the energy distribution of particles follows a Kappa, rather than Maxwellian, distribution and thus the plasma does not have a single thermodynamic temperature. We also discuss the importance of the suprathermal tail as a tool to diagnose the processes responsible for particle energisation in collisionless plasmas. Such energisation is a common feature in collisionless plasmas, especially in terms of the primary science targets for CrossScale: reconnection, shocks and turbulence. Finally we also touch on the value of using CrossScale studies to provide ground truth measurements for a number of astrophysical techniques that exploit the effects of energetic electrons in the distant universe. Throughout the paper, we stress that suprathermal (30 keV-1 MeV) measurements are essential to fully characterise particle distributions. We show that such measurements will benefit greatly from the improved spatial and temporal resolution (compared to Cluster) that is proposed for the HEP instrument on CrossScale.     

On the role of injection in kinetic approaches to non-linear particle acceleration at non-relativistic shock waves

The dynamical reaction of the particles accelerated at a shock front by the first-order Fermi process can be determined within kinetic models that account for both the hydrodynamics of the shocked fluid and the transport of the accelerated particles. These models predict the appearance of multiple solutions, all physically allowed. We discuss here the role of injection in selecting the real solution, in the framework of a simple phenomenological recipe, which is a variation of what is sometimes referred to as thermal leakage. In this context we show that multiple solutions basically disappear and when they are present they are limited to rather peculiar values of the parameters. We also provide a quantitative calculation of the efficiency of particle acceleration at cosmic ray modified shocks and we identify the fraction of energy which is advected downstream and that of particles escaping the system from upstream infinity at the maximum momentum. The consequences of efficient particle acceleration for shock heating are also discussed.     

Contribution of nuclei accelerated by gamma-ray pulsars to cosmic rays in the Galaxy

We consider the galactic population of gamma-ray pulsars as possible sources of cosmic rays at and just above the “knee” in the observed cosmic ray spectrum at 10¹⁵–10¹⁶ eV. We suggest that iron nuclei may be accelerated in the outer gaps of pulsars, and then suffer partial photo-disintegration in the non-thermal radiation fields of the outer gaps. As a result, protons, neutrons, and surviving heavier nuclei are injected into the expanding supernova remnant. We compute the spectra of nuclei escaping from supernova remnants into the interstellar medium, taking into account the observed population of radio pulsars.

Our calculations, which include a realistic model for acceleration and propagation of nuclei in pulsar magnetospheres and supernova remnants, predict that heavy nuclei accelerated directly by gamma-ray pulsars could contribute about 20% of the observed cosmic rays in the knee region. Such a contribution of heavy nuclei to the cosmic ray spectrum at the knee can significantly increase the average value of lnA with increasing energy as is suggested by recent observations.     

Current constraints on the epoch of cosmic acceleration

The cosmographic expansion history of the universe is investigated by using the 557 type Ia supernovae from the Union2 Compilation set along with the current estimates involving the product of the CMB acoustic scale ?_A and the BAO peak at two different redshifts. Using a well-behaved parameterization for the deceleration parameter, q(z) = q₀ + q₁z/(1 + z), we estimate the accelerating redshift z_acc = −q₀/(q₀ + q₁) (at which the universe switches from deceleration to acceleration) and investigate the influence of a non-vanishing spatial curvature on these estimates. We also use the asymptotic value of q(z) at high-z to place more restrictive bounds on the model parameters q₀ and q₁, which results in a more precise determination of the epoch of cosmic acceleration.     

Electron surfing acceleration by electrostatic waves in current sheet

The electron surfing acceleration in the current sheet with perpendicular propagating electrostatic waves is studied using analytical theories and test particle simulations. The trapped electron moving with the phase velocity v _p of wave may be accelerated effectively in the outflow direction by force until the electron is de-trapped from the wave potential. A criterion K>0 for the electron surfing acceleration is obtained. The electron will escape from the boundary of current sheet quickly, if this criterion does not hold. The maximum velocity of surfing acceleration is about the same as the electric drift velocity. Superposed longitudinal magnetic field along the wave propagation is favorable for the electron surfing acceleration in the current sheet.      

High-energy particle acceleration and production of ultra-high-energy cosmic rays in the giant lobes of Centaurus A

The nearby radio galaxy Centaurus A is poorly studied at high frequencies with conventional radio telescopes because of its very large angular size, but is one of a very few extragalactic objects to be detected and resolved by the Wilkinson Microwave Anisotropy Probe ( WMAP ). We have used the five-year WMAP data for Cen A to constrain the high-frequency radio spectra of the 10° giant lobes and to search for spectral changes as a function of position along the lobes. We show that the high-frequency radio spectra of the northern and southern giant lobes are significantly different: the spectrum of the southern lobe steepens monotonically (and is steeper further from the active nucleus) whereas the spectrum of the northern lobe remains consistent with a power law. The inferred differences in the northern and southern giant lobes may be the result of real differences in their high-energy particle acceleration histories, perhaps due to the influence of the northern middle lobe, an intermediate-scale feature which has no detectable southern counterpart. In light of these results, we discuss the prospects for Fermi Gamma-ray Space Telescope detections of inverse-Compton emission from the giant lobes and the lobes' possible role in the production of the ultra-high-energy cosmic rays (UHECR) detected by the Pierre Auger Observatory. We show that the possibility of a Fermi detection depends sensitively on the physical conditions in the giant lobes, with the northern lobe more likely to be detected, and that any emission observed by Fermi is likely to be dominated by photons at the soft end of the Fermi energy band. On the other hand, we argue that the estimated conditions in the giant lobes imply that UHECRs can be accelerated there, with a potentially detectable γ-ray signature at TeV energies.     

Hard electron energy distribution in the relativistic shocks of gamma-ray burst afterglows

Particle acceleration in relativistic shocks is not a very well understood subject. Owing to that difficulty, radiation spectra from relativistic shocks, such as those in gamma-ray burst (GRB) afterglows, have been often modelled by making assumptions about the underlying electron distribution. One such assumption is a relatively soft distribution of the particle energy, which need not be true always, as is obvious from observations of several GRB afterglows. In this paper, we describe modifications to the afterglow standard model to accommodate energy spectra which are 'hard'. We calculate the overall evolution of the synchrotron and Compton flux arising from such a distribution. We also model two afterglows, GRB010222 and GRB020813, under this assumption and estimate the physical parameters.