L. Del Zanna

High-Order Upwind Schemes for Multidimensional Magnetohydrodynamics

A general method for constructing high-order upwind schemes for multidimensional magnetohydrodynamics (MHD), having as a main built-in condition the divergence-free constraint ∇ = 0 for the magnetic field vector , is proposed. The suggested procedure is based on consistency arguments, by taking into account the specific operator structure of MHD equations with respect to the reference Euler equations of gasdynamics. This approach leads in a natural way to a staggered representation of the field numerical data in which the divergence-free condition in the cell-averaged form, corresponding to second-order accurate numerical derivatives, is exactly fulfilled. To extend this property to higher order schemes, we then give general prescriptions to satisfy a (r + 1)th order accurate ∇ = 0 relation for any numerical field having a rth order interpolation accuracy. Consistency arguments lead also to a proper formulation of the upwind procedures needed to integrate the induction equations, assuring the exact conservation in time of the divergence-free condition and the related continuity properties for the vector components. As an application, a third-order code to simulate multidimensional MHD flows of astrophysical interest is developed using essentially nonoscillatory-based reconstruction algorithms. Several test problems to illustrate and validate the proposed approach are finally presented.

An efficient shock-capturing central-type scheme for multidimensional relativistic flows - I. Hydrodynamics

Multidimensional shock-capturing numerical schemes for special relativistic hydrodynamics (RHD) are compu- tationally more expensive than their correspondent Euler versions, due to the nonlinear relations between conservative and primitive variables and to the consequent complexity of the Jacobian matrices (needed for the spectral decomposition in most of the approximate Riemann solvers of common use). Here an ecient and easy-to-implement three-dimensional (3-D) shock- capturing scheme for ideal RHD is presented. Based on the algorithms developed by P. Londrillo & L. Del Zanna (2000, ApJ, 530, 508) for the non-relativistic magnetohydrodynamic (MHD) case, and having in mind its relativistic MHD extension (to ap- pear in a forthcoming paper), the scheme uses high order (third) Convex Essentially Non-Oscillatory (CENO) finite dierence interpolation routines and central-type averaged Riemann solvers, which do not make use of time-consuming characteristic de- composition. The scheme is very ecient and robust, and it gives results comparable to those obtained with more sophisticated algorithms, even in ultrarelativistic multidimensional test problems.

Axially symmetric relativistic MHD simulations of Pulsar Wind Nebulae in Supernova Remnants - On the origin of torus and jet-like features

The structure and the evolution of Pulsar Wind Nebulae (PWNe) are studied by means of two-dimensional axisym- metric relativistic magnetohydrodynamic (RMHD) simulations. After the first imaging of the Crab Nebula with Chandra ,a growing number of objects has been found to show in the X-rays spatial features such as rings and jets, that clearly cannot be accounted for within the standard framework of one-dimensional semi-analytical models. The most promising explanation suggested so far is based on the combined effects of the latitude dependence of the pulsar wind energy flux, shaping the wind termination shock and naturally providing a higher equatorial emission, and of the wind magnetization, likely responsible for the jet collimation by hoop stresses downstream of the shock. This scenario is investigated here by following the evolution of a PWN interacting with the confining Supernova Remnant (SNR), from the free expansion to the beginning of the reverberation phase. Our results confirm the oblate shape of the wind termination shock and the formation of a polar jet with supersonic velocities (v ≈ 0.5−0.7c) for high enough values of the equatorial wind magnetization parameter (σ > 0.01).

On the divergence-free condition in Godunov-type schemes for ideal magnetohydrodynamics: the upwind constrained transport method

We present a general framework to design Godunov-type schemes for multidimensional ideal magnetohydrodynamic (MHD) systems, having the divergence-free relation and the related properties of the magnetic field B as built-in conditions. Our approach mostly relies on the constrained transport (CT) discretization technique for the magnetic field components, originally developed for the linear induction equation, which assures [Δ ċ B]num = 0 and its preservation in time to within machine accuracy in a finite-volume setting. We show that the CT formalism, when fully exploited, can be used as a general guideline to design the reconstruction procedures of the B vector field, to adapt standard upwind procedures for the momentum and energy equations, avoiding the onset of numerical monopoles of O(1) size. and to formulate approximate Riemann solvers for the induction equation. This general framework will be named here upwind constrained transport (UCT). To demonstrate the versatility of our method, we apply it to a variety of schemes, which are finally validated numerically and compared, a novel Implementation for the MHD case of the second-order Roetype positive: scheme by Lin and Lax [J. Comput. Fluid Dyn. 5 (1996) 133], and both the second- and third-order versions of a central-type MHD scheme presented by Londrillo and Del Zanna [Astrophys. J. 530 (2000) 508], where the basic UCT strategies have been first outlined.

Magnetized relativistic jets and long‐duration GRBs from magnetar spin‐down during core‐collapse supernovae

We use ideal axisymmetric relativistic magnetohydrodynamic simulations to calculate the spin-down of a newly formed millisecond, B ∼ 10 15 G, magnetar and its interaction with the surrounding stellar envelope during a core-collapse supernova (SN) explosion. The mass, angular momentum and rotational energy lost by the neutron star are determined self-consistently given the thermal properties of the cooling neutron stars atmosphere and the winds interaction with the surrounding star. The magnetar drives a relativistic magnetized wind into a cavity created by the outgoing SN shock. For high spin-down powers (∼10 51 ―10 52 erg s ―1 ), the magnetar wind is superfast at almost all latitudes, while for lower spin-down powers (∼10 50 erg s ―1 ), the wind is subfast but still super-Alfvenic. In all cases, the rates at which the neutron star loses mass, angular momentum and energy are very similar to the corresponding free wind values (≤30 per cent differences), in spite of the causal contact between the neutron star and the stellar envelope. In addition, in all cases that we consider, the magnetar drives a collimated (∼5―10°) relativistic jet out along the rotation axis of the star. Nearly all of the spin-down power of the neutron star escapes via this polar jet, rather than being transferred to the more spherical SN explosion. The properties of this relativistic jet and its expected late-time evolution in the magnetar model are broadly consistent with observations of long duration gamma-ray bursts (GRBs) and their associated broad-lined Type Ic SN.

Parametric decay of circularly polarized Alfvén waves: Multidimensional simulations in periodic and open domains

The nonlinear evolution of monochromatic large-amplitude circularly polarized Alfven waves subject to the decay instability is studied via numerical simulations in one, two, and three spatial dimensions. The asymptotic value of the cross helicity depends strongly on the plasma beta: in the low beta case multiple decays are observed, with about half of the energy being transferred to waves propagating in the opposite direction at lower wave numbers, for each saturation step. Correspondingly, the other half of the total transverse energy (kinetic and magnetic) goes into energy carried by the daughter compressive waves and to the associated shock heating. In higher beta conditions we find instead that the cross helicity decreases monotonically with time towards zero, implying an asymptotic balance between inward and outward Alfvenic modes, a feature similar to the observed decrease with distance in the solar wind. Although the instability mainly takes place along the propagation direction, in the two and three-dimensional case a turbulent cascade occurs also transverse to the field. The asymptotic state of density fluctuations appears to be rather isotropic, whereas a slight preferential cascade in the transverse direction is seen in magnetic field spectra. Finally, parametric decay is shown to occur also in a non-periodic domain with open boundaries, when the mother wave is continuously injected from one side. In two and three dimensions a strong transverse filamentation is found at long times, reminiscent of density ray-like features observed in the extended solar corona and pressure-balanced structures found in solar wind data.

Relativistic MHD simulations of pulsar bow-shock nebulae

Pulsars out of their parent SNR directly interact with the ISM producing so called Bow-Shock Pulsar Wind Nebulae, the relativistic equivalents of the heliosphere/heliotail system. These have been directly observed from Radio to X-ray, and are found also associated to TeV halos, with a large variety of morphologies. They offer a unique environment where the pulsar wind can be studied by modelling its interaction with the surrounding ambient medium, in a fashion that is different/complementary from the canonical Plerions. These systems have also been suggested as the possible origin of the positron excess detected by AMS and PAMELA, in contrast to dark matter. I will present results from 3D Relativistic MHD simulations of such nebulae. On top of these simulations we computed the expected emission signatures, the properties of high energy particle escape, the role of current sheets in channeling cosmic rays, the level of turbulence and magnetic amplification, and how they depend on the wind structure and magnetisation.

Simulated synchrotron emission from pulsar wind nebulae

Aims. A complete set of diagnostic tools aimed at producing synthetic synchrotron emissivity, polarization, and spectral index maps from relativistic MHD simulations is presented. As a first application we consider here the case of the emission from Pulsar Wind Nebulae (PWNe). Methods. The proposed method is based on the addition, on top of the basic set of MHD equations, of an extra equation describing the evolution of the maximum energy of the emitting particles. This equation takes into account adiabatic and synchrotron losses along streamlines for the distribution of emitting particles and its formulation is such that it is easily implemented in any numerical scheme for relativistic MHD. Results. Application to the axisymmetric simulations of PWNe, analogous to those described by Del Zanna et al. (2004), allows direct comparison between the numerical results and observations of the inner structure of the Crab Nebula, and similar objects, in the optical and X-ray bands. We are able to match most of the observed features typical of PWNe, like the equatorial torus and the polar jets, with velocities in the correct range, as well as finer emission details, like arcs, rings and the bright knot, that turn out to arise mainly from Doppler boosting effects. Spectral properties appear to be well reproduced too: detailed spectral index maps are produced for the first time and show softening towards the PWN outer borders, whereas spectral breaks appear in integrated spectra. The emission details are found to strongly depend on both the average wind magnetization, here σeff ≈ 0.02, and on the magnetic field shape. Conclusions. Our method, in spite of its simplicity, provides a realistic modeling of synchrotron emission properties, and twodimensional axisymmetric relativistic MHD simulations appear to be well suited to explain the main observational features of PWNe.

Non-thermal emission from relativistic MHD simulations of pulsar wind nebulae: from synchrotron to inverse Compton

Aims. We develop a set of diagnostic tools for synchrotron-emitting sources, presented in a previous paper, to include a computation of inverse-Compton radiation from the same relativistic particles that give rise to the synchrotron emission. For the first time, we then study the gamma-ray emission properties of Pulsar Wind Nebulae, in the context of the axisymmetric jet-torus scenario. Methods. We evolve the relativistic MHD equations and the maximum energy of the emitting particles, including adiabatic and synchrotron losses along streamlines. The particle energy distribution function is split into two components: one corresponds to radio-emitting electrons, which are interpreted to be a relic population that is born at the outburst of the supernova, and the other is associated with a wind population that is continuously accelerated at the termination shock and emits up to the gamma-ray band. The inverse Compton emissivity is calculated using the general Klein-Nishina differential cross-section and three different photon targets for the relativistic particles are considered: the nebular synchrotron photons, photons associated with the far-infrared thermal excess, and the cosmic microwave background. Results. When the method is applied to the simulations that match the optical and X-ray morphology of the Crab Nebula, the overall synchrotron spectrum can only be fitted assuming an excess of injected particles and a steeper power law (E −2.7 ) with respect to previous models. The resulting TeV emission has then the correct shape but is in excess of the data. This is related to the magneticfield structure in the nebula, derived using simulations: in particular, the field is strongly compressed close to the termination shock, but with a lower than expected volume average. The jet-torus structure is also found to be visible clearly in high-resolution gamma-ray synthetic maps. We present a preliminary exploration of time variability in X- and gamma-rays. We find variations with timescales of about 2 years in both bands. The variability observed originates in the strongly time-dependent MHD motions inside the nebula.In this paper we complete the set of diagnostic tools for synchrotron emitting sources presented by Del Zanna et al. (Astron. Astrophys. 453, 621, 2006) with the computation of inverse Compton radiation from the same relativistic particles. Moreover we investigate, for the first time, the gamma-ray emission properties of Pulsar Wind Nebulae in the light of the axisymmetric jet-torus scenario. The method consists in evolving the relativistic MHD equations and the maximum energy of the emitting particles. The particle energy distribution function is split in two components: the radio one connected to a relic population born at the outburst of the supernova and the other associated to the wind population continuously accelerated at the termination shock and emitting up to the gamma-ray band. We consider the general Klein-Nishina cross section and three different photon targets: the nebular synchrotron photons, far-infrared thermal ones and the cosmic microwave background. The overall synchrotron spectrum is fitted assuming an excess of injected particles and a steeper power law with respect to previous models. The TeV emission has the correct shape but is in excess of the data. This is due to the nebular magnetic field structure as obtained by the simulations. The jet-torus morphology is visible in high-resolution gamma-ray synthetic maps too. We present a preliminary exploration of time variability in the X and gamma-ray bands.

Magnetic Rayleigh-Taylor instability for Pulsar Wind Nebulae in expanding Supernova Remnants

We present a numerical investigation of the development of Rayleigh-Taylor instability at the interface between an expanding Pulsar Wind Nebula and its surrounding Supernova Remnant. These systems have long been thought to be naturally subject to this kind of instability, given their expansion behavior and the density jump at the contact discontinuity. High resolution images of the Crab Nebula at optical frequencies show the presence of a complex network of line-emitting filaments protruding inside the synchrotron nebula. These structures are interpreted as the observational evidence that Rayleigh-Taylor instability is in fact at work. The development of this instability in the regime appropriate to describe Supernova Remnant-Pulsar Wind Nebula systems is non-trivial. The conditions at the interface are likely close to the stability threshold, and the inclusion of the nebular magnetic field, which might play an important role in stabilizing the system, is essential to the modeling. If Rayleigh-Taylor features can grow efficiently a mixing layer in the outer portion of the nebula might form where most of the supernova material is confined. When a magnetic field close to equipartition is included we find that the interface is stable, and that even a weaker magnetic field affects substantially the growth and shape of the fingers.