We examine limits to the efficiency for particles acceleration in shearing flows, showing that relativistic flow speeds are required for efficient gradual shear acceleration. We estimate maximum achievable particle energies for parameters applicable to relativistic AGN jets. The implications of our estimates is that if large-scale jets are relativistic, then efficient electron acceleration up to several PeV, and proton acceleration up to several EeV energies appears feasible. This suggests that shear particle acceleration could lead to a continued energization of synchrotron X-ray emitting electrons, and be of relevance for the production of ultra-high-energy cosmic-ray particles.

The Bubble Nebula (or NGC 7635) is a parsec-scale seemingly spherical wind-blown bubble around the relatively unevolved O star BD+60$^\circ$2522. The small dynamical age of the nebula and significant space velocity of the star suggest that the Bubble Nebula might be a bow shock. We have run 2D hydrodynamic simulations to model the interaction of the central star's wind with the interstellar medium (ISM). The models cover a range of possible ISM number densities of $n=50-200 {\rm cm}^{-3}$ and stellar velocities of $v_{\star}=20-40$ km s$^{-1}$. Synthetic H$\alpha$ and 24 $\mu$m emission maps predict the same apparent spherical bubble shape with quantitative properties similar to observations. The synthetic maps also predict a maximum brightness similar to that from the observations and agree that the maximum brightness is at the apex of the bow shock. The best-matching simulation had $v_{\star}\approx20$ km s$^{-1}$ into an ISM with $n\sim100 {\rm cm}^{-3}$, at an angle of 60$^\circ$ with respect to the line of sight. Synthetic maps of soft ($0.3-2$ keV) and hard ($2-10$ keV) X-ray emission show that the brightest region is in the wake behind the star and not at the bow shock itself. The unabsorbed soft X-rays have luminosity $\sim10^{32}-10^{33}$ erg s$^{-1}$. The hard X-rays are fainter, luminosity $\sim 10^{30} - 10^{31}$ erg s$^{-1}$, and may be too faint for current X-ray instruments to successfully observe. Our results imply that the O star creates a bow shock as it moves through the ISM and in turn creates an asymmetric bubble visible at optical and infrared wavelengths, and predicted to be visible in X-rays. The Bubble Nebula does not appear to be unique, it could be just a favourably oriented very dense bow shock. The dense ISM surrounding BD+60$^\circ$2522 and its strong wind suggest that it could be a good candidate for detecting non-thermal emission.

Shear flows are naturally expected to occur in astrophysical environments and potential sites of continuous non-thermal Fermi-type particle acceleration. Here we investigate the efficiency of expanding relativistic outflows to facilitate the acceleration of energetic charged particles to higher energies. To this end, the gradual shear acceleration coefficient is derived based on an analytical treatment. The results are applied to the context of the relativistic jets of active galactic nuclei. The inferred acceleration timescale is investigated for a variety of conical flow profiles (i.e., power law, Gaussian, Fermi-Dirac) and compared to the relevant radiative and non-radiative loss timescales. The results exemplify that relativistic shear flows are capable of boosting cosmic-rays to extreme energies. Efficient electron acceleration, on the other hand, requires weak magnetic fields and may thus be accompanied by a delayed onset of particle energization and affect the overall jet appearance (e.g., core, ridge line and limb-brightening).

Stochastic acceleration of particles under a pressure balance condition can accommodate the universal $p^{-5}$ spectra observed under many different conditions in the inner heliosphere. In this model, in order to avoid an infinite build up of particle pressure, a relationship between the momentum diffusion of particles and the adiabatic deceleration in the solar wind must exist. This constrains both the spatial and momentum diffusion coefficients and results in the $p^{-5}$ spectrum in the presence of adiabatic losses in the solar wind. However, this theory cannot explain the presence of such spectra beyond the termination shock, where adiabatic deceleration is negligible. To explain this apparent discrepancy, we include the effect of charge exchange losses, resulting in new forms of both the spatial and momentum diffusion coefficients that have not previously been considered. Assuming that the turbulence is of a large-scale compressible nature, we find that a balance between momentum diffusion and losses can still readily lead to the creation of $p^{-5}$ suprathermal tails, including those found in the outer heliosphere.

We present a new approach to calculate the particle distribution function about relativistic shocks including synchrotron losses using the method of lines with an explicit finite difference scheme. A steady, continuous, one dimensional plasma flow is considered to model thick (modified) shocks, leading to a calculation in three dimensions plus time, the former three being momentum, pitch angle and position. The method accurately reproduces the expected power law behaviour in momentum at the shock for upstream flow speeds ranging from 0.1c to 0.995c (1 < \Gamma < 10). It also reproduces approximate analytical results for the synchrotron cutoff shape for a non-relativistic shock, demonstrating that the loss process is accurately represented. The algorithm has been implemented as a hybrid OpenMP--MPI parallel algorithm to make efficient use of SMP cluster architectures and scales well up to many hundreds of CPUs.

The existing theoretical framework for the energies stored in the synchrotron-emitting lobes of radio galaxies and quasars doesn't properly account for the curved spectral shape that many of them exhibit. We characterise these spectra using parameters that are straightforwardly observable in the era of high-resolution, low-frequency radio astronomy: the spectral curvature and the turnover in the frequency spectrum. This characterisation gives the Lorentz factor at the turnover in the energy distribution (we point out that this is distinctly different from the Lorentz factor corresponding to the turnover frequency in a way that depends on the amount of curvature in the spectrum) and readily gives the equipartition magnetic field strength and the total energy of the radiating plasma obviating the need for any assumed values of the cutoff frequencies to calculate these important physical quantities. This framework readily yields the form of the X-ray emission due to inverse-Compton (IC) scattering of Cosmic Microwave Background (CMB) photons by the electrons in the plasma having Lorentz factors of $\sim$1000. We also present the contribution to CMB anisotropies due to relativistic plasmas such as giant radio galaxy lobes, expressed in terms of the extent to which the lobes have their magnetic field and particle energies are in equipartition with one another.

Stationary solutions to the equations of non-linear diffusive shock acceleration play a fundamental role in the theory of cosmic-ray acceleration. Their existence usually requires that a fraction of the accelerated particles be allowed to escape from the system. Because the scattering mean-free-path is thought to be an increasing function of energy, this condition is conventionally implemented as an upper cut-off in energy space -- particles are then permitted to escape from any part of the system, once their energy exceeds this limit. However, because accelerated particles are responsible for substantial amplification of the ambient magnetic field in a region upstream of the shock front, we examine an alternative approach in which particles escape over a spatial boundary. We use a simple iterative scheme that constructs stationary numerical solutions to the coupled kinetic and hydrodynamic equations. For parameters appropriate for supernova remnants, we find stationary solutions with efficient acceleration when the escape boundary is placed at the point where growth and advection of strongly driven non-resonant waves are in balance. We also present the energy dependence of the distribution function close to the energy where it cuts off - a diagnostic that is in principle accessible to observation.

Observations of minute-scale flares in TeV Blazars place constraints on particle acceleration mechanisms in those objects. The implications for a variety of radiation mechanisms have been addressed in the literature; in this paper we compare four different acceleration mechanisms: diffusive shock acceleration, second-order Fermi, shear acceleration and the converter mechanism. When the acceleration timescales and radiative losses are taken into account, we can exclude shear acceleration and the neutron-based converted mechanism as possible acceleration processes in these systems. The first-order Fermi process and the converter mechanism working via SSC photons are still practically instantaneous, however, provided sufficient turbulence is generated on the timescale of seconds. We propose stochastic acceleration as a promising candidate for the energy-dependent time delays in recent gamma-ray flares of Markarian 501.

The steady-state structure of a disc with a corona is analyzed when the vertical component of the gravitational force due to the self-gravity of the disc is considered. For the energy exchange between the disc and the corona, we assume a fraction f of the dissipated energy inside the accretion disc is transported to the corona via the magnetic tubes. Analytical solutions corresponding to a prescription for f (in which this parameter directly depends on the ratio of the gas pressure to the total pressure) or free f are presented and their physical properties are studied in detail. We show that the existence of the corona not only decreases the temperature of the disc, but also increases the surface density.The vertical component of the gravitational force due to the self-gravity of the disc decreases the self-gravitating radius and the mass of the fragments at this radius. However, as more energy is transported from the disc to the corona, the effect of the vertical component of the gravitational force due to the self-gravity of the disc on the self-gravitating radius becomes weaker, though the mass of the fragments is reduced irrespective of the amount of the energy exchange from the disc to the corona.

The duty-cycle of powerful radio galaxies and quasars such as the prototype Cygnus A is poorly understood. X-ray observations of inverse-Compton scattered Cosmic Microwave Background (ICCMB) photons probe lower Lorentz-factor particles than radio observations of synchrotron emission. Comparative studies of the nearer and further lobes, separated by many 10s of kpc and thus 10s of thousands of years in light-travel time, yield additional temporal resolution in studies of the lifecycles. We have co-added all archival Chandra ACIS-I data and present a deep 200 ks image of Cygnus A. This deep image reveals the presence of X-ray emission from a counterjet i.e. a jet receding from Earth and related to a previous episode of jet activity. The non-thermal X-ray emission, we interpret as ICCMB radiation. There is an absence of any discernible X-ray emission associated with a jet flowing towards Earth. We conclude that: (1) The emission from a relic jet, indicates a previous episode of jet activity, that took place earlier than the current jet activity appearing as synchrotron radio emission. (2) The presence of X-ray emission from a relic counterjet of Cygnus A and the absence of X-ray emission associated with any relic approaching jet constrains the timescale between successive episodes of jet activity to ~10^6 years. (3) Transverse expansion of the jet causes expansion losses which shifts the energy distribution to lower energies. (4) Assuming the electrons cooled due to adiabatic expansion, the required magnetic field strength is substantially smaller than the equipartition magnetic field strength. (5) A high minimum Lorentz factor for the distribution of relativistic particles in the current jet, of a few 10^3, is ejected from the central nucleus of this active galaxy. Abridged.

We investigate the so-called non-resonant 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 non-resonant mode in SN1006 and RX J1713.7-3946 are presented.

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.

We present a 200 ks Chandra ACIS-I image of Cygnus A, and discuss a long linear feature seen in its counterlobe. This feature has a non-thermal spectrum and lies on the line connecting the brighter hotspot on the approaching side and the nucleus. We therefore conclude that this feature is (or was) a jet. However, the outer part of this X-ray jet does not trace the current counterjet observed in radio. No X-ray counterpart is observed on the jet side. Using light-travel time effects we conclude that this X-ray 50 kpc linear feature is a relic jet that contains enough low-energy plasma (gamma ~ 10^3) to inverse-Compton scatter cosmic microwave background photons, producing emission in the X-rays.

We investigate the growth of hydromagnetic waves driven by streaming cosmic rays in the precursor environment of a supernova remnant shock. 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. 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. 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. 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.

Relativistic shocks provide an efficient method for high-energy particle acceleration in many astrophysical sources. Multiple shock systems are even more effective and of importance, for example, in the internal shock model of gamma-ray bursts. We investigate the reacceleration of pre-existing energetic particles at such relativistic internal shocks by the first order Fermi process of pitch angle scattering. We use a well established eigenfunction method to calculate the resulting spectra for infinitely thin shocks. Implications for GRBs and relativistic jets are discussed.

We investigate the acceleration and simultaneous radiative losses of electrons in the vicinity of relativistic shocks. Particles undergo pitch angle diffusion, gaining energy as they cross the shock by the Fermi mechanism and also emitting synchrotron radiation in the ambient magnetic field. A semi-analytic approach is developed which allows us to consider the behaviour of the shape of the spectral cut-off and the variation of that cut-off with the particle pitch angle. The implications for the synchrotron emission of relativistic jets, such as those in gamma ray burst sources and blazars, are discussed.

A microscopic analysis of the viscous energy gain of energetic particles in (gradual) non-relativistic shear flows is presented. We extend previous work and derive the Fokker-Planck coefficients for the average rate of momentum change and dispersion in the general case of a momentum-dependent scattering time $\tau(p) \propto p^{\alpha}$ with $\alpha \geq 0$. We show that in contrast to diffusive shock acceleration the characteristic shear acceleration timescale depends inversely on the particle mean free path which makes the mechanism particularly attractive for high energy seed particles. Based on an analysis of the associated Fokker-Planck equation we show that above the injection momentum $p_0$ power-law differential particle number density spectra $n(p) \propto p^{-(1+ \alpha)}$ are generated for $\alpha >0$ if radiative energy losses are negligible. We discuss the modifications introduced by synchrotron losses and determine the contribution of the accelerated particles to the viscosity of the background flow. Possible implications for the plasma composition in mildly relativistic extragalactic jet sources (WATs) are addressed.

We consider the acceleration of energetic particles by Fermi processes (i.e., diffusive shock acceleration, second order Fermi acceleration, and gradual shear acceleration) in relativistic astrophysical jets, with particular attention given to recent progress in the field of viscous shear acceleration. We analyze the associated acceleration timescales and the resulting particle distributions, and discuss the relevance of these processes for the acceleration of charged particles in the jets of AGNs, GRBs and microquasars, showing that multi-component powerlaw-type particle distributions are likely to occur.

Recently, Bell has reanalysed the problem of wave excitation by cosmic rays propagating in the pre-cursor region of a supernova remnant shock front. He pointed out a strong, non-resonant, current-driven instability that had been overlooked in the kinetic treatments, and suggested that it is responsible for substantial amplification of the ambient magnetic field. Magnetic field amplification is also an important issue in the problem of the formation and structure of relativistic shock fronts, particularly in relation to models of gamma-ray bursts. We have therefore generalised the linear analysis to apply to this case, assuming a relativistic background plasma and a monoenergetic, unidirectional incoming proton beam. We find essentially the same non-resonant instability noticed by Bell, and show that also under GRB conditions, it grows much faster than the resonant waves. We quantify the extent to which thermal effects in the background plasma limit the maximum growth rate.

Gradual shear acceleration of energetic particles in gamma-ray burst (GRB) jets is considered. Special emphasis is given to the analysis of universal structured jets, and characteristic acceleration timescales are determined for a power-law and a Gaussian evolution of the bulk flow Lorentz factor $\gamma_b$ with angle $\phi$ from the jet axis. The results suggest that local power-law particle distributions may be generated and that higher energy particles are generally concentrated closer to the jet axis. Taking several constraints into account we show that efficient electron acceleration in gradual shear flows, with maximum particle energy successively decreasing with time, may be possible on scales larger than $r \sim 10^{15}$ cm, provided the jet magnetic field becomes sufficiently weak and/or decreases rapidly enough with distance, while efficient acceleration of protons to ultra-high energies $> 10^{20}$ eV may be possible under a wide range of conditions.

We review the acceleration of energetic particles in relativistic astrophysical jets characterized by a significant velocity shear. The possible formation of power-law momentum spectra is discussed and typical acceleration timescales are determined for a variety of different conditions such as parallel and azimuthal shear flows. Special attendance is given to the analysis of parallel shear flows with either a linear decreasing or a Gaussian-type velocity profile. It is shown that in the presence of a gradual shear flow and a particle mean free path scaling with the gyroradius, synchrotron radiation losses may no longer be able to stop the acceleration once it has started to work efficiently. Finally, the relevance of shear acceleration in small- and large-scale relativistic jets is addressed.

The acceleration of energetic particles in astrophysical shear flows is analyzed. We show that in the presence of a non-relativistic gradual velocity shear, power law particle momentum distributions $f(p) \propto p^{-(3+\alpha)}$ may be generated, assuming a momentum-dependent scattering time $\tau \propto p^{\alpha}$, with $\alpha > 0$. We consider possible acceleration sites in astrophysical jets and study the conditions for efficient acceleration. It is shown, for example, that in the presence of a gradual shear flow and a gyro-dependent particle mean free path, synchrotron radiation losses no longer stop the acceleration once it has started to work efficiently. This suggests that shear acceleration may naturally account for a second, non-thermal population of energetic particles in addition to a shock-accelerated one. The possible relevance of shear acceleration is briefly discussed with reference to the relativistic jet in the quasar 3C 273.

We consider the acceleration of energetic particles by a velocity shear in the relativistic background flow containing scattering centers. Three possible acceleration sites for astrophysical jets are identified: (1) gradual velocity shear parallel to the jet axis, such as a velocity profile decreasing linearly outward with radial coordinate, (2) gradual velocity shear perpendicular to the jet axis, such as intrinsic jet rotation, and (3) non-gradual/discontinuous, longitudinal velocity shear at the jet-side boundary. We determine the characteristic acceleration timescales, specify the conditions for efficient acceleration and discuss observational features with respect to each process. In particular, it is shown that in the case of (2) the higher energy emission is expected to be concentrated closer to the jet axis, while in the case of (1) and (3) the higher energy particles are likely to be located near the edges of the jet, thus possibly leading to some form of limb-brightening.

We consider the acceleration of particles due to a velocity shear in relativistic astrophysical flows. The basic physical picture and the formation of power law momentum spectra is discussed for a non-relativistic velocity field using a microscopic approach. We identify possible sites for shear acceleration in relativistic astrophysical jets and analyze their associated acceleration timescales. It is shown in particular that for a mean scattering time $\tau$ scaling with the gyro-radius, the acceleration timescale for gradual shear scales in the same manner as the synchrotron cooling timescale, so that losses may no longer be able to stop the acceleration once it has started to work efficiently. Finally, the possible role of shear acceleration is discussed with reference to the relativistic jet in the quasar 3C~273.

Studies of the knotty jets in both quasars and microquasars frequently make use of the ratio of the intensities of corresponding knots on opposite sides of the nucleus in order to infer the product of the intrinsic jet speed (beta) and the cosine of the inclination angle of the jet-axis (cos{theta}), via the formalism I_{a}/I_{r} = ((1+beta cos{theta})/(1-beta cos{theta}))^{3+alpha}, where alpha relates the intensity I_{nu} as a function of frequency nu as I_{nu} propto nu^{-alpha}. Where beta cos{theta} is determined independently, the intensity ratio of a given pair of jet to counter-jet knots is over-predicted by the above formalism compared with the intensity ratio actually measured from radio images. As an example in the case of Cygnus X-3 the original formalism predicts an intensity ratio of about 185, whereas the observed intensity ratio at one single epoch is about 3. Mirabel and Rodriguez (1999) have refined the original formalism, and suggested measuring the intensity ratio of knots when they are at equal angular separations from the nucleus. This method is only applicable where there is sufficient time-sampling with sufficient physical resolution to interpolate the intensities of the knots at equal distances from the nucleus, and can therefore be difficult to apply to microquasars and is impossible to apply to quasars. Accounting for both the light-travel time between the knots and the simple evolution of the knots themselves reconciles this over-prediction and renders the original formalism obsolete.

Multifrequency observations from the VLA, VLBA and OVRO Millimeter Array of a major radio outburst of Cygnus X-3 in 2001 September are presented, measuring the evolution of the spectrum of the source over three decades in frequency, over a period of six days. Following the peak of the flare, as the intensity declines the high-frequency spectrum at frequency nu steepens from nu^{-0.4} to nu^{-0.6}, after which the spectral index remains at this latter terminal value; a trend previously observed but hitherto not satisfactorily explained. VLBA observations, for the first time, track over several days the expansion of a sequence of knots whose initial diameters are approximately 8 milliarcseconds. The light-crossing time within these plasmons is of the same order as the time-scale over which the spectrum is observed to evolve. We contend that properly accounting for light-travel time effects in and between plasmons which are initially optically thick, but which after expansion become optically thin, explains the key features of the spectral evolution, for example the observed timescale. Using the VLBA images, we have directly measured for the first time the proper motions of individual knots, analysis of which shows a two-sided jet whose axis is precessing. The best-fit jet speed is roughly beta = 0.63 and the precession period is about 5 days, significantly lower than fitted for a previous flare. Extrapolation of the positions of the knots measured by the VLBA back to zero-separation shows this to occur approximately 2.5 days after the detection of the rise in flux density of Cygnus X-3.

Relativistic shocks accelerate particles by the first order Fermi mechanism. These particles then emit synchrotron emission in the post shock gas. We have developed a numerical code which integrates the relativistic Euler equations for fluid dynamics with a general equation of state, together with the Liouville equation for the accelerated particles. We present tests of this code and, in addition, we use it to study the gamma ray burst afterglow predicted by the fireball model, along with the hydrodynamics of a relativistic blastwave. We find that, while, broadly speaking, the behaviour of the emission is similar to that already predicted with semi-analytic approaches, the detailed behaviour is somewhat different. The ``breaks'' in the synchrotron spectrum behave differently with time, and the spectrum above the final break is harder than previously expected. These effects are due to the incorporation of the geometry of the (spherical) blastwave, along with relativistic beaming and adiabatic cooling of the energetic particles leading to a mix, in the observed spectrum, between recently injected "uncooled" particles and the older "cooled" population in different parts of the evolving, inhomogeneous flow.

The recent detection of high energy gamma-rays coming from supernova remnants and active galactic nuclei has revived interest in the diffusive shock acceleration of electrons. In the present paper we examine the basis of the so-called ``box'' model for particle acceleration and present a more physical version of it. Using this we determine simple criteria for the conditions under which ``pile-ups'' can occur in shock accelerated electron spectra subject to synchrotron or inverse Compton losses (the latter in the Thompson limit). An extension to include nonlinear effects is proposed.

Observations of both gamma-ray burst sources and certain classes of active galaxy indicate the presence of relativistic shock waves and require the production of high energy particles to explain their emission. In this paper we review the basic theory of shock waves in relativistic hydrodynamics and magneto-hydrodynamics, emphasising the astrophysically interesting cases. This is followed by an overview of the theory of particle acceleration at such shocks. We summarise the applications to the astrophysics of relativistic jets and fireball models of gamma-ray-bursts.

The theory of diffusive acceleration of energetic particles at shock fronts assumes charged particles undergo spatial diffusion in a uniform magnetic field. If, however, the magnetic field is not uniform, but has a stochastic or braided structure, the transport of charged particles across the average direction of the field is more complicated. Assuming quasi-linear behaviour of the field lines, the particles undergo sub-diffusion on short time scales. We derive the propagator for such motion, which differs from the Gaussian form relevant for diffusion, and apply it to a configuration with a plane shock front whose normal is perpendicular to the average field direction. Expressions are given for the acceleration time as a function of the diffusion coefficient of the wandering magnetic field lines and the spatial diffusion coefficient of the charged particles parallel to the local field. In addition we calculate the spatial dependence of the particle density in both the upstream and downstream plasmas. In contrast to the diffusive case, the density of particles at the shock front is lower than it is far downstream. This is a consequence of the partial trapping of particles by structures in the magnetic field. As a result, the spectrum of accelerated particles is a power-law in momentum which is steeper than in the diffusive case. For a phase- space density $f\propto p^{-s}$, we find $s=\sdiff[1+1/(2\rcomp)]$, where $\rcomp$ is the compression ratio of the shock front and $\sdiff$ is the standard result of diffusive acceleration: $\sdiff=3\rcomp/(\rcomp-1)$. A strong shock in a monatomic ideal gas yields a spectrum of $s=4.5$. In the case of electrons, this corresponds to a radio synchrotron spectral index of $\alpha=0.75$.

The limits imposed on diffusive shock acceleration by upstream ion-neutral Alfven wave damping, and by ionisation and Coulomb losses of low energy particles, are calculated. Analytic solutions are given for the steady upstream wave excitation problem with ion-neutral damping and the resulting escaping upstream flux calculated. The time dependent problem is discussed and numerical solutions presented. Finally the significance of these results for possible observational tests of shock acceleration in supernova remnants is discussed.

The theory of first order Fermi acceleration at shocks assumes that particles diffuse due to scattering off slow-moving magnetic irregularities. However, cosmic rays are closely tied to magnetic field lines, and the transport process, particularly across the direction of the field, is likely to be more complicated. To describe cross field transport we employ recent extensions of the Rechester-Rosenbluth theory in localised stochastic regions of magnetic field. During acceleration at a shock and when the motion along the field is diffusive, there is a transition at a critical energy from \lq\lq sub-diffusive\rq\rq\ motion, where the mean square displacement of a particle increases with time as $t^{1/2}$, to compound diffusion, a combined process involving diffusion along a magnetic field which is itself wandering. Requiring this critical energy to be less than the system cut-off in a SNR of radius $R$ places an upper limit on the coherence length, $\lambda$, of the magnetic field for which diffusive shock acceleration can occur at quasi-perpendicular shocks when field line wandering dominates the effective transport of particles: $\lambda<R(U/v)^{1/2}$ where $U$ is the shock speed and $v$ the particle speed. The acceleration rate for particles in the sub-diffusive regime is derived and its implication for the origin of high energy cosmic rays is discussed.

The observations of radio emission from SNR~1987A can be accounted for on the basis of diffusive shock acceleration of electrons by the supernova blast wave. However, with this interpretation the observed spectral index implies that the compression ratio of the gas subshock is roughly $2.7$ rather than the value of $4$ expected of a strong shock front. We propose that in SNR~1987A, ions also undergo diffusive acceleration at the shock, a process that is likely to be rapid. Unlike the electron population, the accelerated ions can have an important effect on the gas dynamics. We calculate this coupled gas and energetic particle dynamics on the basis of the two-fluid model, in which the accelerated ions provide an additional component to the total pressure acting on the fluid. By accelerating and possibly heating the upstream plasma, the initially strong shock is modified and a weaker subshock with an upstream precursor results. The electrons behave as test particles. They are accelerated at the evolving subshock, escape downstream, and emit synchrotron radiation in the swept up magnetic field. Two models are considered for the surroundings of the progenitor: that of a freely expanding wind of number density $n\propto r^{-2}$, and that of a wind confined by a shell of denser material, creating a stagnation zone of roughly constant density beyond the standing shock which terminates the free wind. We model the observed radio light curves and the relatively steep spectrum of SNR~1987A using similar values for the ion acceleration parameters to those used in models of cosmic ray acceleration in older SNRs which can also contain high Mach number shocks, and find a good fit for the case in which the termination shock is located at about $2\times 10^{15}\,$m from the progenitor.

When applied to the blast wave formed by the explosion of a massive star as a supernova (SN), the theory of diffusive particle acceleration at shock fronts predicts a very high energy density in cosmic rays. Almost immediately after particles begin to be injected into the process, the cosmic ray pressure rises until comparable to the ram-pressure encountered by the shock front. Those supernovae which are observed in the radio band i.e., radio supernovae (RSNe), provide direct evidence of particle acceleration in the form of synchrotron emitting electrons. Furthermore, these objects are particularly interesting, since they are usually surrounded by a relatively dense confining medium. The acceleration of cosmic rays can then lead to the production of very high energy (VHE) gamma-rays which arise from collisions between energetic particles and target nuclei. We estimate the cosmic ray energy density assuming a fraction $\phi\lesim1$ of the energy available at the shock front is converted into cosmic rays. Combining this with the parameters describing the environment of the SN progenitor, as deduced from observations, and from more detailed modelling, we compute the flux at Earth $\tevflux$ of photons of energy above $1\,$TeV. For the relatively weak but nearby supernova SN1987A we predict $\tevflux=2\times10^{-13}\,\fluxunits$ before the shock front encounters the ring of dense matter seen by the Hubble Space Telescope. Subsequently, the flux is expected to rise further. The medium around SN1993J in M81 is thought to have a density

When a magnetised, thermal plasma containing a shock wave efficiently accelerates Cosmic Rays (CR) a reaction is exerted on the fluid structure. The simplest description of this process involves three conservation equations for the background fluid, with the CR pressure gradient included as an extra force, plus a hydrodynamic equation for the CR energy density. Since the CR spectrum is not determined in this approach two quantities are introduced, the ratio of CR pressure to energy density and a mean diffusion coefficient, that must be specified to close the system of equations. Physical arguments are presented for how these quantities should evolve and then compared with time dependent numerical solutions of acceleration at a plane parallel, piston driven shock.