-
The early evolution of magnetar rotation -- II. Rapidly rotating magnetars: Implications for Gamma-Ray Bursts and Super Luminous Supernovae
Authors:
Tejas Prasanna,
Matthew S. B. Coleman,
Matthias J. Raives,
Todd A. Thompson
Abstract:
Rapidly rotating magnetars have been associated with gamma-ray bursts (GRBs) and super-luminous supernovae (SLSNe). Using a suite of 2D magnetohydrodynamic simulations at fixed neutrino luminosity and a couple of evolutionary models with evolving neutrino luminosity and magnetar spin period, we show that magnetars are viable central engines for powering GRBs and SLSNe. We also present analytic est…
▽ More
Rapidly rotating magnetars have been associated with gamma-ray bursts (GRBs) and super-luminous supernovae (SLSNe). Using a suite of 2D magnetohydrodynamic simulations at fixed neutrino luminosity and a couple of evolutionary models with evolving neutrino luminosity and magnetar spin period, we show that magnetars are viable central engines for powering GRBs and SLSNe. We also present analytic estimates of the energy outflow rate from the proto-neutron star (PNS) as a function of polar magnetic field strength $B_0$, PNS angular velocity $Ω_{\star}$, PNS radius $R_{\star}$ and mass outflow rate $\dot{M}$. We show that rapidly rotating magnetars with spin periods $P_{\star}\lesssim 4$ ms and polar magnetic field strength $B_0\gtrsim 10^{15}$ G can release $10^{50}-5\times 10^{51}$ ergs of energy during the first $\sim2$ s of the cooling phase. Based on this result, it is plausible that sustained energy injection by magnetars through the relativistic wind phase can power GRBs. We also show that magnetars with moderate field strengths of $B_0\lesssim 5\times 10^{14}$ G do not release a large fraction of their rotational kinetic energy during the cooling phase and hence, are not likely to power GRBs. Although we cannot simulate to times greater than $\sim 3-5$ s after a supernova, we can hypothesize that moderate field strength magnetars can brighten the supernova light curves by releasing their rotational kinetic energy via magnetic dipole radiation on timescales of days to weeks, since these do not expend most of their rotational kinetic energy during the early cooling phase.
△ Less
Submitted 21 October, 2023; v1 submitted 25 May, 2023;
originally announced May 2023.
-
Magnetized Rotating Isothermal Winds
Authors:
Matthias J. Raives,
Matthew S. B. Coleman,
Todd A. Thompson
Abstract:
We consider the general problem of a Parker-type non-relativistic isothermal wind from a rotating and magnetic star. Using the magnetohydrodynamics (MHD) code athena++, we construct an array of simulations in the stellar rotation rate $Ω_\ast$ and the isothermal sound speed $c_T$, and calculate the mass, angular momentum, and energy loss rates across this parameter space. We also briefly consider…
▽ More
We consider the general problem of a Parker-type non-relativistic isothermal wind from a rotating and magnetic star. Using the magnetohydrodynamics (MHD) code athena++, we construct an array of simulations in the stellar rotation rate $Ω_\ast$ and the isothermal sound speed $c_T$, and calculate the mass, angular momentum, and energy loss rates across this parameter space. We also briefly consider the three dimensional case, with misaligned magnetic and rotation axes.
We discuss applications of our results to the spindown of normal stars, highly-irradiated exoplanets, and to nascent highly-magnetic and rapidly-rotating neutron stars born in massive star core collapse.8
△ Less
Submitted 30 October, 2023; v1 submitted 10 February, 2023;
originally announced February 2023.
-
The Early Evolution of Magnetar Rotation I: Slowly Rotating "Normal" Magnetars
Authors:
Tejas Prasanna,
Matthew S. B. Coleman,
Matthias J. Raives,
Todd A. Thompson
Abstract:
In the seconds following their formation in core-collapse supernovae, "proto"-magnetars drive neutrino-heated magneto-centrifugal winds. Using a suite of two-dimensional axisymmetric MHD simulations, we show that relatively slowly rotating magnetars with initial spin periods of $P_{\star0}=50-500$ ms spin down rapidly during the neutrino Kelvin-Helmholtz cooling epoch. These initial spin periods a…
▽ More
In the seconds following their formation in core-collapse supernovae, "proto"-magnetars drive neutrino-heated magneto-centrifugal winds. Using a suite of two-dimensional axisymmetric MHD simulations, we show that relatively slowly rotating magnetars with initial spin periods of $P_{\star0}=50-500$ ms spin down rapidly during the neutrino Kelvin-Helmholtz cooling epoch. These initial spin periods are representative of those inferred for normal Galactic pulsars, and much slower than those invoked for gamma-ray bursts and super-luminous supernovae. Since the flow is non-relativistic at early times, and because the Alfvén radius is much larger than the proto-magnetar radius, spindown is millions of times more efficient than the typically-used dipole formula. Quasi-periodic plasmoid ejections from the closed zone enhance spindown. For polar magnetic field strengths $B_0\gtrsim5\times10^{14}$ G, the spindown timescale can be shorter than than the Kelvin-Helmholtz timescale. For $B_0\gtrsim10^{15}$ G, it is of order seconds in early phases. We compute the spin evolution for cooling proto-magnetars as a function of $B_0$, $P_{\star0}$, and mass ($M$). Proto-magnetars born with $B_0$ greater than $\simeq1.3\times10^{15}\,{\rm\,G}\,(P_{\star0}/{400\,\rm\,ms})^{-1.4}(M/1.4\,{\rm M}_\odot)^{2.2}$ spin down to periods $> 1$ s in just the first few seconds of evolution, well before the end of the cooling epoch and the onset of classic dipole spindown. Spindown is more efficient for lower $M$ and for larger $P_{\star0}$. We discuss the implications for observed magnetars, including the discrepancy between their characteristic ages and supernova remnant ages. Finally, we speculate on the origin of 1E 161348-5055 in the remnant RCW 103, and the potential for other ultra-slowly rotating magnetars.
△ Less
Submitted 18 August, 2022;
originally announced August 2022.
-
The Antesonic Condition for the Explosion of Core-Collapse Supernovae II: Rotation and Turbulence
Authors:
Matthias J. Raives,
Todd A. Thompson,
Sean M. Couch
Abstract:
In the problem of steady free-fall onto a standing shockwave around acentral mass, the "antesonic" condition limits the regime of stable accretion to $c_T^2/v_\mathrm{esc}^2\leq3/16$, where $c_T$ is the isothermal sound speed in the subsonic post-shock flow, and $v_\mathrm{esc}$ is the escape velocity at the shock radius. Above this limit, it is impossible to satisfy both the Euler equation and th…
▽ More
In the problem of steady free-fall onto a standing shockwave around acentral mass, the "antesonic" condition limits the regime of stable accretion to $c_T^2/v_\mathrm{esc}^2\leq3/16$, where $c_T$ is the isothermal sound speed in the subsonic post-shock flow, and $v_\mathrm{esc}$ is the escape velocity at the shock radius. Above this limit, it is impossible to satisfy both the Euler equation and theshock jump conditions, and the system transitions to a wind. This physics explains the existence of a critical neutrino luminosity in steady-state models ofaccretion in the context of core-collapse supernovae. Here, we extend the antesonic condition to flows with rotation and turbulence using a simple one-dimensional formalism. Both effects decrease the critical post-shock sound speed required for explosion. While quite rapid rotation is required for a significant change to the critical condition, we show that the level of turbulence typically achieved in supernova simulations can greatly impact the critical value of $c_T^2/v_\mathrm{esc}^2$. A core angular velocity corresponding to a millisecond rotation period after contraction of the proto-neutron star results in only a $\sim\!5$ per-cent reduction of the critical curve. In contrast, near-sonic turbulence with specific turbulent kinetic energy $K/c_T^2=0.5-1$, leads to a decrease in the critical value of $c_T^2/v_{\rm esc}^2$ by $\sim\!20-40$ per-cent. This analysis provides a framework for understanding the role of post-shock turbulence in instigating explosions in models that would otherwise fail and helps explain why multi-dimensional simulations explode more easily than their one-dimensional counterparts.
△ Less
Submitted 9 September, 2020;
originally announced September 2020.
-
The Antesonic Condition for the Explosion of Core-Collapse Supernovae I: Spherically Symmetric Polytropic Models: Stability & Wind Emergence
Authors:
Matthias J. Raives,
Sean M. Couch,
Johnny P. Greco,
Ondřej Pejcha,
Todd A. Thompson
Abstract:
Shock revival in core-collapse supernovae (CCSNe) may be due to the neutrino mechanism. While it is known that in a neutrino-powered CCSN, explosion begins when the neutrino luminosity of the proto-neutron star exceeds a critical value, the physics of this condition in time-dependent, multidimensional simulations are not fully understood. \citet{Pejcha2012} found that an `antesonic condition' exis…
▽ More
Shock revival in core-collapse supernovae (CCSNe) may be due to the neutrino mechanism. While it is known that in a neutrino-powered CCSN, explosion begins when the neutrino luminosity of the proto-neutron star exceeds a critical value, the physics of this condition in time-dependent, multidimensional simulations are not fully understood. \citet{Pejcha2012} found that an `antesonic condition' exists for time-steady spherically symmetric models, potentially giving a physical explanation for the critical curve observed in simulations. In this paper, we extend that analysis to time-dependent, spherically symmetric polytropic models. We verify the critical antesonic condition in our simulations, showing that models exceeding it drive transonic winds whereas models below it exhibit steady accretion. In addition, we find that (1) high spatial resolution is needed for accurate determination of the antesonic ratio and shock radius at the critical curve, and that low resolution simulations systematically underpredict these quantities, making explosion more difficult at lower resolution; (2) there is an important physical connection between the critical mass accretion rate at explosion and the mass loss rate of the post-explosion wind: the two are directly proportional at criticality, implying that, at criticality, the wind kinetic power is tied directly to the accretion power; (3) the value of the post-shock adiabatic index $Γ$ has a large effect on the length and time scales of the post-bounce evolution of the explosion larger values of $Γ$ result in a longer transition from the accretion to wind phases.
△ Less
Submitted 8 January, 2019; v1 submitted 8 January, 2018;
originally announced January 2018.
-
Accurate, Meshless Methods for Magneto-Hydrodynamics
Authors:
Philip F. Hopkins,
Matthias J. Raives
Abstract:
Recently, we developed a pair of meshless finite-volume Lagrangian methods for hydrodynamics: the 'meshless finite mass' (MFM) and 'meshless finite volume' (MFV) methods. These capture advantages of both smoothed-particle hydrodynamics (SPH) and adaptive mesh-refinement (AMR) schemes. Here, we extend these to include ideal magneto-hydrodynamics (MHD). The MHD equations are second-order consistent…
▽ More
Recently, we developed a pair of meshless finite-volume Lagrangian methods for hydrodynamics: the 'meshless finite mass' (MFM) and 'meshless finite volume' (MFV) methods. These capture advantages of both smoothed-particle hydrodynamics (SPH) and adaptive mesh-refinement (AMR) schemes. Here, we extend these to include ideal magneto-hydrodynamics (MHD). The MHD equations are second-order consistent and conservative. We augment these with a divergence-cleaning scheme, which maintains div*B~0 to high accuracy. We implement these in the code GIZMO, together with a state-of-the-art implementation of SPH MHD. In every one of a large suite of test problems, the new methods are competitive with moving-mesh and AMR schemes using constrained transport (CT) to ensure div*B=0. They are able to correctly capture the growth and structure of the magneto-rotational instability (MRI), MHD turbulence, and the launching of magnetic jets, in some cases converging more rapidly than AMR codes. Compared to SPH, the MFM/MFV methods exhibit proper convergence at fixed neighbor number, sharper shock capturing, and dramatically reduced noise, div*B errors, and diffusion. Still, 'modern' SPH is able to handle most of our tests, at the cost of much larger kernels and 'by hand' adjustment of artificial diffusion parameters. Compared to AMR, the new meshless methods exhibit enhanced 'grid noise' but reduced advection errors and numerical diffusion, velocity-independent errors, and superior angular momentum conservation and coupling to N-body gravity solvers. As a result they converge more slowly on some problems (involving smooth, slowly-moving flows) but more rapidly on others (involving advection or rotation). In all cases, divergence-control beyond the popular Powell 8-wave approach is necessary, or else all methods we consider will systematically converge to unphysical solutions.
△ Less
Submitted 13 December, 2015; v1 submitted 11 May, 2015;
originally announced May 2015.