We consider anisotropic diffusion of Galactic cosmic rays in the Galactic magnetic field, using the Jansson-Farrar model for the field. In this paper we investigate the influence of source position on the cosmic ray flux at Earth in two ways: [1] by considering the contribution from cosmic ray sources located in different intervals in Galacto-centric radius, and [2] by considering the contribution from a number of specific and individual close-by supernova remnants. Our calculation is performed by using a fully three-dimensional stochastic method. This method is based on the numerical solution of a set of stochastic differential equations, equivalent to Ito formulation, that describes the propagation of the Galactic cosmic rays.

The Fermi Collaboration has recently updated their analysis of gamma rays from the center of the Galaxy. They reconfirm the presence of an unexplained emission feature which is most prominent in the region of $1-10$ GeV, known as the Galactic Center GeV excess (GCE). Although the GCE is now firmly detected, an interpretation of this emission as a signal of self-annihilating dark matter (DM) particles is not unambiguously possible due to systematic effects in the gamma-ray modeling estimated in the Galactic Plane. In this paper we build a covariance matrix, collecting different systematic uncertainties investigated in the Fermi Collaboration's paper that affect the GCE spectrum. We show that models where part of the GCE is due to annihilating DM can still be consistent with the new data. We also re-evaluate the parameter space regions of the minimal supersymmetric Standard Model (MSSM) that can contribute dominantly to the GCE via neutralino DM annihilation. All recent constraints from DM direct detection experiments such as PICO, LUX, PandaX and Xenon1T, limits on the annihilation cross section from dwarf spheroidal galaxies and Large Hadron Collider limits are considered in this analysis. Due to a slight shift in the energy spectrum of the GC excess with respect to the previous Fermi analysis, and the recent limits from direct detection experiments, we find a slightly shifted parameter region of the MSSM compared to our previous analysis that is consistent with the GCE. Neutralinos with a mass between $85-220$ GeV can describe the excess via annihilation into a pair of $W$-bosons or top quarks. Remarkably, there are low fine-tuning models among the regions that we have found. The complete set of solutions will be probed by upcoming direct detection experiments and with dedicated searches in the upcoming data of the Large Hadron Collider.

We have carried out a detailed study to understand the observed energy spectrum and composition of cosmic rays with energies up to ~10^18 eV. Our study shows that a single Galactic component with subsequent energy cut-offs in the individual spectra of different elements, optimised to explain the observed spectra below ~10^14 eV and the knee in the all-particle spectrum, cannot explain the observed all-particle spectrum above ~2x10^16 eV. We discuss two approaches for a second component of Galactic cosmic rays -- re-acceleration at a Galactic wind termination shock, and supernova explosions of Wolf-Rayet stars, and show that the latter scenario can explain almost all observed features in the all-particle spectrum and the composition up to ~10^18 eV, when combined with a canonical extra-galactic spectrum expected from strong radio galaxies or a source population with similar cosmological evolution. In this two-component Galactic model, the knee at ~ 3x10^15 eV and the second knee at ~10^17 eV in the all-particle spectrum are due to the cut-offs in the first and second components, respectively. We also discuss several variations of the extra-galactic component, from a minimal contribution to scenarios with a significant component below the ankle (at ~4x10^18 eV), and find that extra-galactic contributions in excess of regular source evolution are neither indicated nor in conflict with the existing data. Our main result is that the second Galactic component predicts a composition of Galactic cosmic rays at and above the second knee that largely consists of helium or a mixture of helium and CNO nuclei, with a weak or essentially vanishing iron fraction, in contrast to most common assumptions. This prediction is in agreement with new measurements from LOFAR and the Pierre Auger Observatory which indicate a strong light component and a rather low iron fraction between ~10^17 and 10^18 eV.

Observations with the Fermi Large Area Telescope (LAT) indicate an excess in gamma rays originating from the center of our Galaxy. A possible explanation for this excess is the annihilation of Dark Matter particles. We have investigated the annihilation of neutralinos as Dark Matter candidates within the phenomenological Minimal Supersymmetric Standard Model (pMSSM). An iterative particle filter approach was used to search for solutions within the pMSSM. We found solutions that are consistent with astroparticle physics and collider experiments, and provide a fit to the energy spectrum of the excess. The neutralino is a Bino/Higgsino or Bino/Wino/Higgsino mixture with a mass in the range $84-92$~GeV or $87-97$~GeV annihilating into W bosons. A third solutions is found for a neutralino of mass $174-187$~GeV annihilating into top quarks. The best solutions yield a Dark Matter relic density $0.06 < \Omega h^2 <0.13$. These pMSSM solutions make clear forecasts for LHC, direct and indirect DM detection experiments. If the MSSM explanation of the excess seen by Fermi-LAT is correct, a DM signal might be discovered soon.

We present results from imaging of the radio filaments in the southern giant lobe of Centaurus A using data from GMRT observations at 325 and 235 MHz, and outcomes from filament modelling. The observations reveal a rich filamentary structure, largely matching the morphology at 1.4 GHz. We find no clear connection of the filaments to the jet. We seek to constrain the nature and origin of the vertex and vortex filaments associated with the lobe and their role in high-energy particle acceleration. We deduce that these filaments are at most mildly overpressured with respect to the global lobe plasma showing no evidence of large-scale efficient Fermi I-type particle acceleration, and persist for ~ 2-3 Myr. We demonstrate that the dwarf galaxy KK 196 (AM 1318-444) cannot account for the features, and that surface plasma instabilities, the internal sausage mode and radiative instabilities are highly unlikely. An internal tearing instability and the kink mode are allowed within the observational and growth time constraints and could develop in parallel on different physical scales. We interpret the origin of the vertex and vortex filaments in terms of weak shocks from transonic MHD turbulence or from a moderately recent jet activity of the parent AGN, or an interplay of both.

Various radio galaxies show signs of having gone through episodic jet outbursts in the past. An example is the class of double-double radio galaxies (DDRGs). However, to follow the evolution of an individual source in real-time is impossible due to the large time scales involved. Numerical studies provide a powerful tool to investigate the temporal behavior of episodic jet outbursts in a (magneto-)hydrodynamical setting. We simulate the injection of two jets from active galactic nuclei (AGN), separated by a short interruption time. Three different jet models are compared. We find that an AGN jet outburst cycle can be divided into four phases. The most prominent phase occurs when the restarted jet is propagating completely inside the hot and inflated cocoon left behind by the initial jet. In that case, the jet-head advance speed of the restarted jet is significantly higher than the initial jet-head. While the head of the initial jet interacts strongly with the ambient medium, the restarted jet propagates almost unimpeded. As a result, the restarted jet maintains a strong radial integrity. Just a very small fraction of the amount of shocked jet material flows back through the cocoon compared to that of the initial jet and much weaker shocks are found at the head of the restarted jet. We find that the features of the restarted jet in this phase closely resemble the observed properties of a typical DDRG.

Observations of the FR I radio galaxy Centaurus A in radio, X-ray and gamma-ray bands provide evidence for lepton acceleration up to several TeV and clues about hadron acceleration to tens of EeV. Synthesising the available observational constraints on the physical conditions and particle content in the jets, inner lobes and giant lobes of Centaurus A, we aim to evaluate its feasibility as an ultra-high-energy cosmic-ray source. We apply several methods of determining jet power and affirm the consistency of various power estimates of ~ 1 x 10^43 erg s^-1. Employing scaling relations based on previous results for 3C 31, we estimate particle number densities in the jets, encompassing available radio through X-ray observations. Our model is compatible with the jets ingesting ~ 3 x 10^21 g s^-1 of matter via external entrainment from hot gas and ~ 7 x 10^22 g s^-1 via internal entrainment from jet-contained stars. This leads to an imbalance between the internal lobe pressure available from radiating particles and magnetic field, and our derived external pressure. Based on knowledge of the external environments of other FR I sources, we estimate the thermal pressure in the giant lobes as 1.5 x 10^-12 dyn cm^-2, from which we deduce a lower limit to the temperature of ~ 1.6 x 10^8 K. Using dynamical and buoyancy arguments, we infer ~ 440-645 Myr and ~ 560 Myr as the sound-crossing and buoyancy ages of the giant lobes respectively, inconsistent with their spectral ages. We re-investigate the feasibility of particle acceleration via stochastic processes in the lobes, placing new constraints on the energetics and on turbulent input to the lobes. The same 'very hot' temperatures that allow self-consistency between the entrainment calculations and the missing pressure also allow stochastic UHECR acceleration models to work.

Current observations have shown that astrophysical jets reveal strong signs of radial structure. They suggest that the inner region of the jet, the jet spine, consists of a low-density, fast-moving gas, while the outer region of the jet consists of a more dense and slower moving gas, called the jet sheath. Moreover, if jets carry angular momentum, the resultant centrifugal forces lead to a radial stratification. Current observations are not able to fully resolve the radial structure, so little is known about its actual profile. We present three AGN jet models in $2.5D$ of which two have been given a radial structure. The first model is a homogeneous jet, the only model that doesn't carry angular momentum; the second model is a spine-sheath jet with an isothermal equation of state; and the third jet model is a (piecewise) isochoric spine-sheath jet, with constant but different densities for jet spine and jet sheath. In this paper, we look at the effects of radial stratification on jet integrity, mixing between the different jet components and global morphology of the jet-head and surrounding cocoon.

We present a more accurate numerical scheme for the calculation of diffusive shock acceleration of cosmic rays using Stochastic Differential Equations. The accuracy of this scheme is demonstrated using a simple analytical flow profile that contains a shock of finite width and a varying diffusivity of the cosmic rays, where the diffusivity decreases across the shock. We compare the results for the slope of the momentum distribution with those obtained from a perturbation analysis valid for finite but small shock width. These calculations show that this scheme, although computationally more expensive, provides a significantly better performance than the Cauchy-Euler type schemes that were proposed earlier in the case where steep gradients in the cosmic ray diffusivity occur. For constant diffusivity the proposed scheme gives similar results as the Cauchy-Euler scheme.

We simulate time-dependent particle acceleration in the blast wave of a young supernova remnant (SNR), using a Monte Carlo approach for the diffusion and acceleration of the particles, coupled to an MHD code. We calculate the distribution function of the cosmic rays concurrently with the hydrodynamic evolution of the SNR, and compare the results with those obtained using simple steady-state models. The surrounding medium into which the supernova remnant evolves turns out to be of great influence on the maximum energy to which particles are accelerated. In particular, a shock going through a $\rho \propto r^{-2}$ density profile causes acceleration to typically much higher energies than a shock going through a medium with a homogeneous density profile. We find systematic differences between steady-state analytical models and our time-dependent calculation in terms of spectral slope, maximum energy, and the shape of the cut-off of the particle spectrum at the highest energies. We also find that, provided that the magnetic field at the reverse shock is sufficiently strong to confine particles, cosmic rays can be easily re-accelerated at the reverse shock.

Observations show that the magnetic field in young supernova remnants (SNRs) is significantly stronger than can be expected from the compression of the circumstellar medium (CSM) by a factor of four expected for strong blast waves. Additionally, the polarization is mainly radial, which is also contrary to expectation from compression of the CSM magnetic field. Cosmic rays (CRs) may help to explain these two observed features. They can increase the compression ratio to factors well over those of regular strong shocks by adding a relativistic plasma component to the pressure, and by draining the shock of energy when CRs escape from the region. The higher compression ratio will also allow for the contact discontinuity, which is subject to the Rayleigh-Taylor (R-T) instability, to reach much further out to the forward shock. This could create a preferred radial polarization of the magnetic field. With an adaptive mesh refinement MHD code (AMRVAC), we simulate the evolution of SNRs with three different configurations of the initial CSM magnetic field, and look at two different equations of state in order to look at the possible influence of a CR plasma component. The spectrum of CRs can be simulated using test particles, of which we also show some preliminary results that agree well with available analytical solutions.

Supernova remnants (SNR) are now widely believed to be a source of cosmic rays (CRs) up to an energy of 1 PeV. The magnetic fields required to accelerate CRs to sufficiently high energies need to be much higher than can result from compression of the circumstellar medium (CSM) by a factor 4, as is the case in strong shocks. Non-thermal synchrotron maps of these regions indicate that indeed the magnetic field is much stronger, and for young SNRs has a dominant radial component while for old SNRs it is mainly toroidal. How these magnetic fields get enhanced, or why the field orientation is mainly radial for young remnants, is not yet fully understood. We use an adaptive mesh refinement MHD code, AMRVAC, to simulate the evolution of supernova remnants and to see if we can reproduce a mainly radial magnetic field in early stages of evolution. We follow the evolution of the SNR with three different configurations of the initial magnetic field in the CSM: an initially mainly toroidal field, a turbulent magnetic field, and a field parallel to the symmetry axis. Although for the latter two topologies a significant radial field component arises at the contact discontinuity due to the Rayleigh-Taylor instability, no radial component can be seen out to the forward shock. Ideal MHD appears not sufficient to explain observations. Possibly a higher compression ratio and additional turbulence due to dominant presence of CRs can help us to better reproduce the observations in future studies.

We present hydrodynamical models for the Cassiopeia A (Cas A) supernova remnant and its observed jet / counter-jet system. We include the evolution of the progenitor's circumstellar medium, which is shaped by a slow red supergiant wind that is followed by a fast Wolf-Rayet (WR) wind. The main parameters of the simulations are the duration of the WR phase and the jet energy. We find that the jet is destroyed if the WR phase is sufficiently long and a massive circumstellar shell has formed. We therefore conclude that the WR phase must have been short (a few thousand yr), if present at all. Since the actual jet length of Cas A is not known we derive a lower limit for the jet energy, which is ~10^{48} erg. We discuss the implications for the progenitor of Cas A and the nature of its explosion.

A search for TeV - PeV muon neutrinos from unresolved sources was performed on AMANDA-II data collected between 2000 and 2003 with an equivalent livetime of 807 days. This diffuse analysis sought to find an extraterrestrial neutrino flux from sources with non-thermal components. The signal is expected to have a harder spectrum than the atmospheric muon and neutrino backgrounds. Since no excess of events was seen in the data over the expected background, an upper limit of E^{2}\Phi_{90% C.L.} < 7.4 x 10^{-8} GeV cm^{-2} s^{-1} sr^{-1} is placed on the diffuse flux of muon neutrinos with a \Phi \propto E^{-2} spectrum in the energy range 16 TeV to 2.5 PeV. This is currently the most sensitive \Phi \propto E^{-2} diffuse astrophysical neutrino limit. We also set upper limits for astrophysical and prompt neutrino models, all of which have spectra different than \Phi \propto E^{-2}.

The IceCube neutrino detector is a cubic kilometer TeV to PeV neutrino detector under construction at the geographic South Pole. The dominant population of neutrinos detected in IceCube is due to meson decay in cosmic-ray air showers. These atmospheric neutrinos are relatively well-understood and serve as a calibration and verification tool for the new detector. In 2006, the detector was approximately 10% completed, and we report on data acquired from the detector in this configuration. We observe an atmospheric neutrino signal consistent with expectations, demonstrating that the IceCube detector is capable of identifying neutrino events. In the first 137.4 days of livetime, 234 neutrino candidates were selected with an expectation of 211 +/- 76.1(syst.) +/- 14.5(stat.) events from atmospheric neutrinos.

We present the results of the analysis of neutrino observations by the Antarctic Muon and Neutrino Detector Array (AMANDA) correlated with photon observations of more than 400 gamma-ray bursts (GRBs) in the Northern Hemisphere from 1997 to 2003. During this time period, AMANDA's effective collection area for muon neutrinos was larger than that of any other existing detector. Based on our observations of zero neutrinos during and immediately prior to the GRBs in the dataset, we set the most stringent upper limit on muon neutrino emission correlated with gamma-ray bursts. Assuming a Waxman-Bahcall spectrum and incorporating all systematic uncertainties, our flux upper limit has a normalization at 1 PeV of E^2{\Phi}_{\nu} {\leq} 6.0 \times 10^{-9} GeV cm^{-2} s^{-1} sr^{-1}, with 90% of the events expected within the energy range of ~10 TeV to ~3 PeV. The impact of this limit on several theoretical models of GRBs is discussed, as well as the future potential for detection of GRBs by next generation neutrino telescopes. Finally, we briefly describe several modifications to this analysis in order to apply it to other types of transient point sources.

Using the neutrino telescope AMANDA-II, we have conducted two analyses searching for neutrino-induced cascades from gamma-ray bursts. No evidence of astrophysical neutrinos was found, and limits are presented for several models. We also present neutrino effective areas which allow the calculation of limits for any neutrino production model. The first analysis looked for a statistical excess of events within a sliding window of 1 or 100 seconds (for short and long burst classes, respectively) during the years 2001-2003. The resulting upper limit on the diffuse flux normalization times E^2 for the Waxman-Bahcall model at 1 PeV is 1.6 x 10^-6 GeV cm^-2 s^-1 sr^-1 (a factor of 120 above the theoretical prediction). For this search 90% of the neutrinos would fall in the energy range 50 TeV to 7 PeV. The second analysis looked for neutrino-induced cascades in coincidence with 73 bursts detected by BATSE in the year 2000. The resulting upper limit on the diffuse flux normalization times E^2, also at 1 PeV, is 1.5 x 10^-6 GeV cm^-2 s^-1 sr^-1 (a factor of 110 above the theoretical prediction) for the same energy range. The neutrino-induced cascade channel is complementary to the up-going muon channel. We comment on its advantages for searches of neutrinos from GRBs and its future use with IceCube.

We present hydrodynamical models of circumstellar medium (CSM) of long gamma-ray burst (GRB) progenitor candidates. These are massive stars that have lost a large amount of mass in the form of stellar wind during their evolution. There are two possible ways to probe the CSM of long GRB progenitors. Firstly, the GRB afterglow consists of synchrotron radiation, emitted when the GRB jet sweeps up the surrounding medium. Therefore, the lightcurve is directly related to the density profile of the CSM. The density can either decrease with the radius squared (as is the case for a freely expanding stellar wind) or be constant (as we would expect for shocked wind or the interstellar medium). Secondly, material between the GRB and the observer will absorb part of the afterglow radiation, causing absorption lines in the afterglow spectrum. In some cases, such absorption lines are blue-shifted relative to the source indicating that the material is moving away from the progenitor star. This can be explained in terms of wind interactions in the CSM. We can use the CSM of these stars to investigate their prior evolutionary stage.

This collection of proceedings to the TeV PA II Conference presents some of the latest results of the IceCube Collaboration.

We report the results of a five-year survey of the northern sky to search for point sources of high energy neutrinos. The search was performed on the data collected with the AMANDA-II neutrino telescope in the years 2000 to 2004, with a live-time of 1001 days. The sample of selected events consists of 4282 upward going muon tracks with high reconstruction quality and an energy larger than about 100 GeV. We found no indication of point sources of neutrinos and set 90% confidence level flux upper limits for an all-sky search and also for a catalog of 32 selected sources. For the all-sky search, our average (over declination and right ascension) experimentally observed upper limit \Phi^{0}=(E/TeV)^\gamma d\Phi/dE to a point source flux of muon and tau neutrino (detected as muons arising from taus) is \Phi_{\nu_\mu}^{0} + \Phi_{\nu_\tau}^{0} = 11.1 x 10^{-11} TeV^-1 cm^-2 s^-1, in the energy range between 1.6 TeV and 2.5 PeV for a flavor ratio \Phi_{\nu_\mu}^{0} / \Phi_{\nu_\tau}^{0}= 1 and assuming a spectral index \gamma=2. It should be noticed that this is the first time we set upper limits to the flux of muon and tau neutrinos. In previous papers we provided muon neutrino upper limits only neglecting the sensitivity to a signal from tau neutrinos, which improves the limits by 10% to 16%. The value of the average upper limit presented in this work corresponds to twice the limit on the muon neutrino flux \Phi_{\nu_\mu}^{0} = 5.5 \cdot 10^{-11} TeV^-1 cm^-2 s^-1. A stacking analysis for preselected active galactic nuclei and a search based on the angular separation of the events were also performed. We report the most stringent flux upper limits to date, including the results of a detailed assessment of systematic uncertainties.

The progenitor stars of long Gamma-Ray Bursts (GRBs) are thought to be Wolf-Rayet stars, which generate a massive and energetic wind. Nevertheless, about 25 percent of all GRB afterglows light curves indicate a constant density medium close to the exploding star. We explore various ways to produce this, by creating situations where the wind termination shock arrives very close to the star, as the shocked wind material has a nearly constant density. Typically, the distance between a Wolf-Rayet star and the wind termination shock is too large to allow afterglow formation in the shocked wind material. Here, we investigate possible causes allowing for a smaller distance: A high density or a high pressure in the surrounding interstellar medium (ISM), a weak Wolf-Rayet star wind, the presence of a binary companion, and fast motion of the Wolf-Rayet star relative to the ISM. We find that all four scenarios are possible in a limited parameter space, but that none of them is by itself likely to explain the large fraction of constant density afterglows. A low GRB progenitor metallicity, and a high GRB energy make the occurrence of a GRB afterglow in a constant density medium more likely. This may be consistent with constant densities beingpreferentially found for energetic, high redshift GRBs.

We discuss magnetic field generation by the proton Weibel instability in relativistic shocks, a situation that applies to the external shocks in the fireball model for Gamma-ray Bursts, and possibly also to internal shocks. Our analytical estimates show that the linear phase of the instability ends well before it has converted a significant fraction of the energy in the proton beam into magnetic energy: the conversion efficiency is much smaller (of order m_e/m_p) in electron-proton plasmas than in pair plasmas. We find this estimate by modelling the plasma in the shock transition zone with a waterbag momentum distribution for the protons and with a background of hot electrons. For ultra-relativistic shocks we find that the wavelength of the most efficient mode for magnetic field generation equals the electron skin depth, that the relevant nonlinear stabilization mechanism is magnetic trapping, and that the presence of the hot electrons limits the typical magnetic field strength generated by this mode so that it does not depend on the energy content of the protons. We conclude that other processes than the linear Weibel instability must convert the free energy of the protons into magnetic fields.

We present an analytical estimate for the magnetic field strength generated by the Weibel instability in ultra-relativistic shocks in a hydrogen plasma. We find that the Weibel instability is, by itself, not capable of converting the kinetic energy of protons penetrating the shock front into magnetic field energy. Other (nonlinear) processes must determine the magnetic field strength in the wake of the shock.

Hydrodynamical simulations are presented of a pulsar wind emitted by a supersonically moving pulsar. The pulsar moves through the interstellar medium or, in the more interesting case, through the supernova remnant created at its birth event. In both cases there exists a three-fold structure consisting of the wind termination shock, contact discontinuity and a bow shock bounding the pulsar wind nebula. Using hydrodynamical simulations we study the behaviour of the pulsar wind nebula inside a supernova remnant, and in particular the interaction with the outer shell of swept up interstellar matter and the blast wave surrounding the remnant. This interaction occurs when the pulsar breaks out of the supernova remnant. We assume the remnant is in the Sedov stage of its evolution. Just before break-through, the Mach number associated with the pulsar motion equals ${\cal M}_{\rm psr} = 7/\sqrt{5}$, {\em independent} of the supernova explosion energy and pulsar velocity. The bow shock structure is shown to survive this break-through event.

We reconsider the rejuvenation mechanism as proposed by Shull, Fesen, & Saken (1989). These authors suggest that an active pulsar can catch up with, and rejuvenate the shell of the associated supernova remnant. The morphology of the SNRs G5.4-1.2 and CTB80 seem to confirm this rejuvenation mechanism. The spindown energy is deposited by the pulsar as a relativistic pulsar wind, and has a sufficient power to explain the observed radio emission observed in these remnants. Shull et al. (1989) did {\it not} explain the observed lengthscales of the rejuvenated parts of the SNR shell. therefore one needs to consider the diffusive transport of the injected electrons by the pulsar wind. We propose to apply a diffusion mechanism as introduced by Jokipii (1987), which makes a distinction between diffusion along the magnetic field lines and perpendicular to the magnetic field lines, parameterised by the gyro factor $\eta$. We show that one has to assume a high value for the gyro factor, $\eta\simeq 10^3-10^4$, i.e. diffusion of the electrons along the magnetic field line is much faster then perpendicular to the magnetic field line, in order for the rejuvenation mechanism to work on the observed lengthscales.

We review recent theoretical developments on pulsar winds, their nebulae and relativistic shock acceleration, and show how they illuminate unsolved problems in plerion spectra, in particular the multiple spectral breaks in the Crab and the low-frequency breaks of plerions such as G 21.5-0.9 and 3C 58. Recent work on Fermi acceleration theory at relativistic shocks shows that a particle spectral index of 2.2-2.3, compatible with the X-ray spectra of plerions, results under a wide variety of assumptions. If pulsar winds contain ions as well as electrons and positrons, the mechanism of Hoshino et al. (1992), which yields harder spectra, would operate at lower energies and may explain the flat radio spectral indices of plerions. This scenario implies wind parameters in the Crab compatible with the pulsar wind acceleration model of Lyubarsky & Kirk (2001). Recent hydrodynamical simulations of plerion evolution inside SNR blast waves demonstrate that the passage of the reverse shock rapidly compresses the plerion. Using a simple isobaric model, we investigate the influence of the resulting magnetic field compression and decrease in shock radius on the evolution of the plerion spectrum. We suggest that the passage of the reverse shock may explain the low-frequency breaks in 3C 58 and G 21.5-0.9, as well as the increase in 3C 58's radio flux.

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.

A spherically symmetric model is presented for the interaction of a pulsar wind with the associated supernova remnant. This results in a pulsar wind nebula whose evolution is coupled to the evolution of the surrounding supernova remnant. This evolution can be divided in three stages. The first stage is characterised by a supersonic expansion of the pulsar wind nebula into the freely expanding ejecta of the progenitor star. In the next stage the pulsar wind nebula is not steady; the pulsar wind nebula oscillates between contraction and expansion due to interaction with the reverse shock of the supernova remnant: reverberations which propagate forward and backward in the remnant. After the reverberations of the reverse shock have almost completely vanished and the supernova remnant has relaxed to a Sedov solution, the expansion of the pulsar wind nebula proceeds subsonically. In this paper we present results from hydrodynamical simulations of a pulsar wind nebula through all these stages in its evolution. The simulations were carried out with the Versatile Advection Code.

The cosmic-ray spectrum up to the knee ($E\sim 10^{15}$ eV) is attributed to acceleration processes taking place at the blastwaves which bound supernova remnants. Theoretical predictions give a similar estimate for the maximum energy which can be reached at supernova remnant shocks by particle acceleration. Electrons with energies of the order $\sim 10^{15}$ eV should give a nonthermal X-ray component in young supernova remnants. Recent observations of SN1006 and G347.3-0.5 confirm this prediction. We present a method which uses hydrodynamical simulations to describe the evolution of a young remnant. These results are combined with an algorithm which simultaneously calculates the associated particle acceleration. We use the test particle approximation, which means that the back-reaction on the dynamics of the remnant by the energetic particles is neglected. We present synchrotron maps in the X-ray domain, and present spectra of the energies of the electrons in the supernova remnant. Some of our results can be compared directly with earlier semi-analytical work on this subject by Reynolds [1].

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.

We discuss the intergalactic propagation of ultra-high-energy cosmic rays (UHECRs) with energies E \geq 10^{18.5} eV. We consider the propagation of UHECRs under the influence of the energy-dependent deflection by a weak random magnetic field in the intergalactic medium and energy losses by photo-pion and pair production. We calculate arrival spectra taking full account of the kinematics of photo-pion production and the Poisson statistics of the photo-pion interaction rate. We give estimates for the deflection of UHECRs from the line of sight to the source, time delays with respect to photons from the same source, arrival spectra and source statistics. These estimates are confirmed by numerical simulations of the propagation in energy evolution of UHECRs. These simulations demonstrate that the often-used continuous approximation in the treatment of energy losses due to photo-pion production on the cosmic microwave background (CMWB) cannot be justified for UHECRs. We discuss the implications of these results for the observed flux of particles above the Greisen-Zatsepin-Kuz'min cut-off in two of the scenarios that have been proposed for the production of these particles: continuous production in the large shock waves associated with powerful radio galaxies, or possibly large-scale structure formation, and the impulsive production at relativistic blast waves associated with cosmological gamma-ray bursts.

We consider the acceleration of charged particles at the ultra-relativistic shocks, with Lorentz factors \Gamma_s >> 1 relative to the upstream medium, arising in relativistic fireball models of gamma-ray bursts (GRBs). We show that for Fermi-type shock acceleration, particles initially isotropic in the upstream medium can gain a factor of order \Gamma_s^2 in energy in the first shock crossing cycle, but that the energy gain factor for subsequent shock crossing cycles is only of order 2, because for realistic deflection processes particles do not have time to re-isotropise upstream before recrossing the shock. We evaluate the maximum energy attainable and the efficiency of this process, and show that for a GRB fireball expanding into a typical interstellar medium, these exclude the production of ultra-high-energy cosmic rays (UHECRs), with energies in the range 10^{18.5} - 10^{20.5} eV, by the blast wave. We propose, however, that in the context of neutron star binaries as the progenitors of GRBs, relativistic ions from the pulsar wind bubbles produced by these systems could be accelerated by the blast wave. We show that if the known binary pulsars are typical, the maximum energy, efficiency, and spectrum in this case can account for the observed population of UHECRs.

In the theory of diffusive acceleration at oblique shock fronts the question of the existence of a discontinuity of energetic particle density is contentious. The resolution of this problem is interesting from a theoretical point of view, and potentially for the interpretation of observations of particle densities at heliospheric shocks and of high-resolution radio observations of the rims of supernova remnants. It can be shown analytically that an isotropic particle distribution at a shock front implies continuity of the particle density -- whether or not the shock is oblique. However, if the obliquity of the shock induces an anisotropy, a jump is permitted. Both semi-analytic computations and Monte-Carlo simulations are used to show that, for interesting parameter ranges, a jump is indeed produced, with accelerated particles concentrated in a precursor ahead of the shock front.