BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

L. Del Zanna

L. Del Zanna appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEgr-qcnucl-thhep-phastro-ph.SRphysics.plasm-phastro-ph.IMhep-thmath-phmath.MPphysics.geo-ph

Catalog footprint

What is connected

23works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

23 published item(s)

preprint2020arXiv

Axisymmetric equilibrium models for magnetised neutron stars in Scalar-Tensor Theories

Among the possible extensions of General Relativity that have been put forward in order to address some long standing issues in our understanding of the Universe, Scalar-Tensor Theories have received a lot of attention for their simplicity. Interestingly, some of these predict a potentially observable non-linear phenomenon, known as \textit{spontaneous scalarisation}, in the presence of highly compact matter distributions, like the case of neutron stars. Neutron stars are ideal laboratories to investigate the properties of matter under extreme conditions, and in particular they are known to harbour the strongest magnetic fields in the Universe. Here, for the first time, we present a detailed study of magnetised neutron stars in Scalar-Tensor Theories. First, we show that the formalism developed for the study of magnetised neutron stars in General Relativity, based on the \textit{eXtended Conformally Flat Condition}, can easily be extended in the presence of a non-minimally coupled scalar field, retaining many of its numerical advantages. We then carry out a study of the parameter space considering the two extreme geometries of purely toroidal and purely poloidal magnetic fields, varying both the strength of the magnetic field and the intensity of scalarisation. We compare our results with magnetised general-relativistic solutions and un-magnetised scalarised solutions, showing how the mutual interplay between magnetic and scalar fields affect the magnetic and the scalarisation properties of neutron stars. In particular, we focus our discussion on magnetic deformability, maximum mass and range of scalarisation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

General Relativistic Mean-field Dynamo Model for Proto-neutron Stars

Neutron stars, and magnetars in particular, are known to host the strongest magnetic fields in the Universe. The origin of these strong fields is a matter of controversy. In this preliminary work, via numerical simulations, we study, for the first time in non-ideal general relativistic magnetohydrodynamic (GRMHD) regime, the growth of the magnetic field due to the action of the mean-field dynamo due to sub-scale, unresolved turbulence. The dynamo process, combined with the differential rotation of the (proto-)star, is able to produce an exponential growth of any initial magnetic seed field up to the values required to explain the observations. By varying the dynamo coefficient we obtain different growth rates. We find a quasi-linear dependence of the growth rates on the intensity of the dynamo. Furthermore, the time interval in which exponential growth occurs and the growth rates also seems to depend on the initial configuration of the magnetic field.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Magnetic deformation of neutron stars in scalar-tensor theories

Scalar-tensor theories are among the most promising alternatives to general relativity that have been developed to account for some long standing issues in our understanding of gravity. Some of these theories predict the existence of a non-linear phenomenon, spontaneous scalarisation, which can lead to the appearance of sizeable modifications to general relativity in the presence of compact matter distributions, namely neutron stars. On the one hand, one of the effects of the scalar field is to modify the emission of gravitational waves, both due to variations in the quadrupolar deformation of the star and to the presence of additional modes of emission. On the other hand, neutron stars are known to harbour extremely powerful magnetic fields which can affect their structure and shape, leading in turn to the emission of gravitational waves, this time due to a magnetic quadrupolar deformation. We investigate here how the presence of spontaneous scalarisation can affect the magnetic deformation of neutron stars and their emission of quadrupolar gravitational waves, both of tensor and scalar nature. We will show that it is possible to provide simple parameterisations of the magnetic deformation and gravitational wave power of neutron stars in terms of their baryonic mass, circumferential radius and scalar charge, and that a universal scaling exists independently of the magnetic field geometry and of the parameters of the scalar-tensor theory. Finally, we comment on the observability of the deviations of the gravitational waves strain from general relativity by current and future observatories.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

General relativistic magnetohydrodynamic dynamo in thick accretion disks: fully nonlinear simulations

The recent imaging of the M87 black hole at millimeter wavelengths by the Event Horizon Telescope (EHT) collaboration has triggered a renewed interest in numerical models for the accretion of magnetized plasma in the regime of general relativistic magnetohydrodynamics (GRMHD). Here non-ideal simulations, including both the resistive effects and, above all, the mean-field dynamo action due to sub-scale, unresolved turbulence, are applied for the first time to such systems in the fully nonlinear regime. Combined with the differential rotation of the disk, the dynamo process is able to produce an exponential growth of any initial seed magnetic field up to the values required to explain the observations, when the instability tends to saturate even in the absence of artificial quenching effects. Before reaching the final saturation stage we observe a secondary regime of exponential growing, where the magnetic field increases more slowly due to accretion, which is modifying the underlying equilibrium. By varying the dynamo coefficient we obtain different growth rates, though the field seems to saturate at approximately the same level, at least for the limited range of parameters explored here, providing substantial values for the MAD parameter for magnetized accretion. For reasonable values of the central mass density and the commonly employed recipes for synchrotron emission by relativistically hot electrons, our model is able to reproduce naturally the observed flux of Sgr A*, the next target for EHT.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Activation of MHD reconnection on ideal timescales

Magnetic reconnection in laboratory, space and astrophysical plasmas is often invoked to explain explosive energy release and particle acceleration. However, the timescales involved in classical models within the macroscopic MHD regime are far too slow to match the observations. Here we revisit the tearing instability by performing visco-resistive two-dimensional numerical simulations of the evolution of thin current sheets, for a variety of initial configurations and of values of the Lunquist number $S$, up to $10^7$. Results confirm that when the critical aspect ratio of $S^{1/3}$ is reached in the reconnecting current sheets, the instability proceeds on ideal (Alfvénic) macroscopic timescales, as required to explain observations. Moreover, the same scaling is seen to apply also to the local, secondary reconnection events triggered during the nonlinear phase of the tearing instability, thus accelerating the cascading process to increasingly smaller spatial and temporal scales. The process appears to be robust, as the predicted scaling is measured both in inviscid simulations and when using a Prandtl number $P=1$ in the viscous regime.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Fast reconnection in relativistic plasmas: the magnetohydrodynamics tearing instability revisited

Fast reconnection operating in magnetically dominated plasmas is often invoked in models for magnetar giant flares, for magnetic dissipation in pulsar winds, or to explain the gamma-ray flares observed in the Crab nebula, hence its investigation is of paramount importance in high-energy astrophysics. Here we study, by means of two dimensional numerical simulations, the linear phase and the subsequent nonlinear evolution of the tearing instability within the framework of relativistic resistive magnetohydrodynamics, as appropriate in situations where the Alfven velocity approaches the speed of light. It is found that the linear phase of the instability closely matches the analysis in classical MHD, where the growth rate scales with the Lundquist number S as S^-1/2, with the only exception of an enhanced inertial term due to the thermal and magnetic energy contributions. In addition, when thin current sheets of inverse aspect ratio scaling as S^-1/3 are considered, the so-called "ideal" tearing regime is retrieved, with modes growing independently on S and extremely fast, on only a few light crossing times of the sheet length. The overall growth of fluctuations is seen to solely depend on the value of the background Alfven velocity. In the fully nonlinear stage we observe an inverse cascade towards the fundamental mode, with Petschek-type supersonic jets propagating at the external Alfven speed from the X-point, and a fast reconnection rate at the predicted value R~(ln S)^-1.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Multi-D magnetohydrodynamic modelling of pulsar wind nebulae: recent progress and open questions

In the last decade, the relativistic magnetohydrodynamic (MHD) modelling of pulsar wind nebulae, and of the Crab nebula in particular, has been highly successful, with many of the observed dynamical and emission properties reproduced down to the finest detail. Here, we critically discuss the results of some of the most recent studies: namely the investigation of the origin of the radio emitting particles and the quest for the acceleration sites of particles of different energies along the termination shock, by using wisps motion as a diagnostic tool; the study of the magnetic dissipation process in high magnetization nebulae by means of new long-term three-dimensional simulations of the pulsar wind nebula evolution; the investigation of the relativistic tearing instability in thinning current sheets, leading to fast reconnection events that might be at the origin of the Crab nebula gamma-ray flares.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Quark deconfinement in the proto-magnetar model of Long Gamma-Ray Bursts

We investigate the possible implications of quark deconfinement on the phenomenology of Long Gamma-Ray Bursts focusing, in particular, on the possibility to describe multiple prompt emission phases in the context of the proto-magnetar model. Starting from numerical models of rotating Hadron Stars and Quark Stars in full general relativity we track the electromagnetic spin-down evolution in both the hadronic and quark phase, linking the two families through conservation of baryon number and angular momentum. We give estimates of the timescales and the energetics involved in the spin-down process deriving, in the relevant spin range, the relation between the initial and the final masses and rotational energies, whenever hadron-quark conversion is possible. We show how the results can be used in relevant astrophysical cases such as the double burst GRB 110709B.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The "ideal" tearing mode: theory and resistive MHD simulations

Classical MHD reconnection theories, both the stationary Sweet-Parker model and the tearing instability, are known to provide rates which are too slow to explain the observations. However, a recent analysis has shown that there exists a critical threshold on current sheet's thickness, namely a/L~S^(-1/3), beyond which the tearing modes evolve on fast macroscopic Alfvenic timescales, provided the Lunquist number S is high enough, as invariably found in solar and astrophysical plasmas. Therefore, the classical Sweet-Parker scenario, for which the diffusive region scales as a/L~S^(-1/2) and thus can be up to ~100 times thinner than the critical value, is likely to be never realized in nature, as the current sheet itself disrupts in the elongation process. We present here two-dimensional, compressible, resistive MHD simulations, with S ranging from 10^5 to 10^7, that fully confirm the linear analysis. Moreover, we show that a secondary plasmoid instability always occurs when the same critical scaling is reached on the local, smaller scale, leading to a cascading explosive process, reminiscent of the flaring activity.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A study of vorticity formation in high energy nuclear collisions

We present a quantitative study of vorticity formation in peripheral ultrarelativistic heavy ion collisions at sqrt(s)NN = 200 GeV by using the ECHO-QGP numerical code, implementing relativistic dissipative hydrodynamics in the causal Israel-Stewart framework in 3+1 dimensions with an initial Bjorken flow profile. We consider and discuss different definitions of vorticity which are relevant in relativistic hydrodynamics. After demonstrating the excellent capabilities of our code, which proves to be able to reproduce Gubser flow up to 8 fm/c, we show that, with the initial conditions needed to reproduce the measured directed flow in peripheral collisions corresponding to an average impact parameter b=11.6 fm and with the Bjorken flow profile for a viscous Quark Gluon Plasma with η/s=0.1 fixed, a vorticity of the order of some 10^{-2} c/fm can develop at freezeout. The ensuing polarization of Lambda baryons does not exceed 1.4% at midrapidity. We show that the amount of developed directed flow is sensitive to both the initial angular momentum of the plasma and its viscosity.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Modeling the structure of magnetic fields in Neutron Stars: from the interior to the magnetosphere

The phenomenology of the emission of pulsars and magnetars depends dramatically on the structure and properties of their magnetic field. In particular it is believed that the outbursting and flaring activity observed in AXPs and SRGs is strongly related to their internal magnetic field. Recent observations have moreover shown that charges are present in their magnetospheres supporting the idea that their magnetic field is tightly twisted in the vicinity of the star. In principle these objects offer a unique opportunity to investigate physics in a regime beyond what can be obtained in the laboratory. We will discuss the properties of equilibrium models of magnetized neutron stars, and we will show how internal and external currents can be related. These magnetic field configurations will be discussed considering also their stability, relevant for their origin and possibly connected to events like SNe and GRBs. We will also show what kind of deformations they induce in the star, that could lead to emission of gravitational waves. In the case of a twisted magnetosphere we will show how the amount of twist regulates their general topology. A general formalism based on the simultaneous numerical solution of the general relativistic Grad-Shafranov equation and Einstein equations will be presented.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Axisymmetric equilibrium models for magnetized neutron stars in General Relativity under the Conformally Flat Condition

Extremely magnetized neutron stars with magnetic fields as strong as $\sim 10^{15-16}$ G, or magnetars, have received considerable attention in the last decade due to their identification as a plausible source for Soft Gamma Repeaters and Anomalous X-ray Pulsars. Moreover, this class of compact objects has been proposed as a possible engine capable of powering both Long and Short Gamma-Ray Bursts, if the rotation period in their formation stage is short enough (~1 ms). Such strong fields are expected to induce substantial deformations of the star and thus to produce the emission of gravitational waves. Here we investigate, by means of numerical modeling, axisymmetric static equilibria of polytropic and strongly magnetized stars in full general relativity, within the ideal magneto-hydrodynamic regime. The eXtended Conformally Flat Condition (XCFC) for the metric is assumed, allowing us to employ the techniques introduced for the X-ECHO code [Bucciantini & Del Zanna, 2011, Astron. Astrophys. 528, A101], proven to be accurate, efficient, and stable. The updated XNS code for magnetized neutron star equilibria is made publicly available for the community (see www.arcetri.astro.it/science/ahead/XNS). Several sequences of models are here retrieved, from the purely toroidal (resolving a controversy in the literature) or poloidal cases, to the so-called twisted torus mixed configurations, expected to be dynamically stable, which are solved for the first time in the non-perturbative regime.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dynamo action in thick disks around Kerr black holes: high-order resistive GRMHD simulations

We present the first kinematic study of an $αΩ$-dynamo in the General Relativistic Magneto-HydroDynamics (GRMHD) regime, applied to thick disks orbiting around Kerr black holes and using a fully covariant mean field dynamo closure for the Ohm law. We show that the $αΩ$-dynamo mechanism leads to a continuous exponential growth of the magnetic field within the disk and to the formation of dynamo waves drifting away or toward the equatorial plane. Since the evolution of the magnetic field occurs qualitatively in the same fashion as in the Sun, we present also butterfly diagrams that characterize our models and show the establishment of an additional timescale, which depends on the microscopic properties of the turbulent motions, possibly providing an alternative explanation to periodicities observed in many high-energy astrophysical sources where accretion onto a rotating black hole is believed to operate.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

General relativistic neutron stars with twisted magnetosphere

Soft Gamma-Ray Repeaters and Anomalous X-Ray Pulsars are extreme manifestations of the most magnetized neutron stars: magnetars. The phenomenology of their emission and spectral properties strongly support the idea that the magnetospheres of these astrophysical objects are tightly twisted in the vicinity of the star. Previous studies on equilibrium configurations have so far focused on either the internal or the external magnetic field configuration, without considering a real coupling between the two fields. Here we investigate numerical equilibrium models of magnetized neutron stars endowed with a confined twisted magnetosphere, solving the general relativistic Grad-Shafranov equation both in the interior and in the exterior of the compact object. A comprehensive study of the parameters space is provided to investigate the effects of different current distributions on the overall magnetic field structure.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

High-order schemes for non-ideal $3+1$ GRMHD: a study of the kinematic dynamo process in accretion tori

We present the first astrophysical application of the ECHO code in its recent version supplemented by a generalized Ohm law, namely a kinematic study of dynamo effects in thick accretion disks. High-order \emph{implicit-explicit} Runge-Kutta time-stepping routines are implemented and validated within $3+1$ \emph{General Relativistic MagnetoHydroDynamics} (GRMHD). The scheme is applied to a differentially rotating torus orbiting a Kerr black hole, where the mean-field dynamo process leads to strong amplification of seed magnetic fields. We show that the interplay between the toroidal and poloidal components occurs qualitatively in the same fashion as in the Sun, butterfly diagrams are reproduced, and a typical time-scale for the field evolution is found, depending on the dynamo and resistivity numbers, which could explain periodicities as observed in several accreting systems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Solving the 3+1 GRMHD equations in the eXtended Conformally Flat Condition: the XNS code for magnetized neutron stars

High-energy phenomena in astrophysics involve quite generally a combination of relativistic motions and strong gravity. The simultaneous solution of Einstein equations and General Relativistic MHD equations is thus necessary to model with accuracy such phenomena. The so-called Conformally Flat Condition (CFC) allows a simplified treatment of Einstein equations, that can be particularly efficient in those contexts where gravitational wave emission is negligible, like core-collapse, or the formation/evolution of neutron stars. We have developed a set of codes to model axisymmetric MHD flows, in General Relativity, where the solution of Einstein equations is achieved with a semi-spectral scheme. Here, we will show how this framework is particularly well suited to investigate neutron star equilibrium models in the presence of strong magnetic fields and we will present the XNS code, that has been recently developed and here updated to treat poloidal and mixed configurations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The role of currents distribution in general relativistic equilibria of magnetized neutron stars

Magnetic fields play a critical role in the phenomenology of neutron stars. There is virtually no observable aspect which is not governed by them. Despite this, only recently efforts have been done to model magnetic fields in the correct general relativistic regime, characteristic of these compact objects. In this work we present, for the first time a comprehensive and detailed parameter study, in general relativity, of the role that the current distribution, and the related magnetic field structure, have in determining the precise structure of neutron stars. In particular, we show how the presence of localized currents can modify the field strength at the stellar surface, and we look for general trends, both in terms of energetic properties, and magnetic field configurations. Here we verify that, among other things, for a large class of different current distributions the resulting magnetic configurations are always dominated by the poloidal component of the current.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

General Relativistic Equilibrium Models of Magnetized Neutron Stars

Magnetic fields play a crucial role in many astrophysical scenarios and, in particular, are of paramount importance in the emission mechanism and evolution of Neutron Stars (NSs). To understand the role of the magnetic field in compact objects it is important to obtain, as a first step, accurate equilibrium models for magnetized NSs. Using the conformally flat approximation we solve the Einstein's equations together with the GRMHD equations in the case of a static axisymmetryc NS taking into account different types of magnetic configuration. This allows us to investigate the effect of the magnetic field on global properties of NSs such as their deformation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Relativistic distribution function for particles with spin at local thermodynamical equilibrium

We present an extension of relativistic single-particle distribution function for weakly interacting particles at local thermodynamical equilibrium including spin degrees of freedom, for massive spin 1/2 particles. We infer, on the basis of the global equilibrium case, that at local thermodynamical equilibrium particles acquire a net polarization proportional to the vorticity of the inverse temperature four-vector field. The obtained formula for polarization also implies that a steady gradient of temperature entails a polarization orthogonal to particle momentum. The single-particle distribution function in momentum space extends the so-called Cooper-Frye formula to particles with spin 1/2 and allows to predict their polarization in relativistic heavy ion collisions at the freeze-out.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Relativistic viscous hydrodynamics for heavy-ion collisions with ECHO-QGP

We present ECHO-QGP, a numerical code for $(3+1)$-dimensional relativistic viscous hydrodynamics designed for the modeling of the space-time evolution of the matter created in high energy nuclear collisions. The code has been built on top of the \emph{Eulerian Conservative High-Order} astrophysical code for general relativistic magneto-hydrodynamics [\emph{Del Zanna et al., Astron. Astrophys. 473, 11, 2007}] and here it has been upgraded to handle the physics of the Quark-Gluon Plasma. ECHO-QGP features second-order treatment of causal relativistic viscosity effects in both Minkowskian or Bjorken coordinates; partial or complete chemical equilibrium of hadronic species before kinetic freeze-out; initial conditions based on the optical Glauber model, including a Monte-Carlo routine for event-by-event fluctuating initial conditions; a freeze-out procedure based on the Cooper-Frye prescription. The code is extensively validated against several test problems and results always appear accurate, as guaranteed by the combination of the conservative (shock-capturing) approach and the high-order methods employed. ECHO-QGP can be extended to include evolution of the electromagnetic fields coupled to the plasma.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A fully covariant mean-field dynamo closure for numerical 3+1 resistive GRMHD

The powerful high-energy phenomena typically encountered in astrophysics invariably involve physical engines, like neutron stars and black hole accretion disks, characterized by a combination of highly magnetized plasmas, strong gravitational fields, and relativistic motions. In recent years numerical schemes for General Relativistic MHD (GRMHD) have been developed to model the multidimensional dynamics of such systems, including the possibility of an evolving spacetime. Such schemes have been also extended beyond the ideal limit including the effects of resistivity, in an attempt to model dissipative physical processes acting on small scales (sub-grid effects) over the global dynamics. Along the same lines, magnetic fields could be amplified by the presence of turbulent dynamo processes, as often invoked to explain the high values of magnetization required in accretion disks and neutron stars. Here we present, for the first time, a further extension to include the possibility of a mean-field dynamo action within the framework of numerical 3+1 (resistive) GRMHD. A fully covariant dynamo closure is proposed, in analogy with the classical theory, assuming a simple alpha-effect in the comoving frame. Its implementation into a finite-difference scheme for GRMHD in dynamical spacetimes [the X-ECHO code: (Bucciantini and Del Zanna 2011)] is described, and a set of numerical test is presented and compared with analytical solutions wherever possible.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

GRMHD in axisymmetric dynamical spacetimes: the X-ECHO code

We present a new numerical code, X-ECHO, for general relativistic magnetohydrodynamics (GRMHD) in dynamical spacetimes. This is aimed at studying astrophysical situations where strong gravity and magnetic fields are both supposed to play an important role, such as for the evolution of magnetized neutron stars or for the gravitational collapse of the magnetized rotating cores of massive stars, which is the astrophysical scenario believed to eventually lead to (long) GRB events. The code is based on the extension of the Eulerian conservative high-order (ECHO) scheme [Del Zanna et al., A&A 473, 11 (2007)] for GRMHD, here coupled to a novel solver for the Einstein equations in the extended conformally flat condition (XCFC). We fully exploit the 3+1 Eulerian formalism, so that all the equations are written in terms of familiar 3D vectors and tensors alone, we adopt spherical coordinates for the conformal background metric, and we consider axisymmetric spacetimes and fluid configurations. The GRMHD conservation laws are solved by means of shock-capturing methods within a finite-difference discretization, whereas, on the same numerical grid, the Einstein elliptic equations are treated by resorting to spherical harmonics decomposition and solved, for each harmonic, by inverting band diagonal matrices. As a side product, we build and make available to the community a code to produce GRMHD axisymmetric equilibria for polytropic relativistic stars in the presence of differential rotation and a purely toroidal magnetic field. This uses the same XCFC metric solver of the main code and has been named XNS. Both XNS and the full X-ECHO codes are validated through several tests of astrophysical interest.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Magnetized Relativistic Jets and Long-Duration GRBs from Magnetar Spindown during Core-Collapse Supernovae

We use ideal axisymmetric relativistic magnetohydrodynamic simulations to calculate the spindown 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 star's atmosphere and the wind's interaction with the surrounding star. The magnetar drives a relativistic magnetized wind into a cavity created by the outgoing SN shock. For high spindown powers (~ 10^{51}-10^{52} ergs/s), the magnetar wind is super-fast at almost all latitudes, while for lower spindown powers (~ 10^{50} erg/s), the wind is sub-fast 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% 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 deg.) relativistic jet out along the rotation axis of the star. Nearly all of the spindown 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.

0 citations0 reviewsNo code linkedOpen access