The physics of cosmic-ray modulation

The theory of the modulation of galactic cosmic rays by the solar wind is reviewed. The basic transport equation is presented, interpreted and then applied to cosmic-ray transport in a model heliosphere immersed in a constant uniform bath of galactic cosmic rays. The results of numerical modelling are presented and the dominant physical effects analyzed. A variety of observational tests of the model which were reported over the last several years are summarized and shown, generally, to support a model in which particle drifts play an important role. Recent measurements which show that the latitudinal gradient of cosmic rays changed sign in the recent sunspot minimum (relative the last sunspot minimum) are shown to provide additional, strong, support for the model. A new picture of the interplanetary magnetic field is presented, which gives promise of improving considerably the agreement between the theory and observations in the few remaining problem areas.     

Particle acceleration by coronal and interplanetary shock waves

Utilizing many years of observation from deep space and near-earth spacecraft a theoretical understanding has evolved on how ions and electrons are accelerated in interplanetary shock waves. This understanding is now being applied to solar flare-induced shock waves propagating through the solar atmosphere. Such solar flare phenomena as γ-ray line and neutron emissions, interplanetary energetic electron and ion events, and Type II and moving Type IV radio bursts appear understandable in terms of particle accleration in shock waves.     

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.     

Cosmic-ray composition and its relation to shock acceleration by supernova remnants

An overview is given on the present status of the understanding of the origin of galactic cosmic rays. Recent measurements of charged cosmic rays and photons are reviewed. Their impact on the contemporary knowledge about the sources and acceleration mechanisms of cosmic rays and their propagation through the Galaxy is discussed. Possible reasons for the knee in the energy spectrum and scenarios for the end of the galactic cosmic-ray component are described.     

Heliospheric implications of structure in the interstellar medium

Equilibrium models of diffuse interstellar material (ISM) near the Sun show a range of cloud densities, ionization, and temperatures which are consistent with data, although the local ISM must be inhomogeneous over ∼2 pc scales. The ISM close to the Sun has properties that are consistent with the sheetlike warm neutral (and partially ionized) gas detected in the Arecibo Millennium Survey. Local interstellar magnetic fields are poorly understood, but data showing weak polarization for nearby stars indicate dust may be trapped in fields or currents in the heliosheath nose region. Implications of this dust capture are widespread, and may impact the interpretation of the cosmic microwave background data. Observations of interstellar H⁰ inside of the solar system between 1975 and 2000 do not suggest any variation in the properties or structure of local interstellar H⁰ over distance scales of ∼750 AU to within the uncertainties.     

Interplanetary acceleration of low-energy protons as observed during the 25 september 1978 shock event

ISEE-3 observations of a long-lasting low-energy proton intensity increase during the 25 September 1978 shock event are presented as an example for interplanetary particle acceleration in association with shock waves. The observations are discussed in the light of current models for particle acceleration. The particular shape of the time intensity behaviour of the particle intensity increase, the existence of a shock spike and the observed particle distributions indicate that the particles are accelerated at the shock by the induced electric field $\text{E}\text{= ?}\text{1}\text{c}\text{V}\text{×}\text{B}$.     

Nonlinear processes in cosmic-ray precursor of strong supernova shock: Maximum energy and average energy spectrum of accelerated particles

The instability in the cosmic-ray precursor of a supernova shock is studied. The level of turbulence in this region determines the maximum energy of accelerated particles. The consideration is not limited by the case of weak turbulence. It is assumed that the Kolmogorov type nonlinear wave interactions together with the ion-neutral collisions restrict the amplitude of random magnetic field. As a result, the maximum energy of accelerated particles strongly depends on the age of a SNR. The average spectrum of cosmic rays injected in the interstellar medium in the course of adiabatic SNR evolution takes the approximate form E⁻² at energies larger than 10–30 GeV/nucleon with the maximum energy that is close to the position of the knee in cosmic-ray spectrum at 4 × 10¹⁵ eV. At an earlier stage of SNR evolution – the ejecta-dominated stage, the particles are accelerated to higher energies and have a rather steep power-law distribution. These results suggest that the knee may mark the transition from the ejecta-dominated to the adiabatic evolution of SNR shocks which accelerate cosmic rays.     

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.     

Propagation of cosmic-ray electrons in the Galaxy

We formulate the global propagation model of cosmic-ray electrons including the source region, which is currently considered to be supernova remnants (SNRs). The model is characterized by the escape rate of electrons from SNRs into the interstellar space. It becomes clear that the energy index of the escape rate influences the high energy side of the interstellar spectrum and makes it possible to explain the observed data up to 2 TeV in the case of source spectral index smaller than 2.2 that is expected from the radio spectrum in SNRs. The escape lifetime of electrons in SNRs is also discussed by using the ratio of the radio flux in two regions: SNRs and the Galaxy. The result shows the mean lifetime in SNRs of ∼10⁴ yr around 1 GeV, which corresponds to the SNR age in the Sedov phase.     

Origin of coronal and interplanetary shock and particle acceleration of a flare/CME event

By using radio data from ground-based telescopes (from 270 MHz to 25 MHz), and from the Radio and Plasma Wave experiment (WAVES) on board the WIND spacecraft (1–14 MHz and several kHz-11 MHz), as well as FY -2 satellite data, the origin of coronal and interplanetary shock and particle acceleration of the 14 July 2000 flare/CME event (the Bastille day event) have been studied. Main conclusions are as follows: (1) We investigate the causal relationship between metric type 11 bursts observed by the digital IZMIRAN radio spectrograph and type II radio emissions in the frequency range from 1–14 MHz and several kHz-11 MHz observed by the WAVES/WIND. The analysis indicate that the fast CME is the origin of both coronal and interplanetary shocks. (2)According to the time profiles of Hard X-ray, and energetic particles (include proton, ³He, and ⁴He) from FY-2 satellite, it is obvious that the Bastille day event is the event, in which both impulsive and gradual phenomena occur. The energetic particles accelerated not only in flare but also in CME.     

Fractal Structure of the Heliospheric Plasma Sheet at the Earth''''s Orbit

An analysis of the data from the Wind and IMP-8 spacecraft revealed that a slow solar wind, flowing in the heliospheric plasma sheet, represents a set of magnetic tubes with plasma of increased density (N > 10cm-3 at the Earth's orbit). They have a fine structure at several spatial scales (fractality), from 2°-3°(at the Earth's orbit, it is equivalent to 3.6-5.4h, or (5.4-8.0)×106km) to the minimum about 0.025°, i.e. the angular size of the nested tubes is changed nearly by two orders of magnitude. The magnetic tubes at each observed spatial scale are diamagnetic, i.e. their surface sustains a flow of diamagnetic (or drift) current that decreases the magnetic field within the tube itself and increases it outside the tube. Furthermore, the value of β= 8π[N(Te + Tp)]/B2 within the tube exceeds the value of βoutside the tube. In many cases total pressure P = N(Te + Tp) + B2/8πis almost constant within and outside the tubes at any one of the aforementioned scales.     

Heliospheric asymmetries due to the action of the interstellar magnetic field

We discuss the asymmetry of the heliospheric discontinuities obtained from the analysis of 3D modeling of the solar wind (SW) interaction with local interstellar medium (LISM). The flow of charged particles is governed by the ideal MHD equations and the flow of neutral particles is described by the Boltzmann equation. The emphasis is made on the asymmetries of the termination shock (TS) and the heliopause under the combined action of the interstellar and interplanetary magnetic fields (ISMF and IMF) in the presence of neutral hydrogen atoms whose transport through the heliosphere is modeled kinetically, using a Monte Carlo approach. We show that the deflection of neutral hydrogen flow from its original direction in the unperturbed LISM is highly anisotropic and evaluate a possible angle between the hydrogen deflection plane measured in the SOHO SWAN experiment and the plane containing the ISMF and LISM velocity vectors for different ISMF strengths. It is shown that the ISMF of a strength greater than 4 μG can account for the 10 AU difference in the TS heliocentric difference observed during its crossing by the Voyager 1 and Voyager 2 spacecraft, which however results in a larger discrepancy between the calculated and observed velocity distributions. The effect of a strong ISMF on the distribution of plasma quantities in the inner heliosheath and on 2–3 kHz radio emission is discussed.     

Interplanetary acceleration of low energy protons in association with shock waves: The 11 February 1979 event

The problem of interplanetary acceleration of low energy protons in association with shock waves is examined in the context of the specific event observed on 11 February 1979 on board the ISEE-3 spacecraft. This event has been selected for special study as it apparently was not associated with a solar flare event. The low energy proton telescope system on ISEE-3 measures the proton distribution function with good spectral, directional and temporal resolution from E_p = 35 keV. The evolution of the anisotropies and of the energy spectrum during the event are consistent with particle acceleration taking place in the vicinity of the shock wave.     

The IGY and beyond: A brief history of ground-based cosmic-ray detectors

The state of art of ground-based cosmic-ray research from its discovery to present is reviewed. After discovery of cosmic rays by Hess in 1912, the nature of the primary and secondary radiation was established from recordings by a variety of instruments, sensitive to various components of cosmic rays and operated at different latitudes, longitudes and altitudes, including instruments carried by balloons. The IGY formalized international co-operation and coordinated study of cosmic rays, which is vital for meaningful interpretation of cosmic-ray data. Data collected at different geographic locations require an effective cutoff rigidity as a data ordering parameter. This parameter is obtained from tracing trajectories of primary cosmic rays in the Earth’s magnetic field. After 50 years the world’s neutron monitor network remains still the backbone for studying intensity variations of primary cosmic rays in the rigidity ranges between 1 and 15 GV, associated with transport and with transient events. Also the penetrating muon and neutrino components of secondary cosmic rays have a long history of recording and fundamental problem investigations. Valuable data about composition and spectrum of primary cosmic rays in ever increasing high-energy regions have been obtained during the years of investigations with various configurations and types of extensive air shower detectors. The culture of personal involvement of the physicist in carrying out experiments and data acquisition characterized the continued vitality of cosmic-ray investigations ranging from its atmospheric, geomagnetic and heliospheric transport through to its solar and astrophysical origins.     

Comparison of the effects of two models for perpendicular diffusion on cosmic-ray latitudinal gradients

We compare the effects of two different models for perpendicular diffusion on the latitudinal gradients of galactic cosmic ray protons during solar minimum conditions. These two models correspond to the newly developed non-linear guiding center theory [Matthaeus, W.H., Qin, G., Bieber, J.W., Zank, G.P. Nonlinear collisionless perpendicular diffusion of charged particles. Astrophys. J. Lett., 590 (1), L53–L56, 2003] and the theory based on a velocity correlation function approach [Bieber, J.W., Matthaeus, W.H. Perpendicular diffusion and drift at intermediate cosmic-ray energies. Astrophys. J., 485 (2) 655–659, 1997]. In this ab initio study a steady-state two-dimensional numerical modulation model is used which incorporates a state-of-the-art turbulence model. We show that the non-linear guiding center theory predicts a mean free path that has a rigidity dependence that better accounts for the latitudinal gradients measured by Ulysses during its first fast latitude scan in 1994/1995.     

Kappa distribution functions and the main properties of auroral particle acceleration

The problem of auroral particle acceleration is analyzed taking into account the filling of the magnetosphere by accelerated particles of ionospheric origin and the trapping of particles between the magnetic and electrostatic barriers. It is taken into account that kappa distributions describe experimentally measured distribution functions better than do maxwellians. The existence of conjugate field-aligned potential drops in the north and south hemispheres is suggested. Field-aligned potential drops are formed if the value of field-aligned current is higher than the threshold determined by the free gazodynamic flow of electrons along field lines. Results of experimental observations showing the possibility of field-aligned potential drop concentration in double layers are summarized. The theory of kinetic double layers is developed taking into account the nonmaxwellian forms of distribution functions. It is shown that the kinetic treatment and the existence of nonmaxwellian distributions lead to real changes in the criteria required for double-layer formation.     

Radio galaxies and the acceleration of the universe beyond a redshift of one

Radio galaxies provide a means to determine the coordinate distance, the luminosity distance, the dimension-less luminosity distance, or the angular size distance to sources with redshifts as large as two. Dimensionless coordinate distances for 55 supernovae and 20 radio galaxies are presented and discussed here. The radio galaxy results are consistent with those obtained using supernovae, suggesting that neither method is plagued by unknown systematic errors. The acceleration parameter q(z) and the expansion rate H(z) or dimensionless expansion rate E(z) can be determined directly from the data without having to make assumptions regarding the nature or evolution of the “dark energy”. The expansion rate E(z) can be determined from the first derivative of the dimensionless coordinate distance, (dy/dz)⁻¹, and the acceleration parameter can be determined from a combination of the first and second derivatives of the dimensionless coordinate distance. A model-independent determination of E(z) will allow the properties and redshift evolution of the “dark energy” to be determined, and a model-independent determination of q(z) will allow the redshift at which the universe transitions from acceleration to deceleration to be determined directly. Determination of E(z) and q(z) may also elucidate possible systematic errors in the determinations of the dimensionless coordinate distances.     

Effects of various dissipation range onset models on the 26-day variations of low-energy galactic cosmic-ray electrons

The effect of various models presented by Leamon et al. (2000) for the dissipation range cutoff wavenumber on the 26-day variations of galactic cosmic-ray electrons in a Fisk-Parker hybrid field is investigated, by means of a three-dimensional steady-state numerical modulation code. Analytical expressions for the mean free paths parallel and perpendicular to the heliospheric magnetic field are adapted from the works of 31 and 28, respectively. Note that only solar minimum conditions are considered, and that only qualitative agreement with data is sought. Effective diffusion for galactic electrons pertaining to 26-day variations is found to be dominated by the ratio of the perpendicular to parallel mean free paths at low energies, and the relationship between changes in cosmic-ray intensities and the modulation parameter postulated by Zhang (1997) is found to no longer hold when this ratio drops below a critical value. Use of ion inertial scale dependent models for the dissipation range cutoff leads to possible second linearities in the relative amplitudes as functions of latitude gradient.     

Origin of complex type III-L events and electron acceleration

Type III-L bursts constitute a class of type III bursts that are intense, complex, and of long duration at hectometric wavelengths. They are often associated with major flares and fast coronal mass ejections. Several observations suggested that the electron beams that produce these complex hectometric emissions could be accelerated and injected in the low or in the middle corona. In this study, we revisit the origin of these bursts by tracing the progression of the events from the low corona to the interplanetary medium. We show that type III-L features are related to sudden changes in the radio emission observed at metric and decametric wavelengths, in particular the onset of new emitting sources at positions that can be at large distances from the flare site.