Constraints from LIGO O3 Data on Gravitational-wave Emission Due to R-modes in the Glitching Pulsar PSR J0537–6910

We present a search for continuous gravitational-wave emission due to r-modes in the pulsar PSR J0537–6910 using data from the LIGO–Virgo Collaboration observing run O3. PSR J0537–6910 is a young energetic X-ray pulsar and is the most frequent glitcher known. The inter-glitch braking index of the pulsar suggests that gravitational-wave emission due to r-mode oscillations may play an important role in the spin evolution of this pulsar. Theoretical models confirm this possibility and predict emission at a level that can be probed by ground-based detectors. In order to explore this scenario, we search for r-mode emission in the epochs between glitches by using a contemporaneous timing ephemeris obtained from NICER data. We do not detect any signals in the theoretically expected band of 86–97 Hz, and report upper limits on the amplitude of the gravitational waves. Our results improve on previous amplitude upper limits from r-modes in J0537-6910 by a factor of up to 3 and place stringent constraints on theoretical models for r-mode-driven spin-down in PSR J0537–6910, especially for higher frequencies at which our results reach below the spin-down limit defined by energy conservation.

Modelling spin evolution of magnetars

The origin and fate of magnetars (young, extremely magnetized neutron stars, NSs) remains unsolved. Probing their evolution is therefore crucial for investigating possible links to other species of isolated NSs, such as the X-ray dim NSs (XDINSs) and rotating radio transients (RRATs). Here we investigate the spin evolution of magnetars. Two avenues of evolution are considered: one with exponentially decaying B-fields, the other with sub- and super-exponential decay. Using Monte Carlo methods, we synthesize magnetar populations using different input distributions and physical parameters, such as for the initial spin period, its time derivative and the B-field decay timescale. Additionally, we introduce a fade-away procedure that can account for the fading of old magnetars, and we briefly discuss the effect of alignment of the B-field and spin axes. Imposing the Galactic core-collapse supernova rate of ∼20 kyr−1 as a strict upper limit on the magnetar birthrate and comparing the synthetic populations to the observed one using both manual and automatic optimization algorithms for our input parameter study, we find that the B-field must decay exponentially or super-exponentially with a characteristic decay timescale of 0.5 − 10 kyr (with a best value of ∼4 kyr). In addition, the initial spin period must be less than 2 sec. If these constraints are kept, we conclude that there are multiple choices of input physics that can reproduce the observed magnetar population reasonably well. We also conclude that magnetars may well be evolutionary linked to the population of XDINSs, whereas they are in general unlikely to evolve into RRATs.

Massive black hole binary inspiral and spin evolution in a cosmological framework

Massive black hole (MBH) binary inspiral time scales are uncertain, and their spins are even more poorly constrained. Spin misalignment introduces asymmetry in the gravitational radiation, which imparts a recoil kick to the merged MBH. Understanding how MBH binary spins evolve is crucial for determining their recoil velocities, their gravitational wave (GW) waveforms detectable with LISA, as well as their retention rate in galaxies. Here we introduce a sub-resolution model for gas- and GW-driven MBH binary spin evolution using accreting MBHs from the Illustris cosmological hydrodynamics simulations. We also model binary inspiral via dynamical friction, stellar scattering, viscous gas drag, and GW emission. Our model assumes that the circumbinary disk always removes angular momentum from the binary. It also assumes differential accretion, which causes greater alignment of the secondary MBH spin in unequal-mass mergers. We find that 47% of the MBHs in our population merge by z = 0. Of these, 19% have misaligned primaries and 10% have misaligned secondaries at the time of merger in our fiducial model with initial eccentricity of 0.6 and accretion rates from Illustris. The MBH misalignment fraction depends strongly on the accretion disc parameters, however. Reducing accretion rates by a factor of 100, in a thicker disc, yields 79% and 42% misalignment for primaries and secondaries, respectively. Even in the more conservative fiducial model, more than 12% of binaries experience recoils of > 500km/s, which could displace them at least temporarily from galactic nuclei. We additionally find that a significant number of systems experience strong precession.

Black hole spin evolution in warped accretion discs

ABSTRACT Massive black holes (BHs) inhabiting galactic nuclei can be described by two parameters only, i.e. mass and spin, that change through cosmic time in response to accretion and merger events. While most numerical simulations accurately track the BH mass, spin evolution is rarely taken into account. In this work, we implement and validate a self-consistent sub-grid model for the evolution of the BH mass and spin via gas accretion in the hydrodynamics code gizmo. The model assumes that accretion from resolved scales does not occur instantaneously but is mediated by a sub-grid geometrically thin α-disc. After validating our model semi-analytically, we test it in an idealized environment consisting of a circumnuclear disc, where gas accretion on to the accretion disc is consistently determined by gizmo. In the absence of any accretion-related feedback, the spin evolution closely traces that observed in the semi-analytical models, and depends on the free parameters of our implementation, such as the initial BH spin, angular momentum of the accretion disc, and radius at which the gas inflow circularizes. In gizmo, we also couple our model with the biconical-outflow model presented in a companion paper, wherein the feedback axis is always aligned with the BH spin. In this last case, the evolution of the central BH differs significantly from the previous cases, since the feedback process modifies the gas dynamics and its inflow rates from resolved scales. Such an interaction cannot be modelled by simple semi-analytical models and should be treated using full N-body hydrodynamical simulations.