-
Magnetic helicity estimations in models and observations of the solar magnetic field. Part IV: Application to solar observations
Authors:
J. K. Thalmann,
M. K. Georgoulis,
Y. Liu,
E. Pariat,
G. Valori,
S. Anfinogentov,
F. Chen,
Y. Guo,
K. Moraitis,
S. Yang,
A. Mastrano
Abstract:
In this ISSI-supported series of studies on magnetic helicity in the Sun, we systematically implement different magnetic helicity calculation methods on high-quality solar magnetogram observations. We apply finite-volume, discrete flux tube (in particular, connectivity-based) and flux-integration methods to data from Hinode's Solar Optical Telescope. The target is NOAA active region 10930 during a…
▽ More
In this ISSI-supported series of studies on magnetic helicity in the Sun, we systematically implement different magnetic helicity calculation methods on high-quality solar magnetogram observations. We apply finite-volume, discrete flux tube (in particular, connectivity-based) and flux-integration methods to data from Hinode's Solar Optical Telescope. The target is NOAA active region 10930 during a ~1.5 day interval in December 2006 that included a major eruptive flare (SOL2006-12-13T02:14X3.4). Finite-volume and connectivity-based methods yield instantaneous budgets of the coronal magnetic helicity, while the flux-integration methods allow an estimate of the accumulated helicity injected through the photosphere. The objectives of our work are twofold: A cross-validation of methods, as well as an interpretation of the complex events leading to the eruption. To the first objective, we find (i) strong agreement among the finite-volume methods, (ii) a moderate agreement between the connectivity-based and finite-volume methods, (iii) an excellent agreement between the flux-integration methods, and (iv) an overall agreement between finite-volume and flux-integration based estimates regarding the predominant sign and magnitude of the helicity. To the second objective, we are confident that the photospheric helicity flux significantly contributed to the coronal helicity budget, and that a right-handed structure erupted from a predominantly left-handed corona during the X-class flare. Overall, we find that the use of different methods to estimate the (accumulated) coronal helicity may be necessary in order to draw a complete picture of an active-region corona, given the careful handling of identified data (preparation) issues, which otherwise would mislead the event analysis and interpretation.
△ Less
Submitted 19 August, 2021;
originally announced August 2021.
-
Reconstructing Highly-twisted Magnetic Fields
Authors:
Victor M. Demcsak,
Michael S. Wheatland,
Alpha Mastrano,
Kai E. Yang
Abstract:
We investigate the ability of a nonlinear force-free code to calculate highly-twisted magnetic field configurations using the Titov and Démoulin (1999) equilibrium field as a test case. The code calculates a force-free field using boundary conditions on the normal component of the field in the lower boundary, and the normal component of the current density over one polarity of the field in the low…
▽ More
We investigate the ability of a nonlinear force-free code to calculate highly-twisted magnetic field configurations using the Titov and Démoulin (1999) equilibrium field as a test case. The code calculates a force-free field using boundary conditions on the normal component of the field in the lower boundary, and the normal component of the current density over one polarity of the field in the lower boundary. The code can also use the current density over both polarities of the field in the lower boundary as a boundary condition. We investigate the accuracy of the reconstructions with increasing flux-rope surface twist number $N_{\textrm{t}}$, achieved by decreasing the sub-surface line current in the model. We find that the code can approximately reconstruct the Titov-Démoulin field for surface twist numbers up to $N_{\textrm{t}} \approx 8.8$. This includes configurations with bald patches. We investigate the ability to recover bald patches, and more generally identify the limitations of our method for highly-twisted fields. The results have implications for our ability to reconstruct coronal magnetic fields from observational data.
△ Less
Submitted 7 August, 2020;
originally announced August 2020.
-
A rotation-driven pulsar radio emission mechanism
Authors:
D. B. Melrose,
M. Z. Rafat,
A. Mastrano
Abstract:
We propose and discuss an alternative pulsar radio emission mechanism that relies on rotation-driven plasma oscillations, rather than on a beam-driven instability, and suggest that it may be the generic radio emission mechanism for pulsars. We identify these oscillations as superluminal longitudinal waves in the pulsar plasma, and point out that these waves can escape directly in the O~mode. We ar…
▽ More
We propose and discuss an alternative pulsar radio emission mechanism that relies on rotation-driven plasma oscillations, rather than on a beam-driven instability, and suggest that it may be the generic radio emission mechanism for pulsars. We identify these oscillations as superluminal longitudinal waves in the pulsar plasma, and point out that these waves can escape directly in the O~mode. We argue that the frequency of the oscillations is $ω_0\approxω_{\rm p}(2\langleγ\rangle)^{1/2}/γ_{\rm s}$, where $γ_{\rm s}$ is the Lorentz factor of bulk streaming motion and $\langleγ\rangle$ is the mean Lorentz factor in the rest frame of the plasma. The dependence of the plasma frequency $ω_{\rm p}$ on radial distance implies a specific frequency-to-radius mapping, $ω_0\propto r^{-3/2}$. Escape of the energy in these oscillations is possible if they are generated in overdense, field-aligned regions that we call fibers; the wave energy is initially refracted into underdense regions between the fibers, which act as ducts. Some implications of the model for the interpretation of pulsar radio emission are discussed.
△ Less
Submitted 11 November, 2020; v1 submitted 27 June, 2020;
originally announced June 2020.
-
Pulsar radio emission mechanisms: a critique
Authors:
D. B. Melrose,
M. Z. Rafat,
A. Mastrano
Abstract:
We consider critically the three most widely favored pulsar radio emission mechanisms: coherent curvature emission (CCE), beam-driven relativistic plasma emission (RPE) and anomalous Doppler emission (ADE). We assume that the pulsar plasma is one dimensional (1D), streaming outward with a bulk Lorentz factor $γ_{\rm s} \gg \langle γ\rangle -1 \gtrsim 1$, where $\langle γ\rangle $ is the intrinsic…
▽ More
We consider critically the three most widely favored pulsar radio emission mechanisms: coherent curvature emission (CCE), beam-driven relativistic plasma emission (RPE) and anomalous Doppler emission (ADE). We assume that the pulsar plasma is one dimensional (1D), streaming outward with a bulk Lorentz factor $γ_{\rm s} \gg \langle γ\rangle -1 \gtrsim 1$, where $\langle γ\rangle $ is the intrinsic spread in the rest frame of the plasma. We argue that the formation of beams in a multi-cloud model is ineffective in the intrinsically relativistic case for plausible parameters, because the overtaking takes too long. We argue that the default choice for the particle distribution in the rest frame is a J{ü}ttner distribution and that relativistic streaming should be included by applying a Lorentz transformation to the rest-frame distribution, rather than the widely assumed relativistically streaming Gaussian distribution. We find that beam-driven wave growth is severely restricted by (a) the wave properties in pulsar plasma, (b) a separation condition between beam and background, and (c) the inhomogeneity of the plasma in the pulsar frame. The growth rate for the kinetic instability is much smaller and the bandwidth of the growing waves is much larger for a J{ü}ttner distribution than for a relativistically streaming Gaussian distribution. No reactive instability occurs at all for a J{ü}ttner distribution. We conclude that none of CCE, RPE and ADE in tenable as the generic pulsar radio emission mechanism for ``plausible'' assumptions about the pulsar plasma.
△ Less
Submitted 1 November, 2020; v1 submitted 26 June, 2020;
originally announced June 2020.
-
Self-consistent Nonlinear Force-free Field Reconstruction from Weighted Boundary Conditions
Authors:
Alpha Mastrano,
Kai E. Yang,
Michael S. Wheatland
Abstract:
Vector magnetogram data are often used as photospheric boundary conditions for force-free coronal magnetic field extrapolations. In general, however, vector magnetogram data are not consistent with the force-free assumption. In this article, we demonstrate a way to deal with inconsistent boundary data, by generalizing the "self-consistency procedure" of Wheatland & Regnier (2009). In that procedur…
▽ More
Vector magnetogram data are often used as photospheric boundary conditions for force-free coronal magnetic field extrapolations. In general, however, vector magnetogram data are not consistent with the force-free assumption. In this article, we demonstrate a way to deal with inconsistent boundary data, by generalizing the "self-consistency procedure" of Wheatland & Regnier (2009). In that procedure, the inconsistency is resolved by an iterative process of constructing two solutions based on the values of the force-free parameter alpha on the two polarities of the field in the boundary (the P and N polarities), and taking uncertainty-weighted averages of the boundary alpha values in the P and N solutions. When the alpha values in the P and N regions are very different, the self-consistent solution may lose high alpha values from the boundary conditions. We show how, by altering the weighting of the uncertainties in the P or N boundary conditions, we can preserve high alpha values in the self-consistent solution. The weighted self-consistent extrapolation method is demonstrated on an analytic bipole field and applied to vector magnetogram data taken by the Helioseismic and Magnetic Imager (HMI) instrument for NOAA active region AR 12017 on 2014 March 29.
△ Less
Submitted 26 April, 2020;
originally announced April 2020.
-
Wave dispersion in pulsar plasma: 3. Beam-driven instabilities
Authors:
M. Z. Rafat,
D. B. Melrose,
A. Mastrano
Abstract:
Beam-driven instabilities are considered in a pulsar plasma assuming that both the background plasma and the beam are relativistic Jüttner distributions. In the rest frame of the background, the only waves that can satisfy the resonance condition are in a tiny range of slightly subluminal phase speeds. The growth rate for the kinetic (or maser) version of the weak-beam instability is much smaller…
▽ More
Beam-driven instabilities are considered in a pulsar plasma assuming that both the background plasma and the beam are relativistic Jüttner distributions. In the rest frame of the background, the only waves that can satisfy the resonance condition are in a tiny range of slightly subluminal phase speeds. The growth rate for the kinetic (or maser) version of the weak-beam instability is much smaller than has been estimated for a relativistically streaming Gaussian distribution, and the reasons for this are discussed. The growth rate for the reactive version of the weak-beam instability is treated in a conventional way. We compare the results with exact calculations, and find that the approximate solutions are not consistent with the exact results. We conclude that, for plausible parameters, there is no reactive version of the instability. The growth rate in the pulsar frame is smaller than that in the rest frame of the background plasma by a factor $2γ_{\rm s}$, where $γ_{\rm s} = 10^2 {\rm -} 10^3$ is the Lorentz factor of the bulk motion of the background plasma, placing a further constraint on effective wave growth. Based on these results, we argue that beam-driven wave growth probably plays no role in pulsar radio emission.
△ Less
Submitted 5 June, 2019;
originally announced June 2019.
-
A weak coronal heating event associated with periodic particle acceleration episodes
Authors:
Atul Mohan,
Patrick I. McCauley,
Alpha Mastrano,
Divya Oberoi
Abstract:
Weak heating events are frequent and ubiquitous in solar corona. They derive their energy from the local magnetic field and form a major source of local heating, signatures of which are seen in EUV and X-ray bands. Their radio emission arise from various plasma instabilities that lead to coherent radiation, making even a weak flare appear very bright. The radio observations hence probe non-equilib…
▽ More
Weak heating events are frequent and ubiquitous in solar corona. They derive their energy from the local magnetic field and form a major source of local heating, signatures of which are seen in EUV and X-ray bands. Their radio emission arise from various plasma instabilities that lead to coherent radiation, making even a weak flare appear very bright. The radio observations hence probe non-equilibrium dynamics providing complementary information on plasma evolution. However, a robust study of radio emission from a weak event among many simultaneous events, requires high dynamic range imaging at sub-second andsub-MHz resolutions owing to their high spectro-temporal variability. Such observations were not possible until recently.This is among the first spatially resolved studies of an active region loop hosting a transient brightening (ARTB) and dynamically linked to a metrewave type-I noise storm. It uses imaging observations at metrewave, EUV and X-ray bands, along with magnetogram data. We believe this is the first spectroscopic radio imaging study of a type-I source, the data for which was obtained using the Murchison Widefield Array. We report the discovery of 30 s quasi-periodic oscillations (QPOs) in the radio light curve, riding on a coherent baseline flux. The strength of the QPOs and the baseline flux enhanced during a mircoflare associated with the ARTB. Our observations suggest a scenario where magnetic stress builds up over an Alfvén timescale (30s) across the typical magnetic field braiding scale and then dissipates via a cascade of weak reconnection events.
△ Less
Submitted 24 July, 2019; v1 submitted 27 February, 2019;
originally announced February 2019.
-
Wave dispersion in pulsar plasma: 1. Plasma rest frame
Authors:
M. Z. Rafat,
D. B. Melrose,
A. Mastrano
Abstract:
Wave dispersion in a pulsar plasma (a 1D, strongly magnetized, pair plasma streaming highly relativistically with a large spread in Lorentz factors in its rest frame) is discussed, motivated by interest in beam-driven wave turbulence and the pulsar radio emission mechanism. In the rest frame of the pulsar plasma there are three wave modes in the low-frequency, non-gyrotropic approximation. For par…
▽ More
Wave dispersion in a pulsar plasma (a 1D, strongly magnetized, pair plasma streaming highly relativistically with a large spread in Lorentz factors in its rest frame) is discussed, motivated by interest in beam-driven wave turbulence and the pulsar radio emission mechanism. In the rest frame of the pulsar plasma there are three wave modes in the low-frequency, non-gyrotropic approximation. For parallel propagation these are referred to as the X, A and L modes, with the X and A modes having dispersion relation $z=z_A\approx1-1/2β_A^2$, where $z=ω/k_\parallel c$ is the phase speed and $β_Ac$ is the Alfven speed. The L mode dispersion relation is determined by a relativistic plasma dispersion function, $z^2W(z)$, which is negative for $ z < z_0 $ and has a sharp maximum at $z=z_m$, with $1-z_m<1-z_0\ll1$. We give numerical estimates for the maximum of $z^2W(z)$ and for $z_m$ and $z_0$ for a 1D Juttner distribution. The L and A modes reconnect, for $z_A>z_0$, to form the O and Alfven modes for oblique propagation ($θ\neq0$). For $z_A<z_0$ the Alfven and O~mode curves reconnect forming a new mode that exists only for $\tan^2θ>z_0^2-z_A^2$.
The L mode is the nearest counterpart to Langmuir waves in a nonrelativistic plasma, but we argue that there are no `Langmuir-like' waves in pulsar plasma, identifying three features of the L~mode (dispersion relation, ratio of electric to total energy and group speed) that are not Langmuir-like. A beam-driven instability requires a beam speed equal to the phase speed of the wave. This resonance condition can be satisfied for the O mode, but only for an implausibly energetic beam and only for a tiny range of angles for the O~mode around $θ\approx0$. The resonance is also possible for the Alfven mode but only near a turnover frequency that has no counterpart for Alfven waves in a nonrelativistic plasma.
△ Less
Submitted 12 May, 2019; v1 submitted 17 December, 2018;
originally announced December 2018.
-
Wave dispersion in pulsar plasma: 2. Pulsar frame
Authors:
M. Z. Rafat,
D. B. Melrose,
A. Mastrano
Abstract:
Wave dispersion in a pulsar plasma is discussed emphasizing the relevance of different inertial frames, notably the plasma rest frame ${\cal K}$ and the pulsar frame ${\cal K}'$ in which the plasma is streaming with speed $β_{\rm s}$. The effect of a Lorentz transformation on both subluminal, $|z|<1$, and superluminal, $|z|>1$, waves is discussed. It is argued that the preferred choice for a relat…
▽ More
Wave dispersion in a pulsar plasma is discussed emphasizing the relevance of different inertial frames, notably the plasma rest frame ${\cal K}$ and the pulsar frame ${\cal K}'$ in which the plasma is streaming with speed $β_{\rm s}$. The effect of a Lorentz transformation on both subluminal, $|z|<1$, and superluminal, $|z|>1$, waves is discussed. It is argued that the preferred choice for a relativistically streaming distribution should be a Lorentz-transformed Jüttner distribution; such a distribution is compared with other choices including a relativistically streaming Gaussian distribution. A Lorentz transformation of the dielectric tensor is written down, and used to derive an explicit relation between the relativistic plasma dispersion functions in ${\cal K}$ and ${\cal K}'$. It is shown that the dispersion equation can be written in an invariant form, implying a one-to-one correspondence between wave modes in any two inertial frames. Although there are only three modes in the plasma rest frame, it is possible for backward-propagating or negative-frequency solutions in ${\cal K}$ to transform into additional forward-propagating, positive-frequency solutions in ${\cal K}'$ that may be regarded as additional modes.
△ Less
Submitted 3 June, 2019; v1 submitted 17 December, 2018;
originally announced December 2018.
-
Nonlinear force-free modeling of flare-related magnetic field changes at the photosphere and chromosphere
Authors:
Lucia Kleint,
Michael S. Wheatland,
Alpha Mastrano,
Patrick I. McCauley
Abstract:
Rapid and stepwise changes of the magnetic field are often observed during flares but cannot be explained by models yet. Using a 45 min sequence of SDO/HMI 135 s fast-cadence vector magnetograms of the X1 flare on 2014-03-29 we construct, at each timestep, nonlinear force-free models for the coronal magnetic field. Observed flare-related changes in the line-of-sight magnetic field $B_{\rm LOS}$ at…
▽ More
Rapid and stepwise changes of the magnetic field are often observed during flares but cannot be explained by models yet. Using a 45 min sequence of SDO/HMI 135 s fast-cadence vector magnetograms of the X1 flare on 2014-03-29 we construct, at each timestep, nonlinear force-free models for the coronal magnetic field. Observed flare-related changes in the line-of-sight magnetic field $B_{\rm LOS}$ at the photosphere and chromosphere are compared with changes in the magnetic fields in the models. We find a moderate agreement at the photospheric layer (the basis for the models), but no agreement at chromospheric layers. The observed changes at the photosphere and chromosphere are surprisingly different, and are unlikely to be reproduced by a force-free model. The observed changes are likely to require a change in the magnitude of the field, not just in its direction.
△ Less
Submitted 21 August, 2018;
originally announced August 2018.
-
Photospheric response to a flare
Authors:
Michael S. Wheatland,
Donald B. Melrose,
Alpha Mastrano
Abstract:
Flares produce sudden and permanent changes in the horizontal photospheric magnetic field. In particular flares generally produce increased magnetic shear in the photospheric field along the neutral line. Recent observations show also that flares can produce sudden photospheric motion. We present a model for the observed changes as the response of the photosphere to a large-amplitude shear Alfvén…
▽ More
Flares produce sudden and permanent changes in the horizontal photospheric magnetic field. In particular flares generally produce increased magnetic shear in the photospheric field along the neutral line. Recent observations show also that flares can produce sudden photospheric motion. We present a model for the observed changes as the response of the photosphere to a large-amplitude shear Alfvén wave propagating down from the corona on either side of the neutral line. The Alfvénic front is assumed to impact the photosphere close to the neutral line first, and then successively further away with time, such that the line of impact coincides with the flare ribbon. The wave introduces magnetic shear and velocity shear. The magnetic shear introduced at the photosphere has the same sign on either side of the neutral line, while the velocity shear has the opposite sign. We discuss the possibility that this process is responsible for particle acceleration in flares.
△ Less
Submitted 9 August, 2018;
originally announced August 2018.
-
Gravitational radiation from neutron stars deformed by crustal Hall drift
Authors:
Arthur George Suvorov,
Alpha Mastrano,
Ulrich Geppert
Abstract:
A precondition for the radio emission of pulsars is the existence of strong, small-scale magnetic field structures (`magnetic spots') in the polar cap region. Their creation can proceed via crustal Hall drift out of two qualitatively and quantitatively different initial magnetic field configurations: a field confined completely to the crust and another which penetrates the whole star. The aim of t…
▽ More
A precondition for the radio emission of pulsars is the existence of strong, small-scale magnetic field structures (`magnetic spots') in the polar cap region. Their creation can proceed via crustal Hall drift out of two qualitatively and quantitatively different initial magnetic field configurations: a field confined completely to the crust and another which penetrates the whole star. The aim of this study is to explore whether these magnetic structures in the crust can deform the star sufficiently to make it an observable source of gravitational waves. We model the evolution of these field configurations, which can develop, within $\sim 10^4$ -- $10^5$ yr, magnetic spots with local surface field strengths $\sim 10^{14}$ G maintained over $\gtrsim 10^6$ yr. Deformations caused by the magnetic forces are calculated. We show that, under favourable initial conditions, a star undergoing crustal Hall drift can have ellipticity $ε\sim 10^{-6}$, even with sub-magnetar polar field strengths, after $\sim 10^5$ yr. A pulsar rotating at $\sim 10^2$ Hz with such $ε$ is a promising gravitational-wave source candidate. Since such large deformations can be caused only by a particular magnetic field configuration that penetrates the whole star and whose maximum magnetic energy is concentrated in the outer core region, gravitational wave emission observed from radio pulsars can thus inform us about the internal field structures of young neutron stars.
△ Less
Submitted 14 April, 2016;
originally announced April 2016.
-
Electromagnetic spin down of a transient accreting millisecond pulsar during quiescence
Authors:
Andrew Melatos,
Alpha Mastrano
Abstract:
The measured spin-down rates in quiescence of the transient accreting millisecond pulsars IGR J00291+5934, XTE J1751-305, SAX J1808.4-3658, and Swift J1756.9-2508 have been used to estimate the magnetic moments of these objects assuming standard magnetic dipole braking. It is shown that this approach leads to an overestimate, if the amount of residual accretion is enough to distort the magnetosphe…
▽ More
The measured spin-down rates in quiescence of the transient accreting millisecond pulsars IGR J00291+5934, XTE J1751-305, SAX J1808.4-3658, and Swift J1756.9-2508 have been used to estimate the magnetic moments of these objects assuming standard magnetic dipole braking. It is shown that this approach leads to an overestimate, if the amount of residual accretion is enough to distort the magnetosphere away from a force-free configuration, through magnetospheric mass loading or crushing, so that the lever arm of the braking torque migrates inside the light cylinder. We derive an alternative spin-down formula and calculate the residual accretion rates where the formula is applicable. As a demonstration, we apply the alternative spin-down formula to produce updated magnetic moment estimates for the four objects above. We note that, based on current uncertain observations of quiescent accretion rates, magnetospheric mass loading and crushing are neither firmly indicated nor ruled out in these four objects. Because quiescent accretion rates are not measured directly (only upper limits are placed), it is impossible to be confident without more data whether the thresholds for magnetospheric mass loading or crushing are reached or not.
△ Less
Submitted 27 December, 2015;
originally announced December 2015.
-
Spin paramagnetic deformation of a neutron star
Authors:
A. G. Suvorov,
A. Mastrano,
A. Melatos
Abstract:
Quantum mechanical corrections to the hydromagnetic force balance equation, derived from the microscopic Schrödinger-Pauli theory of quantum plasmas, modify the equilibrium structure and hence the mass quadrupole moment of a neutron star. It is shown here that the dominant effect --- spin paramagnetism --- is most significant in a magnetar, where one typically has…
▽ More
Quantum mechanical corrections to the hydromagnetic force balance equation, derived from the microscopic Schrödinger-Pauli theory of quantum plasmas, modify the equilibrium structure and hence the mass quadrupole moment of a neutron star. It is shown here that the dominant effect --- spin paramagnetism --- is most significant in a magnetar, where one typically has $μ_{B}|\boldsymbol{B}|\gtrsim k_B T_e$, where $μ_{B}$ is the Bohr magneton, $\boldsymbol{B}$ is the magnetic field, and $T_e$ is the electron temperature. The spin paramagnetic deformation of a nonbarotropic magnetar with a linked poloidal-toroidal magnetic field is calculated to be up to ${\sim 10}$ times greater than the deformation caused solely by the Lorentz force. It depends on the degree of Pauli blocking by conduction electrons and the propensity to form magnetic domains, processes which are incompletely modelled at magnetar field strengths. The star becomes more oblate, as the toroidal field component strengthens. The result implies that existing classical predictions underestimate the maximum strength of the gravitational wave signal from rapidly spinning magnetars at birth. Turning the argument around, future gravitational-wave upper limits of increasing sensitivity will place ever-stricter constraints on the physics of Pauli blocking and magnetic domain formation under magnetar conditions.
△ Less
Submitted 19 November, 2015;
originally announced November 2015.
-
Interpreting the AXP 1E 2259+586 antiglitch as a change in internal magnetization
Authors:
Alpha Mastrano,
Arthur George Suvorov,
Andrew Melatos
Abstract:
The sudden spin-down event ('anti-glitch') observed in AXP 1E 2259+586 on 2012 April 21 was arguably caused by a decay of its internal toroidal magnetic field component, which turns a stable prolate configuration into an unstable one. We refine previous models of this process by modelling the star's magnetic field self-consistently as a 'twisted torus' configuration in non-barotropic equilibrium (…
▽ More
The sudden spin-down event ('anti-glitch') observed in AXP 1E 2259+586 on 2012 April 21 was arguably caused by a decay of its internal toroidal magnetic field component, which turns a stable prolate configuration into an unstable one. We refine previous models of this process by modelling the star's magnetic field self-consistently as a 'twisted torus' configuration in non-barotropic equilibrium (which allows us to explore a greater range of equilibrium configurations). It is shown that, if the star's magnetic field is purely dipolar, the change in the toroidal field strength required to produce an anti-glitch of the observed size can be ~ 10 times larger than previously calculated. If the star has a quadrupolar magnetic field component, then an anti-glitch of similar magnitude can be produced via a decay of the quadrupole component, in addition to a decay of the toroidal component. We show that, if the quadrupole component decays, the minimum initial toroidal field strength and the change in toroidal field strength needed to produce the observed anti-glitch are lower than in the pure dipole twisted torus. In addition, we predict the maximum anti-glitch sizes, assuming that they are caused by a change in ellipticity, in four glitching magnetars and discuss the implications for energetics of accompanying X-ray bursts.
△ Less
Submitted 18 August, 2015; v1 submitted 27 July, 2015;
originally announced July 2015.
-
Neutron star deformation due to poloidal-toroidal magnetic fields of arbitrary multipole order: a new analytic approach
Authors:
Alpha Mastrano,
Arthur George Suvorov,
Andrew Melatos
Abstract:
A recipe is presented to construct an analytic, self-consistent model of a non-barotropic neutron star with a poloidal-toroidal field of arbitrary multipole order, whose toroidal component is confined in a torus around the neutral curve inside the star, as in numerical simulations of twisted tori. The recipe takes advantage of magnetic-field-aligned coordinates to ensure continuity of the mass den…
▽ More
A recipe is presented to construct an analytic, self-consistent model of a non-barotropic neutron star with a poloidal-toroidal field of arbitrary multipole order, whose toroidal component is confined in a torus around the neutral curve inside the star, as in numerical simulations of twisted tori. The recipe takes advantage of magnetic-field-aligned coordinates to ensure continuity of the mass density at the surface of the torus. The density perturbation and ellipticity of such a star are calculated in general and for the special case of a mixed dipole-quadrupole field as a worked example. The calculation generalises previous work restricted to dipolar, poloidal-toroidal and multipolar, poloidal-only configurations. The results are applied, as an example, to magnetars whose observations (e.g., spectral features and pulse modulation) indicate that the internal magnetic fields may be at least one order of magnitude stronger than the external fields, as inferred from their spin downs, and are not purely dipolar.
△ Less
Submitted 6 January, 2015;
originally announced January 2015.
-
Neutron star deformation due to multipolar magnetic fields
Authors:
Alpha Mastrano,
Paul D. Lasky,
Andrew Melatos
Abstract:
Certain multi-wavelength observations of neutron stars, such as intermittent radio emissions from rotation-powered pulsars beyond the pair-cascade death line, the pulse profile of the magnetar SGR 1900+14 after its 1998 August 27 giant flare, and X-ray spectral features of PSR J0821-4300 and SGR 0418+5729, suggest that the magnetic fields of non-accreting neutron stars are not purely dipolar and m…
▽ More
Certain multi-wavelength observations of neutron stars, such as intermittent radio emissions from rotation-powered pulsars beyond the pair-cascade death line, the pulse profile of the magnetar SGR 1900+14 after its 1998 August 27 giant flare, and X-ray spectral features of PSR J0821-4300 and SGR 0418+5729, suggest that the magnetic fields of non-accreting neutron stars are not purely dipolar and may contain higher-order multipoles. Here, we calculate the ellipticity of a non-barotropic neutron star with (i) a quadrupole poloidal-toroidal field, and (ii) a purely poloidal field containing arbitrary multipoles, deriving the relation between the ellipticity and the multipole amplitudes. We present, as a worked example, a purely poloidal field comprising dipole, quadrupole, and octupole components. We show the correlation between field energy and ellipticity for each multipole, that the l=4 multipole has the lowest energy, and that l=5 has the lowest ellipticity. We show how a mixed multipolar field creates an observationally testable mismatch between the principal axes of inertia (to be inferred from gravitational wave data) and the magnetic inclination angle. Strong quadrupole and octupole components (with amplitudes ~ 10^2 times higher than the dipole) in SGR 0418+5729 still yield ellipticity ~ 10^{-8}, consistent with current gravitational wave upper limits. The existence of higher multipoles in fast-rotating objects (e.g., newborn magnetars) has interesting implications for the braking law and hence phase tracking during coherent gravitational wave searches.
△ Less
Submitted 18 July, 2013; v1 submitted 19 June, 2013;
originally announced June 2013.
-
Stability of magnetic fields in non-barotropic stars: an analytic treatment
Authors:
Taner Akgün,
Andreas Reisenegger,
Alpha Mastrano,
Pablo Marchant
Abstract:
Magnetic fields in upper main-sequence stars, white dwarfs, and neutron stars are known to persist for timescales comparable to their lifetimes. From a theoretical perspective this is problematic, as it can be shown that simple magnetic field configurations are always unstable. In non-barotropic stars, stable stratification allows for a much wider range of magnetic field structures than in barotro…
▽ More
Magnetic fields in upper main-sequence stars, white dwarfs, and neutron stars are known to persist for timescales comparable to their lifetimes. From a theoretical perspective this is problematic, as it can be shown that simple magnetic field configurations are always unstable. In non-barotropic stars, stable stratification allows for a much wider range of magnetic field structures than in barotropic stars, and helps stabilize them by making it harder to induce radial displacements. Recent simulations by Braithwaite and collaborators have shown that, in stably stratified stars, random initial magnetic fields evolve into nearly axisymmetric configurations with both poloidal and toroidal components, which then remain stable for some time. It is desirable to provide an analytic study of the stability of such fields. We write an explicit expression for a plausible equilibrium structure of an axially symmetric magnetic field with both poloidal and toroidal components of adjustable strengths, in a non-barotropic static fluid star, and study its stability using the energy principle. We construct a displacement field that should be a reasonable approximation to the most unstable mode of a toroidal field, and confirm Braithwaite's result that a given toroidal field can be stabilized by a poloidal field containing much less energy than the former. This is consistent with the speculation that the toroidal field is the main reservoir powering magnetar activity. The deformation of a neutron star caused by the hidden toroidal field can also cause emission of gravitational waves.
△ Less
Submitted 6 June, 2013; v1 submitted 1 February, 2013;
originally announced February 2013.
-
Updated gravitational-wave upper limits on the internal magnetic field strength of recycled pulsars
Authors:
Alpha Mastrano,
Andrew Melatos
Abstract:
Recent calculations of the hydromagnetic deformation of a stratified, non-barotropic neutron star are generalized to describe objects with superconducting interiors, whose magnetic permeability μis much smaller than the vacuum value μ_0. It is found that the star remains oblate if the poloidal magnetic field energy is \gtrsim 40% of total magnetic field energy, that the toroidal field is confined…
▽ More
Recent calculations of the hydromagnetic deformation of a stratified, non-barotropic neutron star are generalized to describe objects with superconducting interiors, whose magnetic permeability μis much smaller than the vacuum value μ_0. It is found that the star remains oblate if the poloidal magnetic field energy is \gtrsim 40% of total magnetic field energy, that the toroidal field is confined to a torus which shrinks as μdecreases, and that the deformation is much larger (by a factor \sim μ_0/μ) than in a non-superconducting object. The results are applied to the latest direct and indirect upper limits on gravitational-wave emission from Laser Interferometer Gravitational Wave Observatory (LIGO) and radio pulse timing (spin-down) observations of 81 millisecond pulsars, to show how one can use these observations to infer the internal field strength. It is found that the indirect spin-down limits already imply astrophysically interesting constraints on the poloidal-toroidal field ratio and diamagnetic shielding factor (by which accretion reduces the observable external magnetic field, e.g. by burial). These constraints will improve following gravitational-wave detections, with implications for accretion-driven magnetic field evolution in recycled pulsars and the hydromagnetic stability of these objects' interiors.
△ Less
Submitted 7 December, 2011;
originally announced December 2011.
-
Instability-driven interfacial dynamo in protoneutron stars
Authors:
Alpha Mastrano,
Andrew Melatos
Abstract:
The existence of a tachocline in the Sun has been proven by helioseismology. It is unknown whether a similar shear layer, widely regarded as the seat of magnetic dynamo action, also exists in a protoneutron star. Sudden jumps in magnetic diffusivity $η$ and turbulent vorticity $α$, for example at the interface between the neutron-finger and convective zones, are known to be capable of enhancing me…
▽ More
The existence of a tachocline in the Sun has been proven by helioseismology. It is unknown whether a similar shear layer, widely regarded as the seat of magnetic dynamo action, also exists in a protoneutron star. Sudden jumps in magnetic diffusivity $η$ and turbulent vorticity $α$, for example at the interface between the neutron-finger and convective zones, are known to be capable of enhancing mean-field dynamo effects in a protoneutron star. Here we apply the well-known, plane-parallel, MacGregor-Charbonneau analysis of the Solar interfacial dynamo to the protoneutron star problem and calculate the growth rate analytically under a range of conditions. It is shown that, like the Solar dynamo, it is impossible to achieve self-sustained growth if the discontinuities in $α$, $η$, and shear are coincident and the magnetic diffusivity is isotropic. In contrast, when the jumps in $η$ and $α$ are situated away from the shear layer, self-sustained growth is possible for $P\lesssim 49.8$ ms (if the velocity shear is located at $0.3R$) or $P\lesssim 83.6$ ms (if the velocity shear is located at $0.6R$). This translates into stronger shear and/or $α$-effect than in the Sun. Self-sustained growth is also possible if the magnetic diffusivity if anisotropic, through the ${\bfΩ}\times{\bf{J}}$ effect, even when the $α$, $η$, and shear discontinuities are coincident.
△ Less
Submitted 31 July, 2011;
originally announced August 2011.
-
Gravitational wave emission from a magnetically deformed non-barotropic neutron star
Authors:
Alpha Mastrano,
Andrew Melatos,
Andreas Reisenegger,
Taner Akgün
Abstract:
A strong candidate for a source of gravitational waves is a highly magnetised, rapidly rotating neutron star (magnetar) deformed by internal magnetic stresses. We calculate the mass quadrupole moment by perturbing a zeroth-order hydrostatic equilibrium by an axisymmetric magnetic field with a \emph{linked poloidal-toroidal structure}. In this work, we do \emph{not} require the model star to obey a…
▽ More
A strong candidate for a source of gravitational waves is a highly magnetised, rapidly rotating neutron star (magnetar) deformed by internal magnetic stresses. We calculate the mass quadrupole moment by perturbing a zeroth-order hydrostatic equilibrium by an axisymmetric magnetic field with a \emph{linked poloidal-toroidal structure}. In this work, we do \emph{not} require the model star to obey a barotropic equation of state (as a realistic neutron star is not barotropic), allowing us to explore the hydromagnetic equilibria with fewer constraints. We derive the relation between the ratio of poloidal-to-total field energy $Λ$ and ellipticity $ε$ and briefly compare our results to those obtained using the barotropic assumption. Then, we present some examples of how our results can be applied to astrophysical contexts. First, we show how our formulae, in conjunction with current gravitational wave (non-)detections of the Crab pulsar and the Cassiopeia A central compact object (Cas A CCO), can be used to constrain the strength of the internal toroidal fields of those objects. We find that, for the Crab pulsar (whose canonical equatorial dipole field strength, inferred from spin down, is $4\times 10^8$ T) to emit detectable gravitational radiation, the neutron star must have a strong toroidal field component, with maximum internal toroidal field strength $B_{\mathrm{tm}}=7\times 10^{12}$ T; for gravitational waves to be detected from the Cas A CCO at 300 Hz, $B_{\mathrm{tm}}\sim 10^{13}$ T, whereas detection at 100 Hz would require $B_{\mathrm{tm}}\sim 10^{14}$ T. Using our results, we also show how the gravitational wave signal emitted by a magnetar immediately after its birth (assuming it is born rapidly rotating, with $Λ\lesssim 0.2$) makes such a newborn magnetar a stronger candidate for gravitational wave detection than, for example, an SGR giant flare.
△ Less
Submitted 31 July, 2011;
originally announced August 2011.
-
Nonideal evolution of nonaxisymmetric, force-free magnetic fields in a magnetar
Authors:
Alpha Mastrano,
Andrew Melatos
Abstract:
Recent numerical magnetohydrodynamic calculations by Braithwaite and collaborators support the `fossil field' hypothesis regarding the origin of magnetic fields in compact stars and suggest that the resistive evolution of the fossil field can explain the reorganisation and decay of magnetar magnetic fields. Here, these findings are modelled analytically by allowing the stellar magnetic field to…
▽ More
Recent numerical magnetohydrodynamic calculations by Braithwaite and collaborators support the `fossil field' hypothesis regarding the origin of magnetic fields in compact stars and suggest that the resistive evolution of the fossil field can explain the reorganisation and decay of magnetar magnetic fields. Here, these findings are modelled analytically by allowing the stellar magnetic field to relax through a quasistatic sequence of nonaxisymmetric, force-free states, by analogy with spheromak relaxation experiments, starting from a random field. Under the hypothesis that the force-free modes approach energy equipartition in the absence of resistivity, the output of the numerical calculations is semiquantitatively recovered: the field settles down to a linked poloidal-toroidal configuration, which inflates and becomes more toroidal as time passes. A qualitatively similar (but not identical) end state is reached if the magnetic field evolves by exchanging helicity between small and large scales according to an $α$-dynamo-like, mean-field mechanism, arising from the fluctuating electromotive force produced by the initial random field. The impossibility of matching a force-free internal field to a potential exterior field is discussed in the magnetar context.
△ Less
Submitted 25 April, 2008;
originally announced April 2008.
-
Kelvin Helmholtz Instability and Circulation Transfer at an Isotropic-Anisotropic Superfluid Interface in a Neutron Star
Authors:
A. Mastrano,
A. Melatos
Abstract:
A recent laboratory experiment (Blaauwgeers et al. 2003) suggests that a Kelvin-Helmholtz (KH) instability at the interface between two superfluids, one rotating and anisotropic, the other stationary and isotropic, may trigger sudden spin-up of the stationary superfluid. This result suggests that a KH instability at the crust-core ($^1S_0$-$^3P_2$-superfluid) boundary of a neutron star may provi…
▽ More
A recent laboratory experiment (Blaauwgeers et al. 2003) suggests that a Kelvin-Helmholtz (KH) instability at the interface between two superfluids, one rotating and anisotropic, the other stationary and isotropic, may trigger sudden spin-up of the stationary superfluid. This result suggests that a KH instability at the crust-core ($^1S_0$-$^3P_2$-superfluid) boundary of a neutron star may provide a trigger mechanism for pulsar glitches. We calculate the dispersion relation of the KH instability involving two different superfluids including the normal fluid components and their effects on stability, particularly entropy transport. We show that an entropy difference between the core and crust superfluids reduces the threshold differential shear velocity and threshold crust-core density ratio. We evaluate the wavelength of maximum growth of the instability for neutron star parameters and find the resultant circulation transfer to be within the range observed in pulsar glitches.
△ Less
Submitted 26 May, 2005;
originally announced May 2005.