--- title: Magnetar Spin-down Model url: https://www.emergentmind.com/topics/magnetar-spin-down-model type: topic --- # Magnetar Spin-down Model Searching arXiv for recent and foundational papers on magnetar spin-down models. The magnetar spin-down model denotes the family of torque prescriptions used to describe the rotational evolution of highly magnetized neutron stars whose external activity is powered primarily by magnetic-field decay rather than by rotation. In the minimal formulation, spin evolution is set by electromagnetic dipole braking; in contemporary usage, however, the term also encompasses particle-wind braking, twisted-magnetosphere torques, gravitational-wave losses, magnetic-field decay, inclination-angle evolution, and coupling of engine power to nebulae, gamma-ray-burst afterglows, and supernova ejecta [1403.7898][2109.07484][1808.04241]. The model therefore functions both as a timing framework for isolated magnetars and as an engine model for transients powered by newly born millisecond magnetars. ## 1. Canonical torque laws and generalized braking form The baseline description of magnetar spin evolution is the magnetic-dipole law. In the vacuum approximation, the rotational-energy loss is obtained by equating the loss of rotational energy, \(E_{\rm rot}=\tfrac12 I\Omega^2\), to dipole radiation, which gives \[ \dot\Omega = -\frac{2R^6}{3Ic^3} B^2 \Omega^3 \sin^2\alpha, \qquad \dot P(t)=K\,B(t)^2 \sin^2\alpha(t)\,P(t)^{-1}, \] with \(K=8\pi^2R^6/(3Ic^3)\) [2109.07484]. In this limit the braking index is \(n=3\), corresponding to \(\dot\Omega\propto-\Omega^3\), and the evolutionary tracks on the \(\log P\)–\(\log\dot P\) plane have gradient \(2-n=-1\) [2109.07484]. A more general formulation writes the torque as \[ \dot\Omega=-K\Omega^n, \] where \(n\) is the braking index and \(K\) is a constant set by the operative torque mechanism. Integrating gives \[ \Omega(t)=\Omega_0\Bigl[1+\frac{t}{t_{\rm SD}}\Bigr]^{-1/(n-1)}, \qquad L_{\rm SD}(t)=L_0\Bigl(1+\frac{t}{t_{\rm SD}}\Bigr)^{-(n+1)/(n-1)}, \] with \(t_{\rm SD}=\Omega_0^{\,1-n}/[(n-1)K]\) or equivalently \(t_{\rm SD}=P_0/[(n-1)\dot P_0]\) [2308.12997][2008.05745]. This generalization is central because many magnetar applications do not assume \(n=3\). In the supernova context, \(n=5\) gives \(L\propto t^{-3/2}\), whereas \(n<3\) steepens the late-time decay [2308.12997]. Secular evolution is commonly modeled by allowing the dipole field to decay. A general prescription is \[ \frac{dB}{dt}=-a B^{1+\beta}, \] with solution \[ B(t)= \begin{cases} B_0\,[1+(\beta t)/\tau_B]^{-1/\beta}, & \beta\neq 0,\\[4pt] B_0\,e^{-t/\tau_B}, & \beta=0. \end{cases} \] For constant \(\alpha\) and \(\beta=0\), the period evolves as \[ P(t)=\sqrt{P_0^2+K B_0^2\sin^2\alpha_0\,\tau_B\,[1-e^{-2t/\tau_B}]}, \] and field decay causes the characteristic age \(\tau=P/(2\dot P)\) to exceed the true age once \(t\sim\tau_B\) [2109.07484]. In population synthesis, decay timescales of \(0.5\)–\(10\) kyr, with a best value of \(\sim4\) kyr, are required to reproduce the observed Galactic magnetar population while respecting the Galactic core-collapse supernova rate constraint of \(\sim20\,{\rm kyr}^{-1}\) [2109.07484]. ## 2. Wind braking and the open–closed magnetospheric partition A major extension of the dipole picture is wind braking, in which magnetic-energy release generates a quasi-steady particle outflow with total luminosity \(L_p\sim L_X\sim10^{35}\,{\rm erg\,s^{-1}}\). In this framework, the outflow splits into a “wind region” on open field lines and a “trap region” on closed field lines. Particles in the wind region escape to infinity and dominate the torque; particles in the trap region remain confined, heat the magnetosphere or surface, and radiate as X-rays [1403.7898]. The rotational-energy loss is then enhanced from the vacuum-dipole value \(\dot E_d\) to \[ \dot E_w=\dot E_d\Bigl(\frac{L_w}{\dot E_d}\Bigr)^{1/2}, \qquad N_{\rm wind}=N_{\rm dipole}\Bigl(\frac{L_w}{\dot E_d}\Bigr)^{1/2}, \] where the wind luminosity is set by the fraction of particles on open lines, \[ L_w=L_p\,\frac{\theta_{\rm open}^2}{\theta_s^2}. \] The opening angle is determined by the condition that the particle kinetic-energy density equals the magnetic-energy density, leading to \[ r_{\rm open}=7\times10^9\,b_0^{2/3}L_{p,35}^{-1/3}\Bigl(\frac{\theta_s}{0.05}\Bigr)^{2/3}\ {\rm cm}, \qquad \theta_{\rm open}\approx(r_0/r_{\rm open})^{1/2}\propto\theta_s^{-1/3} \] [1403.7898]. In this geometry, the X-ray luminosity comes from trapped particles, \[ L_X=L_p-L_w =L_p\Bigl(1-\frac{\theta_{\rm open}^2}{\theta_s^2}\Bigr). \] For roughly constant \(L_p\), a decrease in the polar-cap half-angle \(\theta_s\) increases the fraction of outflow entering the wind and decreases the trapped fraction. Quantitatively, \[ \dot E_w\propto L_w^{1/2}\propto \theta_s^{-4/3}, \qquad \dot P\propto \theta_s^{-4/3}, \] while \(L_X\) drops as \(\theta_s\) shrinks [1403.7898]. This mechanism was proposed to resolve the puzzling behavior of SGR J1745–2900, whose X-ray luminosity decreased while its spin-down rate increased. Fitting the observed factor-of-2 drop in \(L_X\) and factor-of-2.6 rise in \(\dot P\) yields \[ \theta_{s,f}\simeq0.5\,\theta_{s,i}, \] with \(\theta_{s,i}\approx3.3\times10^{-2}\) and total \(L_p=1.7\times10^{35}\,{\rm erg\,s^{-1}}\). In the extreme limit \(L_w\to L_p\), the model gives a maximum period derivative \[ \dot P_{\rm max}\simeq2.2\times10^{-11} \] for \(L_p=1.7\times10^{35}\,{\rm erg\,s^{-1}}\) and \(B_0=1.8\times10^{14}\,{\rm G}\), roughly twice the then-current \(\dot P\approx1.1\times10^{-11}\), with an estimated timescale of \(\sim240\) days to reach that maximum state [1403.7898]. ## 3. Twisted magnetospheres, outbursts, and low braking indices A second major branch of magnetar spin-down modeling attributes torque variability to a twisted external magnetosphere. In the transient magnetar XTE J1810–197, the measured spin-down displayed two regimes: during outburst decay \(\dot\nu\) varied in the range \(-(2\text{--}4.5)\times10^{-13}\,{\rm Hz\,s^{-1}}\), while during quiescence it was more stable at an average value of \(-1\times10^{-13}\,{\rm Hz\,s^{-1}}\). Over \(\sim3000\) days of quiescence, a phase-connected solution yielded \(\dot\nu=-4.9\times10^{-14}\,{\rm Hz\,s^{-1}}\) and \(\ddot\nu=1.8\times10^{-22}\,{\rm Hz\,s^{-2}}\), behavior interpreted as the decay of a strong magnetospheric twist [1602.03359]. The quantitative statement given in that analysis is not an explicit fitted torque law \(N_{\rm twist}=N_{\rm dipole}f(\Delta\phi)\), but the small-twist estimate from Beloborodov: \[ \frac{\Delta\mu}{\mu}\sim \frac{\psi^2}{4\pi}\ln\!\Bigl(\frac{u_*}{u_{\rm LC}}\Bigr), \qquad \frac{\Delta\dot\nu}{\nu}\simeq 2\,\frac{\Delta\mu}{\mu}, \] valid for \(\psi<1\) rad. Because the observed outburst-to-quiescence change in \(|\dot\nu|\) is a factor of \(\sim4\), the paper concludes that a small twist cannot account for the data and that the twist must have been large, \(\psi\gtrsim1\) [1602.03359]. An important methodological point is that the same paper explicitly does not present a direct fit of \(N_{\rm dipole}\) versus \(N_{\rm twist}(\psi)\) to the timing data. Tong and Huang framed outbursts and spin-down glitches in a broader magnetospheric toy model. They assume exponential decay of magnetic free energy, \[ \frac{dE_{\rm mf}}{dt}=-\,\frac{E_{\rm mf}}{\tau}, \qquad E_{\rm mf}(t)=E_{\rm mf,0}e^{-t/\tau}, \] identify \(|\dot E_{\rm mf}|\) with the X-ray luminosity, and relate the decay of a global twist to a shrinking hot spot and evolving torque [2005.11281]. In their phenomenological field-amplification prescription, \[ \frac{B(t)}{B_0}=1+A(1-e^{-t/\tau_1})e^{-t/\tau_2}, \] so that \(\dot P(t)\propto [B(t)]^2\). This two-timescale form is used to explain delayed torque enhancements relative to the X-ray peak in sources such as 1E 1048.1–5937 and PSR J1119–6127. For PSR J1119–6127 they adopt \(A\simeq2\), \(\tau_1\simeq50\) d, and \(\tau_2\simeq1\) yr to reproduce a torque peak occurring weeks after the X-ray maximum [2005.11281]. The low braking-index problem of Swift J1834.9–0846 has motivated a further unified formulation that combines magnetic dipole radiation, gravitational-wave emission, and wind braking: \[ N_{\rm tot}=N_{\rm dip}+N_{\rm gw}+N_{\rm wind}, \qquad N_{\rm wind}=-\kappa I\Omega. \] Within this framework the wind braking parameter is constrained to \(\kappa\in[13,37]\), and wind braking is found to contribute substantially, \(17\%\text{--}51\%\), to the current spin-down torque. The observed braking index \(n=1.08\pm0.04\) is then interpreted as the result of the combined effect of wind torque, slow field decay, and magnetic-inclination evolution, with the analysis favoring a toroidally-dominated internal magnetic field and constraining the number of precession cycles to \(\xi\sim10^4\text{--}10^5\) [2602.06615]. ## 4. Newly born magnetars as central engines For newborn millisecond magnetars, spin-down models are used as central-engine prescriptions for GRB afterglows, extended emission, and supernova light curves. In the standard magnetic-dipole case, \[ L_{\rm sd}(t)=L_0\,(1+t/\tau_{\rm sd})^{-2}, \] with \[ \tau_{\rm sd}=\frac{6Ic^3}{B^2R_{\rm ns}^6\Omega_0^2} \approx1.3\times10^5\,B_{14}^{-2}P_{0,{\rm ms}}^2\ {\rm s}, \qquad E_{\rm rot}=\tfrac12 I\Omega_0^2 \approx2\times10^{52}I_{45}P_{0,{\rm ms}}^{-2}\ {\rm erg} \] [2606.18842]. This same \(L_{\rm sd}\propto(1+t/\tau_{\rm sd})^{-2}\) form underlies short-GRB plateau modeling and broadband afterglow energy-injection calculations [1302.3643][1411.5477]. If gravitational-wave losses are important, the spin-down channels become \[ L_{\rm md}=\frac{B_p^2R^6\Omega^4}{6c^3}, \qquad L_{\rm gw}=\frac{32GI^2\epsilon^2\Omega^6}{5c^5}. \] The corresponding characteristic timescales are \[ \tau_{\rm md}\simeq2.0\times10^3\,{\rm s}\, \bigl(I_{45}B_{p,15}^{-2}P_{0,-3}^2R_6^{-6}\bigr), \] \[ \tau_{\rm gw}\simeq9.1\times10^3\,{\rm s}\, \bigl(I_{45}^{-1}\epsilon_{-3}^{-2}P_{0,-3}^4\bigr), \] and the electromagnetic luminosity can be written \[ L_{\rm em}(t)=L_0\Bigl(1+\frac{t}{\tau_c}\Bigr)^{\alpha}, \] with \(\alpha=-1\) in the GW-dominated regime and \(\alpha=-2\) in the MD-dominated regime. A smooth break from \(-1\) to \(-2\) is therefore a predicted signature of a transition from GW-dominated to MD-dominated spin-down, while collapse to a black hole before \(\tau_c\) yields a plateau followed by a sharp drop [1808.04241]. Sarin et al. generalized this picture further by allowing an arbitrary constant braking index \(n\) and radiative losses in the external shock: \[ \frac{d\Omega}{dt}=-K\Omega^n, \qquad L_{\rm sd}(t)=L_0\,[1+t/\tau_0]^{(1+n)/(1-n)}, \] \[ \frac{dE}{dt}=L_{\rm sd}(t)-\kappa\,\frac{E}{t}, \qquad L(t)=A t^\Gamma + H(t-t_0)\,\kappa\,\frac{E(t)}{t}. \] Applied to Swift X-ray afterglows, this framework yielded \(n=4.85^{+0.11}_{-0.15}\) for GRB 061121, suggesting that the millisecond magnetar in that burst spins down predominantly through gravitational-wave emission [2008.05745]. The r-mode variant of GW-dominated spin-down has been used for GRB 130831A. In that model the X-ray plateau luminosity is written as a fraction of \(L_{\rm md}\), with a piecewise temporal dependence \(F(t)\): \(F(t)\propto t^0\) for \(tT_c\). Matching \((1+z)T_g\simeq269\) s and \((1+z)T_{\rm col}\simeq10^5\) s gives \(P_i\approx0.8\) ms and \(B_p\approx10^{14}\) G [1604.04817]. In short GRBs with extended emission, the same dipole formalism has been used to infer post-EE and birth spin periods. The relations employed are \[ L(t)=\frac{B^2R^6\Omega^4}{6c^3}, \qquad \tau=\frac{3c^3I}{B^2R^6\Omega^2}, \qquad \Delta E=2\pi^2 I\,(P_i^{-2}-P_0^{-2}). \] For a sample of nine extended-emission bursts, the mean values were roughly \(\langle B\rangle\sim1.4\times10^{16}\) G, \(\langle P_0\rangle\sim20\) ms, and \(\langle P_i\rangle\sim3\) ms [1302.3643]. In superluminous-supernova and luminous-transient modeling, the same engine law is coupled to ejecta dynamics and diffusion. Metzger et al. identified a transition region in the \(B_d\)–\(P_0\) plane in which a magnetar can power both a long GRB and a luminous supernova, with a \(2\) ms magnetar and \(t_{\rm sd}\sim10^4\) s proposed for GRB 111209A/SN 2011kl [1508.02712]. Omand and Sarin generalized supernova fitting to non-dipole spin-down by retaining the \(n\)-dependent luminosity law and coupling it to one-zone ejecta dynamics [2308.12997]. More recently, a time-dependent radiative-diffusion treatment coupled magnetar injection to a pulsar-wind nebula and forward shock, showing that the model naturally produces well-separated double peaks, partially merged peaks, or single broad peaks; in an illustrative fit to LSQ14bdq, the adopted magnetar parameters were \(P_0\approx1.2\) ms, \(B\approx7\times10^{13}\) G, and \(\tau_{\rm sd}\approx30\) d [2606.18842]. ## 5. Evolution from proto-magnetars to old magnetars The earliest spin evolution of a magnetar need not be dipole dominated. Two-dimensional axisymmetric MHD simulations of proto-magnetars with initial spin periods \(P_{\star0}=50\text{--}500\) ms show that neutrino-heated magneto-centrifugal winds can remove angular momentum far more efficiently than the canonical dipole formula during the Kelvin–Helmholtz cooling epoch. The torque is written \[ N(t)=\dot J(t)=\dot M(t)\,\Omega(t)\,\langle R_A(t)^2\rangle, \] so that \[ I\frac{d\Omega}{dt}=-\dot M\,\Omega\,R_A^2, \qquad \tau_{\rm sd}=\frac{I}{\dot M R_A^2}. \] For \(B_0\gtrsim10^{15}\) G, the early spin-down timescale can be of order seconds, and a fit for spinning down to \(P>1\) s within \(\sim5\) s is \[ B_0\gtrsim1.3\times10^{15}\,{\rm G} \Bigl(\frac{P_{\star0}}{400\,{\rm ms}}\Bigr)^{-1.4} \Bigl(\frac{M}{1.4\,M_\odot}\Bigr)^{2.2} \] [2208.09042]. On kiloyear timescales, field decay reshapes the observed \(P\)–\(\dot P\) distribution. Monte Carlo synthesis with exponentially decaying or super-exponentially decaying fields indicates that the dipole field must decay on a characteristic timescale of \(0.5\text{--}10\) kyr, with a best value of \(\sim4\) kyr, and that the initial spin period must be less than \(2\) sec. Under these constraints there are multiple choices of input physics that can reproduce the observed magnetar population reasonably well, and the faded synthetic populations overlap substantially with X-ray dim isolated neutron stars while generally not favoring an evolutionary link to RRATs [2109.07484]. At still later stages, additional torques become relevant. For wind-fed accreting magnetars, the total torque is written \[ I\frac{d\Omega}{dt}=N_{\rm tot} =\zeta\,\dot M\,R_m^2\,\Omega_K(R_m)\,[1-\Omega/\Omega_K(R_m)], \] with \[ R_m=\Bigl(\frac{\mu^4}{2GM\dot M^2}\Bigr)^{1/7}, \qquad R_{\rm co}=(GM/\Omega^2)^{1/3}, \qquad P_{\rm eq}=6.7\times10^3\,\mu_{33}^{6/7}\dot M_{15}^{-3/7}\ {\rm s}. \] In that picture, 4U 2206+54 is placed in early propeller spin-down, AX J1910.7+0917 at the transition point, and 2S 0114+65 in late spin-up associated with magnetic-field decay and transient-disk formation [1912.03839]. A related isolated-neutron-star extension invokes a pulsar-to-propeller transition to explain long-period radio transients. In this framework the early stage is pure dipole braking, but once \(R_{\rm LC}\gtrsim R_m\) the star enters the propeller phase. Population synthesis shows that two propeller prescriptions, Models E and F, can account for most of the observed LPT periods, \(P\sim1\text{--}400\) min, and \(\dot P<10^{-9}\,{\rm s\,s^{-1}}\). The same study finds that a transition from the pulsar to the propeller phase is required to reach the observed LPT period range \(P>10^3\) s [2602.15024]. ## 6. Alternative braking channels, diagnostics, and open issues Not all proposed magnetar spin-down models are wind- or twist-based. One alternative is the joint “magnetic-dipole + quantum vacuum friction” scenario, in which an additional dissipative torque arises from interaction between the strong surface field and the polarized quantum vacuum. In that model, \[ \dot P =\frac{2\pi^2R^6\sin^2\theta}{3c^3I}\frac{B_{\rm dip}^2}{P} +\frac{3\alpha R^4\sin^2\theta}{64\pi B_c^2 c I}\,B_{\rm surf}^4\,P, \] with \(\xi\equiv B_{\rm surf}/B_{\rm dip}\). The ratio of QVF-to-dipole losses is then \[ \frac{\dot E_{\rm QVF}}{\dot E_{\rm dip}} \simeq 7.69\times10^{-24}\,B_{\rm dip}^2 P^2 \xi^4, \] and QVF dominates when \[ P^3\dot P > 0.63\times10^{-16}\,\xi^{-4}\ {\rm s}^2. \] In the pure-dipole limit the braking index tends to \(3\); in the pure-QVF limit it tends to \(1\) [1504.02204]. This provides an explicit alternative interpretation for low braking indices and low inferred dipole fields. A separate observational diagnostic is the magnetar wind nebula. In one broadband nebular model, the stellar torque is parameterized phenomenologically as \[ \dot\Omega=-k\Omega^n, \qquad L_{\rm spin}(t)=L_0\Bigl(1+\frac{t}{t_0}\Bigr)^{-(n+1)/(n-1)}, \] and the spin-down power is mapped into a broken-power-law pair-injection spectrum [1606.01391]. However, the same study emphasizes that early-injected particles suffer severe adiabatic and synchrotron losses, so the present-day nebular spectrum does not allow one to test whether 1E 1547.0–5408 had a millisecond birth period. This suggests that nebular calorimetry alone is not a universal proxy for the birth spin of a magnetar [1606.01391]. Several papers also identify the limits of their own torque models. The XTE J1810–197 timing study explicitly states that it does not work out or fit a full analytic \(N_{\rm dipole}\) versus \(N_{\rm twist}\) derivation [1602.03359]. The outburst model of Tong and Huang is explicitly a toy model and does not derive, from first principles, how a twist injected on closed lines feeds open-field-line currents [2005.11281]. In the transient-light-curve literature, semi-analytic supernova and afterglow models likewise rely on one-zone diffusion, gray opacities, constant braking indices, or phenomenological leakage prescriptions [2308.12997][2606.18842][2008.05745]. A recurrent misconception is that magnetar timing can be reduced to a single \(n=3\) vacuum-dipole law. The literature instead supports a hierarchy of regimes: neutrino-driven wind torques immediately after birth, dipole and gravitational-wave losses in millisecond engines, wind braking and twist-mediated torques in active Galactic magnetars, field-decay-driven secular evolution on kiloyear timescales, and propeller or accretion torques in older systems [2208.09042][1808.04241][1403.7898][2109.07484][2602.15024]. The principal unresolved issue is therefore not whether a magnetar spins down, but which torque channel dominates in a given source, epoch, and observational band. Source: https://www.emergentmind.com/topics/magnetar-spin-down-model