The physical processes involved in diffusion of Galactic cosmic rays in the interstellar medium are addressed. We study the possibility that the nonlinear MHD cascade sets the power-law spectrum of turbulence which scatters charged energetic particles. We find that the dissipation of waves due to the resonant interaction with cosmic ray particles may terminate the Kraichnan-type cascade below wavelengths 10^13 cm. The effect of this wave dissipation has been incorporated in the GALPROP numerical propagation code in order to asses the impact on measurable astrophysical data. The energy-dependence of the cosmic-ray diffusion coefficient found in the resulting self-consistent model may explain the peaks in the secondary to primary nuclei ratios observed at about 1 GeV/nucleon.

Shock jump conditions are obtained for steady-state, plane shocks with oblique magnetic fields and arbitrary flow speeds. For ultrarelativistic and nonrelativistic shocks, the jump conditions may be solved analytically. For mildly relativistic shocks, analytic solutions are obtained for isotropic pressure using an approximation for the adiabatic index that is valid in high sonic Mach number cases. In the more general case of gyrotropic pressure, the jump conditions cannot be solved analytically without additional assumptions, and the effects of gyrotropic pressure are investigated by parameterizing the distribution of pressure parallel and perpendicular to the magnetic field. Our numerical solutions reveal that relatively small departures from isotropy (e.g., about 20%) produce significant changes in the shock compression ratio, r, at all shock Lorentz factors, including ultrarelativistic ones, where an analytic solution with gyrotropic pressure is obtained. In particular, either dynamically important fields or significant pressure anisotropies can incur marked departures from the canonical gas dynamic value of r=3 for a shocked ultrarelativistic flow and this may impact models of particle acceleration in gamma-ray bursts and other environments where relativistic shocks are inferred. The jump conditions presented apply directly to test-particle acceleration, and will facilitate future self-consistent numerical modeling of particle acceleration at oblique, relativistic shocks.

The diffusion of galactic cosmic rays (CR) is considered. It is anticipated that the nonlinear MHD cascade sets the power-law spectrum of turbulence from the principle scale 100 pc to much smaller scales. We found that the dissipation of waves due to the resonant interaction with energetic particles may terminate the cascade at less than 1.e13 cm. The shape of CR diffusion coefficient that was found may explain the observed peaks of secondary to primary nuclei ratios at a few GeV/n.

The origin and evolution of isotopes of the lightest elements H2, He3, Li, Be, B in the universe is a key problem in such fields as astrophysics of CR, Galactic evolution, non-thermal nucleosynthesis, and cosmological studies. One of the major sources of these species is spallation by CR nuclei in the interstellar medium. On the other hand, it is the B/C ratio in CR and Be10 abundance which are used to fix the propagation parameters and thus the spallation rate. We study the production and Galactic propagation of isotopes of elements Z<6 using the numerical propagation code GALPROP and updated production cross sections.

In a recent paper (Moskalenko et al., 2002), it has been shown that the flux of secondary cosmic ray (CR) antiprotons appears to be contradictory to measurements of secondary to primary nuclei ratios in cosmic rays when calculated in the same Galactic propagation model. The contradiction appears as a value of the diffusion coefficient necessary to match the secondary ratios antiproton/proton and B/C. In particlular, it was shown that the reacceleration models designed to match secondary to primary nuclei ratios produce too few antiprotons. It is, however, clear that some reacceleration is unavoidable in the turbulent interstellar medium. Here we discuss an idea of how to improve reacceleration model by allowing for the damping of interstellar turbulences on the small scale by cosmic rays, mostly protons. This would lead to increase in the mean free path lenghts at low energies, the well-known phenomena empirically discovered in the Leaky Box models, thus producing less secondary nuclei. Antiprotons will remain almost non-affected due to their high energy threshold of production cross section.

In an attempt to understand the source and propagation of galactic cosmic rays we have employed the Modified Weighted Slab technique along with recent values of the relevant cross sections to compute primary to secondary ratios including B/C and Sub-Fe/Fe for different galactic propagation models. The models that we have considered are the disk-halo diffusion model, the dynamical halo wind model, the turbulent diffusion model and a model with minimal reacceleration. The modified weighted slab technique will be briefly discussed and a more detailed description of the models will be given. We will also discuss the impact that the various models have on the problem of anisotropy at high energy and discuss what properties of a particular model bear on this issue.

We have modeled the injection and acceleration of pickup ions at the solar wind termination shock and investigated the parameters needed to produce the observed Anomalous Cosmic Ray (ACR) fluxes. A non-linear Monte Carlo technique was employed which in effect solves the Boltzmann equation and is not restricted to near isotropic particle distribution functions. This technique models the injection of thermal and pickup ions, the acceleration of these ions, and the determination of the shock structure under the influence of the accelerated ions. The essential effects of injection are treated in a mostly self-consistent manner, including effects from shock obliquity, cross-field diffusion, and pitch-angle scattering. Using recent determinations of pickup ion densities, we are able to match the absolute flux of hydrogen in the ACRs by assuming that pickup ion scattering mean free paths, at the termination shock, are much less than an AU and that moderately strong cross-field diffusion occurs. Simultaneously, we match the flux ratios He$^{+}$/H$^+$ or O$^{+}$/H$^+$ to within a factor ~ 3. If the conditions of strong scattering apply, no pre-termination-shock injection phase is required and the injection and acceleration of pickup ions at the termination shock is totally analogous to the injection and acceleration of ions at highly oblique interplanetary shocks recently observed by the Ulysses spacecraft.

There has recently been interest in the role of inverse bremsstrahlung, the emission of photons by fast suprathermal ions in collisions with ambient electrons possessing relatively low velocities, in tenuous plasmas in various astrophysical contexts. This follows a long hiatus in the application of suprathermal ion bremsstrahlung to astrophysical models since the early 1970s. The potential importance of inverse bremsstrahlung relative to normal bremsstrahlung, i.e. where ions are at rest, hinges upon the underlying velocity distributions of the interacting species. In this paper, we identify the conditions under which the inverse bremsstrahlung emissivity is significant relative to that for normal bremsstrahlung in shocked astrophysical plasmas. We determine that, since both observational and theoretical evidence favors electron temperatures almost comparable to, and certainly not very deficient relative to proton temperatures in shocked plasmas, these environments generally render inverse bremsstrahlung at best a minor contributor to the overall emission. Hence inverse bremsstrahlung can be safely neglected in most models invoking shock acceleration in discrete sources such as supernova remnants. However, on scales > 100pc distant from these sources, Coulomb collisional losses can deplete the cosmic ray electrons, rendering inverse bremsstrahlung, and perhaps bremsstrahlung from knock-on electrons, possibly detectable.

We have modeled the injection and acceleration of pickup ions at the solar wind termination shock and investigated the parameters needed to produce the observed Anomalous Cosmic Ray (ACR) fluxes. A non-linear Monte Carlo technique was employed, which in effect solves the Boltzmann equation and is not restricted to near-isotropic particle distribution functions. This technique models the injection of thermal and pickup ions, the acceleration of these ions, and the determination of the shock structure under the influence of the accelerated ions. The essential effects of injection are treated in a mostly self-consistent manner, including effects from shock obliquity, cross-field diffusion, and pitch-angle scattering. Using recent determinations of pickup ion densities, we are able to match the absolute flux of hydrogen in the ACRs by assuming that pickup ion scattering mean free paths, at the termination shock, are much less than an AU and that modestly strong cross-field diffusion occurs. Simultaneously, we match the flux ratios He$^{+}$/H$^{+}$ or O$^{+}$/H$^{+}$ to within a factor $\sim 5$. If the conditions of strong scattering apply, no pre-termination-shock injection phase is required and the injection and acceleration of pickup ions at the termination shock is totally analogous to the injection and acceleration of ions at highly oblique interplanetary shocks recently observed by the Ulysses spacecraft. The fact that ACR fluxes can be modeled with standard shock assumptions suggests that the much-discussed "injection problem" for highly oblique shocks stems from incomplete (either mathematical or computer) modeling of these shocks rather than from any actual difficulty shocks may have in injecting and accelerating thermal or quasi-thermal particles.

We give a rigorous derivation of a theorem showing that charged particles in an arbitrary electromagnetic field with at least one ignorable spatial coordinate remain forever tied to a given magnetic-field line. Such a situation contrasts the significant motions normal to the magnetic field that are expected in most real three-dimensional systems. It is pointed out that, while the significance of the theorem has not been widely appreciated, it has important consequences for a number of problems and is of particular relevance for the acceleration of cosmic rays by shocks.

We have developed a Monte Carlo technique for self-consistently calculating the hydrodynamic structure of oblique, steady-state shocks, together with the first-order Fermi acceleration process and associated non-thermal particle distributions. This is the first internally consistent treatment of modified shocks that includes cross-field diffusion of particles. Our method overcomes the injection problem faced by analytic descriptions of shock acceleration, and the lack of adequate dynamic range and artificial suppression of cross-field diffusion faced by plasma simulations; it currently provides the most broad and versatile description of collisionless shocks undergoing efficient particle acceleration. We present solutions for plasma quantities and particle distributions upstream and downstream of shocks, illustrating the strong differences observed between non-linear and test-particle cases. It is found that there are only marginal differences in the injection efficiency and resultant spectra for two extreme scattering modes, namely large-angle scattering and pitch-angle diffusion, for a wide range of shock parameters, i.e., for subluminal shocks with field obliquities less than or equal to 75 degrees and de Hoffmann-Teller frame speeds much less than the speed of light.

The rate at which particles are accelerated by the first-order Fermi mechanism in shocks depends on the angle, \teq{\Tbone}, that the upstream magnetic field makes with the shock normal. The greater the obliquity the greater the rate, and in quasi-perpendicular shocks rates can be hundreds of times higher than those seen in parallel shocks. In many circumstances pertaining to evolving shocks (\eg, supernova blast waves and interplanetary traveling shocks), high acceleration rates imply high maximum particle energies and obliquity effects may have important astrophysical consequences. However, as is demonstrated here, the efficiency for injecting thermal particles into the acceleration mechanism also depends strongly on obliquity and, in general, varies inversely with \teq{\Tbone}. The degree of turbulence and the resulting cross-field diffusion strongly influences both injection efficiency and acceleration rates. The test particle \mc simulation of shock acceleration used here assumes large-angle scattering, computes particle orbits exactly in shocked, laminar, non-relativistic flows, and calculates the injection efficiency as a function of obliquity, Mach number, and degree of turbulence. We find that turbulence must be quite strong for high Mach number, highly oblique shocks to inject significant numbers of thermal particles and that only modest gains in acceleration rates can be expected for strong oblique shocks over parallel ones if the only source of seed particles is the thermal background.