Diffusive shock acceleration (DSA) at relativistic shocks is widely thought to be an important acceleration mechanism in various astrophysical jet sources, including radio-loud active galactic nuclei such as blazars. Such acceleration can produce the non-thermal particles that emit the broadband continuum radiation that is detected from extragalactic jets. An important recent development for blazar science is the ability of Fermi-LAT spectroscopy to pin down the shape of the distribution of the underlying non-thermal particle population. This paper highlights how multi-wavelength spectra spanning optical to X-ray to gamma-ray bands can be used to probe diffusive acceleration in relativistic, oblique, magnetohydrodynamic (MHD) shocks in blazar jets. Diagnostics on the MHD turbulence near such shocks are obtained using thermal and non-thermal particle distributions resulting from detailed Monte Carlo simulations of DSA. These probes are afforded by the characteristic property that the synchrotron $\nu F_{\nu}$ peak energy does not appear in the gamma-ray band above 100 MeV. We investigate self-consistently the radiative synchrotron and inverse Compton signatures of the simulated particle distributions. Important constraints on the diffusive mean free paths of electrons, and the level of electromagnetic field turbulence are identified for three different case study blazars, Mrk 501, BL Lacertae and AO 0235+164. The X-ray excess of AO 0235+164 in a flare state can be modelled as the signature of bulk Compton scattering of external radiation fields, thereby tightly constraining the energy-dependence of the diffusion coefficient for electrons. The concomitant interpretations that turbulence levels decline with remoteness from jet shocks, and the probable significant role for non-gyroresonant diffusion, are posited.

Diffusive shock acceleration (DSA) at relativistic shocks is likely to be an important acceleration mechanism in various astrophysical jet sources, including radio-loud AGN. An important recent development for blazar science is the ability of Fermi-LAT data to pin down the power-law index of the high energy portion of emission in these sources, and therefore also the index of the underlying non-thermal particle population. This diagnostic potential was not possible prior to Fermi launch, when gamma-ray information was dominated by the highly-absorbed TeV band. This paper highlights how multiwavelength spectra including X-ray band and Fermi data can be used to probe diffusive acceleration in relativistic, oblique, MHD shocks in blazar jets. The spectral index of the non-thermal particle distributions resulting from Monte Carlo simulations of DSA, and the fraction of thermal particles accelerated to non-thermal energies, depend sensitively on the particles' mean free path scale, and also on the magnetic field obliquity to the shock normal. We investigate self-consistently the radiative synchrotron/Compton signatures of the resulting thermal and non-thermal particle distributions. Important constraints on the frequency of particle scattering and the level of field turbulence are identified for the blazar AO 0235+164. The possible interpretation that turbulence levels decline with remoteness from jet shocks, and a significant role for non-gyroresonant diffusion, are discussed.

Diffusive shock acceleration (DSA) at relativistic shocks is expected to be an important acceleration mechanism in a variety of astrophysical objects including extragalactic jets in active galactic nuclei and gamma ray bursts. These sources remain good candidate sites for the generation of ultra-high energy cosmic rays. In this paper, key predictions of DSA at relativistic shocks that are germane to production of relativistic electrons and ions are outlined. The technique employed to identify these characteristics is a Monte Carlo simulation of such diffusive acceleration in test-particle, relativistic, oblique, magnetohydrodynamic (MHD) shocks. Using a compact prescription for diffusion of charges in MHD turbulence, this approach generates particle angular and momentum distributions at any position upstream or downstream of the shock. Simulation output is presented for both small angle and large angle scattering scenarios, and a variety of shock obliquities including superluminal regimes when the de Hoffmann-Teller frame does not exist. The distribution function power-law indices compare favorably with results from other techniques. They are found to depend sensitively on the mean magnetic field orientation in the shock, and the nature of MHD turbulence that propagates along fields in shock environs. An interesting regime of flat spectrum generation is addressed; we provide evidence for it being due to shock drift acceleration, a phenomenon well-known in heliospheric shock studies. The impact of these theoretical results on blazar science is outlined. Specifically, Fermi-LAT gamma-ray observations of these relativistic jet sources are providing significant constraints on important environmental quantities for relativistic shocks, namely the field obliquity, the frequency of scattering and the level of field turbulence.

Diffusive shock acceleration (DSA) at relativistic shocks is expected to be an important acceleration mechanism in a variety of astrophysical objects including extragalactic jets in active galactic nuclei and gamma ray bursts. These sources remain strong and interesting candidate sites for the generation of ultra-high energy cosmic rays. In this paper, key predictions of DSA at relativistic shocks that are salient to the issue of cosmic ray ion and electron production are outlined. Results from a Monte Carlo simulation of such diffusive acceleration in test-particle, relativistic, oblique, MHD shocks are presented. Simulation output is described for both large angle and small angle scattering scenarios, and a variety of shock obliquities including superluminal regimes when the de Hoffman-Teller frame does not exist. The distribution function power-law indices compare favorably with results from other techniques. They are found to depend sensitively on the mean magnetic field orientation in the shock, and the nature of MHD turbulence that propagates along fields in shock environs. An interesting regime of flat spectrum generation is addressed, providing evidence for its origin being due to shock drift acceleration. The impact of these theoretical results on gamma-ray burst and blazar science is outlined. Specifically, Fermi gamma-ray observations of these cosmic sources are already providing significant constraints on important environmental quantities for relativistic shocks, namely the frequency of scattering and the level of field turbulence.

The acceleration of interstellar pick-up ions as well as solar wind species has been observed at a multitude of interplanetary (IP) shocks by different spacecraft. This paper expands upon previous work modeling the phase space distributions of accelerated ions associated with the shock event encountered on day 292 of 1991 by the Ulysses mission at 4.5 AU. A kinetic Monte Carlo simulation is employed here to model the diffusive acceleration process. This exposition presents recent developments pertaining to the incorporation into the simulation of the diffusive characteristics incurred by field line wandering (FLW), according to the work of Giacalone and Jokipii. For a pure field-line wandering construct, it is determined that the upstream spatial ramp scales are too short to accommodate the HI-SCALE flux increases for 200 keV protons, and that the distribution function for H+ somewhat underpopulates the combined SWICS/HI-SCALE spectra at the shock. This contrasts our earlier theory/data comparison where it was demonstrated that diffusive transport in highly turbulent fields according to kinetic theory can successfully account for both the proton distributions and upstream ramp scales, using a single turbulence parameter. The principal conclusion here is that, in a FLW scenario, the transport of ions across the mean magnetic field is slightly less efficient than is required to effectively trap energetic ions within a few Larmor radii of the shock layer and thereby precipitate efficient acceleration. This highlights the contrast between ion transport in highly turbulent shock environs and remote, less-disturbed interplanetary regions.

The observed spectra of blazars, their intrinsic emission, and the underlying populations of radiating particles are intimately related. The use of these sources as probes of the extragalactic infrared background, a prospect propelled by recent advances in TeV-band telescopes, soon to be augmented by observations by NASA's upcoming Gamma-Ray Large Area Space Telescope (GLAST), has been a topic of great recent interest. Here, it is demonstrated that if particles in blazar jets are accelerated at relativistic shocks, then gamma-ray spectra with indices less than 1.5 can be produced. This, in turn, loosens the upper limits on the near infrared extragalactic background radiation previously proposed. We also show evidence hinting that TeV blazars with flatter spectra have higher intrinsic TeV gamma-ray luminosities and we indicate that there may be a correlation of flatness and luminosity with redshift.

Recent advances in the understanding of the properties of supernova remnant shocks have been precipitated by the Chandra and XMM X-ray Observatories, and the HESS Atmospheric Cerenkov Telescope in the TeV band. A critical problem for this field is the understanding of the relative degree of dissipative heating/energization of electrons and ions in the shock layer. This impacts the interpretation of X-ray observations, and moreover influences the efficiency of injection into the acceleration process, which in turn feeds back into the thermal shock layer energetics and dynamics. This paper outlines the first stages of our exploration of the role of charge separation potentials in non-relativistic electron-ion shocks where the inertial gyro-scales are widely disparate, using results from a Monte Carlo simulation. Charge density spatial profiles were obtained in the linear regime, sampling the inertial scales for both ions and electrons, for different magnetic field obliquities. These were readily integrated to acquire electric field profiles in the absence of self-consistent, spatial readjustments between the electrons and the ions. It was found that while diffusion plays little role in modulating the linear field structure in highly oblique and perpendicular shocks, in quasi-parallel shocks, where charge separations induced by gyrations are small, and shock-layer electric fields are predominantly generated on diffusive scales.

Heliospheric shocks are excellent systems for testing theories of particle acceleration in their environs. These generally fall into two classes: (1) interplanetary shocks that are linear in their ion acceleration characteristics, with the non-thermal ions serving as test particles, and (2) non-linear systems such as the Earth's bow shock and the solar wind termination shock, where the accelerated ions strongly influence the magnetohydrodynamic structure of the shock. This paper explores the modelling of diffusive acceleration at a particular interplanetary shock, with an emphasis on explaining in situ measurements of ion distribution functions. The observational data for this event was acquired on day 292 of 1991 by the Ulysses mission. The modeling is performed using a well-known kinetic Monte Carlo simulation, which has yielded good agreement with observations at several heliospheric shocks, as have other theoretical techniques, namely hybrid plasma simulations, and numerical solution of the diffusion-convection equation. In this theory/data comparison, it is demonstrated that diffusive acceleration theory can, to first order, successfully account for both the proton distribution data near the shock, and the observation of energetic protons farther upstream of this interplanetary shock than lower energy pick-up protons, using a single turbulence parameter. The principal conclusion is that diffusive acceleration of inflowing upstream ions can model this pick-up ion-rich event without the invoking any seed pre-acceleration mechanism, though this investigation does not rule out the action of such pre-acceleration.

The acceleration of interstellar pick-up ions as well as solar wind species has been observed at a multitude of interplanetary (IP) shocks by different spacecraft. The efficiency of injection of the pick-up ion component differs from that of the solar wind, and is expected to be strongly enhanced at highly oblique and quasi-perpendicular shock events, in accord with inferences from {\it in situ} observations. This paper explores theoretical modeling of the phase space distributions of accelerated ions obtained by the Ulysses mission for the Day 292, 1991 shock associated with a corotating interaction region, encountered before Ulysses' fly-by of Jupiter. A Monte Carlo simulation is used to model the acceleration process, adapting a technique that has been successfully tested on earlier IP shocks possessing minimal pick-up ion presence. Phase space distributions from the simulation technique for various low mass ions are compared with SWICS and HI-SCALE data to deduce values of a ``turbulence parameter'' that controls the efficiency of injection, and the degree of cross-field diffusion. Acceptable fits are obtained for the $H^+$ and $He^+$ populations using standard prescriptions for the pick-up ion distribution; $He^{++}$ spectral data was only fit well for scenarios very close to the Bohm diffusion limit. It is also found that the simulation successfully accounts for the observation of energetic protons farther upstream of the forward shock than lower energy pick-up protons, using the same turbulence parameter that is required to achieve reasonable spectral fits.