uvex_transients.models.arnett.get_magnetar_engine#

uvex_transients.models.arnett.get_magnetar_engine(spin_period: PhysicalInput, magnetic_field: PhysicalInput, ns_mass: PhysicalInput = <Quantity 1.4 solMass>) → Callable[[PhysicalInput], FloatArray][source]#

Build a magnetar spin-down energy source for compute_arnett_luminosity().

\[F_\mathrm{mag}(t) = \frac{E_\mathrm{mag}}{t_\mathrm{mag}}\frac{1}{(1 + t/t_\mathrm{mag})^2},\]

with \(E_\mathrm{mag} = 2.6\times10^{52}\,(M_\mathrm{NS}/1.4\,M_\odot)^{3/2}(P/1\,\mathrm{ms})^{-2}\) erg and \(t_\mathrm{mag} = 1.3\times10^{5}\,(M_\mathrm{NS}/1.4\,M_\odot)^{3/2}(P/1\,\mathrm{ms})^{2} (B_\perp/10^{14}\,\mathrm{G})^{-2}\) s [1].

Parameters:

  • spin_period (Quantity or float) – Initial spin period \(P\) (a float is taken to be in s).

  • magnetic_field (Quantity or float) – Perpendicular dipole field \(B_\perp\) (a float is taken to be in G). Quantities may be in u.G, u.T or the Gaussian-cgs base unit GAUSS_CGS.

  • ns_mass (Quantity or float, optional) – Neutron star mass (a float is taken to be in g). Default 1.4 Msun.

Returns:

engine(t) returning the injected power as a bare float64 array in erg/s. t may be a Quantity (a float is taken to be in s).

Return type:

callable

References