Normalized to: Haskell, B.
[1]
oai:arXiv.org:2007.02748 [pdf] - 2128536
Turbulent, pinned superfluids in neutron stars and pulsar glitch
recoveries
Submitted: 2020-07-06
Pulsar glitches offer an insight into the dynamics of superfluids in the high
density interior of a neutron star. To model these phenomena, however, one
needs to have an understanding of the dynamics of a turbulent array of
superfluid vortices moving through a pinning lattice. In this paper we develop
a theoretical approach to describe vortex mediated mutual friction in a pinned,
turbulent and rotating superfluid. Our model is then applied to the study of
the post glitch rotational evolution in the Vela pulsar and in PSR J0537-6910.
We show that in both cases a turbulent model fits the evolution of the spin
frequency derivative better than a laminar one. We also predict that the second
derivative of the frequency after a glitch should be correlated with the
waiting time since the previous glitch, which we find to be consistent with
observational data for these pulsars. The main conclusion of this paper is that
in the post-glitch rotational evolution of these two pulsars we are most likely
observing the response to the glitch of a pinned turbulent region of the star
(possibly the crust) and not the laminar response of a regular straight vortex
array.
[2]
oai:arXiv.org:2006.09843 [pdf] - 2117136
When the entropy has no maximum: a new perspective on the instability of
the first-order theories of dissipation
Submitted: 2020-06-17
The first-order relativistic fluid theories of dissipation proposed by Eckart
and Landau-Lifshitz have been proved to be unstable. They admit solutions which
start in proximity of equilibrium and depart exponentially from it. We show
that this behaviour is due to the fact that the total entropy of these fluids,
restricted to the dynamically accessible states, has no upper bound. As a
result, these systems have the tendency to constantly change according to the
second law of thermodynamics and the unstable modes represent the directions of
growth of the entropy in state space. We, then, verify that the conditions of
stability of Israel and Stewart's theory are exactly the requirements for the
entropy to have an absolute maximum. Hence, we explain how the instability of
the first-order theories is a direct consequence of the truncation of the
entropy current at the first order, which turns the maximum into a saddle point
of the total entropy. Finally, we show that recently proposed first-order
stable theories, constructed using more general frames, do not solve the
instability problem by providing a maximum for the entropy, but, rather, are
made stable by allowing for small violations of the second law.
[3]
oai:arXiv.org:2001.08951 [pdf] - 2112590
A universal formula for the relativistic correction to the mutual
friction coupling time-scale in neutron stars
Submitted: 2020-01-24, last modified: 2020-05-14
Vortex-mediated mutual friction governs the coupling between the superfluid
and normal components in neutron star interiors. By, for example, comparing
precise timing observations of pulsar glitches with theoretical predictions it
is possible to constrain the physics in the interior of the star, but to do so
an accurate model of the mutual friction coupling in general relativity is
needed. We derive such a model directly from Carter's multifluid formalism, and
study the vortex structure and coupling time-scale between the components in a
relativistic star. We calculate how general relativity modifies the shape and
the density of the quantized vortices and show that, in the quasi-Schwarzschild
coordinates, they can be approximated as straight lines for realistic neutron
star configurations. Finally, we present a simple universal formula (given as a
function of the stellar compactness alone) for the relativistic correction to
the glitch rise-time, which is valid under the assumption that the superfluid
reservoir is in a thin shell in the crust or in the outer core. This universal
relation can be easily employed to correct, a posteriori, any Newtonian
estimate for the coupling time-scale, without any additional computational
expense.
[4]
oai:arXiv.org:2002.10357 [pdf] - 2093206
Magnetic field configurations in neutron stars from MHD simulations
Submitted: 2020-02-24, last modified: 2020-04-29
We have studied numerically the evolution of magnetic fields in barotropic
neutron stars, by performing nonlinear magnetohydrodynamical simulations with
the code PLUTO. For both initially predominantly poloidal and toroidal fields,
with varying strengths, we find that the field settles down to a mixed
poloidal-toroidal configuration, where the toroidal component contributes
between 10% and 20% of the total magnetic energy. This is, however, not a
strict equilibrium, as the instability leads to the development of turbulence,
which in turn gives rise to an inverse helicity cascade, which determines the
final 'twisted torus' setup. The final field configuration is thus dictated by
the non-linear saturation of the instability and is not stationary. The average
energy of the poloidal and toroidal components, however, is approximately
stable in our simulations, and a complex multipolar structure emerges at the
surface, while the magnetic field is dipolar at the exterior boundary, outside
the star.
[5]
oai:arXiv.org:1908.05038 [pdf] - 2075583
Asymmetric accretion and thermal `mountains' in magnetized neutron star
crusts
Submitted: 2019-08-14, last modified: 2020-04-05
Accreting neutron stars are one of the main targets for continuous
gravitational wave searches, as asymmetric accretion may lead to quadrupolar
deformations, or `mountains', on the crust of the star, which source
gravitational wave emission at twice the rotation frequency. The gravitational
wave torque may also impact on the spin evolution of the star, possibly
dictating the currently observed spin periods of neutron stars in Low Mass
X-ray Binaries and leading to the increased spindown rate observed during
accretion in PSR J1023+0038. Previous studies have shown that deformed reaction
layers in the crust of the neutron star lead to thermal and compositional
gradients that can lead to gravitational wave emission. However, there are no
realistic constraints on the level of asymmetry that is expected. In this paper
we consider a natural source of asymmetry, namely the magnetic field, and
calculate the density and pressure perturbations that are expected in the crust
of accreting neutron stars. In general we find that only the outermost reaction
layers of the neutron star are strongly perturbed. The mass quadrupole that we
estimate is generally small and cannot explain the increase of spin-down rate
of PSR J1023+0038. However, if strong shallow heating sources are present at
low densities in the crust, as cooling observations suggest, these layers will
be strongly perturbed and the resulting quadrupole could explain the observed
spindown of PSR J1023+0038, and lead to observable gravitational wave signals
from systems with higher accretion rates.
[6]
oai:arXiv.org:2003.04609 [pdf] - 2061851
Bulk viscosity in relativistic fluids: from thermodynamics to
hydrodynamics
Submitted: 2020-03-10
Bulk viscosity plays a key role in several astrophysical phenomena linked to
neutron stars, from the damping of modes of oscillation in mature stars to the
dynamics of the hot remnant of a binary merger. To quantitatively model these
systems and constrain the imprint of the equation of state of dense matter it
is, however, necessary to resort to a general relativistic description of bulk
viscosity. Several prescriptions exist to modify the standard Navier-Stokes
equations, which are parabolic and lead to non-causal solutions, to obtain
hyperbolic equations that can be used in a relativistic setting. The connection
between microphysics and the parameters introduced is not, however, always
transparent. In this paper we present a relativistic formalism where any
thermodynamic process which contributes to the bulk viscosity is modelled as a
set of chemical reactions, whose reaction coordinates are abstract parameters
describing the displacement from local thermodynamic equilibrium. The result is
a non-equilibrium thermodynamic theory for bulk-viscous fluids which does not
rely on any near-equilibrium assumption, and naturally gives rise to
telegraph-type equations for the reaction processes. We present also two
concrete applications of our theory, to the bulk viscosity in neutron stars and
to a non interacting neutron gas and compare our formulation with that of
Israel and Stewart (1979), which for a mature neutron star emerges as a
perturbative expansion of our theory near equilibrium.
[7]
oai:arXiv.org:2002.04310 [pdf] - 2046571
The effect of non-linear mutual friction on pulsar glitch sizes and rise
times
Submitted: 2020-02-11
Observations of pulsar glitches have the potential to provide constraints on
the dynamics of the high density interior of neutron stars. However, to do so,
realistic glitch models must be constructed and compared to the data. We take a
step towards this goal by testing non-linear models for the mutual friction
force, which is responsible for the exchange of angular momentum between the
neutron superfluid and the observable normal component in a glitch. In
particular, we consider a non-linear dependence of the drag force on the
relative velocity between superfluid vortices and the normal component, in
which the contributions of both kelvon and phonon excitations are included.
This non-linear model produces qualitatively new features, and is able to
reproduce the observed bimodal distribution of glitch sizes in the pulsar
population. The model also suggests that the differences in size distributions
in individual pulsars may be due to the glitches being triggered in regions
with different pinning strengths, as stronger pinning leads to higher vortex
velocities and a qualitatively different mutual friction coupling with respect
to the weak pinning case. Glitches in pulsars that appear to glitch
quasi-periodically with similar sizes may thus be due to the same mechanisms as
smaller events in pulsars that have no preferred glitch size, but simply
originate in stronger pinning regions, possibly in the core of the star.
[8]
oai:arXiv.org:1906.05177 [pdf] - 2011876
Hydrodynamical instabilities in superfluid neutron stars with background
flows between the components
Submitted: 2019-06-12
The interiors of mature neutron stars are expected to be superfluid.
Superfluidity of matter on the microscopic scale can have a number of large
scale, potentially observable consequences, as the superfluid component of the
star can now flow relative to the `normal' component that is tracked by
electromagnetic emission. The most spectacular of such phenomena are pulsar
glitches, sudden spin-ups observed in many pulsars. A background flow of a
normal-fluid component with respect to the superfluid is also known to lead to
a number of instabilities in laboratory superfluids, possibly leading to
turbulence and modifying the nature of the mutual friction coupling between the
two fluids. In this paper we consider modes of oscillation in the crust and
core of a superfluid neutron star, by conducting a plane-wave analysis. We
explicitly account for a background flow between the two components (as would
be expected in the presence of pinning) and the entrainment, and we consider
both standard (Hall-Vinen) and isotropic (Gorter-Mellink) forms of the mutual
friction. We find that for standard mutual friction there are families of
unstable inertial and sound waves both in the case of a counter-flow along the
superfluid vortex axis and for counterflow perpendicular to the vortex axis and
find that entrainment leads to a quantitative difference between instabilities
in the crust and core of the star. For isotropic mutual friction we find no
unstable modes, and speculate that instabilities in a straight vortex array may
be linked to glitching behaviour, which then ceases until the turbulence has
decayed.
[10]
oai:arXiv.org:1905.02877 [pdf] - 1886626
Collective, glitch-like vortex motion in a neutron star with an annular
pinning barrier
Submitted: 2019-05-07
Neutron star glitches are commonly believed to occur, when angular momentum
is transferred suddenly from the star's interior to the crust by the collective
unpinning and repinning of large numbers of superfluid vortices. In general,
the pinning potential associated with nuclei in the crustal lattice varies as a
function of radius. We explore vortex dynamics under these conditions by
solving the three-dimensional Gross-Pitaevskii equation in a rotating, harmonic
trap with an axisymmetric `moat' of deeper pinning sites on an otherwise
uniform, corotating pinning grid. The moat is designed to resemble crudely a
radially dependent pinning profile in a neutron star crust, although the values
of the pinning potential are not astrophysically realistic due to computational
constraints. It is shown that vortices accumulate in the moat, inducing large
differential rotation which can trigger mass unpinning events. It is also shown
that the system self-adjusts, such that the net vortex flux out of the system
is the same with and without a moat, as the trap spins down, but glitches are
less frequent and larger when the moat is present. The results, generated for
an idealized system, represent a first step towards including stratified
pinning in quantum mechanical models of neutron star glitches.
[11]
oai:arXiv.org:1806.02822 [pdf] - 1737524
Evidence for a minimum ellipticity in millisecond pulsars
Submitted: 2018-06-07, last modified: 2018-08-22
Neutron stars spin down over time due to a number of energy-loss processes.
We provide tantalizing population-based evidence that millisecond pulsars
(MSPs) have a minimum ellipticity of $\epsilon\approx10^{-9}$ around their spin
axis and that, consequently, some spin down mostly through gravitational-wave
emission. We discuss the implications of such a minimum ellipticity in terms of
the internal magnetic field strengths and nuclear matter composition of neutron
stars and show it would result in the Advanced LIGO and Virgo
gravitational-wave detectors, or their upgrades, detecting gravitational waves
from some known MSPs in the near future.
[12]
oai:arXiv.org:1806.10168 [pdf] - 1706616
Core and crust contributions in pulsar glitches: constraints from the
slow rise of the largest glitch observed in the Crab pulsar
Submitted: 2018-06-26
Pulsar glitches are attributed to the sudden re-coupling of very weakly
coupled large scale superfluid components in the neutron star interior. This
process leads to rapid exchange of angular momentum and an increase in spin
frequency. The transfer of angular momentum is regulated by a dissipative
mutual friction, whose strength defines the spin-up timescale of a glitch.
Hence, observations of glitch rises can be used to shed light on the dominant
microphysical interactions at work in the high density interior of the star. We
present a simple analytical model, complemented with more detailed numerical
simulations, which produces a fast spin-up followed by a more gradual rise.
Such features are observed in some large glitches of the Crab pulsar, including
the largest recent glitch of 2017. We also use observations to constrain the
mutual friction coefficient of the glitch-driving region for two possible
locations: the inner crust and outer core of the star. We find that the
features of Crab glitches require smaller values of the mutual friction
coefficient than those needed to explain the much faster Vela spin-ups. This
suggests a crustal origin for the former but an outer core contribution for the
latter.
[13]
oai:arXiv.org:1801.01413 [pdf] - 1705181
Modelling pulsar glitches: the hydrodynamics of superfluid vortex
avalanches in neutron stars
Submitted: 2018-01-04, last modified: 2018-06-25
The dynamics of quantised vorticity in neutron star interiors is at the heart
of most pulsar glitch models. However, the large number of vortices (up to
$\approx 10^{13}$) involved in a glitch and the huge disparity in scales
between the femtometer scale of vortex cores and the kilometre scale of the
star makes quantum dynamical simulations of the problem computationally
intractable. In this paper we take a first step towards developing a mean field
prescription to include the dynamics of vortices in large scale hydrodynamical
simulations of superfluid neutron stars. We consider a one dimensional setup
and show that vortex accumulation and differential rotation in the neutron
superfluid lead to propagating waves, or `avalanches', as solutions for the
equations of motion for the superfluid velocities. We introduce an additional
variable, the fraction of free vortices, and test different prescriptions for
its advection with the superfluid flow. We find that the new terms lead to
solutions with a linear component in the rise of a glitch, and that, in
specific setups, they can give rise to glitch precursors and even to decreases
in frequency, or `anti-glitches'.
[14]
oai:arXiv.org:1805.11277 [pdf] - 1790738
Fundamental physics and the absence of sub-millisecond pulsars
Submitted: 2018-05-29
Observations of the spin distribution of rapidly rotating neutron stars show
evidence for a lack of stars spinning at frequencies larger than $f\approx 700$
Hz, well below the predictions of theoretical equations of state. This has
generally been taken as evidence of an additional spin-down torque operating in
these systems and it has been suggested that gravitational wave torques may be
operating and be linked to a potentially observable signal. In this paper we
aim to determine whether additional spin-down torques are necessary, or whether
the observed limit of $f\approx 700$ Hz could correspond to the mass-shedding
frequency for the observed systems and is simply a consequence of the,
currently unknown, state of matter at high densities. Given our ignorance with
regard to the true equation of state of matter above nuclear saturation
densities, we make minimal physical assumption and only demand causality in the
core. We then connect our causally-limited equation of state to a realistic
microphysical crustal equation of state for densities below nuclear saturation
density. This produces a limiting model that will give the lowest possible
maximum frequency, which we compare to observational constraints on neutron
star masses and frequencies. We also compare our findings with the constraints
on the tidal deformability obtained in the observations of the GW170817 event.
We find that the lack of pulsars spinning faster than $f\approx 700$ Hz is not
compatible with our causal limited `minimal' equation of state, for which the
breakup frequency cannot be lower than $f_{\rm max}\approx 1200$ Hz. A low
frequency cutoff, around $f\approx 800$ Hz could only be possible if we assume
that these systems do not contain neutron stars with masses above $M\approx 2
M_\odot$. This would have to be due either to selection effects, or possibly to
a phase transition in the interior of the neutron star.
[15]
oai:arXiv.org:1803.08813 [pdf] - 1739900
Estimating the EOS from the measurement of NS radii with 5% accuracy
Submitted: 2018-03-22, last modified: 2018-05-23
Observations of heavy (${\simeq}2\,M_\odot$) neutron stars in addition to the
recent measurement of tidal deformability from the binary neutron-star merger
GW170817, place interesting constraints on theories of dense matter. Current
and future observatories, such as the NICER and ATHENA are expected to collect
information on the global parameters of neutron stars, namely masses and radii,
with the accuracy of a few percent. Such accuracy will allow for precise
comparisons of measurements to models of compact objects. Here we investigate
how the measurement accuracy of the NICER and ATHENA missions will improve our
understanding of the dense-matter interior of neutron stars. We compare global
parameters of stellar configurations obtained using three different equations
of state: a reference (SLy4 EOS) and two piecewise polytropes manufactured to
produce mass-radius relations indistinguishable from the observational point of
view i.e. within the predicted error of the radius measurement. We assume
observational errors on the radius determination corresponding to the expected
accuracies. The effect of rotation is examined using high-precision numerical
relativity computations. Due to the fact that masses and rotational frequencies
might be determined very precisely in the most optimistic scenario, only the
influence of observational errors on the radius measurements is investigated.
We show that ${\pm}5\%$ errors in radius measurement lead to ${\sim}10\%$ and
${\sim}40\%$ accuracy in central parameter estimation, for low-mass and
high-mass neutron stars, respectively. Global parameters, such as oblateness
and surface area, can be established with $8-10\%$ accuracy, even if only
compactness (instead of mass and radius) is measured. We also report on the
range of tidal deformabilities corresponding to the estimated masses of
GW170817, for the assumed uncertainty in radius.
[16]
oai:arXiv.org:1611.10223 [pdf] - 1635681
Constraints on pulsar masses from the maximum observed glitch
Submitted: 2016-11-30, last modified: 2018-02-15
Pulsar glitches, sudden jumps in frequency in otherwise steadily spinning
down radio pulsars, offer a unique glimpse into the superfluid interior of
neutron stars. The exact trigger of these events remains, however, elusive and
this has hampered attempts to use glitch observations to constrain fundamental
physics. In this paper we propose a new method to measure the mass of glitching
pulsars, using observations of the maximum glitch recorded in a star, together
with state of the art microphysical models of the pinning interaction between
superfluid vortices and ions in the crust. Studying systematically all the
presently observed large glitchers, we find an inverse correlation between size
of the maximum glitch and the pulsar mass. Our procedure will allow current and
future observations of glitching pulsars to constrain not only the physics of
glitch models but also the equation of state of dense matter in neutron star
interiors.
[17]
oai:arXiv.org:1712.01547 [pdf] - 1736245
Probing neutron star interiors with pulsar glitches
Submitted: 2017-12-05
Pulsar glitches are thought to be probes of the superfluid interior of
neutron stars. These sudden jumps in frequency observed in many pulsars are
generally assumed to be the macroscopic manifestation of superfluid vortex
motion on a microscopic scale. Resolving and modelling such phenomena on the
scale of a neutron star is, however, a challenging problem which still remains
open, fifty years after the discovery of pulsars. In this article I will review
recent theoretical progress, both on the microscopic level and on the
macroscopic level, and discuss which constraints on the models can be provided
by observations.
[18]
oai:arXiv.org:1711.05550 [pdf] - 1749239
The enigmatic spin evolution of PSR J0537-6910: r-modes, gravitational
waves and the case for continued timing
Submitted: 2017-11-15
We discuss the unique spin evolution of the young X-ray pulsar PSR
J0537-6910, a system in which the regular spin down is interrupted by glitches
every few months. Drawing on the complete timing data from the Rossi X-ray
Timing Explorer (RXTE, from 1999-2011), we argue that a trend in the
inter-glitch behaviour points to an effective braking index close to $n=7$,
much larger than expected. This value is interesting because it would accord
with the neutron star spinning down due to gravitational waves from an unstable
r-mode. We discuss to what extent this, admittedly speculative, scenario may be
consistent and if the associated gravitational-wave signal would be within
reach of ground based detectors. Our estimates suggest that one may, indeed, be
able to use future observations to test the idea. Further precision timing
would help enhance the achievable sensitivity and we advocate a joint observing
campaign between the Neutron Star Interior Composition ExploreR (NICER) and the
LIGO-Virgo network.
[19]
oai:arXiv.org:1709.10340 [pdf] - 1842290
Superfluidity and Superconductivity in Neutron Stars
Submitted: 2017-09-29
This review focuses on applications of the ideas of superfluidity and
superconductivity in neutron stars in a broader context, ranging from the
microphysics of pairing in nucleonic superfluids to macroscopic manifestations
of superfluidity in pulsars. The exposition of the basics of pairing, vorticity
and mutual friction can serve as an introduction to the subject. We also review
some topics of recent interest, including the various types of pinning of
vortices, glitches, and oscillations in neutron stars containing superfluid
phases of baryonic matter.
[20]
oai:arXiv.org:1705.07669 [pdf] - 1598049
The Spin Distribution of Fast Spinning Neutron Stars in Low Mass X-Ray
Binaries: Evidence for Two Sub-Populations
Submitted: 2017-05-22
We study the current sample of rapidly rotating neutron stars in both
accreting and non-accreting binaries in order to determine whether the spin
distribution of accreting neutron stars in low-mass X-ray binaries can be
reconciled with current accretion torque models. We perform a statistical
analysis of the spin distributions and show that there is evidence for two
sub-populations among low-mass X-ray binaries, one at relatively low spin
frequency, with an average of ~300 Hz and a broad spread, and a peaked
population at higher frequency with average spin frequency of ~575 Hz. We show
that the two sub-populations are separated by a cut-point at a frequency of
~540 Hz. We also show that the spin frequency of radio millisecond pulsars does
not follow a log-normal distribution and shows no evidence for the existence of
distinct sub-populations. We discuss the uncertainties of different accretion
models and speculate that either the accreting neutron star cut-point marks the
onset of gravitational waves as an efficient mechanism to remove angular
momentum or some of the neutron stars in the fast sub-population do not evolve
into radio millisecond pulsars.
[21]
oai:arXiv.org:1703.08374 [pdf] - 1582108
Are gravitational waves spinning down PSR J1023+0038?
Submitted: 2017-03-24
The pulsar J1203+0038 rotates with a frequency $\nu\approx 592$ Hz and has
been observed to transition between a radio state, during which it is visible
as a millisecond radio pulsar, and and a Low Mass X-ray Binary state, during
which accretion powered X-ray pulsations are visible. Timing during the two
phases reveals that during the LMXB phase the neutron star is spinning down at
a rate of $\dot{\nu}\approx -3 \times 10^{-15}$ Hz/s, which is approximately
27\% faster than the rate measured during the radio phase, $\dot{\nu}\approx
-2.4 \times 10^{-15}$ Hz/s, and at odds with the predictions of accretion
models. In this letter we suggest that the increase in spin-down rate is
compatible with gravitational wave emission, and in particular to the creation
of a `mountain' during the accretion phase. We show that asymmetries in
pycno-nuclear reaction rates in the crust can lead to a large enough mass
quadrupole to explain the observed spin-down rate, which so far has no other
self-consistent explanation, and that radio timing at the onset of the next
millisecond radio pulsar phase can test this scenario. Another possibility is
that an unstable $r$-mode with amplitude $\alpha\approx 5\times10^{-8}$ may be
present in the system.
[22]
oai:arXiv.org:1512.07903 [pdf] - 1411264
Hydrodynamic simulations of pulsar glitch recovery
Submitted: 2015-12-08, last modified: 2016-05-03
Glitches are sudden jumps in the spin frequency of pulsars believed to
originate in the superfluid interior of neutron stars. Superfluid flow in a
model neutron star is simulated by solving the equations of motion of a
two-component superfluid consisting of a viscous proton-electron plasma and an
inviscid neutron condensate in a spherical Couette geometry. We examine the
response of our model neutron star to glitches induced in three different ways:
by instantaneous changes of the spin frequency of the inner and outer
boundaries, and by instantaneous recoupling of the fluid components in the
bulk. All simulations are performed with strong and weak mutual friction. It is
found that the maximum size of a glitch that originates in the bulk decreases
as the mutual friction strengthens. It is also found that mutual friction
determines the fraction of the frequency jump which is later recovered, a
quantity known as the 'healing parameter'. These behaviours may explain some of
the diversity in observed glitch recoveries.
[23]
oai:arXiv.org:1603.04304 [pdf] - 1436327
The effect of superfluid hydrodynamics on pulsar glitch sizes and
waiting times
Submitted: 2016-03-14
Pulsar glitches, sudden jumps in frequency observed in many radio pulsars,
may be the macroscopic manifestation of superfluid vortex avalanches on the
microscopic scale. Small scale quantum mechanical simulations of vortex motion
in a decelerating container have shown that such events are possible and
predict power-law distributions for the size of the events, and exponential
distributions for the waiting time. Despite a paucity of data, this prediction
is consistent with the size and waiting time distributions of most glitching
pulsars. Nevertheless a few object appear to glitch quasi-periodically, and
exhibit many large glitches, while a recent study of the Crab pulsar has
suggested a cut-off deviations from a power-law distribution for smaller
glitches. In this paper we incorporate the results of quantum mechanical
simulations in a macroscopic scale superfluid hydrodynamics simulation. We show
that the effect of vortex coupling to the neutron and proton fluids in the
neutron star naturally leads to deviations from power-law distributions for
sizes and from exponential distributions for waiting times. In particular we
predict a cut-off in the size distribution for small glitches.
[24]
oai:arXiv.org:1408.4298 [pdf] - 1307156
Mesoscopic pinning forces in neutron star crusts
Submitted: 2014-08-19, last modified: 2015-11-07
The crust of a neutron star is thought to be comprised of a lattice of nuclei
immersed in a sea of free electrons and neutrons. As the neutrons are
superfluid their angular momentum is carried by an array of quantized vortices.
These vortices can pin to the nuclear lattice and prevent the neutron
superfluid from spinning down, allowing it to store angular momentum which can
then be released catastrophically, giving rise to a pulsar glitch. A crucial
ingredient for this model is the maximum pinning force that the lattice can
exert on the vortices, as this allows us to estimate the angular momentum that
can be exchanged during a glitch. In this paper we perform, for the first time,
a detailed and quantitative calculation of the pinning force \emph{per unit
length} acting on a vortex immersed in the crust and resulting from the
mesoscopic vortex-lattice interaction. We consider realistic vortex tensions,
allow for displacement of the nuclei and average over all possible orientation
of the crystal with respect to the vortex. We find that, as expected, the
mesoscopic pinning force becomes weaker for longer vortices and is generally
much smaller than previous estimates, based on vortices aligned with the
crystal. Nevertheless the forces we obtain still have maximum values of order
$f_{\rm{pin}}\approx 10^{15}$ dyn/cm, which would still allow for enough
angular momentum to be stored in the crust to explain large Vela glitches, if
part of the star is decoupled during the event.
[25]
oai:arXiv.org:1510.03136 [pdf] - 1422171
Pinned vortex hopping in a neutron star crust
Submitted: 2015-10-12
The motion of superfluid vortices in a neutron star crust is at the heart of
most theories of pulsar glitches. Pinning of vortices to ions can decouple the
superfluid from the crust and create a reservoir of angular momentum. Sudden
large scale unpinning can lead to an observable glitch. In this paper we
investigate the scattering of a free vortex off a pinning potential and
calculate its mean free path, in order to assess whether unpinned vortices can
skip multiple pinning sites and come close enough to their neighbours to
trigger avalanches, or whether they simply hop from one pinning site to another
giving rise to a more gradual creep. We find that there is a significant range
of parameter space in which avalanches can be triggered, thus supporting the
hypothesis that they may lie at the origin of pulsar glitches. For realistic
values of the pinning force and superfluid drag parameters we find that
avalanches are more likely in the higher density regions of the crust where
pinning is stronger. Physical differences in stellar parameters, such as mass
and temperature, may lead to a switch between creep-like motion and avalanches,
explaining the different characteristics of glitching pulsars.
[26]
oai:arXiv.org:1509.04370 [pdf] - 1277251
R-modes in neutron stars: Theory and observations
Submitted: 2015-09-14
In this article I will review the theory behind the gravitational wave driven
r-mode instability in rapidly rotating neutron stars and discuss which
constraints can be derived from observations of spins and temperatures in Low
Mass X-ray Binaries. I will discuss how a standard, `minimal' neutron star
model is not consistent with the data, and discuss some of the additional
physical mechanisms that could reconcile theory with observations. In
particular I will focus on additional forms of damping due to exotic cores and
on strong mutual friction due to superfluid vortices cutting through
superconducting flux tubes, and examine the repercussions these effects could
have on the saturation amplitude of the mode. Finally I will also discuss the
possibility that oscillations due to r-modes may have been recently observed in
the X-ray light curves of two Low Mass X-ray Binaries.
[27]
oai:arXiv.org:1506.01445 [pdf] - 1303998
Observational constraints on neutron star crust-core coupling during
glitches
Submitted: 2015-06-03
We demonstrate that observations of glitches in the Vela pulsar can be used
to investigate the strength of the crust-core coupling in a neutron star, and
suggest that recovery from the glitch is dominated by torque exerted by the
re-coupling of superfluid components of the core that were decoupled from the
crust during the glitch. Assuming that the recoupling is mediated by mutual
friction between the superfluid neutrons and the charged components of the
core, we use the observed magnitudes and timescales of the shortest timescale
components of the recoveries from two recent glitches in the Vela pulsar to
infer the fraction of the core that is coupled to the crust during the glitch,
and hence spun up by the glitch event. Within the framework of a two-fluid
hydrodynamic model of glitches, we analyze whether crustal neutrons alone are
sufficient to drive the glitch activity observed in the Vela pulsar. We use two
sets of neutron star equations of state (EOSs), both of which span crust and
core consistently and cover a range of the slope of the symmetry energy at
saturation density $30 < L <120$ MeV. One set produces maximum masses
$\approx$2.0$M_{\odot}$, the second $\approx$2.6$M_{\odot}$. We also include
the effects of entrainment of crustal neutrons by the superfluid lattice. We
find that for medium to stiff EOSs, observations imply $>70\%$ of the moment of
inertia of the core is coupled to the crust during the glitch, though for
softer EOSs $L\approx 30$MeV as little as $5\%$ could be coupled. No EOS is
able to reproduce the observed glitch activity with crust neutrons alone, but
extending the region where superfluid vortices are strongly pinned into the
core by densities as little as 0.016fm$^{-3}$ above the crust-core transition
density restores agreement with the observed glitch activity.
[28]
oai:arXiv.org:1502.07062 [pdf] - 940402
Models of Pulsar Glitches
Submitted: 2015-02-25
Radio pulsars provide us with some of the most stable clocks in the universe.
Nevertheless several pulsars exhibit sudden spin-up events, known as glitches.
More than forty years after their first discovery, the exact origin of these
phenomena is still open to debate. It is generally thought that they an
observational manifestation of a superfluid component in the stellar interior
and provide an insight into the dynamics of matter at extreme densities. In
recent years there have been several advances on both the theoretical and
observational side, that have provided significant steps forward in our
understanding of neutron star interior dynamics and possible glitch mechanisms.
In this article we review the main glitch models that have been proposed and
discuss our understanding, in the light of current observations.
[29]
oai:arXiv.org:1501.06039 [pdf] - 1043096
Detecting gravitational waves from mountains on neutron stars in the
Advanced Detector Era
Submitted: 2015-01-24
Rapidly rotating Neutron Stars (NSs) in Low Mass X-ray Binaries (LMXBs) are
thought to be interesting sources of Gravitational Waves (GWs) for current and
next generation ground based detectors, such as Advanced LIGO and the Einstein
Telescope. The main reason is that many of the NS in these systems appear to be
spinning well below their Keplerian breakup frequency, and it has been
suggested that torques associated with GW emission may be setting the observed
spin period. This assumption has been used extensively in the literature to
assess the strength of the likely gravitational wave signal. There is now,
however, a significant amount of theoretical and observation work that suggests
that this may not be the case, and that GW emission is unlikely to be setting
the spin equilibrium period in many systems. In this paper we take a different
starting point and predict the GW signal strength for two physical mechanisms
that are likely to be at work in LMXBs: crustal mountains due to thermal
asymmetries and magnetically confined mountains. We find that thermal crustal
mountains in transient LMXBs are unlikely to lead to detectable GW emission,
while persistent systems are good candidates for detection by Advanced LIGO and
by the Einstein Telescope. Detection prospects are pessimistic for the magnetic
mountain case, unless the NS has a buried magnetic field of $B\approx 10^{12}$
G, well above the typically inferred exterior dipole fields of these objects.
Nevertheless, if a system were to be detected by a GW observatory, cyclotron
resonant scattering features in the X-ray emission could be used to distinguish
between the two different scenarios.
[30]
oai:arXiv.org:1407.8254 [pdf] - 1216034
Gravitational waves from rapidly rotating neutron stars
Haskell, Brynmor;
Andersson, Nils;
D`Angelo, Caroline;
Degenaar, Nathalie;
Glampedakis, Kostas;
Ho, Wynn C. G.;
Lasky, Paul D.;
Melatos, Andrew;
Oppenoorth, Manuel;
Patruno, Alessandro;
Priymak, Maxim
Submitted: 2014-07-30
Rapidly rotating neutron stars in Low Mass X-ray Binaries have been proposed
as an interesting source of gravitational waves. In this chapter we present
estimates of the gravitational wave emission for various scenarios, given the
(electromagnetically) observed characteristics of these systems. First of all
we focus on the r-mode instability and show that a 'minimal' neutron star model
(which does not incorporate exotica in the core, dynamically important magnetic
fields or superfluid degrees of freedom), is not consistent with observations.
We then present estimates of both thermally induced and magnetically sustained
mountains in the crust. In general magnetic mountains are likely to be
detectable only if the buried magnetic field of the star is of the order of
$B\approx 10^{12}$ G. In the thermal mountain case we find that gravitational
wave emission from persistent systems may be detected by ground based
interferometers. Finally we re-asses the idea that gravitational wave emission
may be balancing the accretion torque in these systems, and show that in most
cases the disc/magnetosphere interaction can account for the observed spin
periods.
[31]
oai:arXiv.org:1311.1352 [pdf] - 791769
Nuclear Equation of State from Observations of Short Gamma-Ray Burst
Remnants
Submitted: 2013-11-06, last modified: 2014-02-11
The favoured progenitor model for short $\gamma$-ray bursts (SGRBs) is the
merger of two neutron stars that triggers an explosion with a burst of
collimated $\gamma$-rays. Following the initial prompt emission, some SGRBs
exhibit a plateau phase in their $X$-ray light curves that indicates additional
energy injection from a central engine, believed to be a rapidly rotating,
highly magnetised neutron star. The collapse of this `protomagnetar' to a black
hole is likely to be responsible for a steep decay in $X$-ray flux observed at
the end of the plateau. In this letter, we show that these observations can be
used to effectively constrain the equation of state of dense matter. In
particular, we show that the known distribution of masses in binary neutron
star systems, together with fits to the $X$-ray light curves, provide
constraints that exclude the softest and stiffest plausible equations of state.
We further illustrate how a future gravitational wave observation with Advanced
LIGO/Virgo can place tight constraints on the equation of state, by adding into
the picture a measurement of the chirp mass of the SGRB progenitor.
[32]
oai:arXiv.org:1309.3885 [pdf] - 843048
On the universality of I-Love-Q relations in magnetized neutron stars
Submitted: 2013-09-16, last modified: 2014-01-28
Recently, general relations among the quadrupole moment (Q), the moment of
inertia (I), and the tidal deformability (Love number) of a neutron star were
shown to exist. They are nearly independent of the nuclear matter equation of
state and would be of great aid in extracting parameters from observed
gravitational-waves and in testing general relativity. These relations,
however, do not account for strong magnetic fields. We consider this problem by
studying the effect of a strong magnetic field on slowly rotating relativistic
neutron stars and show that, for simple magnetic field configurations that are
purely poloidal or purely toroidal, the relation between Q and I is again
nearly universal. However, different magnetic field geometries lead to
different I-Q relations, and, in the case of a more realistic twisted-torus
magnetic field configuration, the relation depends significantly on the
equation of state, losing its universality. I-Love-Q relations must thus be
used with very great care, since universality is lost for stars with long spin
periods, i.e. P > 10 s, and strong magnetic fields, i.e. B > 10^12 G.
[33]
oai:arXiv.org:1307.0985 [pdf] - 1172465
A new mechanism for saturating unstable r-modes in neutron stars
Submitted: 2013-07-03
We consider a new mechanism for damping the oscillations of a mature neutron
star. The new dissipation channel arises if superfluid vortices are forced to
cut through superconducting fluxtubes. This mechanism is interesting because
the oscillation modes need to exceed a critical amplitude in order for it to
operate. Once it acts the effect is very strong (and nonlinear) leading to
efficient damping. The upshot of this is that modes are unlikely to ever evolve
far beyond the critical amplitude. We consider the effect of this new
dissipation channel on the r-modes, that may be driven unstable by the emission
of gravitational waves. Our estimates show that the fluxtube cutting leads to a
saturation threshold for the instability that can be smaller than that of other
proposed mechanisms. This suggests that the idea may be of direct astrophysical
relevance.
[34]
oai:arXiv.org:1306.5214 [pdf] - 1172205
Glitch recoveries in radio-pulsars and magnetars
Submitted: 2013-06-21
Pulsar glitches are sudden increases in the spin frequency of an otherwise
steadily spinning down neutron star. These events are thought to represent a
direct probe of the dynamics of the superfluid interior of the star. However
glitches can differ significantly from one another, not only in size and
frequency, but also in the post-glitch response of the star. Some appear as
simple steps in frequency, while others also display an increase in spin-down
rate after the glitch. Others still show several exponentially relaxing
components in the post-glitch recovery. We show that if glitches are indeed due
to large scale unpinning of superfluid vortices, the different regions in which
this occurs and respective timescales on which they recouple can lead to the
various observed signatures. Furthermore we show that this framework naturally
accounts for the peculiar relaxations of glitches in Anomalous X-ray Pulsars.
[35]
oai:arXiv.org:1303.3282 [pdf] - 639877
The hottest superfluid and superconductor in the Universe: Discovery and
nuclear physics implications
Submitted: 2013-03-13, last modified: 2013-03-15
We present recent work on using astronomical observations of neutron stars to
reveal unique insights into nuclear matter that cannot be obtained from
laboratories on Earth. First, we discuss our measurement of the rapid cooling
of the youngest neutron star in the Galaxy; this provides the first direct
evidence for superfluidity and superconductivity in the supra-nuclear core of
neutron stars. We show that observations of thermonuclear X-ray bursts on
neutron stars can be used to constrain properties of neutron superfluidity and
neutrino emission. We describe the implications of rapid neutron star rotation
rates on aspects of nuclear and superfluid physics. Finally, we show that
entrainment coupling between the neutron superfluid and the nuclear lattice
leads to a less mobile crust superfluid; this result puts into question the
conventional picture of pulsar glitches as being solely due to the crust
superfluid and suggests that the core superfluid also participates.
[36]
oai:arXiv.org:1209.6260 [pdf] - 1151694
Investigating superconductivity in neutron star interiors with glitch
models
Submitted: 2012-09-27
The high density interior of a neutron star is expected to contain
superconducting protons and superfluid neutrons. Theoretical estimates suggest
that the protons will form a type II superconductor in which the stellar
magnetic field is carried by flux tubes. The strong interaction between the
flux tubes and the neutron rotational vortices could lead to strong 'pinning',
i.e. vortex motion could be impeded. This has important implications especially
for pulsar glitch models as it would lead to a large part of the vorticity of
the star being decoupled from the 'normal' component, to which the
electromagnetic emission is locked. In this paper we explore the consequences
of strong pinning in the core on the 'snowplow' model for pulsar glitches
(Pizzochero 2011), making use of realistic equations of state and relativistic
background models for the neutron star. We find that in general a large
fraction of pinned vorticity in the core is not compatible with observations of
giant glitches in the Vela pulsar. The conclusion is thus that either most of
the core is in a type I superconducting state or that the interaction between
vortices and flux tubes is weaker than previously assumed.
[37]
oai:arXiv.org:1205.6647 [pdf] - 1123761
The effect of realistic equations of state and general relativity on the
"snowplow" model for pulsar glitches
Submitted: 2012-05-30
Many pulsars are observed to "glitch", i.e. show sudden jumps in their
rotational frequency $\nu$, some of which can be as large as $\Delta
\nu/\nu\approx 10^{-6}-10^{-5}$ in a subset of pulsars known as giant
glitchers. Recently Pizzochero (2011) has shown that an analytic model based on
realistic values for the pinning forces in the crust and for the angular
momentum transfer in the star can describe the average properties of giant
glitches, such as the inter-glitch waiting time, the step in frequency and that
in frequency derivative. In this paper we extend the model (originally
developed in Newtonian gravity and for a polytropic equation of state) to
realistic backgrounds obtained by integrating the relativistic equations of
stellar structure and using physically motivated equations of state to describe
matter in the neutron star. We find that this more detailed treatment still
reproduces the main features of giant glitches in the Vela pulsar and allows us
to set constraints on the equation of state. In particular we find that stiffer
equations of state are favoured and that it is unlikely that the Vela pulsar
has a high mass (larger than $M\approx 1.5 M_\odot$).
[38]
oai:arXiv.org:1204.2894 [pdf] - 1117948
The dynamics of dissipative multi-fluid neutron star cores
Submitted: 2012-04-13
We present a Newtonian multi-fluid formalism for superfluid neutron star
cores, focussing on the additional dissipative terms that arise when one takes
into account the individual dynamical degrees of freedom associated with the
coupled "fluids". The problem is of direct astrophysical interest as the nature
of the dissipative terms can have significant impact on the damping of the
various oscillation modes of the star and the associated gravitational-wave
signatures. A particularly interesting application concerns the
gravitational-wave driven instability of f- and r-modes. We apply the developed
formalism to two specific three-fluid systems: (i) a hyperon core in which both
Lambda and Sigma^- hyperons are present, and (ii) a core of deconfined quarks
in the colour-flavour-locked phase in which a population of neutral K^0 kaons
is present. The formalism is, however, general and can be applied to other
problems in neutron-star dynamics (such as the effect of thermal excitations
close to the superfluid transition temperature) as well as laboratory
multi-fluid systems.
[39]
oai:arXiv.org:1201.2101 [pdf] - 1092846
Constraining the physics of the r-mode instability in neutron stars with
X-ray and UV observations
Submitted: 2012-01-10
Rapidly rotating Neutron Stars in Low Mass X-ray Binaries (LMXBs) may be an
interesting source of Gravitational Waves (GWs). In particular, several modes
of stellar oscillation may be driven unstable by GW emission, and this can lead
to a detectable signal. Here we illustrate how current X-ray and ultra-violet
(UV) observations can constrain the physics of the r-mode instability. We show
that the core temperatures inferred from the data would place many systems well
inside the unstable region predicted by standard physical models. However, this
is at odds with theoretical expectations. We discuss different mechanisms that
could be at work in the stellar interior, and we show how they can modify the
instability window and make it consistent with the inferred temperatures.
[40]
oai:arXiv.org:1109.0536 [pdf] - 1083790
Gravitational Waves and the Maximum Spin Frequency of Neutron Stars
Submitted: 2011-09-02
In this Letter we re-examine the idea that gravitational waves are required
as a braking mechanism to explain the observed maximum spin-frequency of
neutron stars. We show that for millisecond X-ray pulsars, the existence of
spin equilibrium as set by the disk/magnetosphere interaction is sufficient to
explain the observations. We show as well that no clear correlation exists
between the neutron star magnetic field B and the X-ray outburst luminosity Lx
when considering an enlarged sample size of millisecond X-ray pulsars.
[41]
oai:arXiv.org:1107.5064 [pdf] - 403995
Revealing the physics of r-modes in low-mass X-ray binaries
Submitted: 2011-07-25, last modified: 2011-08-30
We consider the astrophysical constraints on the gravitational-wave driven
r-mode instability in accreting neutron stars in low-mass X-ray binaries. We
use recent results on superfluid and superconducting properties to infer the
core temperature in these neutron stars and show the diversity of the observed
population. Simple theoretical models indicate that many of these systems
reside inside the r-mode instability region. However, this is in clear
disagreement with expectations, especially for the systems containing the most
rapidly rotating neutron stars. The inconsistency highlights the need to
re-evaluate our understanding of the many areas of physics relevant to the
r-mode instability. We summarize the current status of our understanding, and
we discuss directions for future research which could resolve this dilemma.
[42]
oai:arXiv.org:1102.1548 [pdf] - 698656
Implications of burst oscillations from the slowly rotating accreting
pulsar IGR 17480-2446 in the globular cluster Terzan 5
Submitted: 2011-02-08, last modified: 2011-08-22
The recently-discovered accreting X-ray pulsar IGR J17480--2446 spins at a
frequency of ~11 Hz. We show that Type I X-ray bursts from this source display
oscillations at the same frequency as the stellar spin. IGR J17480--2446 is the
first secure case of a slowly rotating neutron star which shows Type I burst
oscillations, all other sources featuring such oscillations spin at hundreds of
Hertz. This means that we can test burst oscillation models in a completely
different regime. We explore the origin of Type I burst oscillations in IGR
J17480--2446 and conclude that they are not caused by global modes in the
neutron star ocean. We also show that the Coriolis force is not able to confine
an oscillation-producing hot-spot on the stellar surface. The most likely
scenario is that the burst oscillations are produced by a hot-spot confined by
hydromagnetic stresses.
[43]
oai:arXiv.org:1107.5295 [pdf] - 1078234
Modelling pulsar glitches with realistic pinning forces: a
hydrodynamical approach
Submitted: 2011-07-26
Although pulsars are one of the most stable clocks in the universe, many of
them are observed to 'glitch', i.e. to suddenly increase their spin frequency
(\nu) with fractional increases that range from \Delta\nu/\nu \approx 10^{-11}
to 10^{-5}. In this paper we focus on the 'giant' glitches, i.e. glitches with
fractional increases in the spin rate of the order of \Delta\nu/{\nu} \approx
10^{-6}, that are observed in a sub class of pulsars including the Vela. We
show that giant glitches can be modelled with a two-fluid hydrodynamical
approach. The model is based on the formalism for superfluid neutron stars of
Andersson and Comer (2006) and on the realistic pinning forces of Grill and
Pizzochero (2011). We show that all stages of Vela glitches, from the rise to
the post-glitch relaxation, can be reproduced with a set of physically
reasonable parameters and that the sizes and waiting times between giant
glitches in other pulsars are also consistent with our model.
[44]
oai:arXiv.org:1106.6264 [pdf] - 1077676
Spin equilibrium with or without gravitational wave emission: the case
of XTE J1814-338 and SAX J1808.4-3658
Submitted: 2011-06-30
In this letter we present a new analysis of the torques acting on the
accreting millisecond X-ray pulsars SAX J1808.4-365 and XTE J1814-338, and show
how our results can be used to constrain theoretical models of the spin
evolution. In particular we find upper limits on any spin-up/down phase of XTE
J1814-338 of $|\dot{\nu}|\lesssim\,1.5\times\,10^{-14}$ Hz/s at 95% confidence
level. We examine the possibility that a gravitational wave torque may be
acting in these systems and suggest that a more likely scenario is that both
systems are close to spin equilibrium, as set by the disc/magnetosphere
interaction.
[45]
oai:arXiv.org:1105.1244 [pdf] - 1076453
Lagrangian perturbation theory for a superfluid immersed in an elastic
neutron star crust
Submitted: 2011-05-06
The inner crust of mature neutron stars, where an elastic lattice of
neutron-rich nuclei coexists with a neutron superfluid, impacts on a range of
astrophysical phenomena. The presence of the superfluid is key to our
understanding of pulsar glitches, and is expected to affect the thermal
conductivity and hence the evolution of the surface temperature. The coupling
between crust and superfluid must also be accounted for in studies of neutron
star dynamics, discussions of global oscillations and associated instabilities.
In this paper we develop Lagrangian perturbation theory for this problem,
paying attention to key issues like superfluid entrainment, potential vortex
pinning, dissipative mutual friction and the star's magnetic field. We also
discuss the nature of the core-crust interface. The results provide a
theoretical foundation for a range of interesting astrophysical applications.
[46]
oai:arXiv.org:1011.1180 [pdf] - 320775
Tkachenko modes in rotating neutron stars: the effect of compressibility
and implications for pulsar timing noise
Submitted: 2010-11-04
Long wavelength oscillations (Tkachenko waves) of the triangular lattice of
quantized vortices in superfluid neutron stars have been suggested as one of
the possible explanations for the timing noise observed in many radio pulsars,
in particular for the 100-1000 day variations in the spin of PSR B1828-11. Most
studies to date have, however, been based on the hydrodynamics developed for
superfluid Helium. In this paper we extend the formulation to a two fluid
neutron and proton system, relevant for neutron star interiors and include the
effect of chemical coupling, compressibility and mutual friction between the
components. In particular we find that chemical coupling and compressibility
can have a drastic effect on the mode structure. However, for the slower
pulsars rotating at 1-10 Hz (such as PSR 1828-11), most choices of parameters
in the equation of state lead to Tkachenko oscillations with frequencies in the
correct range to explain the timing noise. We also investigate the case of more
rapidly rotating pulsars (above 100 Hz) for which we find that there is a vast
portion of parameter space in which there are no Tkachenko modes, but only
modified sound waves at much higher frequencies.
[47]
oai:arXiv.org:1005.1163 [pdf] - 903134
r-modes in low temperature colour-flavour-locked superconducting quark
star
Submitted: 2010-05-07
We present the first multi-fluid analysis of a dense neutron star core with a
deconfined colour-flavour-locked superconducting quark component. Accounting
only for the condensate and (finite temperature) phonons, we make progress by
taking over results for superfluid $^4$He. The resultant two-fluid model
accounts for a number of additional viscosity coefficients (compared to the
Navier-Stokes equations) and we show how they enter the dissipation analysis
for an oscillating star. We provide simple estimates for the gravitational-wave
driven r-mode instability, demonstrating that the various phonon processes that
we consider are not effective damping agents. Even though the results are
likely of little direct astrophysical importance (since we consider an overly
simplistic stellar model) our analysis represents significant technical
progress, laying the foundation for more detailed numerical studies and
preparing the ground for the inclusion of additional aspects (in particular
associated with kaons) of the problem.
[48]
oai:arXiv.org:1003.5849 [pdf] - 1026024
Superfluid hyperon bulk viscosity and the r-mode instability of rotating
neutron stars
Submitted: 2010-03-30
In order to establish whether the unstable r-modes in a rotating neutron star
provide a detectable source of gravitational waves, we need to understand the
details of the many dissipative processes that tend to counteract the
instability. It has been established that the bulk viscosity due to exotic
particles, like hyperons, may be particularly important in this respect.
However, the effects of hyperon superfluidity have so far not been fully
accounted for. While the associated suppression of the reaction rates that give
rise to the bulk viscosity has been estimated, superfluid aspects of the fluid
dynamics have not been considered. In this paper we determine the r-mode
instability window for a neutron star with a $\Sigma^{-}$ hyperon core, using
the appropriate multifluid formalism including, for the first time, the effect
of the "superfluid" bulk viscosity coefficients. We demonstrate that, even
though the extra terms may increase the bulk viscosity damping somewhat, their
presence does not affect the qualitative features of the r-mode instability
window.
[49]
oai:arXiv.org:0812.3569 [pdf] - 19660
Oscillations of Rapidly Rotating Superfluid Stars
Submitted: 2008-12-18, last modified: 2009-03-26
Using time evolutions of the relevant linearised equations we study
non-axisymmetric oscillations of rapidly rotating and superfluid neutron stars.
We consider perturbations of Newtonian axisymmetric background configurations
and account for the presence of superfluid components via the standard
two-fluid model. Within the Cowling approximation, we are able to carry out
evolutions for uniformly rotating stars up to the mass-shedding limit. This
leads to the first detailed analysis of superfluid neutron star oscillations in
the fast rotation regime, where the star is significantly deformed by the
centrifugal force. For simplicity, we focus on background models where the two
fluids (superfluid neutrons and protons) co-rotate, are in beta-equilibrium and
coexist throughout the volume of the star. We construct sequences of rotating
stars for two analytical model equations of state. These models represent
relatively simple generalisations of single fluid, polytropic stars. We study
the effects of entrainment, rotation and symmetry energy on non-radial
oscillations of these models. Our results show that entrainment and symmetry
energy can have a significant effect on the rotational splitting of
non-axisymmetric modes. In particular, the symmetry energy modifies the
inertial mode frequencies considerably in the regime of fast rotation.
[50]
oai:arXiv.org:0807.3457 [pdf] - 14738
Oscillations of rapidly rotating stratified neutron stars
Submitted: 2008-07-22, last modified: 2009-03-26
We use time-evolutions of the linear perturbation equations to study the
oscillations of rapidly rotating neutrons stars. Our models account for the
buoyancy due to composition gradients and we study, for the first time, the
nature of the resultant g-modes in a fast spinning star. We provide detailed
comparisons of non-stratified and stratified models. This leads to an improved
understanding of the relationship between the inertial modes of a
non-stratified star and the g-modes of a stratified system. In particular, we
demonstrate that each g-mode becomes rotation-dominated, i.e. approaches a
particular inertial mode, as the rotation rate of the star is increased. We
also discuss issues relating to the gravitational-wave driven instability of
the various classes of oscillation modes.
[51]
oai:arXiv.org:0902.1149 [pdf] - 1001430
r-modes and mutual friction in rapidly rotating superfluid neutron stars
Submitted: 2009-02-06
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. Firstly, 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 focussing on a
particular mechanism. Thus, we consider not only the weak drag regime, but also
the strong drag regime where the fluid dynamics is significantly different.
[52]
oai:arXiv.org:0812.3023 [pdf] - 19565
On the oscillations of dissipative superfluid neutron stars
Submitted: 2008-12-16
We investigate the oscillations of slowly rotating superfluid stars, taking
into account the vortex mediated mutual friction force that is expected to be
the main damping mechanism in mature neutron star cores. Working to linear
order in the rotation of the star, we consider both the fundamental f-modes and
the inertial r-modes. In the case of the (polar) f-modes, we work out an
analytic approximation of the mode which allows us to write down a closed
expression for the mutual friction damping timescale. The analytic result is in
good agreement with previous numerical results obtained using an energy
integral argument. We extend previous work by considering the full range of
permissible values for the vortex drag, e.g. the friction between each
individual vortex and the electron fluid. This leads to the first ever results
for the f-mode in the strong drag regime. Our estimates provide useful insight
into the dependence on, and relevance of, various equation of state parameters.
In the case of the (axial) r-modes, we confirm the existence of two classes of
modes. However, we demonstrate that only one of these sets remains purely axial
in more realistic neutron star models. Our analysis lays the foundation for
companion studies of the mutual friction damping of the r-modes at second order
in the slow-rotation approximation, the first time evolutions for superfluid
neutron star perturbations and also the first detailed attempt at studying the
dynamics of superfluid neutron stars with both a relative rotation between the
components and mutual friction.
[53]
oai:arXiv.org:0708.2984 [pdf] - 4190
Are neutron stars with crystalline colour superconducting cores
interesting for the LIGO experiment?
Submitted: 2007-08-22, last modified: 2008-05-16
We estimate the maximal deformation that can be sustained by a rotating
neutron star with a crystalline colour superconducting quark core. Our results
suggest that current gravitational-wave data from LIGO have already reached the
level where a detection would have been possible over a wide range of the
poorly constrained QCD parameters. This leads to the non-trivial conclusion
that compact objects \emph{do not} contain maximally strained colour
crystalline cores drawn from this range of parameter space. We discuss the
uncertainties associated with our simple model and how it can be improved in
the future.
[54]
oai:arXiv.org:0705.1780 [pdf] - 1216
Modelling magnetically deformed neutron stars
Submitted: 2007-05-12
Rotating deformed neutron stars are important potential sources for
groundbased gravitational-wave interferometers such as LIGO, GE0600 and VIRGO.
One mechanism that may lead to significant non-asymmetries is the internal
magnetic field. It is well known that a magnetic star will not be spherical
and, if the magnetic axis is not aligned with the spin axis, the deformation
will lead to the emission of gravitational waves. The aim of this paper is to
develop a formalism that would allow us to model magnetically deformed stars,
using both realistic equations of state and field configurations. As a first
step, we consider a set of simplified model problems. Focusing on dipolar
fields, we determine the internal magnetic field which is consistent with a
given neutron star model. We then calculate the associated deformation. We
conclude by discussing the relevance of our results for current
gravitational-wave detectors and future prospects.
[55]
oai:arXiv.org:astro-ph/0609438 [pdf] - 85030
Mountains on Neutron Stars: Accreted vs. Non-Accreted crusts
Submitted: 2006-09-15
The aim of this paper is to compare the two cases of an isolated neutron
star, with a non-accreted crust, and that of an accreting neutron star, with an
accreted crust, and try to estimate which one of the two would make a better
source of gravitational waves. In order to do this we must evaluate the maximum
``mountain'' that the crust can sustain in these two cases. We first do this
using the formalism of Ushomirsky, Cutler and Bildsten (2000) and find that the
maximum quadrupole is very similar in the two cases, with the non-accreted
crust sustaining a slightly larger mountain. We then develop a perturbation
formalism for the problem, that allows us to drop the Cowling approximation and
have more control over the boundaries. The use of this formalism confirms that
there is not much difference between the two cases, but leads to results
approximately one order of magnitude larger than those we obtain with the
formalism of Ushomirsky, Cutler and Bildsten (2000).
[56]
oai:arXiv.org:astro-ph/0411747 [pdf] - 69346
Modelling the spin equilibrium of neutron stars in LMXBs without
gravitational radiation
Submitted: 2004-11-28, last modified: 2005-08-18
In this paper we discuss the spin-equilibrium of accreting neutron stars in
LMXBs. We demonstrate that, when combined with a naive spin-up torque, the
observed data leads to inferred magnetic fields which are at variance with
those of galactic millisecond radiopulsars. This indicates the need for either
additional spin-down torques (eg. gravitational radiation) or an improved
accretion model. We show that a simple consistent accretion model can be
arrived at by accounting for radiation pressure in rapidly accreting systems
(above a few percent of the Eddington accretion rate). In our model the inner
disk region is thick and significantly sub-Keplerian, and the estimated
equilibrium periods are such that the LMXB neutron stars have properties that
accord well with the galactic millisecond radiopulsar sample. The implications
for future gravitational-wave observations are also discussed briefly.
