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.     

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.     

Modeling the acceleration and modulation of anomalous cosmic ray oxygen

After the solar wind termination shock crossings of the Voyager spacecraft, the acceleration of anomalous cosmic rays has become a very contentious subject. In this paper we examine several topics pertinent to anomalous cosmic ray oxygen acceleration and transport using a numerical cosmic ray modulation model. These include the effects of drifts on a purely Fermi I accelerated spectra, the effects of introducing higher charge states of oxygen into the modulation model, examining the viability of momentum diffusion as a re-acceleration process in the heliosheath and examining energy spectra, and intensity gradients, in the inner heliosphere during consecutive drift cycles.     

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.     

Particle acceleration by stochastic fluctuations and dawn-dusk electric field in the Earth’s magnetotail

This work is devoted to investigate the interaction between protons and stochastic time-dependent electromagnetic fields generated by oscillating clouds of finite size, randomly positioned in the x–y  -plane. The geometry of the system is two-dimensional and, beside the time-dependent electromagnetic fluctuations, a steady-state, dawn-dusk electric field, _Ey$E_y$, has been added along the y-direction. The simultaneous presence of the stochastic time-dependent fluctuations and of the constant electric field component in the same system gives rise to two types of acceleration mechanisms operating on test particles: a second order Fermi-like process and a direct acceleration. By performing a parametric study, we extensively study the contribution of the two processes to proton acceleration. The energy values reached by test particles in this simple model are in good agreement with those observed in the Earth’s magnetotail region. Possible applications to the problem of particle acceleration in the terrestrial magnetosphere are widely discussed, and guidelines for future works are drawn.     

Relative ion Fe,C and O abundances in quiet time particle fluxes in the 23 SC

Using ACE and SOHO data the origin of quiet-time low-energy particle fluxes at 1 AU is studied in the 23rd solar cycle. One of the selection criteria of quiet-time periods is to demand that H/He < 10 provided that periods with noticeable contribution of remnants of gradual events have been excluded from consideration. Our results suggest different origin of 0.03–3 MeV/nucleon particles – different seed populations accelerated and different acceleration processes. During the ascending, maximum and descending phases of solar activity quiet-time ions consist of coronal particles accelerated to suprathermal energies in about a half of the quiet periods, the rest of quiet-time fluxes originates from particle acceleration in processes similar to those in small impulsive solar flares rich in Fe. At solar minimum the bulk solar wind particles serve as seed population.     

2008年4月24日行星际激波事件相关联的超热电子90°投掷角的增强

利用测试粒子数值模拟的方法研究了与STEREO-A卫星观测到的2008年4月24日行星际激波事件相关联的超热电子90°投掷角的增强.根据激波到达前给定时刻超热电子的观测分布,拟合得到不同投掷角的初始分布函数;在给定的激波参数下,采用时间向后的方法计算特定能道上激波下游超热电子的投掷角分布.由于超热电子具有较高的共振频率,模拟采用的磁场湍流谱包含了低能电子发生共振的耗散区.对以215.76,151.67,106.63,eV为中心的三个能道进行了模拟.结果表明,不同能道上超热电子在激波下游的投掷角分布均在90°投掷角附近出现峰值,呈现出明显的90°投掷角增强,这与观测结果符合得很好.可以认为在激波对电子的加速过程中,电子与湍流耗散区的共振对90°投掷角的增强具有重要作用.      

Particle acceleration and gamma rays in solar flares: Recent observations and new modeling

Experiments on SMM, GAMMA, Yohkoh, GRANAT, Compton GRO, INTEGRAL, RHESSI and CORONAS-F satellites over the past three decades have provided copious data for fundamental research relating to particle acceleration, transport and energetics of flares and to the ambient abundance of the solar corona, chromosphere and photosphere. We summarize main results of solar gamma-astronomy (including some results of several joint Russian–Chinese projects) and try to appraise critically a real contribution of those results into modern understanding of solar flares, particle acceleration at the Sun and some properties of the solar atmosphere. Recent findings based on the RHESSI, INTEGRAL and CORONAS-F measurements (source locations, spectrum peculiarities, ³He abundance etc.) are especially discussed. Some unusual features of extreme solar events (e.g., 28 October 2003 and 20 January 2005) have been found in gamma-ray production and generation of relativistic particles (solar cosmic rays, or SCR). A number of different plausible assumptions are considered concerning the details of underlying physical processes during large flares: (1) existence of a steeper distribution of surrounding medium density as compared to a standard astrophysical model (HSRA) for the solar atmosphere; (2) enhanced content of the ³He isotope; (3) formation of magnetic trap with specific properties; (4) prevailing non-uniform (e.g., fan-like) velocity (angular) distributions of secondary neutrons, etc. It is emphasized that real progress in this field may be achieved only by combination of gamma-ray data in different energy ranges with multi-wave and energetic particle observations during the same event. We especially note several promising lines for the further studies: (1) resonant acceleration of the ³He ions in the corona; (2) timing of the flare evolution by gamma-ray fluxes in energy range above 90 MeV; (3) separation of gamma-ray fluxes from different sources at/near the Sun (e.g., different acceleration sources/episodes during the same flare, contribution of energetic particles accelerated by the CME-driven shocks etc.); (4) asymmetric magnetic geometry and new magnetic topology models of the near-limb flares; (5) modeling of self-consistent time scenario of the event.     

Positive ion distributions in the morning auroral zone: Local acceleration and drift effects

We report on the typical structure of the large scale ion precipitation in the morning sector of the auroral zone and associated low frequency electromagnetic waves. Data obtained during near radial passes of the AUREOL-3 satellite point to a distinction between two main precipitation regions: 1) In the poleward part of the auroral zone the latitudinal variation of the average energy (or temperature) of the precipitated ions (mainly H⁺) indicate that they are adiabatically accelerated in the outer magnetosphere. This “high energy” (? 3 to > 20 keV) precipitation is usually associated with a low energy (E < 110 eV) upward flowing 0⁺ and H⁺ component, and 2) near the boundary between discrete and diffuse electron aurorae a drastic change in the ion characteristics is observed. The flux of energetic precipitated H⁺ ions is sharply reduced, which suggests the formation of an Alfvén layer. However, intense fluxes of precipitated H⁺, O⁺, and He⁺ ions with energies < 3 keV are observed equatorward of the Alfvén layer, in coincidence with the diffuse aurora and in association with quasi-monochromatic electromagnetic waves with frequencies around the proton gyrofrequency. As the characteristic convection and bounce times of the low energy upward flowing ion component are comparable (τ > 3 hours) we suggest that the precipitation of ionospheric ions inside the diffuse aurora results from convection and corotation of the ions accelerated to suprathermal energies at higher latitudes.     

垂直激波条件下能量粒子加速机制的模拟研究

在具有湍动的磁场和垂直激波条件下对大量测试粒子的轨迹进行了数值计算,研究了激波强度和粒子初始能量对于粒子穿越激波的平均能量变化的影响,分析了漂移加速(SDA)在不同条件下对粒子加速的贡献,并给出了一个与数值结果相符合的漂移加速理论公式△E=amvivup(1-1/s).结果表明,加入磁场湍流后,垂直激波条件下粒子仍主要受到漂移加速作用,而基于粒子引导中心的耗散漂移加速理论在此条件下失效.      

Using gamma-ray burst prompt emission to probe relativistic shock acceleration

It is widely accepted that the prompt transient signal in the 10 keV–10 GeV band from gamma-ray bursts (GRBs) arises from multiple shocks internal to the ultra-relativistic expansion. The detailed understanding of the dissipation and accompanying acceleration at these shocks is a currently topical subject. This paper explores the relationship between GRB prompt emission spectra and the electron (or ion) acceleration properties at the relativistic shocks that pertain to GRB models. The focus is on the array of possible high-energy power-law indices in accelerated populations, highlighting how spectra above 1 MeV can probe the field obliquity in GRB internal shocks, and the character of hydromagnetic turbulence in their environs. It is emphasized that diffusive shock acceleration theory generates no canonical spectrum at relativistic MHD discontinuities. This diversity is commensurate with the significant range of spectral indices discerned in prompt burst emission. Such system diagnostics are now being enhanced by the broad-band spectral coverage of bursts by the Fermi Gamma-Ray Space Telescope; while the Gamma-Ray Burst Monitor (GBM) provides key diagnostics on the lower energy portions of the particle population, the focus here is on constraints in the non-thermal, power-law regime of the particle distribution that are provided by the Large Area Telescope (LAT).     

Detection of suprathermal ionospheric O+ ions inside the plasmasphere

The plasma diagnostic experiments on the AUREOL-3 satellite have revealed flows of low energy 0⁺ ions deep inside the night plasmasphere during a large substorm. Flux gradients of the 0⁺ ions were accompanied by enhancements of ELF electric field noise. The appearance of suprathermal ions at $\text{L ? 2.5 – 3}$ is interpreted within the framework of electrostatic ion-cyclotron acceleration of ionospheric ions in the diffuse auroral zone /12/ followed by a radial displacement of these ions inside the plasmasphere driven by azimuthal electric fields during substorm activity. Electrostatic oscillations observed inside the plasmasphere are apparently associated with gradient instability at the sharp boundaries of suprathermal ion flows.     

Acceleration of charged particles by fluctuating and steady electric fields in a X-type magnetic field

Release of stored magnetic energy via particle acceleration is a characteristic feature of astrophysical plasmas. Magnetic reconnection is one of the mechanisms for releasing energy from magnetized plasmas. Collisionless magnetic reconnection could provide both the energy release mechanism and the particle accelerator in space plasmas. Here we studied particle acceleration when fluctuating (in-time) electric fields are superposed on an static X-type magnetic field in collisionless hot solar plasma. This system is chosen to mimic the reconnective dissipation of a linear MHD disturbance. Our results are compared to particle acceleration from constant electric field superposed on an X-type magnetic field. The constant electric field configuration represents the effects of steady state magnetic reconnection. Time evolution of ion and electron distributions are obtained by numerically integrating particle trajectories. The frequencies of the electric field represent a turbulent range of waves. Depending on the frequency and amplitude of the electric field, electrons and ions are accelerated to different degrees and have energy distributions of bimodal form consisting of a lower energy part and a high energy tail. For frequencies (ω in dimensioless units) in the range 0.5 ? ω ? 1.0 a substantial fraction (20%–30%) of the proton distribution is accelerated to gamma-ray producing energies. For frequencies in the range 1 ? ω ? 100.0 the bulk of the electron distribution is accelerated to hard X-ray producing energies. The acceleration mechanism is important for solar flares and solar noise storms but it could be applicable to all collisionless astrophysical plasmas.     

The plasma regime at comet Giacobini-Zinner

This review of the plasma regime sampled by the encounter of the International Cometary Explorer spacecraft (ICE) with the comet Giacobini-Zinner, discusses the shock, or bow wave, ion pickup, ionization mechanisms, and the cometary plasma tail.

The observations are consistent with the existence of a weak shock, which may be pulsating, but do not exclude the suggestion by Wallis and Dryer that the shock, though present around the sub-solar point, is in process of decaying to a wave on the flanks.

Pickup of cometary ions provokes, by means of several mechanisms, ion cyclotron, mirror, beam and electrostatic instabilities which cause strong turbulence in the inner coma, as indicated in the power spectra of the magnetic field in the coma and the surrounding volume. Heavy mass loading and consequent slowing down of the solar wind is observed. Acceleration of ions by a stochastic mechanism is indicated.

Ionization of cometary neutrals occurs principally by photoionization and charge exchange. Alfvens critical velocity mechanism, likely operates only in the inner coma not visited by ICE. A steep increase of nearly two orders of magnitude in electron density occurs in the tail, where electron velocity distributions show evidence of entry of electrons from the solar wind. The turbulence there is damped by the high ion density and low temperature.

In general, the vicinity of the comet is filled with plasma phenomena and a rich variety of corresponding atomic and molecular processes can be studied there. Comparison between the ICE, Giotto, and Vega observations forms a most valuable future study.     

Particle acceleration and transport at an oblique CME-driven shock

In gradual solar energetic particle (SEP) events, protons and heavy ions are often accelerated to >100 MeV/nucleon at a CME-driven shock. In this work, we study particle acceleration at an oblique shock by extending our earlier particle acceleration and transport in heliosphere (PATH) code to include shocks with arbitrary θ_BN, where θ_BN is the angle between the upstream magnetic field and the shock normal. Instantaneous particle spectra at the shock front are obtained by solving the transport equation using the total diffusion coefficient κ, which is a function of the parallel diffusion coefficient κ_∥ and the perpendicular diffusion coefficient κ_⊥. In computing κ_∥ and κ_⊥, we use analytic expressions derived previously. The particle maximum energy at the shock front as a function of time, the time intensity profiles and particle spectra at 1 AU for five θ_BN’s are calculated for an example shock.     

A model for stochastic acceleration of electrons during geomagnetic storms

The theory of resonant diffusion is extended to fully relativistic plasmas, and we examine resonant interactions between electrons and electromagnetic R mode (whistler) and L-mode (EMIC) waves. Resonant diffusion curves are constructed for plasma parameters representative of the Earth's storm time magnetosphere, both inside and outside the plasmapause. EMIC waves can resonate with electrons > 1 MeV, but the energies remain nearly constant along the diffusion curves. Storm-time EMIC waves can induce rapid pitch—angle scattering, but the waves are ineffective for stochastic acceleration of elections. Substantial energy change can occur along the diffusion curves for interactions between resonant electrons and whistler—mode waves, especially in regions of low plasma density. Specifically, whistlers can accelerate electrons from energies near 100 keV to above 1 MeV outside the plasmapause. A model is proposed comprising energy diffusion by whistler-mode chorus and pitch-angle scattering by EMIC waves to account for the gradual acceleration of electrons over the region 4 ≤ L ≤ 6 during the recovery phase of a geomagnetic storm.     

Shock acceleration of energetic particles in corotating interaction regions

A simple shock model for the acceleration of energetic particles in corotating interaction regions (CIR) in the solar wind is presented. Particles are accelerated at the forward and reverse shocks which bound the CIR by being compressed between the shock fronts and magnetic irregularities upstream from the shocks, or by being compressed between upstream irregularities and those downstream from the shocks. Particles also suffer adiabatic deceleration in the expanding solar wind, an effect not included in previous shock models for acceleration in CIR's. The model is able to account for the observed exponential spectra at earth, the observed behavior of the spectra with radial distance, the observed radial gradients in the intensity, and the observed differences in the intensity and spectra at the forward and reverse shocks.     

Cosmic rays in the dynamic heliosheath

Cosmic ray modulation in the outer heliosphere is discussed from a modeling perspective. Emphasis is on the transport and acceleration of these particles at and beyond the solar wind termination shock in the inner heliosheath region and how this changes over a solar cycle. We will show that by using numerical models, and by comparing results to spacecraft observations, much can be learned about the dependence of cosmic ray modulation on solar cycle changes in the solar wind and heliospheric magnetic field. While the first determines the heliospheric geometry and shock structure, the latter results in a time-dependence of the transport coefficients. Depending on energy, both these effects contribute to cosmic ray intensities in the inner heliosheath changing over a solar cycle.     

Diffusive shock acceleration at interplanetary perpendicular shock waves: Influence of the large scale structure of turbulence on the maximum particle energy

We calculate the maximum energy that a particle can obtain at perpendicular interplanetary shock waves by the mechanism of diffusive shock acceleration. The influence of the energy range spectral index of the two-dimensional modes of the interplanetary turbulence is explored. We show that changes in this parameter lead to energies that differ in at least one order of magnitude. Therefore, the large scale structure of the turbulence is a key input if the maximum particle energy is calculated.