Source author record

Nils Andersson

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

37works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

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

37 published item(s)

preprint2026arXiv

ESO Expanding Horizon White Paper: Revealing the properties of matter at supranuclear densities with gravitational waves

Understanding dense matter under extreme conditions is one of the most fundamental puzzles in modern physics. Complex interactions give rise to emergent, collective phenomena. While nuclear experiments and Earth - based colliders provide valuable insights, much of the quantum chromodynamics phase diagram at high density and low temperature remains accessible only through astrophysical observations of neutron stars, neutron star mergers, and stellar collapse. Astronomical observations thus offer a direct window to the physics on subatomic scales with gravitational waves presenting an especially clean channel. Next-generation gravitational - wave observatories, such as the Einstein Telescope, would serve as unparalleled instruments to transform our understanding of neutron star matter. They will enable the detection of up to tens of thousands of binary neutron star and neutron star - black hole mergers per year, a dramatic increase over the few events accessible with current detectors. They will provide an unprecedented precision in probing cold, dense matter during the binary inspiral, exceeding by at least an order of magnitude what current facilities can achieve. Moreover, these observatories will allow us to explore uncharted regimes of dense matter at finite temperatures produced in a subset of neutron star mergers, areas that remain entirely inaccessible to current instruments. Together with multimessenger observations, these measurements will significantly deepen our knowledge of dense nuclear matter.

preprint2024arXiv

Rotating Love: The dynamical tides of spinning Newtonian stars

We carefully develop the framework required to model the dynamical tidal response of a spinning neutron star in an inspiralling binary system, in the context of Newtonian gravity, making sure to include all relevant details and connections to the existing literature. The tidal perturbation is decomposed in terms of the normal oscillation modes, used to derive an expression for the effective Love number which is valid for any rotation rate. In contrast to previous work on the problem, our analysis highlights subtle issues relating to the orthogonality condition required for the mode-sum representation of the dynamical tide and shows how the prograde and retrograde modes combine to provide the overall tidal response. Utilising a slow-rotation expansion, we show that the dynamical tide (the effective Love number) is corrected at first order in rotation, whereas in the case of the static tide (the static Love number) the rotational corrections do not enter until second order.

preprint2022arXiv

Detecting the impact of nuclear reactions on neutron star mergers through gravitational waves

Nuclear reactions may affect gravitational-wave signals from neutron-star mergers, but the impact is uncertain. In order to quantify the effect, we compare two numerical simulations representing intuitive extremes. In one case reactions happen instantaneously. In the other case, they occur on timescales much slower than the evolutionary timescale. We show that, while the differences in the two gravitational-wave signals are small, they should be detectable by third-generation observatories. To avoid systematic errors in equation of state parameters inferred from observed signals, we need to accurately implement nuclear reactions in future simulations.

preprint2022arXiv

Dynamical tides in superfluid neutron stars

We study the tidal response of a superfluid neutron star in a binary system, focussing on Newtonian models with superfluid neutrons present throughout the star's core and the inner crust. Within the two-fluid formalism, we consider the main aspects that arise from the presence of different regions inside the star, with particular focus on the various interfaces. Having established the relevant theory, we determine the tidal excitation of the most relevant oscillation modes during binary inspiral. Our results suggest that superfluid physics has a negligible impact on the static tidal deformation. The overwhelming contribution to the Love number is given by, as for normal matter stars, the ordinary fundamental mode (f-mode). Strong entrainment, here described by a phenomenological expression which mimics the large effective neutron mass expected at the bottom of the crust, is shown to have significant impact on the superfluid modes, but our results for the dynamical tide are nevertheless similar to the static limit: the fundamental modes are the ones most significantly excited by the tidal interaction, with the ordinary f-mode dominating the superfluid one. We also discuss the strain built up in the star's crust during binary inspiral, showing that the superfluid f-mode may (depending on entrainment) reach the limit where the crust breaks, although it does so after the ordinary f-mode. Overall, our results suggest that the presence of superfluidity may be difficult to establish from binary neutron star gravitational-wave signals.

preprint2021arXiv

Linearizing a Non-linear Formulation for General Relativistic Dissipative Fluids

Fully non-linear equations of motion for dissipative general relativistic multi-fluids can be obtained from an action principle involving the explicit use of lower dimensional matter spaces. More traditional strategies for incorporating dissipation-like the famous Mueller-Israel-Stewart model-are based on expansions away from equilibrium defined, in part, by the laws of thermodynamics. The goal here is to build a formalism to facilitate comparison of the action-based results with those based on the traditional approach. The first step of the process is to use the action-based approach itself to construct self-consistent notions of equilibrium. Next, first-order deviations are developed directly on the matter spaces, which motivates the latter as the natural arena for the underlying thermodynamics. Finally, we identify the dissipation terms of the action-based model with first-order "thermodynamical" fluxes, on which the traditional models are built. The description is developed in a general setting so that the formalism can be used to describe multi-fluid systems, for which causal and stable models are not yet available. As an illustration of the approach, a simple application of a single viscous fluid is considered and, even though the expansion is halted at first order, we sketch how a causal response can be implemented through Cattaneo-type equations.

preprint2020arXiv

Early neutron star evolution in high-mass X-ray binaries

The application of standard accretion theory to observations of X-ray binaries provides valuable insights into neutron star properties, such as their spin period and magnetic field. However, most studies concentrate on relatively old systems, where the neutron star is in its late propeller, accretor, or nearly spin equilibrium phase. Here we use an analytic model from standard accretion theory to illustrate the evolution of high-mass X-ray binaries early in their life. We show that a young neutron star is unlikely to be an accretor because of the long duration of ejector and propeller phases. We apply the model to the recently discovered ~4000 yr old high-mass X-ray binary XMMU J051342.6-672412 and find that the system's neutron star, with a tentative spin period of 4.4 s, cannot be in the accretor phase and has a magnetic field B > (a few)x10^13 G, which is comparable to the magnetic field of many older high-mass X-ray binaries and is much higher than the spin equilibrium inferred value of (a few)x10^11 G. The observed X-ray luminosity could be the result of thermal emission from a young cooling magnetic neutron star or a small amount of accretion that can occur in the propeller phase.

preprint2020arXiv

Gravitational waves from transient neutron star f-mode oscillations

During their most recent observing run, LIGO/Virgo reported the gravitational wave (GW) transient S191110af, a burst signal at a frequency of 1.78 kHz that lasted for 0.104 s. While this signal was later deemed non-astrophysical, genuine detections of uncertain origin will occur in the future. Here we study the potential for detecting GWs from neutron star fluid oscillations, which have mode frequency and duration matching those of S191110af and which can be used to constrain the equation of state of nuclear matter. Assuming that such transient oscillations can be excited to energies typical of a pulsar glitch, we use measured properties of known glitching pulsars to estimate the amplitude of GWs produced by such events. We find that current GW detectors may observe nearby pulsars undergoing large events with energy similar to Vela pulsar glitch energies, while next generation detectors could observe a significant number of events. Finally, we show that it is possible to distinguish between GWs produced by rapidly rotating and slowly rotating pulsars from the imprint of rotation on the f-mode frequency.

preprint2020arXiv

Instabilities in neutron-star postmerger remnants

Using nonlinear, fully relativistic, simulations we investigate the dynamics and gravitational-wave signature associated with instabilities in neutron star postmerger remnants. For simplified models of the remnant we establish the presence of an instability in stars with moderate $T/|W|$, the ratio between the kinetic and the gravitational potential energies. Detailed analysis of the density oscillation pattern reveals a \emph{local} instability in the inner region of the more realistic differential rotation profile. We apply Rayleigh's inflection theorem and Fjørtoft's theorem to analyze the stability criteria concluding that this \emph{inner local} instability originates from a shear instability close to the peak of the angular velocity profile, and that it later evolves into a fast-rotating m=2 oscillation pattern. We discuss the importance of the presence of a corotation point in the fluid, its connection with the shear instability, and comparisons to the Rossby wave and Papaloizou Pringle instabilities considered in the wider literature.

preprint2020arXiv

Tidal deformations of hybrid stars with sharp phase transitions and elastic crusts

Gravitational wave astronomy is expected to provide independent constraints on neutron star properties, such as their dense matter equation of state. This is possible with the measurements of binary components' tidal deformability, which alter the point-particle gravitational waveforms of the late inspiral phase of neutron-star binaries. Although current gravitational wave detectors are not sensitive enough for a precise determination of the individual tidal deformations of the components, a large number of combined observations with future detectors will decrease uncertainties in this quantity. Here we provide a first study of the tidal deformability effects due to the elasticity/solidity of the crust (hadronic phase) in a hybrid neutron star, as well as the influence of a quark-hadronic phase density jump on tidal deformations. We employ the framework of nonradial perturbations with zero frequency and study hadronic phases presenting elastic aspects when perturbed (with the shear modulus approximately $1\%$ of the pressure). We find that the relative tidal deformation change in a hybrid star with a perfect-fluid quark phase and a hadronic phase presenting an elastic part is never larger than about $2-4\%$ (with respect to a perfect-fluid counterpart). These maximum changes occur when the elastic region of a hybrid star is larger than approximately $60\%$ of the star's radius, which may happen when its quark phase is small and the density jump is large enough, or even when a hybrid star has an elastic mixed phase. For other cases, the relative tidal deformation changes due to an elastic crust are negligible ($10^{-5}-10^{-1}\%$), therefore unlikely to be measured even with third generation detectors. Thus, only when the size of the elastic hadronic region of a hybrid star is over half of its radius, the effects of elasticity could have a noticeable impact on tidal deformations.

preprint2020arXiv

Tidal deformations of neutron stars with elastic crusts

With the first detections of binary neutron star mergers by gravitational-wave detectors, it proves timely to consider how the internal structure of neutron stars affects the way in which they can be asymmetrically deformed. Such deformations may leave measurable imprints on gravitational-wave signals and can be sourced through tidal interactions or the formation of mountains. We detail the formalism that describes fully-relativistic neutron star models with elastic crusts undergoing static perturbations. This formalism primes the problem for studies into a variety of mechanisms that can deform a neutron star. We present results for a barotropic equation of state and a realistic model for the elastic crust, which enables us to compute relevant quantities such as the tidal deformability parameter. We find that the inclusion of an elastic crust provides a very small correction to the tidal deformability. The results allow us to demonstrate when and where the crust starts to fail during a binary inspiral and we find that the majority of the crust will remain intact up until merger.

preprint2016arXiv

Ejector and propeller spin-down: How might a superluminous supernova millisecond magnetar become the 6.67 hr pulsar in RCW103

The X-ray source 1E 161348-5055 in the supernova remnant RCW 103 recently exhibited X-ray activity typical of magnetars, i.e. neutron stars with magnetic fields > 10^14-10^15 G. However, 1E 161348-5055 has an observed period of 6.67 hr, in contrast to magnetars which have a spin period of seconds. Here we describe a simple model which can explain the spin evolution of 1E 161348-5055, as well as other magnetars, from an initial period of milliseconds that would be required for dynamo generation of magnetar-strength magnetic fields. We propose that the key difference between 1E 161348-5055 and other magnetars is the persistence of a remnant disk of small total mass. This disk caused 1E 161348-5055 to undergo ejector and propeller phases in its life, during which strong torques caused a rapid increase of its spin period. By matching its observed spin period and ~1-3 kyr age, we find that 1E 161348-5055 has the (slightly) highest magnetic field of all known magnetars, with B~5x10^15 G, and that its disk had a mass of ~10^24 g, comparable to that of the asteroid Ceres.

preprint2016arXiv

Measuring the neutron star equation of state using X-ray timing

One of the primary science goals of the next generation of hard X-ray timing instruments is to determine the equation of state of the matter at supranuclear densities inside neutron stars, by measuring the radius of neutron stars with different masses to accuracies of a few percent. Three main techniques can be used to achieve this goal. The first involves waveform modelling. The flux we observe from a hotspot on the neutron star surface offset from the rotational pole will be modulated by the star's rotation, giving rise to a pulsation. Information about mass and radius is encoded into the pulse profile via relativistic effects, and tight constraints on mass and radius can be obtained. The second technique involves characterising the spin distribution of accreting neutron stars. The most rapidly rotating stars provide a very clean constraint, since the mass-shedding limit is a function of mass and radius. However the overall spin distribution also provides a guide to the torque mechanisms in operation and the moment of inertia, both of which can depend sensitively on dense matter physics. The third technique is to search for quasi-periodic oscillations in X-ray flux associated with global seismic vibrations of magnetars (the most highly magnetized neutron stars), triggered by magnetic explosions. The vibrational frequencies depend on stellar parameters including the dense matter equation of state. We illustrate how these complementary X-ray timing techniques can be used to constrain the dense matter equation of state, and discuss the results that might be expected from a 10m$^2$ instrument. We also discuss how the results from such a facility would compare to other astronomical investigations of neutron star properties. [Modified for arXiv]

preprint2015arXiv

Magnetic Field Evolution in Superconducting Neutron Stars

The presence of superconducting and superfluid components in the core of mature neutron stars calls for the rethinking of a number of key magnetohydrodynamical notions like resistivity, the induction equation, magnetic energy and flux-freezing. Using a multi-fluid magnetohydrodynamics formalism, we investigate how the magnetic field evolution is modified when neutron star matter is composed of superfluid neutrons, type-II superconducting protons and relativistic electrons. As an application of this framework, we derive an induction equation where the resistive coupling originates from the mutual friction between the electrons and the vortex/fluxtube arrays of the neutron and proton condensates. The resulting induction equation allows the identification of two timescales that are significantly different from those of standard magnetohydrodynamics. The astrophysical implications of these results are briefly discussed.

preprint2015arXiv

Pinning down the superfluid and measuring masses using pulsar glitches

Pulsars are known for their superb timing precision, although glitches can interrupt the regular timing behavior when the stars are young. These glitches are thought to be caused by interactions between normal and superfluid matter in the crust of the star. However, glitching pulsars such as Vela have been shown to require a superfluid reservoir that greatly exceeds that available in the crust. We examine a model in which glitches tap the superfluid in the core. We test a variety of theoretical superfluid models against the most recent glitch data and find that only one model can successfully explain up to 45 years of observational data. We develop a new technique for combining radio and X-ray data to measure pulsar masses, thereby demonstrating how current and future telescopes can probe fundamental physics such as superfluidity near nuclear saturation.

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.

preprint2014arXiv

Probing the neutron star interior and the Equation of State of cold dense matter with the SKA

With an average density higher than the nuclear density, neutron stars (NS) provide a unique test-ground for nuclear physics, quantum chromodynamics (QCD), and nuclear superfluidity. Determination of the fundamental interactions that govern matter under such extreme conditions is one of the major unsolved problems of modern physics, and -- since it is impossible to replicate these conditions on Earth -- a major scientific motivation for SKA. The most stringent observational constraints come from measurements of NS bulk properties: each model for the microscopic behaviour of matter predicts a specific density-pressure relation (its `Equation of state', EOS). This generates a unique mass-radius relation which predicts a characteristic radius for a large range of masses and a maximum mass above which NS collapse to black holes. It also uniquely predicts other bulk quantities, like maximum spin frequency and moment of inertia. The SKA, in Phase 1 and particularly in Phase 2 will, thanks to the exquisite timing precision enabled by its raw sensitivity, and surveys that dramatically increase the number of sources: 1) Provide many more precise NS mass measurements (high mass NS measurements are particularly important for ruling out EOS models); 2) Allow the measurement of the NS moment of inertia in highly relativistic binaries such as the Double Pulsar; 3) Greatly increase the number of fast-spinning NS, with the potential discovery of spin frequencies above those allowed by some EOS models; 4) Improve our knowledge of new classes of binary pulsars such as black widows and redbacks (which may be massive as a class) through sensitive broad-band radio observations; and 5) Improve our understanding of dense matter superfluidity and the state of matter in the interior through the study of rotational glitches, provided that an ad-hoc campaign is developed.

preprint2014arXiv

The intimate relation between the low T/W instability and the co-rotation point

We study the low T/W instability associated with the f-mode of differentially rotating stars. Our stellar models are described by a polytropic equation of state and the rotation profile is given by the standard j-constant law. The properties of the relevant oscillation modes, including the instability growth time, are determined from time evolutions of the linearised dynamical equations in Newtonian gravity. In order to analyse the instability we monitor also the canonical energy and angular momentum. Our results demonstrate that the l=m=2 f-mode becomes unstable as soon as a co-rotation point develops inside the star (i.e. whenever there is a point where the mode's pattern speed matches the bulk angular velocity). Considering various degrees of differential rotation, we show that the instability grows faster deep inside the co-rotation region and deduce an empirical relation that correlates the mode frequency and the star's parameters, which captures the main features of the l=m=2 f-mode growth time. This function is proportional to the product of the kinetic to gravitational energy ratio and the gradient of the star's spin, strengthening further the relationship between the co-rotation point and the low T/W instability. We briefly consider also the l=m=2 r-mode and demonstrate that it never moves far inside the co-rotation region even for significant differential rotation.

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.

preprint2013arXiv

Equilibrium spin pulsars unite neutron star populations

Many pulsars are formed with a binary companion from which they can accrete matter. Torque exerted by accreting matter can cause the pulsar spin to increase or decrease, and over long times, an equilibrium spin rate is achieved. Application of accretion theory to these systems provides a probe of the pulsar magnetic field. We compare the large number of recent torque measurements of accreting pulsars with a high-mass companion to the standard model for how accretion affects the pulsar spin period. We find that many long spin period (P > 100 s) pulsars must possess either extremely weak (B < 10^10 G) or extremely strong (B > 10^14 G) magnetic fields. We argue that the strong-field solution is more compelling, in which case these pulsars are near spin equilibrium. Our results provide evidence for a fundamental link between pulsars with the slowest spin periods and strong magnetic fields around high-mass companions and pulsars with the fastest spin periods and weak fields around low-mass companions. The strong magnetic fields also connect our pulsars to magnetars and strong-field isolated radio/X-ray pulsars. The strong field and old age of our sources suggests their magnetic field penetrates into the superconducting core of the neutron star.

preprint2013arXiv

Rotational evolution of young pulsars due to superfluid decoupling

Pulsars are rotating neutron stars that are seen to slow down, and the spin-down rate is thought to be due to magnetic dipole radiation. This leads to a prediction for the braking index n, which is a combination of spin period and its first and second time derivatives. However, all observed values of n are below the predicted value of 3. Here we provide a simple model that can explain the rotational evolution of young pulsars, including the n=2.51 of the 958-year-old pulsar in the Crab nebula. The model is based on a decrease in effective moment of inertia due to an increase in the fraction of the stellar core that becomes superfluid as the star cools via neutrino emission. The results suggest that future large radio monitoring campaigns of pulsars will yield measurements of the neutron star mass, nuclear equation of state, and superfluid properties.

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.

preprint2013arXiv

The Transient Gravitational-Wave Sky

Interferometric detectors will very soon give us an unprecedented view of the gravitational-wave sky, and in particular of the explosive and transient Universe. Now is the time to challenge our theoretical understanding of short-duration gravitational-wave signatures from cataclysmic events, their connection to more traditional electromagnetic and particle astrophysics, and the data analysis techniques that will make the observations a reality. This paper summarizes the state of the art, future science opportunities, and current challenges in understanding gravitational-wave transients.

preprint2012arXiv

Hydromagnetic equilibrium in non-barotropic multifluid neutron stars

Traditionally, the subject of hydromagnetic equilibrium in neutron stars has been addressed in the context of standard magnetohydrodynamics, with matter obeying a barotropic equation of state. In this paper we take a step towards a more realistic treatment of the problem by considering neutron stars with interior superfluid components. In this multifluid model stratification associated with a varying matter composition (the relative proton to neutron density fraction) enters as a natural ingredient, leading to a non-barotropic system. After formulating the hydromagnetic equilibrium of superfluid/superconducting neutron stars as a perturbation problem, we focus on the particular case of a three-fluid system consisting of superfluid neutrons and normal protons and electrons. We determine the equilibrium structure of dipolar magnetic fields with a mixed poloidal-toroidal composition. We find that, with respect to barotropic models, stratification has the generic effect of leading to equilibria with a higher fraction of magnetic energy stored in the toroidal component. However, even in models with strong stratification the poloidal and toroidal components are comparable, with the former contributing the bulk of the magnetic energy.

preprint2012arXiv

Magnetars: super(ficially) hot and super(fluid) cool

We examine to what extent the inferred surface temperature of magnetars in quiescence can constrain the presence of a superfluid in the neutron star core and the role of magnetic field decay in the core. By performing detailed simulations of neutron star cooling, we show that extremely strong heating from field decay in the core cannot produce the high observed surface temperatures nor delay the onset of neutron superfluidity in the core. We verify the results of Kaminker et al., namely that the high magnetar surface temperatures require heating in the neutron star crust, and crust heating is decoupled from cooling/heating in the core. Therefore, because crust heating masks core heating, it is not possible to conclude that magnetar cores are in a non-superfluid state purely from high surface temperatures. From our interior temperature evolutions and after accounting for proton superconductivity in the core, we find that neutron superfluidity in the core occurs less than a few hundred years after neutron star formation. This onset time is unaffected by heating due to core field decay at fields < 10^16 G. Thus all known neutron stars, including magnetars, without a core containing exotic particles, should have a core of superfluid neutrons and superconducting protons.

preprint2011arXiv

Cosmic recycling of millisecond pulsars

We compare the rotation rate of neutron stars in low-mass X-ray binaries (LMXBs) with the orbital period of the binaries. We find that, while short orbital period LMXBs span a range of neutron star rotation rates, all the long period LMXBs have fast rotators. We also find that the rotation rates are highest for the systems with the highest mean mass accretion rates, as can be expected if the accretion rate correlates with the orbital period. We show that these properties can be understood by a balance between spin-up due to accretion and spin-down due to gravitational radiation. Our scenario indicates that the gravitational radiation emitted by these systems may be detectable by future ground-based gravitational wave detectors.

preprint2011arXiv

Magneto-rotational neutron star evolution: the role of core vortex pinning

We consider the pinning of superfluid (neutron) vortices to magnetic fluxtubes associated with a type II (proton) superconductor in neutron star cores. We demonstrate that core pinning affects the spin-down of the system significantly, and discuss implications for regular radio pulsars and magnetars. We find that magnetars are likely to be in the pinning regime, while most radio pulsars are not. This suggests that the currently inferred magnetic field for magnetars may be overestimated. We also obtain a new timescale for the magnetic field evolution which could be associated with the observed activity in magnetars, provided that the field has a strong toroidal component.

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.

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.

preprint2005arXiv

A Flux-Conservative Formalism for Convective and Dissipative Multi-Fluid Systems, with Application to Newtonian Superfluid Neutron Stars

We develop a flux-conservative formalism for a Newtonian multi-fluid system, including dissipation and entrainment (i.e. allowing the momentum of one fluid to be a linear combination of the velocities of all fluids). Maximum use is made of mass, energy, and linear and angular momentum conservation to specify the equations of motion. Also used extensively are insights gleaned from a convective variational action principle, key being the distinction between each velocity and its canonically conjugate momentum. Dissipation is incorporated to second order in the ``thermodynamic forces'' via the approach pioneered by Onsager. An immediate goal of the investigation is to understand better the number, and form, of independent dissipation terms required for a consistent set of equations of motion in the multi-fluid context. A significant, but seemingly innocuous detail, is that one must be careful to isolate ``forces'' that can be written as total gradients, otherwise errors can be made in relating the net internal force to the net externally applied force. Our long-range aim is to provide a formalism that can be used to model dynamical multi-fluid systems both perturbatively and via fully nonlinear 3D numerical evolutions. To elucidate the formalism we consider the standard model for a heat-conducting, superfluid neutron star, which is believed to be dominated by superfluid neutrons, superconducting protons, and a highly degenerate, ultra-relativistic gas of normal fluid electrons. We determine that in this case there are, in principle, 19 dissipation coefficients in the final set of equations.

preprint2004arXiv

Gravitational-wave astronomy: the high-frequency window

This contribution is divided in two parts. The first part provides a text-book level introduction to gravitational radiation. The key concepts required for a discussion of gravitational-wave physics are introduced. In particular, the quadrupole formula is applied to the anticipated ``bread-and-butter'' source for detectors like LIGO, GEO600, EGO and TAMA300: inspiralling compact binaries. The second part provides a brief review of high frequency gravitational waves. In the frequency range above (say) 100Hz, gravitational collapse, rotational instabilities and oscillations of the remnant compact objects are potentially important sources of gravitational waves. Significant and unique information concerning the various stages of collapse, the evolution of protoneutron stars and the details of the supranuclear equation of state of such objects can be drawn from careful study of the gravitational-wave signal. As the amount of exciting physics one may be able to study via the detections of gravitational waves from these sources is truly inspiring, there is strong motivation for the development of future generations of ground based detectors sensitive in the range from hundreds of Hz to several kHz.

preprint1999arXiv

Construction of Astrophysical Initial Data for Perturbations of Relativistic Stars

We develop a framework for constructing initial data sets for perturbations about spherically symmetric matter distributions. This framework facilitates setting initial data representing astrophysical sources of gravitational radiation involving relativistic stars. The procedure is based on Lichnerowicz-York's conformal approach to solve the constraints in Einstein's equations. The correspondence of these initial data sets in terms of the standard gauge perturbation variables in the Regge-Wheeler perturbation variables is established, and examples of initial data sets of merging neutron stars under the close-limit approximation are presented.

preprint1998arXiv

Dynamics of scalar fields in the background of rotating black holes II: A note on superradiance

We analyze the amplification due to so-called superradiance from the scattering of pulses off rotating black holes as a numerical time evolution problem. We consider the "worst possible case" of scalar field pulses for which superradiance effects yield amplifications $< 1%$. We show that this small effect can be isolated by numerically evolving quasi-monochromatic, modulated pulses with a recently developed Teukolsky code. The results show that it is possible to study superradiance in the time domain, but only if the initial data is carefully tuned. This illustrates the intrinsic difficulties of detecting superradiance in more general evolution scenarios.

preprint1997arXiv

Gravitational waves from pulsating stars: Evolving the perturbation equations for a relativistic star

We consider the perturbations of a relativistic star as an initial-value problem. Having discussed the formulation of the problem (the perturbation equations and the appropriate boundary conditions at the centre and the surface of the star) in detail we evolve the equations numerically from several different sets of initial data. In all the considered cases we find that the resulting gravitational waves carry the signature of several of the star's pulsation modes. Typically, the fluid $f$-mode, the first two $p$-modes and the slowest damped gravitational $w$-mode are present in the signal. This indicates that the pulsation modes may be an interesting source for detectable gravitational waves from colliding neutron stars or supernovae. We also survey the literature and find several indications of mode presence in numerical simulations of rotating core collapse and coalescing neutron stars. If such mode-signals can be detected by future gravitational-wave antennae one can hope to infer detailed information about neutron stars. Since a perturbation evolution should adequately describe the late time behaviour of a dynamically excited neutron star, the present work can also be used as a bench-mark test for future fully nonlinear simulations.

preprint1996arXiv

Evolving test-fields in a black-hole geometry

We consider the initial value problem for a massless scalar field in the Schwarzschild geometry. When constructed using a complex-frequency approach the necessary Green's function splits into three components. We discuss all of these in some detail: 1) The contribution from the singularities (the quasinormal modes of the black hole) is approximated and the mode-sum is demonstrated to converge after a certain well defined time in the evolution. A dynamic description of the mode-excitation is introduced and tested. 2) It is shown how a straightforward low-frequency approximation to the integral along the branch cut in the black-hole Green's function leads to the anticipated power-law fall off at very late times. We also calculate higher order corrections to this tail and show that they provide an important complement to the leading order. 3) The high-frequency problem is also considered. We demonstrate that the combination of the obtained approximations for the quasinormal modes and the power-law tail provide a complete description of the evolution at late times. Problems that arise (in the complex-frequency picture) for early times are also discussed, as is the fact that many of the presented results generalize to, for example, Kerr black holes.

preprint1996arXiv

Spacetime modes of relativistic stars

The problem of relativistic stellar pulsations is studied in a somewhat ad hoc approximation that ignores all fluid motions. This {\em Inverse Cowling Approximation} (ICA) is motivated by two observations: 1) For highly damped ($w$-mode) oscillations the fluid plays very little role. 2) If the fluid motion is neglected the problem for polar oscillation modes becomes similar to that for axial modes. Using the ICA we find a polar mode spectrum that has all features of the $w$-mode spectrum of the full problem. Moreover, in the limit of superdense stars, we find the ICA spectrum to be qualitatively similar to that of the axial modes. These results clearly show the importance of general relativity for the pulsation modes of compact stars, and that there are modes whose existence do not depend on motions of the fluid at all; pure ``spacetime'' modes.

preprint1995arXiv

A NEW NUMERICAL APPROACH TO THE OSCILLATION MODES OF RELATIVISTIC STARS

The oscillation modes of a simple polytropic stellar model are studied. Using a new numerical approach (based on integration for complex coordinates) to the problem for the stellar exterior we have computed the eigenfrequencies of the highly damped w-modes. The results obtained agree well with recent ones of Leins, Nollert and Soffel (1993) Specifically, we are able to explain why several modes in this regime of the complex frequency plane could not be identified within the WKB approach of Kokkotas and Schutz (1992). Furthermore, we have established that the ``kink'' that was a prominent feature of the spectra of Kokkotas and Schutz, but did not appear in the results of Leins {\em et al.}, was a numerical artefact. Using our new numerical code we are also able to compute, for the first time, several of the slowly damped (p) modes for the considered stellar models. For very compact stars we find, somewhat surprisingly, that the damping of these modes does not decrease monotonically as one proceeds to higher oscillation frequencies. The existence of low-order modes that damp away much faster than anticipated may have implications for questions regarding stellar stability and the lifetime of gravitational-wave sources. The present results illustrate the accuracy and reliability of the complex-coordinate method and indicate that the method could prove to be of great use also in problems involving rotating stars. There is no apparent reason why the complex-coordinate approach should not extend to rotating stars, whereas it is accepted that all previous methods will fail to do so.

preprint1995arXiv

ON THE OSCILLATION SPECTRA OF ULTRA COMPACT STARS

Quasinormal modes of ultra compact stars with uniform energy density have been calculated. For less compact stars, there is only one very slowly damped polar mode (corresponding to the Kelvin f-mode) for each spherical harmonic index $l$. Further long-lived modes become possible for a sufficiently compact star (roughly when $M/R \ge 1/3$). We compare the characteristic frequencies of these resonant polar modes to the axial modes first found by Chandrasekhar and Ferrari [{\em Proc. Roy. Soc. London A} {\bf 434} 449 (1991)]. We find that the two spectra approach each other as the star is made more compact. The oscillation frequencies of the corresponding polar and axial modes agree to within a percent for stars more compact than $M/R = 0.42$. At the same time, the damping times are slightly different. The results illustrate that there is no real difference between the origin of these axial and polar modes: They are essentially spacetime modes.