Mutual friction in superfluid neutron stars

We discuss vortex-mediated mutual friction in the two-fluid model for superfluid neutron star cores. Our discussion is based on the general formalism developed by Carter and collaborators, which makes due distinction between transport velocity and momentum for each fluid. This is essential for an implementation of the so-called entrainment effect, whereby the flow of one fluid imparts momentum in the other and vice versa. The mutual friction follows by balancing the Magnus effect that acts on the quantized neutron vortices with resistivity due to the scattering of electrons off of the magnetic field with which each vortex core is endowed. We derive the form of the macroscopic mutual friction force which is relevant for a model based on smooth-averaging over a collection of vortices. We discuss the coefficients that enter the expression for this force, and the time-scale on which the two interpenetrating fluids in a neutron star core are coupled. This discussion confirms that our new formulation accords well with previous work in this area.     

Differential rotation of relativistic superfluid in neutron stars

It is shown how to set up a mathematically elegant and fully relativistic superfluid model that can provide a realistic approximation (neglecting small anisotropies due to crust solidity, magnetic fields, etc., but allowing for the regions with vortex pinning) of the global structure of a rotating neutron star, in terms of just two independently moving constituents. One of these represents the differentially rotating neutron superfluid, while the other part represents the combination of all the other ingredients, including the degenerate electrons, the superfluid protons in the core, and the ions in the crust, the electromagnetic interactions of which will tend to keep them locked together in a state of approximately rigid rotation. Order of magnitude estimates are provided for relevant parameters such as the resistive drag coefficient.     

Numerical simulations of aligned neutron star magnetospheres

We present detailed numerical simulations of the magnetosphere of an isolated neutron star in which the spin and magnetic dipole axes of the star are aligned. We demonstrate that stable charge distributions are always found, rather than particle outflows. A stable magnetosphere consists of a dome above the polar cap containing plasma of one charge and an equatorial belt containing plasma of the other sign: E · B =0 inside both of these. These are separated by a vacuum gap in which E · B ≠0 ( ρ =0 instead). We show that the charge distribution used in the 'standard' Goldreich–Julian pulsar model is inherently unstable: it collapses to a stable configuration that is very similar to the others illustrated here. An instructive video of this collapse is available at http://spacsun.rice.edu/~ian/. For typical pulsars, the stable solution has no particles near to the light cylinder, and if there were any there then their loss from the system would not lead to a replacement from the star (in contradiction to the explicit assumption used in the Goldreich–Julian model). We discuss the generic effects of pair creation, in particular as an additional source of ionization in the vacuum gap. The overall effect is simply to reduce the value of E · B in the vacuum gap so that the pair-production rate drops towards zero. A dome, disc and gap geometry is still the resulting solution. In conclusion, we confirm previous studies that the aligned rotator cannot make an active pulsar.     

Type I and two-gap superconductivity in neutron star magnetism

Neutron star inner cores with several charged baryonic components are likely to be analogues of the two-gap superconductor which is of current interest in condensed-matter physics. Consequently, type I superconductivity is less probable than type II but may nevertheless be present in some intervals of matter density. The intermediate-state structure formed at finite magnetic flux densities after the superconducting transitions is subject to buoyancy, frictional and neutron vortex interaction forces. These are estimated and it is shown that the most important frictional force is that produced by the stable stratification of neutron star matter, the irreversible process being diffusion in the normal, finite magnetic flux density, parts of the structure. The length-scale of the structure, in directions perpendicular to the local magnetic field is of crucial importance. For small scales, the flux comoves with the neutron vortices, as do the proton vortices of a type II superconductor. But for much larger length-scales, flux movement tends to that expected for normal charged Fermi systems.     

First-principles point-defect calculations for solid neutron star matter

Formation enthalpies are calculated for a number of point-defect structures in solid neutron star matter at densities above the neutron-drip threshold. The enthalpies obtained show that an amorphous heterogeneous solid phase is formed at temperatures in a glass transition region, and is likely to persist as the star cools. Its structural differences from the homogeneous body-centred cubic lattice previously assumed make it necessary to reconsider predictions of neutron-star magnetic field evolution, and severely limit the role of conventional superfluid neutron vortex pinning in the interpretation of pulsar glitch phenomena.     

On the dynamics of superfluid neutron star cores

We discuss the nature of the various modes of pulsation of superfluid neutron stars using comparatively simple Newtonian models and the Cowling approximation. The matter in these stars is described in terms of a two-fluid model, where one fluid is the neutron superfluid, which is believed to exist in the core and inner crust of mature neutron stars, and the other fluid represents a conglomerate of all other constituents (crust nuclei, protons, electrons, etc.). In our model, we incorporate the non-dissipative interaction known as the entrainment effect, whereby the momentum of one constituent (e.g. the neutrons) carries along part of the mass of the other constituent. We show that there is no independent set of pulsating g-modes in a non-rotating superfluid neutron star core, even though the linearized superfluid equations contain a well-defined (and real-valued) analogue to the so-called Brunt–Väisälä frequency. Instead, what we find are two sets of spheroidal perturbations whose nature is predominately acoustic. In addition, an analysis of the zero-frequency subspace (i.e. the space of time-independent perturbations) reveals two sets of degenerate spheroidal perturbations, which we interpret to be the missing g-modes, and two sets of toroidal perturbations. We anticipate that the degeneracy of all these zero-frequency modes will be broken by the Coriolis force in the case of rotating stars. To illustrate this we consider the toroidal pulsation modes of a slowly rotating superfluid star. This analysis shows that the superfluid equations support a new class of r-modes, in addition to those familiar from, for example, geophysical fluid dynamics. Finally, the role of the entrainment effect on the superfluid mode frequencies is shown explicitly via solutions to dispersion relations that follow from a 'local' analysis of the linearized superfluid equations.     

r modes and mutual friction in rapidly rotating superfluid neutron stars

We develop a new perturbative framework for studying the r modes of rotating superfluid neutron stars. Our analysis accounts for the centrifugal deformation of the star, and considers the two-fluid dynamics at linear order in the perturbed velocities. Our main focus is on a simple model system where the total density profile is that of an   n = 1  polytrope. We derive a partially analytic solution for the superfluid analogue of the classical r mode. This solution is used to analyse the relevance of the vortex-mediated mutual friction damping, confirming that this dissipation mechanism is unlikely to suppress the gravitational-wave-driven instability in rapidly spinning superfluid neutron stars. Our calculation of the superfluid r modes is significantly simpler than previous approaches, because it decouples the r mode from all other inertial modes of the system. This leads to the results being clearer, but it also means that we cannot comment on the relevance of potential avoided crossings (and associated 'resonances') that may occur for particular parameter values. Our analysis of the mutual friction damping differs from previous studies in two important ways. First, we incorporate realistic pairing gaps which means that the regions of superfluidity in the star's core vary with temperature. Secondly, we allow the mutual friction parameters to take the whole range of permissible values rather than focusing on a particular mechanism. Thus, we consider not only the weak drag regime, but also the strong drag regime where the fluid dynamics are significantly different.     

Kelvin–Helmholtz instability and circulation transfer at an isotropic–anisotropic superfluid interface in a neutron star

A recent laboratory experiment 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 (  ¹ S₀–³ P₂  –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.     

The bottom magnetic field and magnetosphere evolution of neutron star in low-mass X-ray binary

The accretion-induced neutron star (NS) magnetic field evolution is studied through considering the accretion flow to drag the field lines aside and dilute the polar-field strength, and as a result the equatorial field strength increases, which is buried inside the crust on account of the accretion-induced global compression of star crust. The main conclusions of model are as follows: (i) the polar field decays with increase in the accreted mass; (ii) the bottom magnetic field strength of about 10⁸ G can occur when the NS magnetosphere radius approaches the star radius, and it depends on the accretion rate as     ; and (iii) the NS magnetosphere radius decreases with accretion until it reaches the star radius, and its evolution is little influenced by the initial field and the accretion rate after accreting  ∼0.01 M_⊙  , which implies that the magnetosphere radii of NSs in low-mass X-ray binaries would be homogeneous if they accreted the comparable masses. As an extension, the physical effects of the possible strong magnetic zone in the X-ray NSs and recycled pulsars are discussed. Moreover, the strong magnetic fields in the binary pulsars PSR 1831−00 and PSR 1718−19 after accreting about  0.5 M_⊙  in the binary-accretion phase,  8.7 × 10¹⁰  and  1.28 × 10¹² G  , respectively, can be explained through considering the incomplete frozen flow in the polar zone. As an expectation of the model, the existence of the low magnetic field  (∼3 × 10⁷ G)  NSs or millisecond pulsars is suggested.     

Pulsar spin–velocity alignment from single and binary neutron star progenitors

The role of binary progenitors of neutron stars (NSs) in the apparent distribution of space velocities and spin–velocity alignment observed in young pulsars is studied. We performed a Monte Carlo synthesis of pulsar populations originated from single and binary stars with different assumptions about the NS natal kick (kick–spin alignment, kick amplitude and kick reduction in electron-capture supernovae in binary progenitors with initial main-sequence masses from the range  8–11 M_⊙  which experienced mass exchange due to Roche lobe overflow). The calculated spin–velocity alignment in pulsars is compared with data inferred from radio polarization measurements. The observed space velocity of pulsars is found to be mostly affected by the natal kick velocity form and its amplitude; the fraction of binaries is not important here for reasonably large kicks. The natal kick–spin alignment is found to strongly affect the spin–velocity correlation of pulsars. Comparison with the observed pulsar spin–velocity angles favours a sizeable fraction of binary progenitors and kick–spin angles  ∼5°–20°  .     

Can a slowly rotating neutron star be a radio pulsar?

It is shown that the radius of curvature of magnetic field lines in the polar region of a rotating magnetized neutron star can be significantly less than the usual radius of curvature of the dipole magnetic field. The magnetic field in the polar cap is distorted by toroidal electric currents flowing in the neutron star crust. These currents close up the magnetospheric currents driven by the electron–positron plasma generation process in the pulsar magnetosphere. Owing to the decrease in the radius of curvature, electron–positron plasma generation becomes possible even for slowly rotating neutron stars, with   PB ^−2/3₁₂ < 10 s  , where P is the period of star rotation and   B ₁₂= B /10¹² G  is the magnitude of the magnetic field on the star surface.     

Resistive relaxation of a magnetically confined mountain on an accreting neutron star

Three-dimensional numerical magnetohydrodynamic (MHD) simulations are performed to investigate how a magnetically confined mountain on an accreting neutron star relaxes resistively. No evidence is found for non-ideal MHD instabilities on a short time-scale, such as the resistive ballooning mode or the tearing mode. Instead, the mountain relaxes gradually as matter is transported across magnetic surfaces on the diffusion time-scale, which evaluates to  τ_I∼ 10⁵–10⁸ yr  (depending on the conductivity of the neutron star crust) for an accreted mass of   M _a= 1.2 × 10⁻⁴ M_⊙  . The magnetic dipole moment simultaneously re-emerges as the screening currents dissipate over  τ_I  . For non-axisymmetric mountains, ohmic dissipation tends to restore axisymmetry by magnetic reconnection at a filamentary neutral sheet in the equatorial plane. Ideal-MHD oscillations on the Alfvén time-scale, which can be excited by external influences, such as variations in the accretion torque, compress the magnetic field and hence decrease  τ_I  by one order of magnitude relative to its standard value (as computed for the static configuration). The implications of long-lived mountains for gravitational wave emission from low-mass X-ray binaries are briefly explored.