Neutron star deformation due to poloidal-toroidal magnetic fields of arbitrary multipole order: a new analytic approach

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.Comment: 14 pages, 6 figures, 1 table. Accepted for publication in the Monthly Notices of the Royal Astronomical Societ

Stability of magnetic fields in non-barotropic stars: an analytic treatment

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.Comment: 23 pages, 6 figures; further discussion and an appendix added; corrected typos; accepted by MNRA

(3+2) Neutrino Scheme From A Singular Double See-Saw Mechanism

We obtain a 3+2 neutrino spectrum within a left-right symmetric framework by invoking a singular double see-saw mechanism. Higgs doublets are employed to break $SU_{R}(2)$ and three additional fermions, singlets under the left-right symmetric gauge group, are included. The introduction of a singularity into the singlet fermion Majorana mass matrix results in a light neutrino sector of three neutrinos containing predominantly $\nu_{\alpha L}$, $\alpha=e,\mu,\tau$, separated from two neutrinos containing a small $\nu_{\alpha L}$ component. The resulting active-sterile mixing in the $5\times 5$ mixing matrix is specified once the mass eigenvalues and the $3\times3$ submatrix corresponding to the MNS mixing matrix are known.Comment: 5 pages, matches published versio

Magnetar Oscillations II: spectral method

The seismological dynamics of magnetars is largely determined by a strong hydro-magnetic coupling between the solid crust and the fluid core. In this paper we set up a "spectral" computational framework in which the magnetar's motion is decomposed into a series of basis functions which are associated with the crust and core vibrational eigenmodes. A general-relativistic formalism is presented for evaluation of the core Alfven modes in the magnetic-flux coordinates, as well for eigenmode computation of a strongly magnetized crust of finite thickness. By considering coupling of the crustal modes to the continuum of Alfven modes in the core, we construct a fully relativistic dynamical model of the magnetar which allows: i) Fast and long simulations without numerical dissipation. ii) Very fine sampling of the stellar structure. We find that the presence of strong magnetic field in the crust results in localizing of some high-frequency crustal elasto-magnetic modes with the radial number n>1 to the regions of the crust where the field is nearly horizontal. While the hydro-magnetic coupling of these localized modes to the Alfven continuum in the core is reduced, their energy is drained on a time-scale much less than 1 second. Therefore the puzzle of the observed QPOs with frequencies larger than 600 Hz still stands.Comment: 15 pages, 11 figures, submitted to MNRA

Kelvin Helmholtz Instability and Circulation Transfer at an Isotropic-Anisotropic Superfluid Interface in a Neutron Star

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.Comment: 17 pages, 8 figures, accepted for publication in MNRA