Evolution of magnetic fields and cosmic ray acceleration in supernova remnants

Observations show that the magnetic field in young supernova remnants (SNRs) is significantly stronger than can be expected from the compression of the circumstellar medium (CSM) by a factor of four expected for strong blast waves. Additionally, the polarization is mainly radial, which is also contrary to expectation from compression of the CSM magnetic field. Cosmic rays (CRs) may help to explain these two observed features. They can increase the compression ratio to factors well over those of regular strong shocks by adding a relativistic plasma component to the pressure, and by draining the shock of energy when CRs escape from the region. The higher compression ratio will also allow for the contact discontinuity, which is subject to the Rayleigh–Taylor (R–T) instability, to reach much further out to the forward shock. This could create a preferred radial polarization of the magnetic field. With an Adaptive Mesh Refinement MHD code (AMRVAC), we simulate the evolution of SNRs with three different configurations of the initial CSM magnetic field, and look at two different equations of state in order to look at the possible influence of a CR plasma component. The spectrum of CRs can be simulated using test particles, of which we also show some preliminary results that agree well with available analytical solutions.     

Gamma-ray astronomy and cosmic ray origin problem

We analyse the results of the nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova remnants (SNRs) in order to describe their relevant properties: the remnant dynamics and the characteristics of nonthermal emission produced by CRs. It is shown that the theory fits the existing data in a satisfactory way and that the magnetic field in SNRs is significantly amplified due to efficient acceleration of the nuclear CR component. From the fact that magnetic field amplification occurs in all the young SNRs for which relevant data exist, and given the strong theoretical connection between magnetic field amplification and efficiently accelerating the nuclear CRs, we tentatively conclude that the Galactic SNRs are the source population of the Galactic CRs. Due to high interior magnetic fields in young SNRs the π⁰-decay γ-rays generated by the nuclear CR component as a rule dominate over γ-rays generated by the electron CR component, and the calculated γ-ray flux fits existing data.     

Propagation of cylindrical and spherical dust–ion acoustic solitary waves in a relativistic dusty plasma

The properties of cylindrical and spherical dust–ion acoustic solitary waves (DIASW) in an unmagnetized dusty plasma comprising of relativistic ions, Boltzmann electrons, and stationary dusty particles are investigated. Under a suitable coordinate transformation, the cylindrical KdV equation can be solved analytically. The change of the DIASW structure due to the effect of geometry, relativistic streaming factor, ion density and electron temperature is studied by numerical calculation of the cylindrical/spherical Kdv equation. It is noted that with ion pressure the effect of relativistic streaming factor to solitary waves structure is different. Without ion pressure, as the relativistic streaming factor decreases, the amplitude of the solitary wave decreases. However, when the ion pressure is taken into account, the amplitude decreases as the relativistic streaming factor increases and is highly sensitive to relativistic streaming factor. Our results may have relevance in the understanding of astrophysical plasmas.     

Alfvén waves propagation in homogeneous and dusty astrophysical plasmas

Astrophysical plasma coexist with dust particles in many situations. These particles are charged either negatively or positively depending on their surrounding plasma environments. This system of such charged dust, electrons, and ions forms a so-called dusty plasma. We discuss the effects of the dust particles on the propagation and absorption of the Alfvén waves in (i) stellar winds and (ii) in star formation regions. In both cases, we have shown the importance of a strong damping of Alfvén waves due to the dust and the consequences for wind acceleration and the changes in the Jeans length related to the star formation process.     

Injection mechanism for cosmic ray acceleration

The quasilinear theory of MHD waves excitation by cosmic rays accelerated at a front of supernova shock has been constructed. It is shown that the energetic particles excite the waves propagating from the shock front, the intensity and the spectrum of these waves is obtained. The role of nonlinear Landau damping in the formation of such spectrum has been analysed. The diffusive scattering length of the high energy particles in the preshock region has been calculated and it is shown with the help of these formulae that the effective Fermi acceleration at the shock front is possible upto the values of the relativistic factor = 10⁴ - 10⁵. The injection mechanism for cosmic rays acceleration has been proposed. It is based on stochastic Fermi acceleration of the thermal plasma by MHD waves excited in the preshock region. Different possibilities for wave phase velocity dispersion needed for stochastic Fermi acceleration are analysed, those are the excitation of the oblique magnetosonic waves as well as the excitation of parallel Alfven waves propagating in opposite directions. The distribution function of the suprathermal particles accelerated by MHD waves is obtained, the cosmic rays density as well as the lower boundary of their energy spectrum realised in the proposed mechanism are also calculated.     

The acceleration of Anomalous Cosmic Rays by stochastic acceleration in the heliosheath

Stochastic acceleration in the heliosheath appears to be a likely mechanism by which Anomalous Cosmic Rays (ACRs) are accelerated. However, most stochastic acceleration mechanisms are not appropriate. The energy density in the ACRs and in the interstellar pickup ions out of which the ACRs are accelerated greatly exceeds the energy density in the turbulence in the heliosheath. Thus, a traditional stochastic acceleration mechanism in which particles are accelerated by damping the turbulence will not work. A stochastic acceleration mechanism has been developed in which the total energy of the pickup ions and the ACRs is conserved. Energy is redistributed from the core pickup ions into a suprathermal tail to create the ACRs. A model for the acceleration of the ACRs in the heliosheath, based on this stochastic acceleration mechanism, is presented. The model provides reasonable fits to the spectra of suprathermal particles and ACRs observed by Voyager.     

Depth dependency of neutron density produced by cosmic rays in the lunar subsurface

Depth dependency of neutrons produced by cosmic rays (CRs) in the lunar subsurface was estimated using the three-dimensional Monte Carlo particle and heavy ion transport simulation code, PHITS, incorporating the latest high energy nuclear data, JENDL/HE-2007. The PHITS simulations of equilibrium neutron density profiles in the lunar subsurface were compared with the measurement by Apollo 17 Lunar Neutron Probe Experiment (LNPE). Our calculations reproduced the LNPE data except for the 350–400 mg/cm² region under the improved condition using the CR spectra model based on the latest observations, well-tested nuclear interaction models with systematic cross section data, and JENDL/HE-2007.     

Stochastic particle acceleration by the forced interaction of relativistic current sheets

We address the problem of interacting relativistic current sheets in self-consistent kinetic plasma simulations within the framework of the Particle-In-Cell model. The interaction is enforced in head-on collisions of up to 10 current sheets at relativistic bulk speeds. The simulations are motivated by the general problem of Poynting flux dissipation in ‘striped wind’ configurations presumably governing the relativistic outflows pervasive in pulsar winds and gamma-ray bursts. We identify the generation of non-thermal particles and formation of a stable power-law shape in the particle energy distributions f(γ) dγ ∝ γ^−s dγ. In 1D, a spectral index s ∼ 2 is observed and attributed to a stochastic Fermi-type acceleration mechanism. In 2D, the generic index of s ∼ 3–4 is retained as in previous simulations of individual current sheets. Whereas in 2D the high energy cut-off is constrained by the limited dissipation of magnetic energy, in 1D the process converts the bulk motion of current sheets towards directed particle momentum of an exclusive class of non-thermal particles.     

Acceleration of interplanetary pick-up ions and anomalous cosmic rays

It may not be doubted anymore that anomalous cosmic rays (ACRs) are produced in the heliosphere from interplanetary pick-up ions through their acceleration at the solar wind termination shock. However, there is no general agreement in the community of heliospheric researchers concerning the mechanism of injection of the pick-up ions into the shock acceleration. We discuss here three possible ways for pick-up ions to be involved into the acceleration process at the termination shock: (1) preacceleration of pick-up ions in the whole region from the Sun up to the termination shock by solar wind turbulences and interplanetary shock waves, (2) local preacceleration of pick-up ions in a vicinity of the termination shock by shock surfing, and (3) formation of high-velocity tails in pick-up ion spectra consisting of secondary pick-up ions which are produced in the supersonic solar wind due to ionization of energetic neutral atoms entering from the inner heliosheath.     

Occurrence features of simultaneous H+- and He+-band EMIC emissions in the outer radiation belt

As an important loss mechanism of radiation belt electrons, electromagnetic ion cyclotron (EMIC) waves show up as three distinct frequency bands below the hydrogen (H⁺), helium (He⁺), and oxygen (O⁺) ion gyrofrequencies. Compared to O⁺-band EMIC waves, H⁺- and He⁺-band emissions generally occur more frequently and result in more efficient scattering removal of <～5?MeV relativistic electrons. Therefore, knowledge about the occurrence of these two bands is important for understanding the evolution of the relativistic electron population. To evaluate the occurrence pattern and wave properties of H⁺- and He⁺-band EMIC waves when they occur concurrently, we investigate 64 events of multi-band EMIC emissions identified from high quality Van Allen Probes wave data. Our quantitative results demonstrate a strong occurrence dependence of the multi-band EMIC emissions on magnetic local time (MLT) and L-shell to mainly concentrate on the dayside region of L?=?～4–6. We also find that the average magnetic field amplitude of H⁺-band waves is larger than that of He⁺-band waves only when L?<?4.5 and AE¹?<?300?nT, and He⁺-band emissions are more intense under all other conditions. In contrast to 5 events that have average H⁺-band amplitude over 2 nT, 19 events exhibit >2 nT He⁺-band amplitude, indicating that the He⁺-band waves can be more easily amplified than the H⁺-band waves under the same circumstances. For simultaneous occurrences of the two EMIC wave bands, their frequencies vary with L-shell and geomagnetic activity: the peak wave frequency of H⁺-band emissions varies between 0.25 and 0.8 f_cp with the average between 0.25 and 0.6 f_cp, while that of He⁺-band emissions varies between 0.03 and 0.23 f_cp with the average between 0.05 and 0.15 f_cp. These newly observed occurrence features of simultaneous H⁺- and He⁺-band EMIC emissions provide improved information to quantify the overall contribution of multi-band EMIC waves to the loss processes of radiation belt electrons.     

Energization of charged particle in a time-dependent chaotic magnetic field with an implication of the production of seed particles in solar energetic particle events

We investigate the acceleration of charged particles in a time-dependent chaotic magnetic field in this work. In earlier works, it has been demonstrated that in an asymmetric wire-loop current systems (WLCSs), the magnetic field is of chaotic in nature. Furthermore, observations also showed that there exist time-varying current loops and current filaments in solar corona. It is therefore natural to conceive that the magnetic field on the solar surface is chaotic and time-dependent. Here, we develop a numerical model to study the acceleration process of charged particles in a time-varying chaotic magnetic field that is generated by an ensemble of 8 WLCSs. We found that the motion of energetic particles in the system is of diffusive in nature and a power law spectrum can quickly develop. The mechanism examined here may serve as an efficient pre-acceleration mechanism that generates the so-called seed particles for diffusive shock acceleration at a coronal mass ejection (CME) driven shock in large solar energetic particle (SEP) events.     

The ground-based large-area wide-angle γ-ray and cosmic-ray experiment HiSCORE

The question of the origin of cosmic rays and other questions of astroparticle and particle physics can be addressed with indirect air-shower observations above 10 TeV primary energy. We propose to explore the cosmic ray and γ-ray sky (accelerator sky) in the energy range from 10 TeV to 1 EeV with the new ground-based large-area wide angle (ΔΩ ∼ 0.85 sterad) air-shower detector HiSCORE (Hundred^∗i Square-km Cosmic ORigin Explorer). The HiSCORE detector is based on non-imaging air-shower Cherenkov light-front sampling using an array of light-collecting stations. A full detector simulation and basic reconstruction algorithms have been used to assess the performance of HiSCORE. First prototype studies for different hardware components of the detector array have been carried out. The resulting sensitivity of HiSCORE to γ-rays will be comparable to CTA at 50 TeV and will extend the sensitive energy range for γ-rays up to the PeV regime. HiSCORE will also be sensitive to charged cosmic rays between 100 TeV and 1 EeV.     

The beginning of cosmic ray astrophysics

A short history of the beginning of cosmic ray (CR) astrophysics is considered: from the hypothesis on CR origin as a result of Supernova explosions in the Metagalaxy, to a model of solar origin of CR, galactic origin based on the stochastic mechanism of charged particle acceleration in interstellar space, to extragalactic and hierarchical models of CR origin, as well as galactic CR origin taking into account radio-astronomical data. We consider also the first balloon results on the chemical contents of primary CR (especially of the contents Li, Be, B), important for any model of CR origin. Investigations of the injection problem, CR drift and diffusion acceleration by shock waves, and CR generation in Supernova remnants were also important steps in the beginning of CR astrophysics.     

Fluctuations of cosmic rays and IMF in the vicinity of interplanetary shocks

Fluctuations of cosmic rays and interplanetary magnetic field upstream of interplanetary shocks are studied using data of ground-based polar neutron monitors as well as measurements of energetic particles and solar wind plasma parameters aboard the ACE spacecraft. It is shown that coherent cosmic ray fluctuations in the energy range from 10 keV to 1 GeV are often observed at the Earth’s orbit before the arrival of interplanetary shocks. This corresponds to an increase of solar wind turbulence level by more than the order of magnitude upstream of the shock. We suggest a scenario where the cosmic ray fluctuation spectrum is modulated by fast magnetosonic waves generated by flux of low-energy cosmic rays which are reflected and/or accelerated by an interplanetary shock.     

Formation of relativistic jets by collapsing stars to black holes

Formation of relativistic jets in the magnetosphere of collapsing stars is considered. These jets will be formed in the polar caps of magnetosphere of collapsing star, where the stellar magnetic field increases during the collapse and the charged particles are accelerated. The jets will generate non-thermal radiation. The analysis of dynamics and emission of particles in the stellar magnetosphere under collapse shows that collapsing stars can be powerful sources of relativistic jets.     

The role of radial perpendicular diffusion and latitude dependent acceleration along the solar wind termination shock

A numerical model, based on Parker’s transport equation, describing the modulation of anomalous cosmic rays and containing diffusive shock acceleration is applied. The role of radial perpendicular diffusion at the solar wind termination shock, and as the dominant diffusion coefficient in the outer heliosphere, is studied, in particular the role it plays in the effectiveness of the acceleration of anomalous protons and helium when its latitude dependence is changed. It is found that the latitudinal enhancement of radial perpendicular diffusion towards the heliospheric poles and along the termination shock has a prominent effect on the acceleration of these particles. It results in a ‘break’ in the energy spectrum for anomalous protons at ∼6.0 MeV, causing the spectral index to change from E^−1.38 to E^−2.23, but for anomalous helium at ∼3.0 MeV, changing the spectral index from E^−1.38 to E^−2.30. When approaching the simulated TS, the changes in the modulated spectra as they unfold to a ‘steady’ power law shape at energies below 50 MeV are much less prominent as a function of radial distances when radial perpendicular diffusion is increased with heliolatitude.     

Acceleration and transport of energeticparticles at CME-driven shocks

Historically, solar energetic particle (SEP) events are classified in two classes as “impulsive” and “gradual”. Whether there is a clear distinction between the two classes is still a matter of debate, but it is now commonly accepted that in large “gradual” SEP events, Fermi acceleration, also known as diffusive shock acceleration, is the underlying acceleration mechanism. At shock waves driven by coronal mass ejections (CMEs), particles are accelerated diffusively at the shock and often reach > MeV energies (and perhaps up to GeV energies). As a CME-driven shock propagates, expands and weakens, the accelerated particles can escape ahead of the shock into the interplanetary medium. These escaping energized particles then propagate along the interplanetary magnetic field, experiencing only weak scattering from fluctuations in the interplanetary magnetic field (IMF). In this paper, we use a Monte-Carlo approach to study the transport of energetic particles escaping from a CME-driven shock. We present particle spectra observed at 1 AU. We also discuss the particle “crossing number” at 1AU and its implication to particle anisotropy. Based on previous models of particle acceleration at CME-driven shocks, our simulation allows us to investigate various characteristics of energetic particles arriving at various distances from the sun. This provides us an excellent basis for understanding the observations of high-energy particles made at 1 AU by ACE and WIND.     

Relationship of solar flare accelerated particles to solar energetic particles (SEPs) observed in the interplanetary medium

Observations of hard X-ray (HXR)/γ-ray continuum and γ-ray lines produced by energetic electrons and ions, respectively, colliding with the solar atmosphere, have shown that large solar flares can accelerate ions up to many GeV and electrons up to hundreds of MeV. Solar energetic particles (SEPs) are observed by spacecraft near 1 AU and by ground-based instrumentation to extend up to similar energies as in large SEP events, but it appears that a different acceleration process, one associated with fast coronal mass ejections is responsible. Much weaker SEP events are observed that are generally rich in electrons, ³He, and heavy elements. The energetic particles in these events appear to be similar to those accelerated in flares. The Ramaty high energy solar spectroscopic imager (RHESSI) mission provides high-resolution spectroscopy and imaging of flare HXRs and γ-rays. Such observations can provide information on the location, energy spectra, and composition of the flare accelerated energetic particles at the Sun. Here, preliminary comparisons of the RHESSI observations with observations of both energetic electron and ion near 1 AU are reviewed, and the implications for the particle acceleration and escape processes are discussed.     

Analytical approach to cosmic ray ionization by nuclei with charge Z in the middle atmosphere – Distribution of galactic CR effects

The effects of galactic and solar cosmic rays (CR) in the middle atmosphere are considered in this work. A new analytical approach for CR ionization by protons and nuclei with charge Z in the lower ionosphere and middle atmosphere is developed in this paper. For this purpose the ionization losses (dE/dh) according to the Bohr–Bethe–Bloch formula for the energetic charged particles are approximated in three different energy intervals. More accurate expressions for energy decrease E(h) and electron production rate profiles q(h) are derived. The obtained formulas allow comparatively easy computer programming. The integrand in q(h) gives the possibility for application of adequate numerical methods – such as Romberg method or Gauss quadrature, for the solution of the mathematical problem. On this way the process of interaction of cosmic ray particles with the upper, middle and lower atmosphere will be described much more realistically. Computations for cosmic ray ionization in the middle atmosphere are made. The full CR composition is taken into account: protons, Helium (α-particles), light L, medium M, heavy H and very heavy VH group of nuclei.