We consider the acceleration of charged particles near ultra-relativistic shocks, with Lorentz factor Gamma_s >> 1. We present simulations of the acceleration process and compare these with results from semi-analytical calculations. We show that the spectrum that results from acceleration near ultra-relativistic shocks is a power law, N(E) \propto E^{-s}, with a nearly universal value s \approx 2.2 - 2.3 for the slope of this power law. We confirm that the ultra-relativistic equivalent of Fermi acceleration at a shock differs from its non-relativistic counterpart by the occurence of large anisotropies in the distribution of the accelerated particles near the shock. In the rest frame of the upstream fluid, particles can only outrun the shock when their direction of motion lies within a small loss cone of opening angle theta_c \approx 1/Gamma_s around the shock normal. We also show that all physically plausible deflection or scattering mechanisms can change the upstream flight direction of relativistic particles originating from downstream by only a small amount: Delta theta ~ 1/Gamma_s. This limits the energy change per shock crossing cycle to Delta E ~ E, except for the first cycle where particles originate upstream. In that case the upstream energy is boosted by a factor ~ Gamma_s^2 for those particles that are scattered back across the shock into the upstream region.

We extend the eigenfunction method of computing the power-law spectrum of particles accelerated at a relativistic shock fronts to apply to shocks of arbitrarily high Lorentz factor. In agreement with the findings of Monte-Carlo simulations, we find the index of the power-law distribution of accelerated particles which undergo isotropic diffusion in angle at an ultrarelativistic, unmagnetized shock is s=4.23 (where s=-d(ln f)/dp with f the Lorentz invariant phase-space density and p the momentum). This corresponds to a synchrotron index for uncooled electrons of a=0.62 (taking cooling into account a=1.12), where a=-d(ln F)/dn, F is the radiation flux and n the frequency. We also present an approximate analytic expression for the angular distribution of accelerated particles, which displays the effect of particle trapping by the shock: compared with the non-relativistic case the angular distribution is weighted more towards the plane of the shock and away from its normal. We investigate the sensitivity of our results to the transport properties of the particles and the presence of a magnetic field. Shocks in which the ratio of Poynting to kinetic energy flux upstream is not small are less compressive and lead to larger values of $s$.

We adapt and modify the eigenfunction method of computing the power-law spectrum of particles accelerated at a relativistic shock front via the first-order Fermi process (Kirk, J.G., Schneider, P., Astrophysical Journal 315, 425 (1987)) to apply to shocks of arbitrarily high Lorentz factor. The power-law index of accelerated particles undergoing isotropic small-angle scattering at an ultrarelativistic, unmagnetized shock is found to be s=4.23 +/- 0.2 (where s=d\ln f/ d\ln p, with f the Lorentz-invariant phase-space density and p the momentum), in agreement with the results of Monte-Carlo simulations. We present results for shocks in plasmas with different equations of state and for Lorentz factors ranging from 5 to infinity.

We examine Fermi-type acceleration at relativistic shocks, and distinguish between the initial boost of the first shock crossing cycle, where the energy gain per particle can be very large, and the Fermi process proper with repeated shock crossings, in which the typical energy gain is of order unity. We calculate by means of numerical simulations the spectrum and angular distribution of particles accelerated by this Fermi process, in particular in the case where particle dynamics can be approximated as small-angle scattering. We show that synchrotron emission from electrons or positrons accelerated by this process can account remarkably well for the observed power-law spectra of GRB afterglows and Crab-like supernova remnants. In the context of a decelerating relativistic fireball, we calculate the maximum particle energy attainable by acceleration at the external blast wave, and discuss the minimum energy for this acceleration process and its consequences for the observed spectrum.