BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Brynmor Haskell

Brynmor Haskell contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.HEgr-qcphysics.flu-dynastro-ph.SRastro-ph.IMcond-mat.quant-gashep-phhep-thnucl-th

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

12 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