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.

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.

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.

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.

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.

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.

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.

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.

On the maximum energy of shock-accelerated cosmic rays at ultra-relativistic shocks

The maximum energy to which cosmic rays can be accelerated at weakly magnetised ultra-relativistic shocks is investigated. We demonstrate that for such shocks, in which the scattering of energetic particles is mediated exclusively by ion skin-depth scale structures, as might be expected for a Weibel-mediated shock, there is an intrinsic limit on the maximum energy to which particles can be accelerated. This maximum energy is determined from the requirement that particles must be isotropized in the downstream plasma frame before the mean field transports them far downstream, and falls considerably short of what is required to produce ultra-high-energy cosmic rays. To circumvent this limit, a highly disorganized field is required on larger scales. The growth of cosmic ray-induced instabilities on wavelengths much longer than the ion-plasma skin depth, both upstream and downstream of the shock, is considered. While these instabilities may play an important role in magnetic field amplification at relativistic shocks, on scales comparable to the gyroradius of the most energetic particles, the calculated growth rates have insufficient time to modify the scattering. Since strong modification is a necessary condition for particles in the downstream region to re-cross the shock, in the absence of an alternative scattering mechanism, these results imply that acceleration to higher energies is ruled out. If weakly magnetized ultra-relativistic shocks are disfavoured as high-energy particle accelerators in general, the search for potential sources of ultra-high-energy cosmic rays can be narrowed.

Galactic Edge Clouds I. Molecular Line Observations and Chemical Modelling of Edge Cloud 2

Edge Cloud 2 (EC2) is a molecular cloud, about 35 pc in size, with one of the largest galactocentric distances known to exist in the Milky Way. We present observations of a peak CO emission region in the cloud and use these to determine its physical characteristics. We calculate a gas temperature of 20 K and a density of n(H2)~10^4 cm-3. Based on our CO maps, we estimate the mass of EC2 at around 10^4 Msolar and continuum observations suggest a dust-to-gas mass ratio as low as 0.001. Chemical models have been developed to reproduce the abundances in EC2, and they indicate that heavy element abundances may be reduced by a factor of 5 relative to the solar neighborhood (similar to dwarf irregular galaxies and damped Lya systems), very low extinction (A_V <4 mag) due to a very low dust-to-gas mass ratio, an enhanced cosmic-ray ionization rate, and a higher UV field compared to local interstellar values. The reduced abundances may be attributed to the low level of star formation in this region and are probably also related to the continuing infall of primordial (or low-metallicity) halo gas since the Milky Way formed. Finally, we note that shocks from the old supernova remnant GSH 138-01-94 may have determined the morphology and dynamics of EC2.

The 90-degree problem of scattering theory revisited: dynamical turbulence effects versus nonlinear effects

The 90° problem of cosmic-ray transport theory is revisited in this paper. By using standard forms of the wave spectrum in the solar wind, the pitch-angle Fokker–Planck coefficient and the parallel mean free path are computed for different resonance functions. A critical comparison is made of the strength of 90° scattering due to plasmawave effects, dynamical turbulence effects and nonlinear effects. It is demonstrated that, only for low-energy cosmic particles, dynamical effects are usually dominant. The novel results presented here are essential for an effective comparison of heliospheric observations for the parallel mean free path with the theoretical model results.

Chemical Processes in Protoplanetary Disks

We have developed a high-resolution combined physical and chemical model of a protoplanetary disk surrounding a typical T Tauri star. Our aims were to use our model to calculate the chemical structure of disks on small scales (submilliarcsecond in the inner disk for objects at the distance of Taurus, ~140 pc) to investigate the various chemical processes thought to be important in disks and to determine potential molecular tracers of each process. Our gas-phase network was extracted from the UMIST Database for Astrochemistry to which we added gas–grain interactions including freezeout and thermal and non-thermal desorption (cosmic-ray-induced desorption, photodesorption, and X-ray desorption), and a grain-surface network. We find that cosmic-ray-induced desorption has the least effect on our disk chemical structure while photodesorption has a significant effect, enhancing the abundances of most gas-phase molecules throughout the disk and affecting the abundances and distribution of HCN, CN, and CS, in particular. In the outer disk, we also see enhancements in the abundances of H2O and CO2. X-ray desorption is a potentially powerful mechanism in disks, acting to homogenize the fractional abundances of gas-phase species across the depth and increasing the column densities of most molecules, although there remain significant uncertainties in the rates adopted for this process. The addition of grain-surface chemistry enhances the fractional abundances of several small complex organic molecules including CH3OH, HCOOCH3, and CH3OCH3 to potentially observable values (i.e., a fractional abundance of greater than 10-11).

Chemical evolution in collapsing cores

We study the chemical evolution in the central core of contracting interstellar clouds. The chemical rate equations and the hydrodynamic equations are integrated simultaneously. The. contraction is followed from very low density (n = 10 cm(-3)) to a high-density core with n > 10(7) cm(-3). The chemical evolution is studied for various physical and chemical conditions, including the effects of varying the cosmic ray ionization rate, in order to understand the observed structures in TMC-1 and the extended ridge cloud in Orion. Our results give good agreement with the observations for models with fast ion-dipole reaction rates, low cosmic ray ionization rates and low depletion of N and S. It is also found that there should be different stages of evolution with different densities in these sources.