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.

We examine the problem of particle acceleration at a relativistic shocks assuming pitch-angle scattering and using a Hartree-Fock method to approximate the associated eigenfunctions. This leads to a simple transcendental equation determining the power-law index, $s$, given the up and downstream velocities. We compare our results with accurate numerical solutions obtained using the eigenfunction method. In addition to the power-law index this method yields the angular and spatial distributions upstream of the shock.

We present both numerical and semi-analytical results on test-particle acceleration in multiple parallel shocks. We apply a kinetic Monte Carlo code and an eigenfunction expansion method to calculate the distribution functions for electron populations accelerated in subsequent parallel shocks with speeds ranging from non- to fully-relativistic. We examine the levels of particle anisotropy at the shocks and discuss the implications for AGN and microquasar jets.

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.