uvex_transients.models.arnett#

Utilities for computing Arnett-style (diffusion-limited) bolometric lightcurves.

Implements the semi-analytic model of Nicholl et al.[1] (Sec. 2; the MOSFiT magnetar model), building on Arnett (1982) and the gamma-ray leakage treatment of Wang et al.[2]: an internal energy-injection rate \(F_\mathrm{in}(t)\) is diffused through a homologously expanding, grey-opacity ejecta,

\[L(t) = e^{-(t/t_d)^2}\left(1 - e^{-A/t^2}\right) \int_0^t 2 F_\mathrm{in}(t')\,\frac{t'}{t_d}\,e^{(t'/t_d)^2}\,\frac{dt'}{t_d},\]

with diffusion time \(t_d = \sqrt{2\kappa M_\mathrm{ej} / (\beta c v_\mathrm{ej})}\) and high-energy leakage parameter \(A = 3\kappa_\gamma M_\mathrm{ej} / (4\pi v_\mathrm{ej}^2)\).

Everything private here works on bare CGS floats/arrays. The two public functions, get_magnetar_engine() and compute_arnett_luminosity(), accept unit-aware Quantity inputs (unitless inputs are assumed to already be in CGS) and convert at the boundary.

References

Attributes

GAUSS_CGS

The Gaussian-cgs gauss (numerically 1 G).

Functions

compute_arnett_luminosity(t_eval, ...)

Bolometric luminosity of an energy source diffusing through expanding ejecta.

get_magnetar_engine(spin_period, ...)

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

get_nickel_engine(nickel_mass, nickel_yield, ...)

Construct a radioactive Ni-56/Co-56 heating engine.

Classes

ArnettDecaySED(**overrides)

Radioactively-powered (Ni-56/Co-56 decay) SED, with the light curve from Arnett-style diffusion.

ArnettMagnetarSpindownSED(**overrides)

Magnetar spin-down powered SED, with the light curve from Arnett-style diffusion.