Neutron stars typically have magnetic field strengths of $\sim 10^{12} \rm \ G$. However, if a neutron star’s rotation period is comparable to the convective overturn time, magnetic fields can be amplified by helical motion in a mean field dynamo. These highly magnetized neutron stars, magnetars, are born with short periods of $\sim 1 \ \rm ms$, which allow them to support an efficient α − Ω dynamo, resulting in large magnetic dipole fields of $10^{14}-10^{15} \rm \ G$ . A magnetar of mass 1.4M⊙, $R = 10 \ \rm km$, and $P = 1 \ \rm ms$ has a rotational energy of E_{\rm rot} = I \omega^2 /2 = {5}MR^2 \left({P}\right)^2 \simeq 4 \times 10^{51} \ \rm erg The strong magnetic stresses and torques damped the rotation and release a large fraction of $E_{\rm rot}$ on the time scale of \tau \simeq 0.6 B_{15}^{-2} (P/1 \ {\rm ms})^2 \ \rm hr . This amount of energy is similar to that released in γ-ray bursts and superluminous supernovae (SLSNe), thus making magnetar an attractive candidate for the central engine that power these explosions. In the past two decades, a connection has been made between supernovae and at least some subclass of long-duration γ-ray bursts (LGRBs), supporting the case that they are powered by the same mechanism. Around the same time that magnetars were proposed, a realization had emerged that GRBs could have cosmological origins. This would require a total energy of $\sim 10^{51} \ \rm erg$ to explain the observed flux assuming isotropic emission (e.g. ). Interestingly, the first proposal of highly magnetized neutron stars as central engines of these bursts, , was published only a day after . The proposed scenario is such that a magnetar forms from a white dwarf via accretion induced collapse (AIC). The WD magnetic field of $\sim 10^9 \ \rm G$ is amplified to $10^{15} \rm \ G$ by magnetic flux conservation. The rotational period of $\sim 1 \ \rm ms$ is a result of angular momentum conservation. [1] The newly formed neutron star then loses its rotational energy quickly due to electromagnetic torque, generating electric fields that accelerate particles to ultra-relativistic energies, which eventually give out γ-ray. The timescale of energy release due to magnetic dipole luminosity and gravitational wave emission of $\sim 20 \ \rm s$ for a typical magnetar is consistent with the timescale of long-duration GRBs. Near the magnetar’s surface out to the light cylinder, the optical depth to this radiation due to Compton scattering, absorption, and pair-cration is large. The radiation has to propagate out to a photosphere radius of $\sim 10^8 \ \rm cm$ before it is released. The typical radiated γ-ray energies of $0.1-1 \ \rm MeV$ is also consistent with those observed from GRBs . Recent works have shown that while magnetars formed via AIC may be responsible for some GRBs, magnetars formed in the core collapse of massive stars like those proposed by are probably more prevalent. The modern view of the mechanism for core-collapse GRB as laid out by is as following. Shortly after the core bounce, a non-relativistic wind heated by neutrino blows through the cavity carved out by the supernova (SN) shock into a _bipolar jet_. The relativistic jet from the newly formed magnetar follows, and emerges as a GRB prompt emission. After $30-100 \ \rm s$, the maximum Lorantz factor increases to σ₀ ≫ 1 rendering magnetic dissipation and jet acceleration ineffective. This ended the prompt GRB within the observed timescale of $\sim 20 \ \rm s$. After the prompt emission ended, the spin-down of the central magnetar continues to power the GRB into its X-ray plateau phase with a correlation between the plateau luminosity and duration (LT correlation). The observed correlation is given by $ \log L_{\rm X} = a + b \log T_a$ where $L_{\rm X}$ is the plateau luminosity and Ta is the rest frame plateau end time. The magnetar model predicts this correlation with b = 1 and a = log(10⁵²I₄₅−1P0, −3−2) which match observations . also showed that the magnetar model is able to produce the evolution of σ₀ that matches observations with no need of fine-tuning, unlike models in which GRBs are powered by rapidly accreting BHs [2]. The requirement of collimated relativistic flows leading to bi-polar jets in GRBs limit the central engines to only magnetars having periods of $1 \ \rm ms$ and magnetic fields of $\sim 10^{15} \ \rm G$. However, the less extreme population of magnetars can still power quite fantastic cosmic fireworks. noted that the spin down timescale given in [eqn:timescale] is shorter than the SN shock breakout time, making SNe that create magnetars brighter than usual. These subclasses of brighter Type II SNe are indeed observed (e.g. ). Recently, a number of rare superluminous SNe (SLSNe) are discovered, emitting the total radiation energy of $\sim 10^{51} \ \rm erg$ (e.g. SN2005ap, ; SN2008es, ). The radioactive decay of $\rm ^{56}Ni$ alone cannot output this amount of energy. showed that magnetars with $B \sim 10^{14} \ \rm G$ and initial period $P_i \sim 2-20 \ \rm ms$ can release rotational energy via magnetic dipole radiation on the spin-down timescale comparable to the effective diffusion time of the ejecta. They showed that these magnetars can enhance the peak luminosity to what given by L_{\rm peak} \sim E_{\rm p} t_{\rm p}/t_{\rm d}^2 \sim 5 \times 10^{43} B_{14}^{-2} ^{-1} M_{5}^{-3/2}E_{51}^{1/2} \ \rm erg s^{-1} where $$ is the electron scattering opacity, M₅ is the ejecta mass in the 5 M⊙ unit. This exceeds $10^{43} \rm \ erg \ s^{-1} $, making these events brighter than normal core-collapse Type II-P SNe. Figure 4 and 5 in show the dependence of the peak luminosity and the time to peak on B and Pi for ejecta mass of 5 and 20 M⊙ respectively. If magnetars with different magnetic field strengths and birth periods are responsible for both LGRBs and SLSNe, at least some LGRBs should be accompanied by SNe. These associations between GRBs and SNe are indeed observed, with the first pair of events being SN1998bw and GRB 980425 . A number of other events have since been observed (see e.g. for a review). showed that magnetars in the transition region in the B − Pi plane with $P_i \sim 2 \ \rm ms$ and the spin-down timescale of $\sim 10^4 \ \rm s$ can explain GRB 111209 and SN2011kl pair and the SLSN ASASSN-15lh, the most luminous SN ever observed. This association between GRBs and SNe is another evidence that these two classes of events share magnetar as their central engine. As modern transients surveys in all wavelengths continue to discover various types of peculiar events, we will have more data to constrain theoretical models, leading to a better understanding of this messy ending chapter of the stellar evolution. [1] also briefly mentioned this scenario as a possible explanation for cosmological GRBs. [2] A lot can be said about the rivalry between these two competing models for the central engine of GRBs. and references therein provide some introduction to both.