Pulsar Wind Nebulae

Rotation powered pulsars lose most of their energy to the production of a relativistic magnetized wind mostly made of electrons and positrons. PWNe are bubbles of relativistic plasma that arise from the confinement of such a wind by the surrounding medium (either the remnant of the parent supernova or the ISM). These sources generally show a broad non-thermal emission spectrum, shining from radio to γ-ray frequencies. The emission is interpreted as mostly due to synchrotron and Inverse Compton radiation from relativistic particles interacting with the ambient magnetic and photon field. The highly relativistic outflow emanating from the pulsar magnetosphere is slowed down at a termination shock where particles are accelerated with extreme efficiency. A number of properties put PWNe among the most intriguing and interesting sources in the sky: they are relativistic sources, likely hosting the most relativistic outflows in Nature (the pulsar winds, with Lorentz factors possibly reaching 107-108); they are the most efficient accelerators in the Galaxy, with efficiencies reaching 20-30%; the only firmly detected Galactic PeVatrons; and the only well established reservoirs of leptonic anti-matter in the Galaxy. In summary, PWNe are sources of chief interest for High Energy Astrophysics, being the nearest and brightest relativistic sources, excellent cosmic accelerators and potentially the primary contributors of cosmic ray positrons, at least above some energy.

Our group has been deeply involved in the study of PWNe, touching virtually all aspects of their physics. We were one of the first two groups to perform 2D relativistic MHD simulations of a pulsar wind interacting with the surrounding SNR. These simulations allowed us to show how the anisotropic energy flux of the pulsar wind, higher at the pulsar rotational equator than at the poles, could explain the findings by Chandra, whose high-resolution X-ray observations had just highlighted the existence of an axisymmetric "jet-torus" morphology in the Crab Nebula and then several other young PWNe. The assumed wind structure was that predicted by theoretical modeling of the wind launching by the pulsar magnetosphere, so that our results also provided independent confirmation of the correct description of the pulsar wind energy flow, including vanishing magnetization in the equatorial region. The overall magnetization of the wind is found instead to be larger than estimated based on 1D modeling: only sufficiently magnetized outflows can lead to substantial hoop stresses downstream of the termination shock so as to divert the post-shock flow towards the axis and give rise to jets. 

Our further studies including radiation diagnostics opened then the way to some other puzzles: within axisymmetric MHD simulations, the morphology and emission spectrum of the Crab Nebula, the class prototype and one of the astrophysical objects for which data are most abundant, could not be simultaneously reproduced with a given value of the wind magnetization. More specifically the magnetization that could best reproduce the X-ray emission morphology would give origin to a nebular magnetic field well below the value inferred from multi-wavelength spectral modeling.

A possible solution to this puzzle, that is currently being investigated, is that the discrepancy be due to the reduced dimensionality of our simulations. Axisymmetry suppresses the development of kink instabilities that could instead be very important in these systems hosting mainly toroidal fields: if kinks can tangle the field lines without substantial magnetic dissipation then a larger magnetic field strength can be accommodated in the nebula without increasing hoop stresses beyond the level that best fits the morphology. This idea is currently being investigated with the help of 3-D relativistic MHD simulations, whose preliminary results are very promising. 

Following up on our early work on these objects, we are currently performing new 2 and 3D relativistic MHD simulations aimed at clarifying the impact of instabilities and the physics of particle escape from these systems. Our final goal is that of assessing more reliably the PWN contribution to the measured leptonic antimatter in Cosmic Rays and their relevance for the production of high energy gamma-ray haloes such as those recently observed by HAWC. This is again very important in view of the advent of CTA.