A New Theory for Perpendicular Transport of Cosmic Rays

We present an improved nonlinear theory for the perpendicular transport of charged particles. This approach is based on an improved nonlinear treatment of field-line random walk in combination with a generalized compound diffusion model. The generalized compound diffusion model employed is more systematic and reliable, in comparison with previous theories. Furthermore, the theory shows remarkably good agreement with test-particle simulations and solar wind observations.

THE EFFECT OF VARYING COSMIC-RAY IONIZATION RATES ON DARK CLOUD CHEMISTRY

We investigate the effects of varying the cosmic ray ionization rate in chemical models of dense interstellar clouds. In the absence of such ionization, a scenario which may be applicable to dark cloud cores, we find that chemi-ionization is able to drive a limited ion-neutral chemistry. Models of clouds in starburst galaxies, which may have enhanced cosmic ray fluxes, are also investigated and enable an upper limit to be derived for the cosmic ray ionization rate in M82. The derived value, which is about 700 times the typical value for Galactic molecular clouds, is in good agreement with that necessary to explain the recent observations of C I in this galaxy.

A filamentation instability for streaming cosmic rays

We demonstrate that cosmic rays form filamentary structures in the precursors of supernova remnant shocks due to their self-generated magnetic fields. The cosmic ray filamentation results in the growth of a long-wavelength instability, and naturally couples the rapid non-linear amplification on small scales to larger length-scales. Hybrid magnetohydrodynamics-particle simulations are performed to confirm the effect. The resulting large-scale magnetic field may facilitate the scattering of high-energy cosmic rays as required to accelerate protons beyond the knee in the cosmic ray spectrum at supernova remnant shocks. Filamentation far upstream of the shock may also assist in the escape of cosmic rays from the accelerator.

Cosmic-ray acceleration and escape from supernova remnants

Galactic cosmic-ray (CR) acceleration to the knee in the spectrum at a few PeV is only possible if the magnetic field ahead of a supernova remnant (SNR) shock is strongly amplified by CRs escaping the SNR. A model formulated in terms of the electric charge carried by escaping CRs predicts the maximum CR energy and the energy spectrum of CRs released into the surrounding medium. We find that historical SNRs such as Cas A, Tycho and Kepler may be expanding too slowly to accelerate CRs to the knee at the present time.

Universal behaviour of shock precursors in the presence of efficient cosmic ray acceleration

The self-consistent interaction between energetic particles and self-generated hydromagnetic waves in a cosmic ray pressure dominated plasma is considered. Using a three-dimensional hybrid magnetohydrodynamics (MHD)-kinetic code, which utilizes a spherical harmonic expansion of the Vlasov-Fokker-Planck equation, high-resolution simulations of the magnetic field growth including feedback on the cosmic rays are carried out. It is found that for shocks with high cosmic ray acceleration efficiency, the magnetic fields become highly disorganized, resulting in near isotropic diffusion, independent of the initial orientation of the ambient magnetic field. The possibility of sub-Bohm diffusion is demonstrated for parallel shocks, while the diffusion coefficient approaches the Bohm limit from below for oblique shocks. This universal behaviour suggests that Bohm diffusion in the root-mean-squared field inferred from observation may provide a realistic estimate for the maximum energy acceleration time-scale in young supernova remnants. Although disordered, the magnetic field is not self-similar suggesting a non-uniform energy-dependent behaviour of the energetic particle transport in the precursor. Possible indirect radiative signatures of cosmic ray driven magnetic field amplification are discussed.

Diffusive shock acceleration at laser-driven shocks: studying cosmic-ray accelerators in the laboratory

The non-thermal particle spectra responsible for the emission from many astrophysical systems are thought to originate from shocks via a first order Fermi process otherwise known as diffusive shock acceleration. The same mechanism is also widely believed to be responsible for the production of high energy cosmic rays. With the growing interest in collisionless shock physics in laser produced plasmas, the possibility of reproducing and detecting shock acceleration in controlled laboratory experiments should be considered. The various experimental constraints that must be satisfied are reviewed. It is demonstrated that several currently operating laser facilities may fulfil the necessary criteria to confirm the occurrence of diffusive shock acceleration of electrons at laser produced shocks. Successful reproduction of Fermi acceleration in the laboratory could open a range of possibilities, providing insight into the complex plasma processes that occur near astrophysical sources of cosmic rays.

Cosmic ray acceleration at oblique shocks

We show that the diffusion approximation breaks down for particle acceleration at oblique shocks with velocities typical of young supernova remnants. Higher order anisotropies flatten the spectral index at quasi-parallel shocks and steepen the spectral index at quasi-perpendicular shocks. We compare the theory with observed spectral indices.

Stochastic particle acceleration in the lobes of giant radio galaxies

We investigate the acceleration of particles by Alfven waves via the second-order Fermi process in the lobes of giant radio galaxies. Such sites are candidates for the accelerators of ultra-high-energy cosmic rays (UHECR). We focus on the nearby Fanaroff-Riley type I radio galaxy Centaurus A. This is motivated by the coincidence of its position with the arrival direction of several of the highest energy Auger events. The conditions necessary for consistency with the acceleration time-scales predicted by quasi-linear theory are reviewed. Test particle calculations are performed in fields which guarantee electric fields with no component parallel to the local magnetic field. The results of quasi-linear theory are, to an order of magnitude, found to be accurate at low turbulence levels for non-relativistic Alfven waves and at both low and high turbulence levels in the mildly relativistic case. We conclude that for pure stochastic acceleration via Alfven waves to be plausible as the generator of UHECR in Cen A, the baryon number density would need to be several orders of magnitude below currently held upper limits.

TWO-DIMENSIONAL PARTICLE-IN-CELL SIMULATIONS OF THE NONRESONANT, COSMIC-RAY-DRIVEN INSTABILITY IN SUPERNOVA REMNANT SHOCKS

In supernova remnants, the nonlinear amplification of magnetic fields upstream of collisionless shocks is essential for the acceleration of cosmic rays to the energy of the "knee" at 10(15.5) eV. A nonresonant instability driven by the cosmic ray current is thought to be responsible for this effect. We perform two-dimensional, particle-in-cell simulations of this instability. We observe an initial growth of circularly polarized nonpropagating magnetic waves as predicted in linear theory. It is demonstrated that in some cases the magnetic energy density in the growing waves can grow to at least 10 times its initial value. We find no evidence of competing modes, nor of significant modification by thermal effects. At late times, we observe saturation of the instability in the simulation, but the mechanism responsible is an artifact of the periodic boundary conditions and has no counterpart in the supernova-shock scenario.

The transport of cosmic rays in self-excited magnetic turbulence

The process of diffusive shock acceleration relies on the efficacy with which hydromagnetic waves can scatter charged particles in the precursor of a shock. The growth of self-generated waves is driven by both resonant and non-resonant processes. We perform high-resolution magnetohydrodynamic simulations of the non-resonant cosmic ray driven instability, in which the unstable waves are excited beyond the linear regime. In a snapshot of the resultant field, particle transport simulations are carried out. The use of a static snapshot of the field is reasonable given that the Larmor period for particles is typically very short relative to the instability growth time. The diffusion rate is found to be close to, or below, the Bohm limit for a range of energies. This provides the first explicit demonstration that self-excited turbulence reduces the diffusion coefficient and has important implications for cosmic-ray transport and acceleration in supernova remnants.

A cosmic ray current-driven instability in partially ionised media

Context. We investigate the growth of hydromagnetic waves driven by streaming cosmic rays in the precursor environment of a supernova remnant shock.

Aims. It is known that transverse waves propagating parallel to the mean magnetic field are unstable to anisotropies in the cosmic ray distribution, and may provide a mechanism to substantially amplify the ambient magnetic field. We quantify the extent to which temperature and ionisation fractions modify this picture.

Methods. Using a kinetic description of the plasma we derive the dispersion relation for a collisionless thermal plasma with a streaming cosmic ray current. Fluid equations are then used to discuss the effects of neutral-ion collisions.

Results. We calculate the extent to which the environment into which the cosmic rays propagate influences the growth of the magnetic field, and determines the range of possible growth rates.

Conclusions. If the cosmic ray acceleration is efficient, we find that very large neutral fractions are required to stabilise the growth of the non-resonant mode. For typical supernova parameters in our Galaxy, thermal effects do not significantly alter the growth rates. For weakly driven modes, ion-neutral damping can dominate over the instability at more modest ionisation fractions. In the case of a supernova shock interacting with a molecular clouds, such as in RX J1713.7-3946, with high density and low ionisation, the modes can be rapidly damped.

ENVIRONMENTAL LIMITS ON THE NONRESONANT COSMIC-RAY CURRENT-DRIVEN INSTABILITY

We investigate the so-called nonresonant cosmic-ray streaming instability, first discussed by Bell (2004). The extent to which thermal damping and ion-neutral collisions reduce the growth of this instability is calculated. Limits on the growth of the nonresonant mode in SN1006 and RX J1713.7-3946 are presented.