The "Mouse" (PWN G359.23-0.82) is a spectacular bow shock pulsar wind nebula, powered by the radio pulsar J1747-2958. The pulsar and its nebula are presumed to have a high space velocity, but their proper motions have not been directly measured. Here we present 8.5 GHz interferometric observations of the Mouse nebula with the Very Large Array, spanning a time baseline of 12 yr. We measure eastward proper motion for PWN G359.23-0.82 (and hence indirectly for PSR J1747-2958) of 12.9+/-1.8 mas/yr, which at an assumed distance of 5 kpc corresponds to a transverse space velocity of 306+/-43 km/s. Considering pressure balance at the apex of the bow shock, we calculate an in situ hydrogen number density of approximately 1.0(-0.2)(+0.4) cm^(-3) for the interstellar medium through which the system is traveling. A lower age limit for PSR J1747-2958 of 163(-20)(+28) kyr is calculated by considering its potential birth site. The large discrepancy with the pulsar's spin-down age of 25 kyr is possibly explained by surface dipole magnetic field growth on a timescale ~15 kyr, suggesting possible future evolution of PSR J1747-2958 to a different class of neutron star. We also argue that the adjacent supernova remnant G359.1-0.5 is not physically associated with the Mouse system but is rather an unrelated object along the line of sight.

We present new X-ray imaging and spectroscopy of a composite supernova remnant G327.1-1.1 using the Chandra and XMM-Newton X-ray observatories. G327.1-1.1 has an unusual morphology consisting of a symmetric radio shell and an off center nonthermal component that indicates the presence of a pulsar wind nebula (PWN). Radio observations show a narrow finger of emission extending from the PWN structure towards the northwest. X-ray studies with ASCA, ROSAT, and BeppoSAX revealed elongated extended emission and a compact source at the tip of the finger that may be coincident with the actual pulsar. The high resolution Chandra observations provide new insight into the structure of the inner region of the remnant. The images show a compact source embedded in a cometary structure, from which a trail of X-ray emission extends in the southeast direction. The Chandra images also reveal two prong-like structures that appear to originate from the vicinity of the compact source and extend into a large bubble that is oriented in the north-west direction, opposite from the bright radio PWN. The emission from the entire radio shell is detected in the XMM data and can be characterized by a thermal plasma model with a temperature of 0.3 keV, which we use to estimate the physical properties of the remnant. The peculiar morphology of G327.1-1.1 may be explained by the emission from a moving pulsar and a relic PWN that has been disrupted by the reverse shock.

"The Duck" is a complicated non-thermal radio system, consisting of the energetic radio pulsar B1757-24, its surrounding pulsar wind nebula G5.27-0.90 and the adjacent supernova remnant (SNR) G5.4-1.2. PSR B1757-24 was originally claimed to be a young (~15 000 yr) and extreme velocity (>~1500 km/s) pulsar which had penetrated and emerged from the shell of the associated SNR G5.4-1.2, but recent upper limits on the pulsar's motion have raised serious difficulties with this interpretation. We here present 8.5 GHz interferometric observations of the nebula G5.27-0.90 over a 12-year baseline, doubling the time-span of previous measurements. These data correspondingly allow us to halve the previous upper limit on the nebula's westward motion to 14 milliarcseconds/yr (5-sigma), allowing a substantive reevaluation of this puzzling object. We rule out the possibility that the pulsar and SNR were formed from a common supernova explosion ~15 000 yrs ago as implied by the pulsar's characteristic age, but conclude that an old (>~70 000 yr) pulsar / SNR association, or a situation in which the pulsar and SNR are physically unrelated, are both still viable explanations.

We present a deep observation with the Chandra X-ray Observatory of the neutron star bow shock G189.22+2.90 in the supernova remnant (SNR) IC 443. Our data confirm the cometary morphology and central point source seen previously, but also reveal considerable new structure. Specifically, we find that the X-ray nebula consists of two distinct components: a "tongue" of bright emission close to the neutron star, enveloped by a larger, fainter "tail". We interpret the tongue and tail as delineating the termination shock and the post-shock flow, respectively, as previously identified also in the pulsar bow shock G359.23-0.82 ("the Mouse"). However, for G189.22+2.90 the tongue is much less elongated than for the Mouse, while the tail is much broader. These differences are consistent with the low Mach number, M >~ 2, expected for a neutron star moving through the hot gas in a SNR's interior, supporting the case for a physical association between G189.22+2.90 and IC 443. We resolve the stand-off distance between the star and the head of the bow shock, which allows us to estimate a space velocity for the neutron star of ~230 km/s, independent of distance. We detect thermal emission from the neutron star surface at a temperature of 102 +/- 22 eV, which is consistent with the age of SNR IC 443 for standard neutron star cooling models. We also identify two compact knots of hard emission located 1-2 arcsec north and south of the neutron star.

The nature of the radio-quiet X-ray halo around the plerionic SNR G21.5-0.9 is under debate. On the basis of spatial and spectral analysis of a large Chandra and XMM-Newton dataset of this source, we have developed a self-consistent scenario which explains all the observational features. We found that the halo is composed by diffuse extended emission due to dust scattering of X-rays from the plerion, by a bright limb which traces particle acceleration in the fast forward shock of the remnant, and by a bright spot (the ``North Spur'') which may be a knot of ejecta in adiabatic expansion. By applying a model of interaction between the PWN, the SNR and supernova environment, we argue that G21.5-0.9 progenitor may be of Type IIP or Ib/Ic, and that the remnant may be young (200-1000 yr).

We present analytical and numerical studies of magnetorotational instabilities occuring in magnetized accretion disks. In these studies we make use of the linearised compressible MHD equations. These calculations are performed for general radially stratified disks in the cylindrical limit. In particular, we investigate the influence of nonvanishing toroidal magnetic field component on the growth rate and oscillation frequency of magnetorotational instabilities in Keplerian disks. We find the persistence of these instabilities in accretion disks close to equipartition. Our calculations show that these eigenmodes become overstable (complex eigenvalue), due to the presence of a toroidal magnetic field component, while their growth rate reduces slightly. Furthermore, we demonstrate the presence of magneto-rotational overstabilities in weakly magnetized sub-Keplerian rotating disks. We show that the growth rate scales with the rotation frequency of the disk. These eigenmodes also have a nonzero oscillation frequency, due to the presence of the dominant toroidal magnetic field component. The overstable character of the MRI increases as the rotation frequency of the disk decreases.

We present a study of instabilities occuring in thick magnetized accretion disks. We calculate the growth rates of these instabilities and characterise precisely the contribution of the magneto-rotational and the convective mechanism. All our calculations are performed in radially stratified disks in the cylindrical limit. The numerical calculations are performed using the appropriate local dispersion equation solver discussed in Blokland et al. (2005). A comparison with recent results by Narayan et al. (2002) shows excellent agreement with their approximate growth rates only if the disks are weakly magnetized. However, for disks close to equipartition, the dispersion equation from Narayan et al. (2002) loses its validity. Our calculations allow for a quantitative determination of the increase of the growth rate due to the magneto-rotational mechanism. We find that the increase of the growth rate for long wavelength convective modes caused by this mechanism is almost neglible. On the other hand, the growth rate of short wavelength instabilities can be significantly increased by this mechanism, reaching values up to 60%.

In the late 1960s the discovery of the Crab pulsar in its associated supernova remnant, launched a new field in supernova remnant research: the study of pulsar-driven or plerionic supernova remnants. In these type of remnants, the relativistic wind emitted by the pulsar, blows a pulsar wind nebula into the interior of its supernova remnant. Now, more then forty years after the discovery of the Crab pulsar, there are more then fifty plerionic supernova remnants known, due to the ever-increasing capacity of observational facilities. I will review our current understanding of the different evolutionary stages of a pulsar wind nebula as it is interacting with its associated supernova remnant.Therefore I will discuss both analytical and more recent numerical (M)HD models.The four main stages of a pulsar wind nebula are: the supersonic expansion stage, the reverse shock interaction stage, the subsonic expansion stage and ultimatelythe stage when the head of the bubble is bounded by a bow shock, due to the supersonic motion of the pulsar. Ultimately this pulsar wind nebula bow shock will break through its associated remnant, after which the pulsar-powered bow shock will interact directly with the interstellar medium. I will discuss recent numerical models from these type of pulsar wind nebulae and their morphology.

We present an investigation of the spectral and spatial structure of the X-ray emission from 3C 58 based on a 350 ks observation with the Chandra X-ray Observatory. This deep image, obtained as part of the Chandra Large Project program, reveals new information on nearly all spatial scales in the pulsar wind nebula (PWN). On the smallest scales we derive an improved limit of T < 1.02 X 10^6 K for blackbody emission from the entire surface of the central neutron star (NS), confirming the need for rapid, nonstandard cooling in the stellar interior. Furthermore, we show that the data are consistent with emission from a light element atmosphere with a similar temperature. Surrounding the NS, a toroidal structure with a jet is resolved, consistent with earlier measurements and indicative of an east-west orientation for the projected rotation axis of the pulsar. A complex of loop-like X-ray filaments fills the nebula interior, and corresponds well with structures seen in the radio band. Several of the structures coincide with optical filaments as well. The emission from the interior of the PWN, including the pulsar, jet, and filaments, is primarily nonthermal in nature. The power law index steepens with radius, but appears to also show small azimuthal variations. The outermost regions of the nebula require a thermal emission component, confirming the presence of an ejecta-rich swept up shell.

We present an observation with the Chandra X-ray Observatory of the unusual radio source G359.23-0.82 ("the Mouse"), along with updated radio timing data from the Parkes radio telescope on the coincident young pulsar J1747-2958. We find that G359.23-0.82 is a very luminous X-ray source (L_X [0.5-8.0 keV] = 5e34 ergs/s for a distance of 5 kpc), whose morphology consists of a bright head coincident with PSR J1747-2958, plus a 45"-long narrow tail whose power-law spectrum steepens with distance from the pulsar. We thus confirm that G359.23-0.82 is a bow-shock pulsar wind nebula powered by PSR J1747-2958; the nebular stand-off distance implies that the pulsar is moving with a Mach number of ~60, suggesting a space velocity ~600 km/s through gas of density ~0.3 cm^-3. We combine the theory of ion-dominated pulsar winds with hydrodynamic simulations of pulsar bow shocks to show that a bright elongated X-ray and radio feature extending 10'' behind the pulsar represents the surface of the wind termination shock. The X-ray and radio "trails" seen in other pulsar bow shocks may similarly represent the surface of the termination shock, rather than particles in the postshock flow as is usually argued. The tail of the Mouse contains two components: a relatively broad region seen only at radio wavelengths, and a narrow region seen in both radio and X-rays. We propose that the former represents material flowing from the wind shock ahead of the pulsar's motion, while the latter corresponds to more weakly magnetized material streaming from the backward termination shock. This study represents the first consistent attempt to apply our understanding of "Crab-like" nebulae to the growing group of bow shocks around high-velocity pulsars.

Observations indicate that most massive stars in the Galaxy appear in groups, called OB associations, where their strong wind activity generates large structures known as superbubbles, inside which the subsequent supernovae (SNe) explode, in tight space and time correlation. Acknowledging this fact, we investigate four main questions: 1) does the clustering of massive stars and SN explosions influence the particle acceleration process usually associated with SNe, and induce collective effects which would not manifest around isolated supernova remnants?; 2) does it make a difference for the general phenomenology of Galactic Cosmic Rays (GCRs), notably for their energy spectrum and composition?; 3) can this help alleviate some of the problems encountered within the standard GCR source model?; and 4) Is the link between superbubbles and energetic particles supported by observational data, and can it be further tested and constrained? We argue for a positive answer to all these questions. Theoretical, phenomenological and observational aspects are treated in separate papers. Here, we discuss the interaction of massive stellar winds and SN shocks inside superbubbles and indicate how this leads to specific acceleration effects. We also show that due to the high SN explosion rate and low diffusion coefficient, low-energy particles experience repeated shock acceleration inside superbubbles.

We present a model of a pulsar wind nebula evolving inside its associated supernova remnant. The model uses a hydrodynamics code to simulate the evolution of this system when the pulsar has a high velocity. The simulation distinguishes four different stages of pulsar wind nebula evolution: the supersonic expansion stage, the reverse shock interaction stage, the subsonic expansion stage and ultimately the bow shock stage. The simulation bears out that, due to the high velocity of the pulsar, the position of the pulsar is off-centered with respect to its pulsar wind nebula, after the passage of the reverse shock. Subsequently the pulsar wind nebula expands subsonically untill the event of the bow shock formation, when the motion of the pulsar becomes supersonic. The bow shock formation event occurs at roughly half the crossing time, when the pulsar is positioned at 0.677 times the radius of the supernova remnant blastwave, in complete agreement with analytical predictions. The crossing time is defined by the age of the supernova remnant, when the pulsar overtakes the blastwave bounding the supernova remnant. The results of the model are applied to three supernova remnants: N157B, G327.1-1.1 and W44. We argue that the head of the pulsar wind nebula, containing the active pulsar, inside the first two systems are not bounded by a bow shock. However, in the case of W44 we argue for a scenario in which the pulsar wind nebula is bounded by a bow shock, due to the supersonic motion of the pulsar.

We present a model of a pulsar-driven supernova remnant, by using a hydrodynamics code, which simulates the evolution of a pulsar wind nebula when the pulsar is moving at a high velocity through its expanding supernova remnant. The simulation shows four different stages of the pulsar wind nebula: the supersonic expansion stage, the reverse shock interaction stage, the subsonic expansion stage and ultimately the bow shock stage. Due to the high velocity of the pulsar, the position of the pulsar is located at the head of the pulsar wind bubble, after the passage of the reverse shock. The resulting morphology of the pulsar wind bubble is therefore similar to the morphology of a bow shock pulsar wind nebula. We show how to distinguish these two different stages, and apply this method to the SNR G327.1-1.1, for which we argue that there is no bow shock around its pulsar wind nebula.

I present a model of a pulsar wind interacting with its associated supernova remnant. I will use the model to argue that one can explain the morphology of the pulsar wind nebula inside N157B, a supernova remnant in the Large Magellanic Cloud, without the need for a bow shock interpretation. The model uses a hydrodynamics code which simulates the evolution of a pulsar wind nebula, when the pulsar is moving at a high velocity (1000 km/sec) through the expanding supernova remnant. The simulation shows three different stages of the pulsar wind nebula: the supersonic expansion stage, the subsonic expansion stage and the bow shock stage. The last stage occurs when the head of the pulsar wind nebula, containing the active pulsar, deforms into a bow shock, due to the motion of the pulsar becoming supersonic. I will argue that the pulsar wind nebula inside N157B is currently in the second stage of its evolution, i.e. the expansion of the pulsar wind nebula is subsonic and there is no bow shock around the pulsar wind bubble. The strongly off-centered position of the pulsar with respect to its pulsar wind nebula is naturally explained by the result of the interaction of the reverse shock with the pulsar wind nebula, as the simulation bears out.

Young supernova remnants are thought to be the sites where cosmic ray acceleration occurs by the mechanism of diffusive shock acceleration. The maximum energy gained in this process is conventionally extimated to have a value close to, but distinctly below, the ``knee'' ($\sim 10^{15}$ eV) of the cosmic-ray spectrum. Bell & Lucek (2001) have suggested that the generated cosmic rays simultaneously amplify the magnetic field around the supernova remnant shock to many times its pre-shock value. In this case the acceleration rate may be significantly increased and protons can easily reach energies up to $10^{17}$ eV. We use a ``simplified'' box model incorporating the magnetic field amplification suggested by Bell & Lucek to investigate the resulting modifications of the cosmic-ray spectrum. The model predicts a spectral break at high energies, close to the ``knee'' region, and in good accordance with observations.

Magnetohydrodynamical simulations are presented of a magnetized pulsar wind interacting directly with the interstellar medium, or, in the case of a surrounding supernova remnant, with the associated freely expanding ejecta of the progenitor star. In both cases the simulations show that the pulsar wind nebula will be elongated due to the dynamical influence of the toroidal magnetic fields, which confirm predictions from a semi-analytical model presented by Begelman & Li. The simulations follow the expansion of the pulsar wind nebula when the latter is bounded by a strong shock and show that the expansion can be modeled with a standard power-law expansion rate. By performing different simulations with different magnetization parameters, I show that the latter weakly correlates with the elongation of the pulsar wind nebula. The results from the simulations are applied to determine the nature of the expansion rate of the pulsar wind nebula 3C58. It is shown that there is both observational and theoretical evidence which supports the scenario in which the pulsar wind nebula 3C58 has caught up with the reverse shock of the associated (but undetected) supernova remnant.

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.

Using the Very Large Array and the Pie Town antenna, we have measured the position of the radio pulsar B1951+32 relative to nearby background radio sources at four epochs between 1989 and 2000. These data show a clear motion for the pulsar of (25 +/- 4) milliarcsec/yr at a position angle (252 +/- 7) degrees (north through east), corresponding to a transverse velocity (240 +/- 40) km/s for a distance to the source of 2 kpc. The measured direction of motion confirms that the pulsar is moving away from the center of its associated supernova remnant, the first time that such a result has been demonstrated. Independent of assumptions made about the pulsar birth-place, we show that the measured proper motion implies an age for the pulsar of (64 +/- 18) kyr, somewhat less than its characteristic age of 107 kyr. This discrepancy can be explained if the initial spin period of the pulsar was (27 +/- 6) ms.

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 propose a method to infer the initial spin period of pulsars residing in composite supernova remnants. Such a remnant consists of both a plerionic and a shell type component, corresponding respectively to the pulsar wind nebula driven by the spindown luminosity of the central pulsar, and the blastwave bounding the supernova remnant. Theoretical investigations including hydrodynamical simulations have shown that at late times (~ 1,000 - 10,000 years), a simple scaling law connects the radius of the supernova shell to the radius of the plerion. The energy content of the plerion and the total mechanical energy of the supernova remnant enter into this scaling law. One can use this scaling law to estimate the initial spin period of pulsars residing in composite remnants. We discuss potential pitfalls of this method, including the effect of a small remnant age and of strong radiative losses in the plerion.

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].