BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Brynmor Haskell

Brynmor Haskell appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEastro-ph.SRgr-qchep-phnucl-thphysics.flu-dynastro-ph.IMcond-mat.quant-gashep-th

Catalog footprint

What is connected

34works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

34 published item(s)

preprint2024arXiv

Glitching pulsars as gravitational wave sources

Spinning neutron stars, when observed as pulsars, are seen to undergo occasional spin-up events known as glitches. Despite several decades of study, the physical mechanisms responsible for glitches are still not well understood, but probably involve an interplay between the star's outer elastic crust, and the superfluid and superconducting core that lies within. Glitches will be accompanied by some level of gravitational wave emission. In this article we review proposed models that link gravitational wave emission to glitches, exploring both short duration burst-like emission, and longer-lived signals. We illustrate how detections (and in some cases, non-detections) of gravitational signals probe both the glitch mechanism, and, by extension, the behaviour of matter at high densities.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Extending Israel and Stewart hydrodynamics to relativistic superfluids via Carter's multifluid approach

We construct a relativistic model for bulk viscosity and heat conduction in a superfluid. Building on the principles of Unified Extended Irreversible Thermodynamics, the model is derived from Carter's multifluid approach for a theory with 3 four-currents: particles, entropy, and quasi-particles. Dissipation arises directly from the fact that the quasi-particle four-current is an independent degree of freedom that does not necessarily comove with the entropy. For small deviations from local thermodynamic equilibrium, the model provides an extension of the Israel-Stewart theory to superfluid systems. It can, therefore, be made hyperbolic, causal and stable if the microscopic input is accurate. The non-dissipative limit of the model is the relativistic two-fluid model of Carter, Khalatnikov and Gusakov. The Newtonian limit of the model is an Extended-Irreversible-Thermodynamic extension of Landau's two-fluid model. The model predicts the existence of four bulk viscosity coefficients and accounts for their microscopic origin, providing their exact formulas in terms of the quasi-particle creation rate. Furthermore, when fast oscillations of small amplitude around the equilibrium are considered, the relaxation-time term in the telegraph-type equations for the bulk viscosities accounts directly for their expected dependence on the frequency.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Long-term GRMHD simulation of magnetic field in isolated neutron stars

Strong magnetic fields play an important role in powering the emission of neutron stars. Nevertheless a full understanding of the interior configuration of the field remains elusive. In this work, we present General Relativistic MagnetoHydroDynamics simulations of the magnetic field evolution in neutron stars lasting 500 ms (5 Alfven crossing times) and up to resolutions of 0.231 km using Athena++. We explore two different initial conditions, one with purely poloidal magnetic field and the other with a dominant toroidal component, and study the poloidal and toroidal field energies, the growth times of the various instability-driven oscillation modes and turbulence. We find that the purely poloidal setup generates a toroidal field which later decays exponentially reaching 1% of the total magnetic energy, showing no evidence of reaching equilibrium. The initially stronger toroidal field setup, on the other hand, loses up to 20% of toroidal energy and maintains this state till the end of our simulation. We also explore the hypothesis, drawn from previous MHD simulations, that turbulence plays an important role in the quasi equilibrium state. An analysis of the spectra in our higher resolution setups reveal, however, that in most cases we are not observing turbulence at small scales, but rather a noisy velocity field inside the star. We also observe that the majority of the magnetic energy gets dissipated as heat increasing the internal energy of the star, while a small fraction gets radiated away as electromagnetic radiation.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Superfluid dynamics in neutron star crusts: the Iordanskii force and chemical gauge covariance

We present a geometrical derivation of the relativistic dynamics of the superfluid inner crust of a neutron star. The resulting model is analogous to the Hall-Vinen-Bekarevich-Khalatnikov hydrodynamics for a single-component superfluid at finite temperature, but particular attention should be paid to the fact that some fraction of the neutrons are locked to the motion of the protons in nuclei. This gives rise to an ambiguity in the definition of the two currents (the normal and the superfluid one) on which the model is built, a problem that manifests itself as a chemical gauge freedom of the theory. To ensure chemical gauge covariance of the hydrodynamic model, the phenomenological equation of motion for a quantized vortex should contain an extra transverse force, that is the relativistic version of the Iordanskii force discussed in the context of superfluid Helium. Hence, we extend the mutual friction model of Langlois et al. (1998) to account for the possible presence of this Iordanskii-like force. Furthermore, we propose that a better understanding of the (still not completely settled) controversy around the presence of the Iordanskii force in superfluid Helium, as well as in neutron stars, may be achieved by considering that the different incompatible results present in the literature pertain to two, opposite, dynamical regimes of the fluid system.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The impact of superconductivity and the Hall effect in models of magnetized neutron stars

Equilibrium configurations of the internal magnetic field of a pulsar play a key role in modeling astrophysical phenomena, from glitches to gravitational wave emission. In this paper we present a numerical scheme for solving the Grad-Shafranov equation and calculating equilibrium configurations of pulsars, accounting for superconductivity in the core of the neutron star, and for the Hall effect in the crust of the star. Our numerical code uses a finite-difference method in which the source term appearing in the Grad-Shafranov equation, used to model the magnetic equilibrium is nonlinear. We obtain solutions by linearizing the source and applying an under-relaxation scheme at each step of computation to improve the solver's convergence. We have developed our code in both C++ and Python, and our numerical algorithm can further be adapted to solve any nonlinear PDEs appearing in other areas of computational astrophysics. We produce mixed toroidal-poloidal field configurations, and extend the portion of parameter space that can be investigated with respect to previous studies. We find that even in the more extreme cases the magnetic energy in the toroidal component does not exceed approximately 5\% of the total. We also find that if the core of the star is superconducting, the toroidal component is entirely confined to the crust of the star, which has important implications for pulsar glitch models which rely on the presence of a strong toroidal field region in the core of the star, where superfluid vortices pin to superconducting fluxtubes.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Thermodynamic stability implies causality

The stability conditions of a relativistic hydrodynamic theory can be derived directly from the requirement that the entropy should be maximised in equilibrium. Here we use a simple geometrical argument to prove that, if the hydrodynamic theory is stable according to this entropic criterion, then localised perturbations to the equilibrium state cannot propagate outside their future light-cone. In other words, within relativistic hydrodynamics, acausal theories must be thermodynamically unstable, at least close to equilibrium. We show that the physical origin of this deep connection between stability and causality lies in the relationship between entropy and information. Our result may be interpreted as an ``equilibrium conservation theorem'', which generalizes the Hawking-Ellis vacuum conservation theorem to finite temperature and chemical potential.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A universal formula for the relativistic correction to the mutual friction coupling time-scale in neutron stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Asymmetric accretion and thermal `mountains' in magnetized neutron star crusts

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Magnetic field configurations in neutron stars from MHD simulations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The effect of non-linear mutual friction on pulsar glitch sizes and rise times

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Turbulent, pinned superfluids in neutron stars and pulsar glitch recoveries

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

When the entropy has no maximum: A new perspective on the instability of the first-order theories of dissipation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The effect of superfluid hydrodynamics on pulsar glitch sizes and waiting times

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Detecting gravitational waves from mountains on neutron stars in the Advanced Detector Era

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Mesoscopic pinning forces in neutron star crusts

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Models of Pulsar Glitches

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Observational constraints on neutron star crust-core coupling during glitches

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Pinned vortex hopping in a neutron star crust

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

R-modes in neutron stars: Theory and observations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gravitational waves from rapidly rotating neutron stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nuclear Equation of State from Observations of Short Gamma-Ray Burst Remnants

The favoured progenitor model for short $γ$-ray bursts (SGRBs) is the merger of two neutron stars that triggers an explosion with a burst of collimated $γ$-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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the universality of I-Love-Q relations in magnetized neutron stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A new mechanism for saturating unstable r-modes in neutron stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Glitch recoveries in radio-pulsars and magnetars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The hottest superfluid and superconductor in the Universe: Discovery and nuclear physics implications

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Constraining the physics of the r-mode instability in neutron stars with X-ray and UV observations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Investigating superconductivity in neutron star interiors with glitch models

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The effect of realistic equations of state and general relativity on the "snowplow" model for pulsar glitches

Many pulsars are observed to "glitch", i.e. show sudden jumps in their rotational frequency $ν$, some of which can be as large as $Δν/ν\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$).

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Gravitational Waves and the Maximum Spin Frequency of Neutron Stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Modelling pulsar glitches with realistic pinning forces: a hydrodynamical approach

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 (ν) with fractional increases that range from Δν/ν\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 Δν/ν \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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Revealing the physics of r-modes in low-mass X-ray binaries

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Spin equilibrium with or without gravitational wave emission: the case of XTE J1814-338 and SAX J1808.4-3658

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ν|\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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Superfluid hyperon bulk viscosity and the r-mode instability of rotating neutron stars

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 $Σ^{-}$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Tkachenko modes in rotating neutron stars: the effect of compressibility and implications for pulsar timing noise

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.

0 citations0 reviewsNo code linkedOpen access