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

Andrew I. MacFadyen

Andrew I. MacFadyen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEastro-phastro-ph.EPgr-qcphysics.flu-dynastro-ph.COastro-ph.GAphysics.comp-phastro-ph.IM

Catalog footprint

What is connected

21works
9topics
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

21 published item(s)

preprint2016arXiv

A transition in circumbinary accretion discs at a binary mass ratio of 1:25

We study circumbinary accretion discs in the framework of the restricted three-body problem (R3Bp) and via numerically solving the height-integrated equations of viscous hydrodynamics. Varying the mass ratio of the binary, we find a pronounced change in the behaviour of the disc near mass ratio $q \equiv M_s/M_p \sim 0.04$. For mass ratios above $q=0.04$, solutions for the hydrodynamic flow transition from steady, to strongly-fluctuating; a narrow annular gap in the surface density around the secondary's orbit changes to a hollow central cavity; and a spatial symmetry is lost, resulting in a lopsided disc. This phase transition is coincident with the mass ratio above which stable orbits do not exist around the L4 and L5 equilibrium points of the R3B problem. Using the DISCO code, we find that for thin discs, for which a gap or cavity can remain open, the mass ratio of the transition is relatively insensitive to disc viscosity and pressure. The $q=0.04$ transition has relevance for the evolution of massive black hole binary+disc systems at the centers of galactic nuclei, as well as for young stellar binaries and possibly planets around brown dwarfs.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Narrow Short-Duration GRB Jet from a Wide Central Engine

We use two-dimensional relativistic hydrodynamic numerical calculations to show that highly collimated relativistic jets can be produced in neutron star merger models of short-duration gamma ray bursts without the need for a highly directed engine or a large net magnetic flux. Even a hydrodynamic engine generating a very wide sustained outflow on small scales can in principle produce a highly collimated relativistic jet, facilitated by a dense surrounding medium which provides a cocoon surrounding the jet core. An oblate geometry to the surrounding gas significantly enhances the collimation process. Previous numerical simulations have shown that the merger of two neutron stars produces an oblate, expanding cloud of dynamical ejecta. We show that this gas can efficiently collimate the central engine power much as the surrounding star does in long-duration GRB models. For typical short-duration GRB central engine parameters, we find jets with opening angles of order 10 degrees in which a large fraction of the total outflow power of the central engine resides in highly relativistic material. These results predict large differences in the opening angles of outflows from binary neutron star mergers versus neutron star-black hole mergers.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A reduced orbital period for the supermassive black hole binary candidate in the quasar PG 1302-102?

Graham et al. (2015) have detected a 5.2 year periodic optical variability of the quasar PG 1302-102 at redshift $z=0.3$, which they interpret as the redshifted orbital period $(1+z)t_{\rm bin}$ of a putative supermassive black hole binary (SMBHB). Here we consider the implications of a $3-8$ times shorter orbital period, suggested by hydrodynamical simulations of circumbinary discs (CBDs) with nearly equal--mass SMBHBs ($q\equiv M_2/M_1\gtrsim 0.3$). With the corresponding $2-4$ times tighter binary separation, PG 1302 would be undergoing gravitational wave dominated inspiral, and serve as a proof that the BHs can be fueled and produce bright emission even in this late stage of the merger. The expected fraction of binaries with the shorter $t_{\rm bin}$, among bright quasars, would be reduced by 1-2 orders of magnitude, compared to the 5.2 year period, in better agreement with the rarity of candidates reported by Graham et al. (2015). Finally, shorter periods would imply higher binary speeds, possibly imprinting periodicity on the light curves from relativistic beaming, as well as measurable relativistic effects on the Fe K$α$ line. The CBD model predicts additional periodic variability on timescales of $t_{\rm bin}$ and $\approx 0.5 t_{\rm bin}$, as well as periodic variation of broad line widths and offsets relative to the narrow lines, which are consistent with the observations. Future observations will be able to test these predictions and hence the binary+CBD hypothesis for PG 1302.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

From Engine to Afterglow: Collapsars Naturally Produce Top-Heavy Jets and Early-Time Plateaus in Gamma Ray Burst Afterglows

We demonstrate that the steep decay and long plateau in the early phases of gamma ray burst (GRB) X-ray afterglows are naturally produced in the collapsar model, by a means ultimately related to the dynamics of relativistic jet propagation through a massive star. We present two-dimensional axisymmetric hydrodynamical simulations which start from a collapsar engine and evolve all the way through the late afterglow phase. The resultant outflow includes a jet core which is highly relativistic after breaking out of the star, but becomes baryon-loaded after colliding with a massive outer shell, corresponding to mass from the stellar atmosphere of the progenitor star which became trapped in front of the jet core at breakout. The prompt emission produced before or during this collision would then have the signature of a high Lorentz factor jet, but the afterglow is produced by the amalgamated post-collision ejecta which has more inertia than the original highly relativistic jet core and thus has a delayed deceleration. This naturally explains the early light curve behavior discovered by Swift, including a steep decay and a long plateau, without invoking late-time energy injection from the central engine. The numerical simulation is performed continuously from engine to afterglow, covering a dynamic range of over ten orders of magnitude in radius. Light curves calculated from the numerical output demonstrate that this mechanism reproduces basic features seen in early afterglow data. Initial steep decays are produced by internal shocks, and the plateau corresponds to the coasting phase of the outflow.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

High-Frequency Voronoi Noise Reduced by Smoothed Mesh Motion

We describe a technique for improving the performance of hydrodynamics codes which employ a moving Voronoi mesh. Currently, such codes are susceptible to high-frequency noise produced by rapid adjustments in the grid topology on the smallest scales. The treatment for this grid noise is simple; instead of moving the mesh-generating marker points with the local fluid velocity, this velocity field is smoothed on small scales, so that neighboring marker points generally have similar velocities. We demonstrate significant improvement gained by this adjustment in several code tests relevant to the physics which moving-mesh codes are designed to capture.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Producing Magnetar Magnetic Fields in the Merger of Binary Neutron Stars

The merger of binary neutron stars (BNSs) can lead to large amplifications of the magnetic field due to the development of turbulence and instabilities in the fluid, such as the Kelvin-Helmholtz shear instability, which drive small-scale dynamo activity. In order to properly resolve such instabilities and obtain the correct magnetic field amplification, one would need to employ resolutions that are currently unfeasible in global general relativistic magnetohydrodynamic (GRMHD) simulations of BNS mergers. Here, we present a subgrid model that allows global simulations to take into account the small-scale amplification of the magnetic field which is caused by the development of turbulence during BNS mergers. Assuming dynamo saturation, we show that magnetar-level fields ($\sim 10^{16}\,{\rm G}$) can be easily reached, and should therefore be expected from the merger of magnetized BNSs. The total magnetic energy can reach values up to $\sim 10^{51}\,{\rm erg}$ and the post-merger remnant can therefore emit strong electromagnetic signals and possibly produce short gamma-ray bursts.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Binary Black Hole Accretion During Inspiral and Merger

We present the results of 2D, moving mesh, viscous hydrodynamical simulations of accretion onto merging supermassive black hole (SMBH) binaries. We include viscous heating, shock heating, and radiative cooling, and simulate the transition from the "pre-decoupling" epoch, where the inspiral timescale is longer than the viscous timescale, to the "post-decoupling" epoch, where the inspiral timescale is shorter than the viscous timescale. We find that there is no abrupt halt to the accretion at decoupling, but rather the accretion shows a slow decay, with significant accretion well after the expected decoupling. Moreover, we find that the luminosity in X-rays is significantly higher prior to the merger, as orbital energy from the SMBH binary is converted to heat via strong shocks inside the cavity, and radiated away. Following the merger, the cavity refills viscously and the accretion rate relaxes to the Shakura-Sunyaev value, while the X-ray luminosity drops as the shocks quickly dissipate.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Characteristic Signatures in the Thermal Emission from Accreting Binary Black Holes

We present the results of a calculation of the thermal spectrum from a 2D, moving mesh, high-accuracy, viscous hydrodynamical simulation of an accreting supermassive black hole binary. We include viscous heating, shock heating, and radiative cooling, evolving for longer than a viscous time so that we reach a quasi-steady accretion state. In agreement with previous work, we find that gas is efficiently stripped from the inner edge of the circumbinary disk and enters the cavity along accretion streams, which feed persistent "mini-disks" surrounding each black hole. We also find that emission from the shock-heated mini-disks and accretion streams prevents any deficit in high-energy emission that may be expected inside the circumbinary cavity, and instead leads to a characteristic brightening of the spectrum beginning in soft X-rays.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Migration of Gap-Opening Planets is not Locked to Viscous Disk Evolution

Most standard descriptions of Type II migration state that massive, gap-opening planets must migrate at the viscous drift rate. This is based on the idea that the disk is separated into an inner and outer region and gas is considered unable to cross the gap. In fact, gas easily crosses the gap on horseshoe orbits, nullifying this necessary premise which would set the migration rate. In this work, it is demonstrated using highly accurate numerical calculations that the actual migration rate is dependent on disk and planet parameters, and can be significantly larger or smaller than the viscous drift rate. In the limiting case of a disk much more massive than the secondary, the migration rate saturates to a constant which is sensitive to disk parameters and is not necessarily of order viscous rate. In the opposite limit of a low-mass disk, the migration rate decreases linearly with disk mass. Steady-state solutions in the low disk mass limit show no pile-up outside the secondary's orbit, and no corresponding drainage of the inner disk.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A "Boosted Fireball" Model for Structured Relativistic Jets

We present a model for relativistic jets which generates a particular angular distribution of Lorentz factor and energy per solid angle. We consider a fireball with specific internal energy E/M launched with bulk Lorentz factor γ_B. This "boosted fireball" model is motivated by the phenomenology of collapsar jets, but is applicable to a wide variety of relativistic flows. In its center-of-momentum frame the fireball expands isotropically, converting its internal energy into radially expanding flow with asymptotic Lorentz factor η_0 ~ E/M. In the lab frame the flow is beamed, expanding with Lorentz factor Γ= 2 η_0 γ_B in the direction of its initial bulk motion and with characteristic opening angle θ_0 ~ 1/γ_B. The flow is jet-like with Γθ_0 ~ 2 η_0 such that jets with Γ> 1/θ_0 are naturally produced. The choice η_0 ~ γ_B ~ 10 yields a jet with Γ~ 200 on-axis and angular structure characterized by opening angle θ_0 ~ 0.1 of relevance for cosmological GRBs, while γ_B >~ 1 may be relevant for low-luminosity GRBs. The model produces a family of outflows, of relevance for different relativistic phenomena with structures completely determined by η_0 and γ_B. We calculate the energy per unit solid angle for the model and use it to compute light curves for comparison with the widely used top-hat model. The jet break in the boosted fireball light curve is greatly subdued when compared to the top-hat model because the edge of the jet is smoother than for a top-hat. This may explain missing jet breaks in afterglow light curves.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Binary Black Hole Accretion From a Circumbinary Disk: Gas Dynamics Inside the Central Cavity

We present the results of 2D hydrodynamical simulations of circumbinary disk accretion using the finite-volume code DISCO. This code solves the 2D viscous Navier-Stokes equations on a high-resolution moving mesh which shears with the fluid flow, greatly reducing advection errors in comparison with a fixed grid. We perform a series of simulations for binary mass ratios in the range 0.026 < q < 1.0, each lasting longer than a viscous time so that we reach a quasi-steady accretion state. In each case, we find that gas is efficiently stripped from the inner edge of the circumbinary disk and enters the cavity along accretion streams, which feed persistent "mini-disks" surrounding each black hole. We find that for q > 0.1, the binary excites eccentricity in the inner region of the circumbinary disk, creating an overdense lump which gives rise to enhanced periodicity in the accretion rate. The dependence of the periodicity on mass ratio may provide a method for observationally inferring mass ratios from measurements of the accretion rate. We also find that for all mass ratios studied, the magnitude of the accretion onto the secondary is sufficient to drive the binary toward larger mass ratio. This suggests a mechanism for biasing mass ratio distributions toward equal mass.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Hyperbolicity of Force-Free Electrodynamics

We analyze the equations of relativistic magnetized plasma dynamics in the limiting case that electromagnetic stress-energy is dominant over pressure and rest mass energy density. The naive formulation of these equations is shown to be not hyperbolic. Modifying the equations by terms that vanish for all physical solutions, we obtain a symmetric hyperbolic evolution system, which should exhibit improved numerical behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Magnetic energy production by turbulence in binary neutron star mergers

The simultaneous detection of electromagnetic and gravitational wave emission from merging neutron star binaries would aid greatly in their discovery and interpretation. By studying turbulent amplification of magnetic fields in local high-resolution simulations of neutron star merger conditions, we demonstrate that magnetar-level (~10^16) G fields are present throughout the merger duration. We find that the small-scale turbulent dynamo converts 60% of the randomized kinetic energy into magnetic fields on a merger time scale. Since turbulent magnetic energy dissipates through reconnection events which accelerate relativistic electrons, turbulence may facilitate the conversion of orbital kinetic energy into radiation. If 10^-4 of the ~ 10^53 erg of orbital kinetic available gets processed through reconnection, and creates radiation in the 15-150 keV band, then the fluence at 200 Mpc would be 10^-7 erg/cm^2, potentially rendering most merging neutron stars in the advanced LIGO and Virgo detection volumes detectable by Swift BAT.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Rayleigh-Taylor Instability in a Relativistic Fireball on a Moving Computational Grid

We numerically calculate the growth and saturation of the Rayleigh-Taylor instability caused by the deceleration of relativistic outflows with Lorentz factor Γ = 10, 30, and 100. The instability generates turbulence whose scale exhibits strong dependence on Lorentz factor, as only modes with angular size smaller than 1/Γ can grow. We develop a simple diagnostic to measure the kinetic energy in turbulent fluctuations, and calculate a ratio of turbulent kinetic energy to thermal energy of .03 in the region affected by the instability. Although our numerical calculation does not include magnetic fields, we argue that small scale turbulent dynamo amplifies magnetic fields to nearly this same fraction, giving a ratio of magnetic to thermal energy of ~ .01, to within a factor of two. The instability completely disrupts the contact discontinuity between the ejecta and the swept up circumburst medium. The reverse shock is stable, but is impacted by the Rayleigh-Taylor instability, which strengthens the reverse shock and pushes it away from the forward shock. The forward shock front is unaffected by the instability, but Rayleigh-Taylor fingers can penetrate of order 10% of the way into the energetic region behind the shock during the two-shock phase of the explosion. We calculate afterglow emission from the explosion and find the reverse shock emission peaks at a later time due to its reduced Lorentz factor and modified density and pressure at the shock front. These calculations are performed using a novel numerical technique that includes a moving computational grid. The moving grid is essential as it maintains contact discontinuities to high precision and can easily evolve flows with extremely large Lorentz factors.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Global Calculations of Density Waves and Gap Formation in Protoplanetary Disks using a Moving Mesh

We calculate the global quasi-steady state of a thin disk perturbed by a low-mass protoplanet orbiting at a fixed radius using extremely high-resolution numerical integrations of Euler's equations in two dimensions. The calculations are carried out using a moving computational domain, which greatly reduces advection errors and allows for much longer time-steps than a fixed grid. We calculate the angular momentum flux and the torque density as a function of radius and compare them with analytical predictions. We discuss the quasi-steady state after 100 orbits and the prospects for gap formation by low mass planets.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Observational implications of gamma-ray burst afterglow jet simulations and numerical light curve calculations

We discuss jet dynamics for narrow and wide gamma-ray burst (GRB) afterglow jets and the observational implications of numerical simulations of relativistic jets in two dimensions. We confirm earlier numerical results that sideways expansion of relativistic jets during the bulk of the afterglow emission phase is logarithmic in time and find that this also applies to narrow jets with half opening angle of 0.05 radians. As a result, afterglow jets remain highly nonspherical until after they have become nonrelativistic. Although sideways expansion steepens the afterglow light curve after the jet break, the jet edges becoming visible dominates the jet break, which means that the jet break is sensitive to the observer angle even for narrow jets. Failure to take the observer angle into account can lead to an overestimation of the jet energy by up to a factor 4. This weakens the challenge posed to the magneter energy limit by extreme events such as GRB090926A. Late time radio calorimetry based on a spherical nonrelativistic outflow model remains relevant when the observer is approximately on-axis and where differences of a few in flux level between the model and the simulation are acceptable. However, this does not imply sphericity of the outflow and therefore does not translate to high observer angles relevant to orphan afterglows. For more accurate calorimetry and in order to model significant late time features such as the rise of the counterjet, detailed jet simulations remain indispensable.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

TESS: A Relativistic Hydrodynamics Code on a Moving Voronoi Mesh

We have generalized a method for the numerical solution of hyperbolic systems of equations using a dynamic Voronoi tessellation of the computational domain. The Voronoi tessellation is used to generate moving computational meshes for the solution of multi-dimensional systems of conservation laws in finite-volume form. The mesh generating points are free to move with arbitrary velocity, with the choice of zero velocity resulting in an Eulerian formulation. Moving the points at the local fluid velocity makes the formulation effectively Lagrangian. We have written the TESS code to solve the equations of compressible hydrodynamics and magnetohydrodynamics for both relativistic and non-relativistic fluids on a dynamic Voronoi mesh. When run in Lagrangian mode, TESS is significantly less diffusive than fixed mesh codes and thus preserves contact discontinuities to high precision while also accurately capturing strong shock waves. TESS is written for Cartesian, spherical and cylindrical coordinates and is modular so that auxilliary physics solvers are readily integrated into the TESS framework and so that the TESS framework can be readily adapted to solve general systems of equations. We present results from a series of test problems to demonstrate the performance of TESS and to highlight some of the advantages of the dynamic tessellation method for solving challenging problems in astrophysical fluid dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

The Hydrodynamics of Gamma-Ray Burst Remnants

This paper reports on the results of a numerical investigation designed to address how the initially anisotropic appearance of a GRB remnant is modified by the character of the circumburst medium and by the possible presence of an accompanying supernova (SN). Axisymmetric hydrodynamical calculations of light, impulsive jets propagating in both uniform and inhomogeneous external media are presented, which show that the resulting dynamics of their remnants since the onset of the non-relativistic phase is different from the standard self-similar solutions. Because massive star progenitors are expected to have their close-in surroundings modified by the progenitor winds, we consider both free winds and shocked winds as possible external media for GRB remnant evolution. Abundant confirmation is provided here of the important notion that the morphology and visibility of GRB remnants are determined largely by their circumstellar environments. For this reason, their detectability is highly biased in favor of those with massive star progenitors; although, in this class of models, the beamed component may be difficult to identify because the GRB ejecta is eventually swept up by the accompanying SN. The number density of asymmetric GRB remnants in the local Universe could be, however, far larger if they expand in a tenuous interstellar medium, as expected for some short GRB progenitor models. In these sources, the late size of the observable, asymmetric remnant could extend over a wide, possibly resolvable angle and may be easier to constrain directly.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

An Eccentric Circumbinary Accretion Disk and the Detection of Binary Massive Black Holes

We present a two-dimensional grid-based hydrodynamic simulation of a thin, viscous, locally-isothermal corotating disk orbiting an equal-mass Newtonian binary point mass on a fixed circular orbit. We study the structure of the disk after multiple viscous times. The binary maintains a central hole in the viscously-relaxed disk with radius equal to about twice the binary semimajor axis. Disk surface density within the hole is reduced by orders of magnitude relative to the density in the disk bulk. The inner truncation of the disk resembles the clearing of a gap in a protoplanetary disk. An initially circular disk becomes elliptical and then eccentric. Disturbances in the disk contain a component that is stationary in the rotating frame in which the binary is at rest; this component is a two-armed spiral density wave. We measure the distribution of the binary torque in the disk and find that the strongest positive torque is exerted inside the central low-density hole. We make connection with the linear theory of disk forcing at outer Lindblad resonances (OLRs) and find that the measured torque density distribution is consistent with forcing at the 3:2 (m=2) OLR, well within the central hole. We also measure the time dependence of the rate at which gas accretes across the hole and find quasi-periodic structure. We discuss implications for variability and detection of active galactic nuclei containing a binary massive black hole.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Ultra-relativistic geometrical shock dynamics and vorticity

Geometrical shock dynamics, also called CCW theory, yields approximate equations for shock propagation in which only the conditions at the shock appear explicitly; the post-shock flow is presumed approximately uniform and enters implicitly via a Riemann invariant. The nonrelativistic theory, formulated by G. B. Whitham and others, matches many experimental results surprisingly well. Motivated by astrophysical applications, we adapt the theory to ultra-relativistic shocks advancing into an ideal fluid whose pressure is negligible ahead of the shock, but one third of its proper energy density behind the shock. Exact results are recovered for some self-similar cylindrical and spherical shocks with power-law pre-shock density profiles. Comparison is made with numerical solutions of the full hydrodynamic equations. We review relativistic vorticity and circulation. In an ultrarelativistic ideal fluid, circulation can be defined so that it changes only at shocks, notwithstanding entropy gradients in smooth parts of the flow.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Stellar Explosions by Magnetic Towers

We propose a magnetic mechanism for the collimated explosion of a massive star relevant for GRBs, XRFs and asymmetric supernovae. We apply Lynden-Bell's magnetic tower scenario to the interior of a massive rotating star after the core has collapsed to form a black hole with an accretion disk or a millisecond magnetar acting as a central engine. We solve the force-free Grad-Shafranov equation to calculate the magnetic structure and growth of a tower embedded in a stellar environment. The pressure of the toroidal magnetic field, continuously generated by differential rotation of the central engine, drives a rapid expansion which becomes vertically collimated after lateral force balance with the surrounding gas pressure is reached. The collimation naturally occurs because hoop stress concentrates magnetic field toward the rotation axis and inhibits lateral expansion. This leads to the growth of a self-collimated magnetic tower. When embedded in a massive star, the supersonic expansion of the tower drives a strong bow shock behind which an over-pressured cocoon forms. The cocoon confines the tower by supplying collimating pressure and provides stabilization against disruption due to MHD instabilities. Because the tower consists of closed field lines starting and ending on the central engine, mixing of baryons from the cocoon into the tower is suppressed. The channel cleared by the growing tower is thus plausibly free of baryons and allows the escape of magnetic energy from the central engine through the star. While propagating down the stellar density gradient, the tower accelerates and becomes relativistic. During the expansion, fast collisionless reconnection becomes possible resulting in dissipation of magnetic energy which may be responsible for GRB prompt emission.

0 citations0 reviewsNo code linkedOpen access