Erratum: The runaway instability of thick discs around black holes - II. Non-constant angular momentum discs

GW170817: Implications for the Stochastic Gravitational-Wave Background from Compact Binary Coalescences

The LIGO Scientific and Virgo Collaborations have announced the first detection of gravitational waves from the coalescence of two neutron stars. The merger rate of binary neutron stars estimated from this event suggests that distant, unresolvable binary neutron stars create a significant astrophysical stochastic gravitational-wave background. The binary neutron star background will add to the background from binary black holes, increasing the amplitude of the total astrophysical background relative to previous expectations. In the Advanced LIGO-Virgo frequency band most sensitive to stochastic backgrounds (near 25 Hz), we predict a total astrophysical background with amplitude $\Omega_{\rm…

Towards modelling the central engine of short GRBs

Numerical relativity simulations of non-vacuum spacetimes have reached a status where a complete description of the inspiral, merger and post-merger stages of the late evolution of close binary neutron systems is possible. Determining the properties of the black-hole-torus system produced in such an event is a key aspect to understand the central engine of short-hard gamma-ray bursts (sGRBs). Of the many properties characterizing the torus, the total rest-mass is the most important one, since it is the torus' binding energy which can be tapped to extract the large amount of energy necessary to power the sGRB emission. In addition, the rest-mass density and angular momentum distribution in t…

Are pulsars born with a hidden magnetic field?

The observation of several neutron stars in the center of supernova remnants and with significantly lower values of the dipolar magnetic field than the average radio-pulsar population has motivated a lively debate about their formation and origin, with controversial interpretations. A possible explanation requires the slow rotation of the proto-neutron star at birth, which is unable to amplify its magnetic field to typical pulsar levels. An alternative possibility, the hidden magnetic field scenario, considers the accretion of the fallback of the supernova debris onto the neutron star as responsible for the submergence (or screening) of the field and its apparently low value. In this paper …

Spanish Relativity Meeting (ERE 2014): almost 100 years after Einstein's revolution

This volume presents the proceedings of the international scientific conference ''Spanish Relativity Meeting (ERE 2014): almost 100 years after Einstein's revolution''. The conference was devoted to discussing the current state-of-the-art of a wide variety of topics of research in the fields of Gravitation and General Relativity in the ''pre-centennial'' year of General Relativity. The name of the conference was chosen to highlight the importance of the upcoming one hundredth anniversary of Einstein's theory of General Relativity, officially established by the Internal Society on General Relativity and Gravitation in November 25th, 2015. In particular, the conference was organized along thr…

Gravitational Waves from the Papaloizou-Pringle Instability in Black-Hole-Torus Systems

Black hole (BH)--torus systems are promising candidates for the central engine of gamma-ray bursts (GRBs), and also possible outcomes of the collapse of supermassive stars to supermassive black holes (SMBHs). By three-dimensional general relativistic numerical simulations, we show that an $m=1$ nonaxisymmetric instability grows for a wide range of self-gravitating tori orbiting BHs. The resulting nonaxisymmetric structure persists for a timescale much longer than the dynamical one, becoming a strong emitter of large amplitude, quasiperiodic gravitational waves. Our results indicate that both, the central engine of GRBs and newly formed SMBHs, can be strong gravitational wave sources observa…

NADA: A new code for studying self-gravitating tori around black holes

We present a new two-dimensional numerical code called Nada designed to solve the full Einstein equations coupled to the general relativistic hydrodynamics equations. The code is mainly intended for studies of self-gravitating accretion disks (or tori) around black holes, although it is also suitable for regular spacetimes. Concerning technical aspects the Einstein equations are formulated and solved in the code using a formulation of the standard 3+1 (ADM) system, the so-called BSSN approach. A key feature of the code is that derivative terms in the spacetime evolution equations are computed using a fourth-order centered finite difference approximation in conjunction with the Cartoon metho…

Estimate of the gravitational-wave background from the observed cosmological distribution of quasars

We study the gravitational-wave background from the observed cosmological quasar distribution. Using the DR9Q quasar catalogue from the ninth data release of the Sloan Digital Sky Survey (SDSS), we create a complete, statistically consistent sample of quasars from $z=0.3$ to $5.4$. Employing the spectroscopic information from the catalogue we estimate the masses of the supermassive black holes hosted by the quasars in the sample, resulting in a log-normal distribution of mean $10^{8.32\pm0.33}M_{\odot}$. The computation of the individual gravitational-wave strains relies on specific functional forms derived from simulations of gravitational collapse and mergers of massive black hole binarie…

Modulating magnetar emission by magneto-elastic oscillations

We present a new numerical tool to calculate the emission of highly magnetized neutron stars (magnetars) and apply it to describe the quasi-periodic oscillations (QPOs) observed in magnetar giant flares. In previous work we have developed a model of magneto-elastic oscillations of magnetars that allows to reproduce the observed frequencies. These QPOs can couple to the star's exterior through the magnetic field and induce currents in the magnetosphere that provide scattering targets for resonant cyclotron scattering of the photons. The scattering is calculated with a Monte-Carlo approach and it is coupled to a code that calculates the momentum distribution of the charge carriers as an one-d…

Dynamical formation of a hairy black hole in a cavity from the decay of unstable solitons

Recent numerical relativity simulations within the Einstein--Maxwell--(charged-)Klein-Gordon (EMcKG) system have shown that the non-linear evolution of a superradiantly unstable Reissner-Nordstr\"om black hole (BH) enclosed in a cavity, leads to the formation of a BH with scalar hair. Perturbative evidence for the stability of such hairy BHs has been independently established, confirming they are the true endpoints of the superradiant instability. The same EMcKG system admits also charged scalar soliton-type solutions, which can be either stable or unstable. Using numerical relativity techniques, we provide evidence that the time evolution of some of these $\textit{unstable}$ solitons leads…

Non-linear axisymmetric pulsations of rotating relativistic stars in the conformal flatness approximation

We study non-linear axisymmetric pulsations of rotating relativistic stars using a general relativistic hydrodynamics code under the assumption of a conformal flatness. We compare our results to previous simulations where the spacetime dynamics was neglected. The pulsations are studied along various sequences of both uniformly and differentially rotating relativistic polytropes with index N = 1. We identify several modes, including the lowest-order l = 0, 2, and 4 axisymmetric modes, as well as several axisymmetric inertial modes. Differential rotation significantly lowers mode frequencies, increasing prospects for detection by current gravitational wave interferometers. We observe an exten…

Imprints of superfluidity on magnetoelastic quasiperiodic oscillations of soft gamma-ray repeaters.

Our numerical simulations show that axisymmetric, torsional, magnetoelastic oscillations of magnetars with a superfluid core can explain the whole range of observed quasiperiodic oscillations (QPOs) in the giant flares of soft gamma-ray repeaters. There exist constant phase QPOs at $f\ensuremath{\lesssim}150\text{ }\text{ }\mathrm{Hz}$ and resonantly excited high-frequency QPOs ($fg500\text{ }\text{ }\mathrm{Hz}$), in good agreement with observations. The range of magnetic field strengths required to match the observed QPO frequencies agrees with that from spin-down estimates. These results suggest that there is at least one superfluid species in magnetar cores.

A Flux-Split Algorithm Applied to Relativistic Flows

The equations of RFD can be written as a hyperbolic system of conservation laws by choosing an appropriate vector of unknowns. We give an explicit formulation of the full spectral decomposition of the Jacobian matrices associated with the fluxes in each spatial direction, which is the essential ingredient of the techniques we propose in this paper. These techniques are based on the recently derived flux formula of Marquina, a new way to compute the numerical flux at a cell interface which leads to a conservative, upwind numerical scheme. Using the spectral decompositions in a fundamental way, we construct high order versions of the basic first-order scheme described by R. Donat and A. Marqu…

All-sky search in early O3 LIGO data for continuous gravitational-wave signals from unknown neutron stars in binary systems

Rapidly spinning neutron stars are promising sources of continuous gravitational waves. Detecting such a signal would allow probing of the physical properties of matter under extreme conditions. A significant fraction of the known pulsar population belongs to binary systems. Searching for unknown neutron stars in binary systems requires specialized algorithms to address unknown orbital frequency modulations. We present a search for continuous gravitational waves emitted by neutron stars in binary systems in early data from the third observing run of the Advanced LIGO and Advanced Virgo detectors using the semicoherent, GPU-accelerated, binaryskyhough pipeline. The search analyzes the most s…

Numerical Simulations of Relativistic Wind Accretion onto Black Holes Using Godunov-Type Methods

We have studied numerically the so-called Bondi-Hoyle (wind) accretion onto a rotating black hole in general relativity. We have used the Kerr-Schild form of the Kerr metric, free of coordinate singularities at the black hole horizon. The ‘test-fluid’ approximation has been adopted, assuming no dynamical evolution of the gravitational field. We have used a formulation of the relativistic hydrodynamic equations which casts them into a first-order hyperbolic system of conservation laws. Our studies were performed using a Godunov-type scheme based on Marquina’s flux-formula.

Universal Relations for Gravitational-Wave Asteroseismology of Protoneutron Stars

State-of-the-art numerical simulations of core-collapse supernovae reveal that the main source of gravitational waves is the excitation of proto-neutron star modes during post-bounce evolution. In this work we derive universal relations that relate the frequencies of the most common oscillation modes observed, i.e. g-modes, p-modes and the f-mode, with fundamental properties of the system, such as the surface gravity of the proto-neutron star or the mean density in the region enclosed by the shock. These relations are independent of the equation of state, the neutrino treatment, and the progenitor mass and hence can be used to build methods to infer proto-neutron star properties from gravit…

Classification of gravitational-wave glitches via dictionary learning

We present a new method for the classification of transient noise signals (or glitches) in advanced gravitational-wave interferometers. The method uses learned dictionaries (a supervised machine learning algorithm) for signal denoising, and untrained dictionaries for the final sparse reconstruction and classification. We use a data set of 3000 simulated glitches of three different waveform morphologies, comprising 1000 glitches per morphology. These data are embedded in non-white Gaussian noise to simulate the background noise of advanced LIGO in its broadband configuration. Our classification method yields a 96% accuracy for a large range of initial parameters, showing that learned diction…

Search for Tensor, Vector, and Scalar Polarizations in the Stochastic Gravitational-Wave Background

The detection of gravitational waves with Advanced LIGO and Advanced Virgo has enabled novel tests of general relativity, including direct study of the polarization of gravitational waves. While general relativity allows for only two tensor gravitational-wave polarizations, general metric theories can additionally predict two vector and two scalar polarizations. The polarization of gravitational waves is encoded in the spectral shape of the stochastic gravitational-wave background, formed by the superposition of cosmological and individually-unresolved astrophysical sources. Using data recorded by Advanced LIGO during its first observing run, we search for a stochastic background of generic…

3-D collapse of rotating stars to Kerr black holes

We study gravitational collapse of uniformly rotating neutron stars to Kerr black holes, using a new three-dimensional, fully general relativistic hydrodynamics code, which uses high-resolution shock-capturing techniques and a conformal traceless formulation of the Einstein equations. We investigate the gravitational collapse by carefully studying not only the dynamics of the matter, but also that of the trapped surfaces, i.e. of both the apparent and event horizons formed during the collapse. The use of these surfaces, together with the dynamical horizon framework, allows for a precise measurement of the black-hole mass and spin. The ability to successfully perform these simulations for su…

General relativistic simulations of pasive-magneto-rotational core collapse with microphysics

This paper presents results from axisymmetric simulations of magneto-rotational stellar core collapse to neutron stars in general relativity using the passive field approximation for the magnetic field. These simulations are performed using a new general relativistic numerical code specifically designed to study this astrophysical scenario. The code is based on the conformally-flat approximation of Einstein's field equations and conservative formulations of the magneto-hydrodynamics equations. The code has been recently upgraded to incorporate a tabulated, microphysical equation of state and an approximate deleptonization scheme. This allows us to perform the most realistic simulations of m…

Robustness of a high-resolution central scheme for hydrodynamic simulations in full general relativity

A recent paper by Lucas-Serrano et al. indicates that a high-resolution central (HRC) scheme is robust enough to yield accurate hydrodynamical simulations of special relativistic flows in the presence of ultrarelativistic speeds and strong shock waves. In this paper we apply this scheme in full general relativity (involving {\it dynamical} spacetimes), and assess its suitability by performing test simulations for oscillations of rapidly rotating neutron stars and merger of binary neutron stars. It is demonstrated that this HRC scheme can yield results as accurate as those by the so-called high-resolution shock-capturing (HRSC) schemes based upon Riemann solvers. Furthermore, the adopted HRC…

General Relativistic Simulations of Accretion Disks Around Tilted Kerr Black Holes

We simulate the dynamics of self-gravitating accretion disks around tilted Kerr black holes (BH) in full 3D general relativity. For this purpose we employ the EinsteinToolkit, using the thorn McLachlan for the evolution of the spacetime via the BSSN formalism of the Einstein equations and the thorn GRHydro for the evolution of the hydrodynamics, using a 3D Cartesian mesh with adaptive mesh refinement. We investigate the effects of the tilt angle between the disk angular momentum and BH spin vector on the dynamics of these systems as the disk evolves in the tilted spacetime. By evolving the spacetime and matter fields, we are able to observe how both BH and disk react and evolve in the tilte…

Rotational properties of hypermassive neutron stars from binary mergers

Determining the differential-rotation law of compact stellar objects produced in binary neutron stars mergers or core-collapse supernovae is an old problem in relativistic astrophysics. Addressing this problem is important because it impacts directly on the maximum mass these objects can attain and hence on the threshold to black-hole formation under realistic conditions. Using the results from a large number of numerical simulations in full general relativity of binary neutron star mergers described with various equations of state and masses, we study the rotational properties of the resulting hypermassive neutron stars. We find that the angular-velocity distribution shows only a modest de…

Search for anisotropic gravitational-wave backgrounds using data from Advanced LIGO and Advanced Virgo's first three observing runs

We report results from searches for anisotropic stochastic gravitational-wave backgrounds using data from the first three observing runs of the Advanced LIGO and Advanced Virgo detectors. For the first time, we include Virgo data in our analysis and run our search with a new efficient pipeline called {\tt PyStoch} on data folded over one sidereal day. We use gravitational-wave radiometry (broadband and narrow band) to produce sky maps of stochastic gravitational-wave backgrounds and to search for gravitational waves from point sources. A spherical harmonic decomposition method is employed to look for gravitational-wave emission from spatially-extended sources. Neither technique found eviden…

Numerical relativity simulations of thick accretion disks around tilted Kerr black holes

In this work we present 3D numerical relativity simulations of thick accretion disks around tilted Kerr BH. We investigate the evolution of three different initial disk models with a range of initial black hole spin magnitudes and tilt angles. For all the disk-to-black hole mass ratios considered (0.044-0.16) we observe significant black hole precession and nutation during the evolution. This indicates that for such mass ratios, neglecting the self-gravity of the disks by evolving them in a fixed background black hole spacetime is not justified. We find that the two more massive models are unstable against the Papaloizou-Pringle (PP) instability and that those PP-unstable models remain unst…

All-sky search for long-duration gravitational wave transients in the first Advanced LIGO observing run

Made available in DSpace on 2018-11-26T17:45:14Z (GMT). No. of bitstreams: 0 Previous issue date: 2018-03-22 Australian Research Council Council of Scientific and Industrial Research of India Department of Science and Technology, India Science AMP; Engineering Research Board (SERB), India Ministry of Human Resource Development, India Spanish Agencia Estatal de Investigacion Vicepresidencia i Conselleria d'Innovacio, Recerca i Turisme Conselleria d'Educacio i Universitat del Govern de les Illes Balears Conselleria d'Educacio, Investigacio, Cultura i Esport de la Generalitat Valenciana National Science Centre of Poland Swiss National Science Foundation (SNSF) Russian Foundation for Basic Rese…

Gravitational-wave parameter inference using Deep Learning

We explore machine learning methods to detect gravitational waves (GW) from binary black hole (BBH) mergers using deep learning (DL) algorithms. The DL networks are trained with gravitational waveforms obtained from BBH mergers with component masses randomly sampled in the range from 5 to 100 solar masses and luminosity distances from 100 Mpc to, at least, 2000 Mpc. The GW signal waveforms are injected in public data from the O2 run of the Advanced LIGO and Advanced Virgo detectors, in time windows that do not coincide with those of known detected signals, and the data from each detector in the Advanced LIGO and Advanced Virgo network is combined into a unique RGB image. We show that a clas…

Total-variation methods for gravitational-wave denoising: Performance tests on Advanced LIGO data

We assess total-variation methods to denoise gravitational-wave signals in real noise conditions, by injecting numerical-relativity waveforms from core-collapse supernovae and binary black hole mergers in data from the first observing run of Advanced LIGO. This work is an extension of our previous investigation where only Gaussian noise was used. Since the quality of the results depends on the regularization parameter of the model, we perform an heuristic search for the value that produces the best results. We discuss various approaches for the selection of this parameter, either based on the optimal, mean, or multiple values, and compare the results of the denoising upon these choices. Mor…

Towards asteroseismology of core-collapse supernovae with gravitational-wave observations – I. Cowling approximation

Gravitational waves from core-collapse supernovae are produced by the excitation of different oscillation modes in the protoneutron star (PNS) and its surroundings, including the shock. In this work we study the relationship between the post-bounce oscillation spectrum of the PNS–shock system and the characteristic frequencies observed in gravitational-wave signals from core-collapse simulations. This is a fundamental first step in order to develop a procedure to infer astrophysical parameters of the PNS formed in core-collapse supernovae. Our method combines information from the oscillation spectrum of the PNS, obtained through linear perturbation analysis in general relativity of a backgr…

On the dynamics of tilted black hole-torus systems

We present results from three-dimensional, numerical relativity simulations of a {\it tilted} black hole-thick accretion disc system. The simulations are analysed using tracer particles in the disc which are advected with the flow. Such tracers, which we employ in these new simulations for the first time, provide a powerful means to analyse in detail the complex dynamics of tilted black hole-torus systems. We show how its use helps to gain insight in the overall dynamics of the system, discussing the origin of the observed black hole precession and the development of a global non-axisymmetric $m=1$ mode in the disc. Our three-dimensional simulations show the presence of quasi-periodic oscil…

Quasi-stationary solutions of self-gravitating scalar fields around black holes

Recent perturbative studies have shown the existence of long-lived, quasi-stationary configurations of scalar fields around black holes. In particular, such configurations have been found to survive for cosmological timescales, which is a requirement for viable dark matter halo models in galaxies based on such type of structures. In this paper we perform a series of numerical relativity simulations of dynamical non-rotating black holes surrounded by self-gravitating scalar fields. We solve numerically the coupled system of equations formed by the Einstein and the Klein-Gordon equations under the assumption of spherical symmetry using spherical coordinates. Our results confirm the existence …

A Numerical Study of Relativistic Jets

Numerical simulations of supersonic jets are able to explain the structures observed in many VLA images of radio sources. The improvements achieved in classical simulations (see Hardee, these proceedings) are in contrast with the almost complete lack of relativistic simulations the reason being that numerical difficulties arise from the highly relativistic flows typical of extragalactic jets. For our study, we have developed a two-dimensional code which is based on (i) an explicit conservative differencing of the special relativistic hydrodynamics (SRH) equations and (ii) the use of an approximate Riemann solver (see Marti et al. 1995a,b and references therein).

The runaway instability of thick discs around black holes. II. Non constant angular momentum discs

We present results from a comprehensive number of relativistic, time-dependent, axisymmetric simulations of the runaway instability of non-constant angular momentum thick discs around black holes. This second paper extends earlier results where only constant angular momentum discs were considered. All relevant aspects of the theory of stationary thick discs around rotating black holes, necessary to build the initial state in our simulations, are presented in great detail. The angular momentum of the discs is assumed to increase outwards with the radial distance according to a power law. The main simplifying assumptions of our approach are not to include magnetic fields and self-gravity in t…

GW190521 as a Merger of Proca Stars: A Potential New Vector Boson of 8.7×10−13  eV

Advanced LIGO-Virgo have reported a short gravitational-wave signal (GW190521) interpreted as a quasicircular merger of black holes, one at least populating the pair-instability supernova gap, that formed a remnant black hole of ${M}_{f}\ensuremath{\sim}142\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$ at a luminosity distance of ${d}_{L}\ensuremath{\sim}5.3\text{ }\text{ }\mathrm{Gpc}$. With barely visible pre-merger emission, however, GW190521 merits further investigation of the pre-merger dynamics and even of the very nature of the colliding objects. We show that GW190521 is consistent with numerically simulated signals from head-on collisions of two (equal mass and spin) horizonless vecto…

Quasistationary solutions of self-gravitating scalar fields around black holes

Recent perturbative studies have shown the existence of long-lived, quasistationary configurations of scalar fields around black holes. In particular, such configurations have been found to survive for cosmological time scales, which is a requirement for viable dark matter halo models in galaxies based on such types of structures. In this paper we perform a series of numerical relativity simulations of dynamical nonrotating black holes surrounded by self-gravitating scalar fields. We solve numerically the coupled system of equations formed by the Einstein and the Klein-Gordon equations under the assumption of spherical symmetry using spherical coordinates. Our results confirm the existence …

Implementation of a simplified approach to radiative transfer in general relativity

We describe in detail the implementation of a simplified approach to radiative transfer in general relativity by means of the well-known neutrino leakage scheme (NLS). In particular, we carry out an extensive investigation of the properties and limitations of the NLS for isolated relativistic stars to a level of detail that has not been discussed before in a general-relativistic context. Although the numerous tests considered here are rather idealized, they provide a well-controlled environment in which to understand the relationship between the matter dynamics and the neutrino emission, which is important in order to model the neutrino signals from more complicated scenarios, such as binar…

Gravitational waves from binary black hole mergers surrounded by scalar field clouds: Numerical simulations and observational implications

We show how gravitational-wave observations of binary black hole (BBH) mergers can constrain the physical characteristics of a scalar field cloud parameterized by mass $\tilde{\mu}$ and strength $\phi_0$ that may surround them. We numerically study the inspiraling equal-mass, non-spinning BBH systems dressed in such clouds, focusing especially on the gravitational-wave signals emitted by their merger-ringdown phase. These waveforms clearly reveal that larger values of $\tilde{\mu}$ or $\phi_0$ cause bigger changes in the amplitude and frequency of the scalar-field-BBH ringdown signals. We show that the numerical waveforms of scalar-field-BBHs can be modelled as chirping sine-Gaussians, with…

Tests of General Relativity with GW170817

The recent discovery by Advanced LIGO and Advanced Virgo of a gravitational wave signal from a binary neutron star inspiral has enabled tests of general relativity (GR) with this new type of source. This source, for the first time, permits tests of strong-field dynamics of compact binaries in presence of matter. In this paper, we place constraints on the dipole radiation and possible deviations from GR in the post-Newtonian coefficients that govern the inspiral regime. Bounds on modified dispersion of gravitational waves are obtained; in combination with information from the observed electromagnetic counterpart we can also constrain effects due to large extra dimensions. Finally, the polari…

Accurate evolutions of unequal-mass neutron-star binaries: properties of the torus and short GRB engines

We present new results from accurate and fully general-relativistic simulations of the coalescence of unmagnetized binary neutron stars with various mass ratios. The evolution of the stars is followed through the inspiral phase, the merger and prompt collapse to a black hole, up until the appearance of a thick accretion disk, which is studied as it enters and remains in a regime of quasi-steady accretion. Although a simple ideal-fluid equation of state with \Gamma=2 is used, this work presents a systematic study within a fully general relativistic framework of the properties of the resulting black-hole--torus system produced by the merger of unequal-mass binaries. More specifically, we show…

Accretion-induced quasinormal mode excitation of a Schwarzschild black hole

By combining the numerical solution of the nonlinear hydrodynamics equations with the solution of the linear inhomogeneous Zerilli-Moncrief and Regge-Wheeler equations we investigate the properties of the gravitational radiation emitted during the axisymmetric accretion of matter onto a Schwarzschild black hole. The matter models considered include quadrupolar dust shells and thick accretion disks, permitting us to simulate situations which may be encountered at the end stages of stellar gravitational collapse or binary neutron star merger. We focus on the interference pattern appearing in the energy spectra of the emitted gravitational waves and on the amount of excitation of the quasi-nor…

Magnetized accretion disks around Kerr black holes with scalar hair - I. Constant angular momentum disks

Testing the true nature of black holes - the no-hair hypothesis - will become increasingly more precise in the next few years as new observational data is collected in both the gravitational wave channel and the electromagnetic channel. In this paper we consider numerically generated spacetimes of Kerr black holes with synchronised scalar hair and build stationary models of magnetized thick disks (or tori) around them. Our approach assumes that the disks are not self-gravitating, they obey a polytropic equation of state, the distribution of their specific angular momentum is constant, and they are marginally stable, i.e. the disks completely fill their Roche lobe. Moreover, contrary to exis…

On the runaway instability of self-gravitating torus around black holes

Black holes surrounded by self-gravitating tori are astrophysical systems which may naturally form following the core collapse of a massive star or the merger of two neutron stars. We present here results from fully general relativistic numerical simulations of such systems in order to assess the influence of the torus self-gravity on the onset of the so-called runaway instability. This instability, which might drive the rapid accretion of the disk on shorter timescales than those required to power a relativistic fireball, potentially challenges current models of gamma-ray bursts. Our simulations indicate that the self-gravity of the torus does not actually favour the onset of the instabili…

Erratum: “Searches for Gravitational Waves from Known Pulsars at Two Harmonics in 2015–2017 LIGO Data” (2019, ApJ, 879, 10)

Two analysis errors have been identified that affect the results for a handful of the high-value pulsars given in Table 1 of Abbott et al. (2019). One affects the Bayesian analysis for the five pulsars that glitched during the analysis period, and the other affects the 5n-vector analysis for J0711-6830. Updated results after correcting the errors are shown in Table 1, which now supersedes the results given for those pulsars in Table 1 of Abbott et al. (2019). Updated versions of figures can be seen in Figures 1-4. Bayesian analysis.-For the glitching pulsars, the signal phase evolution caused by the glitch was wrongly applied twice and was therefore not consistent with our expected model of…

Dynamics of oscillating relativistic tori around Kerr black holes

We present a comprehensive numerical study of the dynamics of relativistic axisymmetric accretion tori with a power-law distribution of specific angular momentum orbiting in the background spacetime of a Kerr black hole. By combining general relativistic hydrodynamics simulations with a linear perturbative approach we investigate the main dynamical properties of these objects over a large parameter space. The astrophysical implications of our results extend and improve two interesting results that have been recently reported in the literature. Firstly, the induced quasi-periodic variation of the mass quadrupole moment makes relativistic tori of nuclear matter densities, as those formed duri…

Status of Advanced Virgo

The LIGO and the Virgo collaborations have recently announced the first detections of Gravitational Waves. Due to their weak amplitude, Gravitational Waves are expected to produce a very small effect on free-falling masses, which undergo a displacement of the order of 10-18 m for a Km-scale mutual distance. This discovery showed that interferometric detectors are suitable to reveal such a feeble effect, and therefore represent a new tool for astronomy, astrophysics and cosmology in the understanding of the Universe. To better reconstruct the position of the Gravitational Wave source and increase the signal-to-noise ratio of the events by means of multiple coincidence, a network of detectors…

Denoising of gravitational wave signals via dictionary learning algorithms

Gravitational wave astronomy has become a reality after the historical detections accomplished during the first observing run of the two advanced LIGO detectors. In the following years, the number of detections is expected to increase significantly with the full commissioning of the advanced LIGO, advanced Virgo and KAGRA detectors. The development of sophisticated data analysis techniques to improve the opportunities of detection for low signal-to-noise-ratio events is, hence, a most crucial effort. In this paper, we present one such technique, dictionary-learning algorithms, which have been extensively developed in the last few years and successfully applied mostly in the context of image…

Dynamical formation of a Reissner-Nordström black hole with scalar hair in a cavity

In a recent Letter [Sanchis-Gual et al., Phys. Rev. Lett. 116, 141101 (2016)], we presented numerical relativity simulations, solving the full Einstein--Maxwell--Klein-Gordon equations, of superradiantly unstable Reissner-Nordstr\"om black holes (BHs), enclosed in a cavity. Low frequency, spherical perturbations of a charged scalar field trigger this instability. The system's evolution was followed into the nonlinear regime, until it relaxed into an equilibrium configuration, found to be a hairy BH: a charged horizon in equilibrium with a scalar field condensate, whose phase is oscillating at the (final) critical frequency. Here, we investigate the impact of adding self-interactions to the …

General Relativistic Hydrodynamics and Magnetohydrodynamics: Hyperbolic Systems in Relativistic Astrophysics

Synchronised gravitational atoms from mergers of bosonic stars

If ultralight bosonic fields exist in Nature as dark matter, superradiance spins down rotating black holes (BHs), dynamically endowing them with equilibrium bosonic clouds, here dubbed synchronised gravitational atoms (SGAs). The self-gravity of these same fields, on the other hand, can lump them into (scalar or vector) horizonless solitons known as bosonic stars (BSs). We show that the dynamics of BSs yields a new channel forming SGAs. We study BS binaries that merge to form spinning BHs. After horizon formation, the BH spins up by accreting the bosonic field, but a remnant lingers around the horizon. If just enough angular momentum is present, the BH spin up stalls precisely as the remnan…

Imprints of superfluidity on magneto-elastic QPOs of SGRs

Our numerical simulations show that axisymmetric, torsional, magneto-elastic oscillations of magnetars with a superfluid core can explain the whole range of observed quasi-periodic oscillations (QPOs) in the giant flares of soft gamma-ray repeaters. There exist constant phase, magneto-elastic QPOs at both low (f<150 Hz) and high frequencies (f>500 Hz), in full agreement with observations. The range of magnetic field strengths required to match the observed QPO frequencies agrees with that from spin-down estimates. These results strongly suggest that neutrons in magnetar cores are superfluid.

Confusing Head-On Collisions with Precessing Intermediate-Mass Binary Black Hole Mergers

We report a degeneracy between the gravitational-wave signals from quasi-circular precessing black-hole mergers and those from extremely eccentric mergers, namely head-on collisions. Performing model selection on numerically simulated signals of head-on collisions using models for quasi-circular binaries we find that, for signal-to-noise ratios of 15 and 25, typical of Advanced LIGO observations, head-on mergers with respective total masses of $M\in (125,300)M_\odot$ and $M\in (200,440)M_\odot$ would be identified as precessing quasi-circular intermediate-mass black hole binaries, located at a much larger distance. Ruling out the head-on scenario would require to perform model selection usi…

Magneto-elastic oscillations of neutron stars: exploring different magnetic field configurations

We study magneto-elastic oscillations of highly magnetized neutron stars (magnetars) which have been proposed as an explanation for the quasi-periodic oscillations (QPOs) appearing in the decaying tail of the giant flares of soft gamma-ray repeaters (SGRs). We extend previous studies by investigating various magnetic field configurations, computing the Alfv��n spectrum in each case and performing magneto-elastic simulations for a selected number of models. By identifying the observed frequencies of 28 Hz (SGR 1900+14) and 30 Hz (SGR 1806-20) with the fundamental Alfv��n QPOs, we estimate the required surface magnetic field strength. For the magnetic field configurations investigated (dipole…

An introduction to relativistic hydrodynamics

We review formulations of the equations of (inviscid) general relativistic hydrodynamics and (ideal) magnetohydrodynamics, along with methods for their numerical solution. Both systems can be cast as first-order, hyperbolic systems of conservation laws, following the explicit choice of an Eulerian observer and suitable fluid and magnetic field variables. During the last fifteen years, the so-called (upwind) high-resolution shock-capturing schemes based on Riemann solvers have been successfully extended from classical to relativistic fluid dynamics, both special and general. Nowadays, general relativistic hydrodynamical simulations in relativistic astrophysics are routinely performed, partic…

Quasistationary solutions of scalar fields around collapsing self-interacting boson stars

There is increasing numerical evidence that scalar fields can form long-lived quasibound states around black holes. Recent perturbative and numerical relativity calculations have provided further confirmation in a variety of physical systems, including both static and accreting black holes, and collapsing fermionic stars. In this work, we investigate this issue yet again in the context of gravitationally unstable boson stars leading to black-hole formation. We build a large sample of spherically symmetric initial models, both stable and unstable, incorporating a self-interaction potential with a quartic term. The three different outcomes of unstable models, namely, migration to the stable b…

Split Bregman Method for Gravitational Wave Denoising

This paper presents a progress report in our aim to develop a Total Variation algorithm for denoising of gravitational waves. These algorithms, are routinely employed in the context of image processing and they do not need any a priori information on the signals. We apply our method to two different types of numerically-simulated gravitational wave signals, namely burst produced from the core collapse of rotating stars and waveforms from binary black hole mergers, and present a preliminary assessment of its capabilities.

Classification of the core-collapse supernova explosion mechanism with learned dictionaries

Core-collapse supernovae (CCSN) are a prime source of gravitational waves. Estimations of their typical frequencies make them perfect targets for the current network of advanced, ground-based detectors. A successful detection could potentially reveal the underlying explosion mechanism through the analysis of the waveform. This has been illustrated using the SupernovaModel Evidence Extractor (SMEE; Logue et al. (2012)), an algorithm based on principal-component analysis and Bayesian model selection. Here, we present a complementary approach to SMEE based on (supervised) dictionary-learning and show that it is able to reconstruct and classify CCSN signals according to their morphology. Our wa…

Magneto‐elastic oscillations modulating the emission of magnetars

Magneto-elastic oscillations of neutron stars are believed to explain observed quasi-periodic oscillations (QPOs) in the decaying tail of the giant flares of highly magnetized neutron stars (magnetars). Strong efforts of the theoretical modelling from different groups have increased our understanding of this phenomenon significantly. Here, we discuss some constraints on the matter in neutron stars that arise if the interpretation of the observations in terms of superfluid, magneto-elastic oscillations is correct. To explain the observed modulation of the light curve of the giant flare, we describe a model that allows the QPOs to couple to the stellar exterior through the magnetic field. In …

Preface: Amazonia in the route of General Relativity

Multi-field, multi-frequency bosonic stars and a stabilization mechanism

Scalar bosonic stars (BSs) stand out as a multi-purpose model of exotic compact objects. We enlarge the landscape of such (asymptotically flat, stationary, everywhere regular) objects by considering multiple fields (possibly) with different frequencies. This allows for new morphologies ${\it and}$ a stabilization mechanism for different sorts of unstable BSs. First, any odd number of complex fields, yields a continuous family of BSs departing from the spherical, equal frequency, $\ell-$BSs. As the simplest illustration, we construct the $\ell$ = ${\it 1}$ ${\it BSs}$ ${\it family}$, that includes several single frequency solutions, including even parity (such as spinning BSs and a toroidal,…

Three-dimensional relativistic simulations of rotating neutron-star collapse to a Kerr black hole

We present a new three-dimensional fully general-relativistic hydrodynamics code using high-resolution shock-capturing techniques and a conformal traceless formulation of the Einstein equations. Besides presenting a thorough set of tests which the code passes with very high accuracy, we discuss its application to the study of the gravitational collapse of uniformly rotating neutron stars to Kerr black holes. The initial stellar models are modelled as relativistic polytropes which are either secularly or dynamically unstable and with angular velocities which range from slow rotation to the mass-shedding limit. We investigate the gravitational collapse by carefully studying not only the dynam…

GW190521: A Binary Black Hole Merger with a Total Mass of 150  M⊙

LIGO Scientific Collaboration and Virgo Collaboration: et al.

Prospects for the inference of inertial modes from hypermassive neutron stars with future gravitational-wave detectors

Some recent, long-term numerical simulations of binary neutron star mergers have shown that the long-lived remnants produced in such mergers might be affected by convective instabilities. Those would trigger the excitation of inertial modes, providing a potential method to improve our understanding of the rotational and thermal properties of neutron stars through the analysis of the modes' imprint in the late post-merger gravitational-wave signal. In this paper we assess the detectability of those modes by injecting numerically generated post-merger waveforms into colored Gaussian noise of second-generation and future detectors. Signals are recovered using BayesWave, a Bayesian data-analysi…

Numerical-relativity simulations of long-lived remnants of binary neutron star mergers

We analyze the properties of the gravitational wave signal emitted after the merger of a binary neutron star system when the remnant survives for more than a 80 ms (and up to 140ms). We employ four different piecewise polytropic equations of state supplemented by an ideal fluid thermal component. We find that the post-merger phase can be subdivided into three phases: an early post-merger phase (where the quadrupole mode and a few subdominant features are active), the intermediate post-merger phase (where only the quadrupole mode is active) and the late post-merger phase (where convective instabilities trigger inertial modes). The inertial modes have frequencies somewhat smaller than the qua…

Relativistic MHD simulations of stellar core collapse and magnetars

We present results from simulations of magneto-rotational stellar core collapse along with Alfven oscillations in magnetars. These simulations are performed with the CoCoA/CoCoNuT code, which is able to handle ideal MHD flows in dynamical spacetimes in general relativity. Our core collapse simulations highlight the importance of genuine magnetic effects, like the magneto-rotational instability, for the dynamics of the flow. For the modelling of magnetars we use the anelastic approximation to general relativistic MHD, which allows for an effective suppression of fluid modes and an accurate description of Alfven waves. We further compute Alfven oscillation frequencies along individual magneti…

Convective Excitation of Inertial Modes in Binary Neutron Star Mergers

We present the first very long-term simulations (extending up to ~140 ms after merger) of binary neutron star mergers with piecewise polytropic equations of state and in full general relativity. Our simulations reveal that at a time of 30-50 ms after merger, parts of the star become convectively unstable, which triggers the excitation of inertial modes. The excited inertial modes are sustained up to several tens of milliseconds and are potentially observable by the planned third-generation gravitational-wave detectors at frequencies of a few kilohertz. Since inertial modes depend on the rotation rate of the star and they are triggered by a convective instability in the postmerger remnant, t…

Measuring the black hole spin direction in 3D Cartesian numerical relativity simulations

We show that the so-called flat-space rotational Killing vector method for measuring the Cartesian components of a black hole spin can be derived from the surface integral of Weinberg's pseudotensor over the apparent horizon surface when using Gaussian normal coordinates in the integration. Moreover, the integration of the pseudotensor in this gauge yields the Komar angular momentum integral in a foliation adapted to the axisymmetry of the spacetime. As a result, the method does not explicitly depend on the evolved lapse $\ensuremath{\alpha}$ and shift ${\ensuremath{\beta}}^{i}$ on the respective time slice, as they are fixed to Gaussian normal coordinates while leaving the coordinate label…

On the black hole from merging binary neutron stars: how fast can it spin?

The merger of two neutron stars will in general lead to the formation of a torus surrounding a black hole whose rotational energy can be tapped to potentially power a short gamma-ray burst. We have studied the merger of equal-mass binaries with spins aligned with the orbital angular momentum to determine the maximum spin the black hole can reach. Our initial data consists of irrotational binaries to which we add various amounts of rotation to increase the total angular momentum. Although the initial data violates the constraint equations, the use of the constraint-damping CCZ4 formulation yields evolutions with violations smaller than those with irrotational initial data and standard formul…

Advanced Virgo Status

Abstract The detection of a gravitational wave signal in September 2015 by LIGO interferometers, announced jointly by LIGO collaboration and Virgo collaboration in February 2016, opened a new era in Astrophysics and brought to the whole community a new way to look at - or “listen” to - the Universe. In this regard, the next big step was the joint observation with at least three detectors at the same time. This configuration provides a twofold benefit: it increases the signal-to-noise ratio of the events by means of triple coincidence and allows a narrower pinpointing of GW sources, and, in turn, the search for Electromagnetic counterparts to GW signals. Advanced Virgo (AdV) is the second ge…

Influence of self-gravity on the runaway instability of black-hole-torus systems.

Results from the first fully general relativistic numerical simulations in axisymmetry of a system formed by a black hole surrounded by a self-gravitating torus in equilibrium are presented, aiming to assess the influence of the torus self-gravity on the onset of the runaway instability. We consider several models with varying torus-to-black hole mass ratio and angular momentum distribution orbiting in equilibrium around a non-rotating black hole. The tori are perturbed to induce the mass transfer towards the black hole. Our numerical simulations show that all models exhibit a persistent phase of axisymmetric oscillations around their equilibria for several dynamical timescales without the …

Spontaneous creation of circularly polarized photons in chiral astrophysical systems

This work establishes a relation between chiral anomalies in curved spacetimes and the radiative content of the gravitational field. In particular, we show that a flux of circularly polarized gravitational waves triggers the spontaneous creation of photons with net circular polarization from the quantum vacuum. Using waveform catalogues we identify precessing binary black holes as astrophysical configurations that emit such gravitational radiation, and then solve the fully non-linear Einstein's equations with numerical relativity to evaluate the net effect. The quantum amplitude for a merger is comparable to the Hawking emission rate of the final black hole, and small to be directly observe…

Neutron star collapse and gravitational waves with a non-convex equation of state

The thermodynamical properties of the equation of state (EoS) of high-density matter (above nuclear saturation density) and the possible existence of exotic states such as phase transitions from nuclear/hadronic matter into quark-gluon plasma, or the appearance of hyperons, may critically influence the stability and dynamics of compact relativistic stars. From a theoretical point of view, establishing the existence of those states requires the analysis of the `convexity' of the EoS. We show indications of the existence of regions in the dense-matter EoS where the thermodynamics may be non-convex as a result of a non-monotonic dependence of the sound speed with the rest-mass density. When th…

DYNAMICAL BAR-MODE INSTABILITY IN DIFFERENTIALLY ROTATING MAGNETIZED NEUTRON STARS

This paper presents a numerical study over a wide parameter space of the likelihood of the dynamical bar-mode instability in differentially rotating magnetized neutron stars. The innovative aspect of this study is the incorporation of magnetic fields in such a context, which have thus far been neglected in the purely hydrodynamical simulations available in the literature. The investigation uses the Cosmos++ code which allows us to perform three dimensional simulations on a cylindrical grid at high resolution. A sample of Newtonian magneto-hydrodynamical simulations starting from a set of models previously analyzed by other authors without magnetic fields has been performed, providing estima…

Magneto-elastic oscillations and the damping of crustal shear modes in magnetars

In a realistic model of magneto-elastic oscillations in magnetars, we find that crustal shear oscillations, often invoked as an explanation of quasi-periodic oscillations (QPOs) seen after giant flares in soft gamma-ray repeaters (SGRs), are damped by resonant absorption on timescales of at most 0.2s, for a lower limit on the dipole magnetic field strength of 5 10 13 G. At higher magnetic field strengths (typical in magnetars) the damping timescale is even shorter, as anticipated by earlier toy-models. We have investigated a range of equations of state and masses and if magnetars are dominated by a dipole magnetic field, our findings exclude torsional shear oscillations of the crust from ex…

First narrow-band search for continuous gravitational waves from known pulsars in advanced detector data

Spinning neutron stars asymmetric with respect to their rotation axis are potential sources of continuous gravitational waves for ground-based interferometric detectors. In the case of known pulsars a fully coherent search, based on matched filtering, which uses the position and rotational parameters obtained from electromagnetic observations, can be carried out. Matched filtering maximizes the signal-to-noise (SNR) ratio, but a large sensitivity loss is expected in case of even a very small mismatch between the assumed and the true signal parameters. For this reason, {\it narrow-band} analyses methods have been developed, allowing a fully coherent search for gravitational waves from known …

Boson stars in Palatini f(R) gravity

Nonlinear dynamics of spinning bosonic stars: formation and stability

We perform numerical evolutions of the fully non-linear Einstein-(complex, massive)Klein-Gordon and Einstein-(complex)Proca systems, to assess the formation and stability of spinning bosonic stars. In the scalar/vector case these are known as boson/Proca stars. Firstly, we consider the formation scenario. Starting with constraint-obeying initial data, describing a dilute, axisymmetric cloud of spinning scalar/Proca field, gravitational collapse towards a spinning star occurs, via gravitational cooling. In the scalar case the formation is transient, even for a non-perturbed initial cloud; a non-axisymmetric instability always develops ejecting all the angular momentum from the scalar star. I…

Numerical evolution of matter in dynamical axisymmetric black hole spacetimes

We have developed a numerical code to study the evolution of self-gravitating matter in dynamic black hole axisymmetric spacetimes in general relativity. The matter fields are evolved with a high-resolution shock-capturing scheme that uses the characteristic information of the general relativistic hydrodynamic equations to build up a linearized Riemann solver. The spacetime is evolved with an axisymmetric ADM code designed to evolve a wormhole in full general relativity. We discuss the numerical and algorithmic issues related to the effective coupling of the hydrodynamical and spacetime pieces of the code, as well as the numerical methods and gauge conditions we use to evolve such spacetime…

Completion of the universal I-Love-Q relations in compact stars including the mass

In a recent paper we applied a rigorous perturbed matching framework to show the amendment of the mass of rotating stars in Hartle's model. Here, we apply this framework to the tidal problem in binary systems. Our approach fully accounts for the correction to the Love numbers needed to obtain the universal $I$-Love-$Q$ relations. We compute the corrected mass vs radius configurations of rotating quark stars, revisiting a classical paper on the subject. These corrections allow us to find a universal relation involving the second-order contribution to the mass $\delta M$. We thus complete the set of universal relations for the tidal problem in binary systems, involving four perturbation param…

GW170814: A Three-Detector Observation of Gravitational Waves from a Binary Black Hole Coalescence

On August 14, 2017 at 10 30:43 UTC, the Advanced Virgo detector and the two Advanced LIGO detectors coherently observed a transient gravitational-wave signal produced by the coalescence of two stellar mass black holes, with a false-alarm rate of 1 in 27 000 years. The signal was observed with a three-detector network matched-filter signal-to-noise ratio of 18. The inferred masses of the initial black holes are 30.5-3.0+5.7M and 25.3-4.2+2.8M (at the 90% credible level). The luminosity distance of the source is 540-210+130 Mpc, corresponding to a redshift of z=0.11-0.04+0.03. A network of three detectors improves the sky localization of the source, reducing the area of the 90% credible regio…

Coherent magneto-elastic oscillations in superfluid magnetars

We study the effect of superfluidity on torsional oscillations of highly magnetised neutron stars (magnetars) with a microphysical equation of state by means of two-dimensional, magnetohydrodynamical- elastic simulations. The superfluid properties of the neutrons in the neutron star core are treated in a parametric way in which we effectively decouple part of the core matter from the oscillations. Our simulations confirm the existence of two groups of oscillations, namely continuum oscillations that are confined to the neutron star core and are of Alfv\'enic character, and global oscillations with constant phase and that are of mixed magneto-elastic type. The latter might explain the quasi-…

Constraining the p -Mode– g -Mode Tidal Instability with GW170817

We analyze the impact of a proposed tidal instability coupling p modes and g modes within neutron stars on GW170817. This nonresonant instability transfers energy from the orbit of the binary to internal modes of the stars, accelerating the gravitational-wave driven inspiral. We model the impact of this instability on the phasing of the gravitational wave signal using three parameters per star: An overall amplitude, a saturation frequency, and a spectral index. Incorporating these additional parameters, we compute the Bayes factor (lnB!pgpg) comparing our p-g model to a standard one. We find that the observed signal is consistent with waveform models that neglect p-g effects, with lnB!pgpg=…

Towards asteroseismology of core-collapse supernovae with gravitational-wave observations - II. Spacetime perturbations

Improvements in ground-based, advanced gravitational wave (GW) detectors may allow in the near future to observe the GW signal of a nearby core-collapse supernova. For the most common type of progenitors, likely with slowly rotating cores, the dominant GW emission mechanisms are the post-bounce oscillations of the proto-neutron star (PNS) before the explosion. We present a new procedure to compute the eigenmodes of the system formed by the PNS and the stalled accretion shock in general relativity including spacetime perturbations. The new method improves on previous results by accounting for perturbations of both the lapse function and the conformal factor. We apply our analysis to two nume…

The advanced Virgo longitudinal control system for the O2 observing run

Following a successful period of data-taking between 2006 and 2011, the Virgo gravitational-wave detector was taken offline for a major upgrade. The changes made to the instrument significantly increased the complexity of the control systems and meant that an extended period of commissioning was required to reach a sensitivity appropriate for science data-taking. This commissioning period was completed in July of 2017 and the second-generation Advanced Virgo detector went on to join the Advanced LIGO detectors in the O2 science run in August of the same year. The upgraded detector was approximately twice as sensitive to binary neutron star mergers as the first-generation instrument. During …

Quasistationary solutions of scalar fields around accreting black holes

Massive scalar fields can form long-lived configurations around black holes. These configurations, dubbed quasi-bound states, have been studied both in the linear and nonlinear regimes. In this paper we show that quasi-bound states can form in a dynamical scenario in which the mass of the black hole grows significantly due to the capture of infalling matter. We solve the Klein-Gordon equation numerically in spherical symmetry, mimicking the evolution of the spacetime through a sequence of analytic Schwarzschild black hole solutions of increasing mass. It is found that the frequency of oscillation of the quasi-bound states decreases as the mass of the black hole increases. In addition, accre…

Magnetoelastic oscillations of neutron stars with dipolar magnetic fields

A new general relativistic magnetohydrodynamics code for dynamical spacetimes

We present a new numerical code which solves the general relativistic magneto-hydrodynamics (GRMHD) equations coupled to the Einstein equations for the evolution of a dynamical spacetime within the conformally-flat approximation. This code has been developed with the main objective of studying astrophysical scenarios in which both, high magnetic fields and strong gravitational fields appear, such as the magneto-rotational collapse of stellar cores, the collapsar model of GRBs, and the evolution of neutron stars. The code is based on an existing and thoroughly tested purely hydrodynamics code and on its extension to accommodate weakly magnetized fluids (passive magnetic field approximation).…

Inference of proto-neutron star properties from gravitational-wave data in core-collapse supernovae

The eventual detection of gravitational waves from core-collapse supernovae (CCSN) will help improve our current understanding of the explosion mechanism of massive stars. The stochastic nature of the late post-bounce gravitational wave signal due to the non-linear dynamics of the matter involved and the large number of degrees of freedom of the phenomenon make the source parameter inference problem very challenging. In this paper we take a step towards that goal and present a parameter estimation approach which is based on the gravitational waves associated with oscillations of proto-neutron stars (PNS). Numerical simulations of CCSN have shown that buoyancy-driven g-modes are responsible …

Relativistic simulations of rotational core collapse : II. Collapse dynamics and gravitational radiation

We have performed hydrodynamic simulations of relativistic rotational supernova core collapse in axisymmetry and have computed the gravitational radiation emitted by such an event. Details of the methodology and of the numerical code have been given in an accompanying paper. We have simulated the evolution of 26 models in both Newtonian and relativistic gravity. Our simulations show that the three different types of rotational supernova core collapse and gravitational waveforms identified in previous Newtonian simulations (regular collapse, multiple bounce collapse, and rapid collapse) are also present in relativistic gravity. However, rotational core collapse with multiple bounces is only …

Quasi-periodic accretion and gravitational waves from oscillating "toroidal neutron stars" around a Schwarzschild black hole

We present general relativistic hydrodynamics simulations of constant specific angular momentum tori orbiting a Schwarzschild black hole. These tori are expected to form as a result of stellar gravitational collapse, binary neutron star merger or disruption, can reach very high rest-mass densities and behave effectively as neutron stars but with a toroidal topology (i.e. ``toroidal neutron stars''). Our attention is here focussed on the dynamical response of these objects to axisymmetric perturbations. We show that, upon the introduction of perturbations, these systems either become unstable to the runaway instability or exhibit a regular oscillatory behaviour resulting in a quasi-periodic …

Boson stars in Palatini $f(\mathcal{R})$ gravity

We explore equilibrium solutions of spherically symmetric boson stars in the Palatini formulation of $f(\mathcal{R})$ gravity. We account for the modifications introduced in the gravitational sector by using a recently established correspondence between modified gravity with scalar matter and general relativity with modified scalar matter. We focus on the quadratic theory $f(\mathcal{R})=R+\xi R^2$ and compare its solutions with those found in general relativity, exploring both positive and negative values of the coupling parameter $\xi$. As matter source, a complex, massive scalar field with and without self-interaction terms is considered. Our results show that the existence curves of bos…

Dynamics of magnetized relativistic tori oscillating around black holes

We present a numerical study of the dynamics of magnetized, relativistic, non-self-gravitating, axisymmetric tori orbiting in the background spacetimes of Schwarzschild and Kerr black holes. The initial models have a constant specific angular momentum and are built with a non-zero toroidal magnetic field component, for which equilibrium configurations have recently been obtained. In this work we extend our previous investigations which dealt with purely hydrodynamical thick discs, and study the dynamics of magnetized tori subject to perturbations which, for the values of the magnetic field strength considered here, trigger quasi-periodic oscillations lasting for tens of orbital periods. Ove…

Hydromagnetic instabilities and magnetic field amplification in core collapse supernovae

Some of the most violent events in the universe, the gamma ray burst, could be related to the gravitational collapse of massive stellar cores. The recent association of long GRBs to some class of type Ic supernova seems to support this view. In such scenario fast rotation, strong magnetic fields and general relativistic effects are key ingredients. It is thus important to understand the mechanism that amplifies the magnetic field under that conditions. I present global simulations of the magneto-rotational collapse of stellar cores in general relativity and semi-global simulations of hydromagnetic instabilities under core collapse conditions. I discuss effect of the magneto-rotational insta…

Total-variation-based methods for gravitational wave denoising

We describe new methods for denoising and detection of gravitational waves embedded in additive Gaussian noise. The methods are based on Total Variation denoising algorithms. These algorithms, which do not need any a priori information about the signals, have been originally developed and fully tested in the context of image processing. To illustrate the capabilities of our methods we apply them to two different types of numerically-simulated gravitational wave signals, namely bursts produced from the core collapse of rotating stars and waveforms from binary black hole mergers. We explore the parameter space of the methods to find the set of values best suited for denoising gravitational wa…

Numerical Hydrodynamics and Magnetohydrodynamics in General Relativity.

This article presents a comprehensive overview of numerical hydrodynamics and magneto-hydrodynamics (MHD) in general relativity. Some significant additions have been incorporated with respect to the previous two versions of this review (2000, 2003), most notably the coverage of general-relativistic MHD, a field in which remarkable activity and progress has occurred in the last few years. Correspondingly, the discussion of astrophysical simulations in general-relativistic hydrodynamics is enlarged to account for recent relevant advances, while those dealing with general-relativistic MHD are amply covered in this review for the first time. The basic outline of this article is nevertheless sim…

Head-on collisions and orbital mergers of Proca stars

Proca stars are self-gravitating Bose-Einstein condensates obtained as numerical stationary solutions of the Einstein-(complex)-Proca system. These solitonic can be both stable and form dynamically from generic initial data by the mechanism of gravitational cooling. In this paper we further explore the dynamical properties of these solitonic objects by performing both head-on collisions and orbital mergers of equal mass Proca stars, using fully non-linear numerical evolutions. For the head-on collisions, we show that the end point and the gravitational waveform from these collisions depends on the compactness of the Proca star. Proca stars with sufficiently small compactness collide leaving…

Exploring gravitational-wave detection and parameter inference using deep learning methods

The data that support the findings of this study are openly available at the following URL/DOI: https://arxiv.org/abs/2011.10425.

Dynamics of oscillating magnetized relativistic tori around a Schwarzschild black hole

We present a comprehensive numerical study of the dynamics of magnetized relativistic axisymmetric tori orbiting in the background spacetime of a Schwarzschild black hole. The tori are modeled as having a purely toroidal magnetic field and a constant distribution of the specific angular momentum. Following previous investigations of tori in a purely hydrodynamical context, the dynamics of these objects has been studied upon the introduction of a perturbation which, for the values of the magnetic field considered here, triggers quasi-periodic oscillations (QPOs) lasting tens of orbital periods. As in the hydrodynamical case, the spectral distribution of the eigenfrequencies shows the presenc…

All-sky search for continuous gravitational waves from isolated neutron stars using Advanced LIGO O2 data

We present results of an all-sky search for continuous gravitational waves (CWs), which can be produced by fast-spinning neutron stars with an asymmetry around their rotation axis, using data from the second observing run of the Advanced LIGO detectors. We employ three different semi-coherent methods ($\textit{FrequencyHough}$, $\textit{SkyHough}$, and $\textit{Time-Domain $\mathcal{F}$-statistic}$) to search in a gravitational-wave frequency band from 20 to 1922 Hz and a first frequency derivative from $-1\times10^{-8}$ to $2\times10^{-9}$ Hz/s. None of these searches has found clear evidence for a CW signal, so we present upper limits on the gravitational-wave strain amplitude $h_0$ (the …

GW190412: Observation of a binary-black-hole coalescence with asymmetric masses

LIGO Scientific Collaboration and Virgo Collaboration: et al.

Time-domain effective-one-body gravitational waveforms for coalescing compact binaries with nonprecessing spins, tides and self-spin effects

We present TEOBResumS, a new effective-one-body (EOB) waveform model for nonprecessing (spin-aligned) and tidally interacting compact binaries.Spin-orbit and spin-spin effects are blended together by making use of the concept of centrifugal EOB radius. The point-mass sector through merger and ringdown is informed by numerical relativity (NR) simulations of binary black holes (BBH) computed with the SpEC and BAM codes. An improved, NR-based phenomenological description of the postmerger waveform is developed.The tidal sector of TEOBResumS describes the dynamics of neutron star binaries up to merger and incorporates a resummed attractive potential motivated by recent advances in the post-Newt…

Magnetized discs and photon rings around Yukawa-like black holes

We present stationary solutions of geometrically thick discs (or tori) endowed with a self-consistent toroidal magnetic field distribution surrounding a non-rotating black hole in an analytical, static, spherically-symmetric $f(R)$-gravity background. These $f(R)$-gravity models introduce a Yukawa-like modification to the Newtonian potential, encoded in a single parameter $\delta$ which controls the strength of the modified potential and whose specific values affect the disc configurations when compared to the general relativistic case. Our models span different magnetic field strengths, from purely hydrodynamical discs to highly magnetized tori. The characteristics of the solutions are ide…

Comparison between the fCCZ4 and BSSN formulations of Einstein equations in spherical polar coordinates

Recently, we generalized a covariant and conformal version of the Z4 system of the Einstein equations using a reference metric approach, that we denote as fCCZ4. We successfully implemented and tested this approach in a 1D code that uses spherical coordinates and assumes spherical symmetry, obtaining from one to three orders of magnitude reduction of the Hamiltonian constraint violations with respect to the BSSN formulation in tests involving neutron star spacetimes. In this work, we show preliminary results obtained with the 3D implementation of the fCCZ4 formulation in a fully 3D code using spherical polar coordinates.

Magnetised Polish doughnuts revisited

We discuss a procedure to build new sequences of magnetised, equilibrium tori around Kerr black holes which combines two approaches previously considered in the literature. For simplicity we assume that the test-fluid approximation holds, and hence we neglect the self-gravity of the fluid. The models are built assuming a particular form of the angular momentum distribution from which the location and morphology of equipotential surfaces can be computed. This ansatz includes, in particular, the constant angular momentum case originally employed in the construction of thick tori - or Polish doughnuts - and it has already been used to build equilibrium sequences of purely hydrodynamical models…

Explosion and Final State of an Unstable Reissner-Nordström Black Hole

A Reissner-Nordstr\"om black hole (BH) is superradiantly unstable against spherical perturbations of a charged scalar field, enclosed in a cavity, with frequency lower than a critical value. We use numerical relativity techniques to follow the development of this unstable system -- dubbed a charged BH bomb -- into the non-linear regime, solving the full Einstein--Maxwell--Klein-Gordon equations, in spherical symmetry. We show that: $i)$ the process stops before all the charge is extracted from the BH; $ii)$ the system settles down into a hairy BH: a charged horizon in equilibrium with a scalar field condensate, whose phase is oscillating at the (final) critical frequency. For low scalar fie…

Gravitational-wave Constraints on the Equatorial Ellipticity of Millisecond Pulsars

We present a search for continuous gravitational waves from five radio pulsars, comprising three recycled pulsars (PSR J0437-4715, PSR J0711-6830, and PSR J0737-3039A) and two young pulsars: the Crab pulsar (J0534+2200) and the Vela pulsar (J0835-4510). We use data from the third observing run of Advanced LIGO and Virgo combined with data from their first and second observing runs. For the first time, we are able to match (for PSR J0437-4715) or surpass (for PSR J0711-6830) the indirect limits on gravitational-wave emission from recycled pulsars inferred from their observed spin-downs, and constrain their equatorial ellipticities to be less than 10-8. For each of the five pulsars, we perfor…

Relativistic Low Angular Momentum Accretion: Long Time Evolution of Hydrodynamical Inviscid Flows

We investigate relativistic low angular momentum accretion of inviscid perfect fluid onto a Schwarzschild black hole. The simulations are performed with a general-relativistic, high-resolution (second-order), shock-capturing, hydrodynamical numerical code. We use horizon-penetrating Eddington-Finkelstein coordinates to remove inaccuracies in regions of strong gravity near the black hole horizon and show the expected convergence of the code with the Michel solution and stationary Fishbone-Moncrief toroids. We recover, in the framework of relativistic hydrodynamics, the qualitative behavior known from previous Newtonian studies that used a Bondi background flow in a pseudo-relativistic gravit…

Combining spectral and shock-capturing methods: A new numerical approach for 3D relativistic core collapse simulations

We present a new three-dimensional general relativistic hydrodynamics code which is intended for simulations of stellar core collapse to a neutron star, as well as pulsations and instabilities of rotating relativistic stars. Contrary to the common approach followed in most existing three-dimensional numerical relativity codes which are based in Cartesian coordinates, in this code both the metric and the hydrodynamics equations are formulated and solved numerically using spherical polar coordinates. A distinctive feature of this new code is the combination of two types of accurate numerical schemes specifically designed to solve each system of equations. More precisely, the code uses spectra…

Stationary models of magnetized viscous tori around a Schwarzschild black hole

We present stationary solutions of magnetized, viscous thick accretion disks around a Schwarzschild black hole. We assume that the tori are not self-gravitating, are endowed with a toroidal magnetic field, and obey a constant angular momentum law. Our study focuses on the role of the black hole curvature in the shear viscosity tensor and in their potential combined effect on the stationary solutions. Those are built in the framework of a causality-preserving, second-order gradient expansion scheme of relativistic hydrodynamics in the Eckart frame description which gives rise to hyperbolic equations of motion. The stationary models are constructed by numerically solving the general relativis…

Fallback accretion onto magnetized neutron stars and the hidden magnetic field model

The observation of several neutron stars with relatively low values of the surface magnetic field found in supernova remnants has led in recent years to controversial interpretations. A possible explanation is the slow rotation of the proto-neutron star at birth which is unable to amplify its magnetic field to typical pulsar levels. An alternative possibility, the hidden magnetic field scenario, seems to be favoured over the previous one due to the observation of three low magnetic field magnetars. This scenario considers the accretion of the fallback of the supernova debris onto the neutron star as the responsible for the observed low magnetic field. In this work, we have studied under whi…

Numerical 3+1 general relativistic magnetohydrodynamics: a local characteristic approach

We present a general procedure to solve numerically the general relativistic magnetohydrodynamics (GRMHD) equations within the framework of the 3+1 formalism. The work reported here extends our previous investigation in general relativistic hydrodynamics (Banyuls et al. 1997) where magnetic fields were not considered. The GRMHD equations are written in conservative form to exploit their hyperbolic character in the solution procedure. All theoretical ingredients necessary to build up high-resolution shock-capturing schemes based on the solution of local Riemann problems (i.e. Godunov-type schemes) are described. In particular, we use a renormalized set of regular eigenvectors of the flux Jac…

GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral

On August 17, 2017 at 12-41:04 UTC the Advanced LIGO and Advanced Virgo gravitational-wave detectors made their first observation of a binary neutron star inspiral. The signal, GW170817, was detected with a combined signal-to-noise ratio of 32.4 and a false-alarm-rate estimate of less than one per 8.0×104 years. We infer the component masses of the binary to be between 0.86 and 2.26 M, in agreement with masses of known neutron stars. Restricting the component spins to the range inferred in binary neutron stars, we find the component masses to be in the range 1.17-1.60 M, with the total mass of the system 2.74-0.01+0.04M. The source was localized within a sky region of 28 deg2 (90% probabili…

Magneto-elastic torsional oscillations of magnetars

We extend a general-relativistic ideal magneto-hydrodynamical code to include the effects of elasticity. Using this numerical tool we analyse the magneto-elastic oscillations of highly magnetised neutron stars (magnetars). In simulations without magnetic field we are able to recover the purely crustal shear oscillations within an accuracy of about a few per cent. For dipole magnetic fields between 5 x 10^13 and 10^15 G the Alfv\'en oscillations become modified substantially by the presence of the crust. Those quasi-periodic oscillations (QPOs) split into three families: Lower QPOs near the equator, Edge QPOs related to the last open field line and Upper QPOs at larger distance from the equa…

Numerical Hydrodynamics in General Relativity

The current status of numerical solutions for the equations of ideal general relativistic hydrodynamics is reviewed. With respect to an earlier version of the article the present update provides additional information on numerical schemes and extends the discussion of astrophysical simulations in general relativistic hydrodynamics. Different formulations of the equations are presented, with special mention of conservative and hyperbolic formulations well-adapted to advanced numerical methods. A large sample of available numerical schemes is discussed, paying particular attention to solution procedures based on schemes exploiting the characteristic structure of the equations through lineariz…

A Gravitational-wave Measurement of the Hubble Constant Following the Second Observing Run of Advanced LIGO and Virgo

This paper presents the gravitational-wave measurement of the Hubble constant (H 0) using the detections from the first and second observing runs of the Advanced LIGO and Virgo detector network. The presence of the transient electromagnetic counterpart of the binary neutron star GW170817 led to the first standard-siren measurement of H 0. Here we additionally use binary black hole detections in conjunction with galaxy catalogs and report a joint measurement. Our updated measurement is H 0 = km s-1 Mpc-1 (68.3% of the highest density posterior interval with a flat-in-log prior) which is an improvement by a factor of 1.04 (about 4%) over the GW170817-only value of km s-1 Mpc-1. A significant …

A "horizon adapted" approach to the study of relativistic accretion flows onto rotating black holes

We present a new geometrical approach to the study of accretion flows onto rotating (Kerr) black holes. Instead of Boyer-Lindquist coordinates, the standard choice in all existing numerical simulations in the literature, we employ the simplest example of a horizon adapted coordinate system, the Kerr-Schild coordinates. This choice eliminates boundary ambiguities and unphysical divergent behavior at the event horizon. Computations of Bondi-Hoyle accretion onto extreme Kerr black holes, performed here for the first time, demonstrate the key advantages of this procedure. We argue it offers the best approach to the numerical study of the, observationally, increasingly more accesible relativisti…

Fully Covariant and Conformal Formulation of the Z4 System Compared to the BSSN Formulation in Spherical Symmetry

We have generalized a covariant and conformal version of the Z4 system of the Einstein equations by adopting a reference metric approach, that we denote as fCCZ4, well suited for curvilinear as well as Cartesian coordinates. We implement this formalism in spherical polar coordinates under the assumption of spherical symmetry using a partially-implicit Runge-Kutta (PIRK) method, without using any regularization scheme, and show that our code can evolve both vacuum and non-vacuum spacetimes without encountering instabilities. We have performed several tests and compared the Hamiltonian constraint violations of the fCCZ4 system, for different choices of certain free parameters, with these of B…

Classification methods for noise transients in advanced gravitational-wave detectors II: performance tests on Advanced LIGO data

The data taken by the advanced LIGO and Virgo gravitational-wave detectors contains short duration noise transients that limit the significance of astrophysical detections and reduce the duty cycle of the instruments. As the advanced detectors are reaching sensitivity levels that allow for multiple detections of astrophysical gravitational-wave sources it is crucial to achieve a fast and accurate characterization of non-astrophysical transient noise shortly after it occurs in the detectors. Previously we presented three methods for the classification of transient noise sources. They are Principal Component Analysis for Transients (PCAT), Principal Component LALInference Burst (PC-LIB) and W…

Quasistationary solutions of self-gravitating scalar fields around collapsing stars

Recent work has shown that scalar fields around black holes can form long-lived, quasistationary configurations surviving for cosmological time scales. Scalar fields thus cannot be discarded as viable candidates for dark matter halo models in galaxies around central supermassive black holes (SMBHs). One hypothesized formation scenario of most SMBHs at high redshift is the gravitational collapse of supermassive stars (SMSs) with masses of $\ensuremath{\sim}{10}^{5}\text{ }\text{ }{\mathrm{M}}_{\ensuremath{\bigodot}}$. Any such scalar field configurations must survive the gravitational collapse of a SMS in order to be a viable model of physical reality. To check for the postcollapse survival …

Non-linear evolutions of magnetized thick discs around black holes: dependence on the initial data

We build equilibrium solutions of magnetised thick discs around a highly spinning Kerr black hole and evolve these initial data up to a final time of about 100 orbital periods. The numerical simulations reported in this paper solve the general relativistic magnetohydrodynamics equations using the BHAC code and are performed in axisymmetry. Our study assumes non-self-gravitating, polytropic, constant angular momentum discs endowed with a purely toroidal magnetic field. In order to build the initial data we consider three approaches, two of which incorporate the magnetic field in a self-consistent way and a third approach in which the magnetic field is included as a perturbation on to an othe…

Dynamical bar-mode instability in spinning bosonic stars

Spinning bosonic stars (SBSs) can form from the gravitational collapse of a dilute cloud of scalar/Proca particles with nonzero angular momentum, via gravitational cooling. The scalar stars are, however, transient due to a nonaxisymmetric instability which triggers the loss of angular momentum. By contrast, no such instability was observed for the fundamental ( m = 1 ) Proca stars. In [N. Sanchis-Gual et al., Phys. Rev. Lett. 123, 221101 (2019)] we tentatively related the different stability properties to the different toroidal/spheroidal morphology of the scalar/Proca models. Here, we continue this investigation, using three-dimensional numerical-relativity simulations of the Einstein-(mas…

All-sky search for short gravitational-wave bursts in the second Advanced LIGO and Advanced Virgo run

We present the results of a search for short-duration gravitational-wave transients in the data from the second observing run of Advanced LIGO and Advanced Virgo. We search for gravitational-wave transients with a duration of milliseconds to approximately one second in the 32-4096 Hz frequency band with minimal assumptions about the signal properties, thus targeting a wide variety of sources. We also perform a matched-filter search for gravitational-wave transients from cosmic string cusps for which the waveform is well-modeled. The unmodeled search detected gravitational waves from several binary black hole mergers which have been identified by previous analyses. No other significant event…

Gravitational waves from galaxy encounters

We discuss the emission of gravitational radiation produced in encounters of dark matter galactic halos. To this aim we perform a number of numerical simulations of typical galaxy mergers, computing the associated gravitational radiation waveforms as well as the energy released in the processes. Our simulations yield dimensionless gravitational wave amplitudes of the order of $10^{-13}$ and gravitational wave frequencies of the order of $10^{-16}$ Hz, when the galaxies are located at a distance of 10 Mpc. These values are of the same order as those arising in the gravitational radiation originated by strong variations of the gravitational field in the early Universe, and therefore, such gra…

Self-gravitating magnetized tori around black holes in general relativity

We investigate stationary, self-gravitating, magnetised disks (or tori) around black holes. The models are obtained by numerically solving the coupled system of the Einstein equations and the equations of ideal general-relativistic magnetohydrodynamics. The mathematical formulation and numerical aspects of our approach are similar to those reported in previous works modeling stationary self-gravitating perfect-fluid tori, but the inclusion of magnetic fields represents a new ingredient. Following previous studies of purely hydrodynamical configurations, we construct our models assuming Keplerian rotation in the disks and both spinning and spinless black holes. We focus on the case of a toro…

Search for Subsolar-Mass Ultracompact Binaries in Advanced LIGO's First Observing Run

We present the first Advanced LIGO and Advanced Virgo search for ultracompact binary systems with component masses between 0.2 $M_\odot$ - 1.0 $M_\odot$ using data taken between September 12, 2015 and January 19, 2016. We find no viable gravitational wave candidates. Our null result constrains the coalescence rate of monochromatic (delta function) distributions of non-spinning (0.2 $M_\odot$, 0.2 $M_\odot$) ultracompact binaries to be less than $1.0 \times 10^6 \text{Gpc}^{-3} \text{yr}^{-1}$ and the coalescence rate of a similar distribution of (1.0 $M_\odot$, 1.0 $M_\odot$) ultracompact binaries to be less than $1.9 \times 10^4 \text{Gpc}^{-3} \text{yr}^{-1}$ (at 90 percent confidence). N…

Spontaneous Scalarization of Charged Black Holes

Extended scalar-tensor-Gauss-Bonnet (eSTGB) gravity has been recently argued to exhibit spontaneous scalarisation of vacuum black holes (BHs). A similar phenomenon can be expected in a larger class of models, which includes e.g. Einstein-Maxwell-scalar (EMS) models, where spontaneous scalarisation of electrovacuum BHs should occur. EMS models have no higher curvature corrections, a technical simplification over eSTGB models that allows us to investigate, fully non-linearly, BH scalarisation in two novel directions. Firstly, numerical simulations in spherical symmetry show, dynamically, that Reissner-Nordstr\"om (RN) BHs evolve into a perturbatively stable scalarised BH. Secondly, we compute…

THREE-DIMENSIONAL RELATIVISTIC SIMULATIONS OF ROTATING NEUTRON-STAR COLLAPSE TO A KERR BLACK HOLE

We present a new three-dimensional fully general-relativistic hydrodynamics code using high-resolution shock-capturing techniques and a conformal traceless formulation of the Einstein equations. Besides presenting a thorough set of tests which the code passes with very high accuracy, we discuss its application to the study of the gravitational collapse of uniformly rotating neutron stars to Kerr black holes. The initial stellar models are modeled as relativistic polytropes which are either secularly or dynamically unstable and with angular velocities which range from slow rotation to the mass-shedding limit. We investigate the gravitational collapse by carefully studying not only the dynami…

Interstellar travel guide: Chronicles of a violent universe

In contrast to observations from our local viewpoint, the universe is extremely violent. Explosions of massive stars, X-ray and gamma-ray binary stars, and active galaxies are related to the presence of compact objects such as neutron stars or black holes. Processes such as matter accretion or the formation and release of relativistic plasma jets are responsible for the generation of large amounts of high-energy radiation, caused by the acceleration of elementary particles, and have a significant impact in their surroundings. Moreover, they are related to the generation of gravitational waves: vibrations in spacetime produced by the acceleration of compact massive objects that produce stron…

Constraining properties of high-density matter in neutron stars with magneto-elastic oscillations

We discuss torsional oscillations of highly magnetised neutron stars (magnetars) using two-dimensional, magneto-elastic-hydrodynamical simulations. Our model is able to explain both the low- and high-frequency quasi-periodic oscillations (QPOs) observed in magnetars. The analysis of these oscillations provides constraints on the breakout magnetic-field strength, on the fundamental QPO frequency, and on the frequency of a particularly excited overtone. More importantly, we show how to use this information to generically constraint properties of high-density matter in neutron stars, employing Bayesian analysis. In spite of current uncertainties and computational approximations, our model-depe…

Dynamical formation and stability of fermion-boson stars

Gravitationally bound structures composed by fermions and scalar particles known as fermion-boson stars are regular and static configurations obtained by solving the coupled Einstein-Klein-Gordon-Euler (EKGE) system. In this work, we discuss one possible scenario through which these fermion-boson stars may form by solving numerically the EKGE system under the simplifying assumption of spherical symmetry. Our initial configurations assume an already existing neutron star surrounded by an accreting cloud of a massive and complex scalar field. The results of our simulations show that once part of the initial scalar field is expelled via gravitational cooling the system gradually oscillates aro…

Magnetised tori with magnetic polarisation around Kerr black holes: variable angular momentum discs

Analytical models of magnetised, geometrically thick discs are relevant to understand the physical conditions of plasma around compact objects and to explore its emitting properties. This has become increasingly important in recent years in the light of the Event Horizon Telescope observations of Sgr A* and M87. Models of thick discs around black holes usually consider constant angular momentum distributions and do not take into account the magnetic response of the fluid to applied magnetic fields. We present a generalisation of our previous work on stationary models of magnetised accretion discs with magnetic polarisation (Pimentel et al. 2018). This extension is achieved by accounting for…

Relativistic simulations of rotational core collapse : I. Methods, initial models, and code tests

We describe an axisymmetric general relativistic code for rotational core collapse. The code evolves the coupled system of metric and fluid equations using the ADM 3+1 formalism and a conformally flat metric approximation of the Einstein equations. The relativistic hydrodynamics equations are formulated as a first-order flux-conservative hyperbolic system and are integrated using high-resolution shock-capturing schemes based on Riemann solvers. We assess the quality of the conformally flat metric approximation for relativistic core collapse and present a comprehensive set of tests which the code successfully passed. The tests include relativistic shock tubes, the preservation of the rotatio…

Dynamical formation of Proca stars and quasistationary solitonic objects

We perform fully non-linear numerical simulations within the spherically symmetric Einstein-(complex)Proca system. Starting with Proca field distributions that obey the Hamiltonian, momentum and Gaussian constraints, we show that the self-gravity of the system induces the formation of compact objects, which, for appropriate initial conditions, asymptotically approach stationary soliton-like solutions known as Proca stars. The excess energy of the system is dissipated by the mechanism of \textit{gravitational cooling} in analogy to what occurs in the dynamical formation of scalar boson stars. We investigate the dependence of this process on the phase difference between the real and imaginary…

Towards relativistic simulations of magneto-rotational core collapse

We present a new general relativistic hydrodynamics code specifically designed to study magneto-rotational, relativistic, stellar core collapse. The code is an extension of an existing (and thoroughly tested) hydrodynamics code, which has been applied in the recent past to study relativistic rotational core collapse. It is based on the conformally-flat approximation of Einstein's field equations and conservative formulations for the magneto-hydrodynamics equations. As a first step towards magneto-rotational core collapse simulations the code assumes a passive (test) magnetic field. The paper is focused on the description of the technical details of the numerical implementation, with emphasi…

Gravitational waves from oscillating accretion tori: Comparison between different approaches

Quasi-periodic oscillations of high density thick accretion disks orbiting a Schwarzschild black hole have been recently addressed as interesting sources of gravitational waves. The aim of this paper is to compare the gravitational waveforms emitted from these sources when computed using (variations of) the standard quadrupole formula and gauge-invariant metric perturbation theory. To this goal we evolve representative disk models using an existing general relativistic hydrodynamics code which has been previously employed in investigations of such astrophysical systems. Two are the main results of this work: First, for stable and marginally stable disks, no excitation of the black hole quas…

Modulating the magnetosphere of magnetars by internal magneto-elastic oscillations

We couple internal torsional, magneto-elastic oscillations of highly magnetized neutron stars (magnetars) to their magnetospheres. The corresponding axisymmetric perturbations of the external magnetic field configuration evolve as a sequence of linear, force-free equilibria that are completely determined by the background magnetic field configuration and by the perturbations of the magnetic field at the surface. The perturbations are obtained from simulations of magneto-elastic oscillations in the interior of the magnetar. While such oscillations can excite travelling Alfv\'en waves in the exterior of the star only in a very limited region close to the poles, they still modulate the near ma…

General Relativistic Simulations of the Collapsar Scenario

We are exploring the viability of the collapsar model for long-soft gamma-ray bursts. For this we perform state-of-the-art general relativistic hydrodynamic simulations in a dynamically evolving space-time with the CoCoNuT code. We start from massive low metallicity stellar models evolved up to core gravitational instability, and then follow the subsequent evolution until the system collapses forming a compact remnant. A preliminary study of the collapse outcome is performed by varying the typical parameters of the scenario, such as the initial stellar mass, metallicity, and rotational profile of the stellar progenitor. 1D models (without rotation) have been used to test our newly developed…

AMR simulations of the low bar-mode instability of neutron stars

It has been recently argued through numerical work that rotating stars with a high degree of differential rotation are dynamically unstable against bar-mode deformation, even for values of the ratio of rotational kinetic energy to gravitational potential energy as low as O(0.01). This may have implications for gravitational wave astronomy in high-frequency sources such as core collapse supernovae. In this paper we present high-resolution simulations, performed with an adaptive mesh refinement hydrodynamics code, of such low T/|W| bar-mode instability. The complex morphological features involved in the nonlinear dynamics of the instability are revealed in our simulations, which show that the…

Morphology and Dynamics of Relativistic Jets

We present a comprehensive analysis of the morphology and dynamics of relativistic pressure-matched axisymmetric jets. The numerical simulations have been carried out with a high-resolution shock-capturing hydrocode based on an approximate relativistic Riemann solver derived from the spectral decomposition of the Jacobian matrices of relativistic hydrodynamics. We discuss the dependence of the jet morphology on several parameters, paying special attention to the relativistic effects caused by high Lorentz factors and large internal energies of the beam flow. The parameter space of our analysis is spanned by the ratio of the beam and ambient medium rest mass density (η), the beam Mach number…

Application of dictionary learning to denoise LIGO’s blip noise transients

Data streams of gravitational-wave detectors are polluted by transient noise features, or ``glitches,'' of instrumental and environmental origin. In this work we investigate the use of total variation methods and learned dictionaries to mitigate the effect of those transients in the data. We focus on a specific type of transient, ``blip" glitches, as this is the most common type of glitch present in the LIGO detectors and their waveforms are easy to identify. We randomly select 100 blip glitches scattered in the data from advanced LIGO's O1 run, as provided by the citizen-science project Gravity Spy. Our results show that dictionary-learning methods are a valid approach to model and subtrac…

CFC+: Improved dynamics and gravitational waveforms from relativistic core collapse simulations

Core collapse supernovae are a promising source of detectable gravitational waves. Most of the existing (multidimensional) numerical simulations of core collapse in general relativity have been done using approximations of the Einstein field equations. As recently shown by Dimmelmeier et al (2002a,b), one of the most interesting such approximation is the so-called conformal flatness condition (CFC) of Isenberg, Wilson and Mathews. Building on this previous work we present here new results from numerical simulations of relativistic rotational core collapse in axisymmetry, aiming at improving the dynamics and the gravitational waveforms. The computer code used for these simulations evolves th…

Numerical evolutions of spherical Proca stars

Vector boson stars, or $\textit{Proca stars}$, have been recently obtained as fully non-linear numerical solutions of the Einstein-(complex)-Proca system. These are self-gravitating, everywhere non-singular, horizonless Bose-Einstein condensates of a massive vector field, which resemble in many ways, but not all, their scalar cousins, the well-known (scalar) $\textit{boson stars}$. In this paper we report fully-non linear numerical evolutions of Proca stars, focusing on the spherically symmetric case, with the goal of assessing their stability and the end-point of the evolution of the unstable stars. Previous results from linear perturbation theory indicate the separation between stable and…

GRAVITATIONAL WAVE SIGNATURES IN BLACK HOLE FORMING CORE COLLAPSE

We present numerical simulations in general relativity of collapsing stellar cores. Our initial model consists of a low metallicity rapidly-rotating progenitor which is evolved in axisymmetry with the latest version of our general relativistic code CoCoNuT, which allows for black hole formation and includes the effects of a microphysical equation of state (LS220) and a neutrino leakage scheme to account for radiative losses. The motivation of our study is to analyze in detail the emission of gravitational waves in the collapsar scenario of long gamma-ray bursts. Our simulations show that the phase during which the proto-neutron star (PNS) survives before ultimately collapsing to a black hol…

Lensing and dynamics of ultracompact bosonic stars

Spherically symmetric bosonic stars are one of the few examples of gravitating solitons that are known to form dynamically, via a classical process of (incomplete) gravitational collapse. As stationary solutions of the Einstein--Klein-Gordon or the Einstein--Proca theory, bosonic stars may also become sufficiently compact to develop light rings and hence mimic, in principle, gravitational-wave observational signatures of black holes (BHs). In this paper, we discuss how these horizonless ultra-compact objects (UCOs) are actually distinct from BHs, both phenomenologically and dynamically. In the electromagnetic channel, the light ring associated phenomenology reveals remarkable lensing patter…

On the dynamical bar-mode instability in spinning bosonic stars

Spinning bosonic stars (SBSs) can form from the gravitational collapse of a dilute cloud of scalar/Proca particles with non-zero angular momentum. In a recent work we found that the scalar stars are transient due to a non-axisymmetric instability which triggers the loss of angular momentum. We further study the dynamical formation of SBSs using 3-dimensional numerical-relativity simulations of the Einstein-(massive, complex)Klein-Gordon system and of the Einstein-(complex)Proca system. We incorporate a quartic self-interaction potential in the scalar case to gauge its effect on the instability; we investigate (m=2) Proca stars to assess their stability; we attempt to relate the instability …

Towards asteroseismology of core-collapse supernovae with gravitational wave observations – II. Inclusion of space–time perturbations

Improvements in ground-based, advanced gravitational wave (GW) detectors may allow in the near future to observe the GW signal of a nearby core-collapse supernova. For the most common type of progenitors, likely with slowly rotating cores, the dominant GW emission mechanisms are the post-bounce oscillations of the proto-neutron star (PNS) before the explosion. We present a new procedure to compute the eigenmodes of the system formed by the PNS and the stalled accretion shock in general relativity including spacetime perturbations. The new method improves on previous results by accounting for perturbations of both the lapse function and the conformal factor. We apply our analysis to two nume…

Alfven QPOs in magnetars in the anelastic approximation

We perform two-dimensional simulations of Alfven oscillations in magnetars, modeled as relativistic stars with a dipolar magnetic field. We use the anelastic approximation to general relativistic magnetohydrodynamics, which allows for an effective suppression of fluid modes and an accurate description of Alfven waves. In addition, we compute Alfven oscillation frequencies along individual magnetic field lines with a semi-analytic approach, employing a short-wavelength approximation. Our main findings are as follows: a) we confirm the existence of two families of quasi-periodic oscillations (QPOs), with harmonics at integer multiples of the fundamental frequency, as was found in the linear s…

Assessment of a high-resolution central scheme for the solution of the relativistic hydrodynamics equations

We assess the suitability of a recent high-resolution central scheme developed by Kurganov & Tadmor (2000) for the solution of the relativistic hydrodynamics equations. The novelty of this approach relies on the absence of Riemann solvers in the solution procedure. The computations we present are performed in one and two spatial dimensions in Minkowski spacetime. Standard numerical experiments such as shock tubes and the relativistic flat-faced step test are performed. As an astrophysical application the article includes two-dimensional simulations of the propagation of relativistic jets using both Cartesian and cylindrical coordinates. The simulations reported clearly show the capabili…

Fully covariant and conformal formulation of the Z4 system in a reference-metric approach: Comparison with the BSSN formulation in spherical symmetry

We adopt a reference-metric approach to generalize a covariant and conformal version of the Z4 system of the Einstein equations. We refer to the resulting system as ``fully covariant and conformal", or fCCZ4 for short, since it is well suited for curvilinear as well as Cartesian coordinates. We implement this fCCZ4 formalism in spherical polar coordinates under the assumption of spherical symmetry using a partially-implicit Runge-Kutta (PIRK) method and show that our code can evolve both vacuum and non-vacuum spacetimes without encountering instabilities. Our method does not require regularization of the equations to handle coordinate singularities, nor does it depend on constraint-preservi…

Accretion-driven gravitational radiation from nonrotating compact objects. Infalling quadrupolar shells

This paper reports results from numerical simulations of the gravitational radiation emitted from non--rotating compact objects(both neutron stars and Schwarzschild black holes) as a result of the accretion of matter. A hybrid procedure is adopted: we evolve, in axisymmetry, the linearized equations describing metric and fluid perturbations, coupled with a nonlinear hydrodynamics code that calculates the motion of the accreting matter. The initial matter distribution is shaped in the form of extended quadrupolar shells of dust or perfect fluid. Self--gravity and radiation reaction effects of the accreting fluid are neglected. This idealized setup is used to understand the qualitative featur…

Convolutional Neural Networks for the classification of glitches in gravitational-wave data streams

We investigate the use of Convolutional Neural Networks (including the modern ConvNeXt network family) to classify transient noise signals (i.e.~glitches) and gravitational waves in data from the Advanced LIGO detectors. First, we use models with a supervised learning approach, both trained from scratch using the Gravity Spy dataset and employing transfer learning by fine-tuning pre-trained models in this dataset. Second, we also explore a self-supervised approach, pre-training models with automatically generated pseudo-labels. Our findings are very close to existing results for the same dataset, reaching values for the F1 score of 97.18% (94.15%) for the best supervised (self-supervised) m…

Axisymmetric core collapse simulations using characteristic numerical relativity

We present results from axisymmetric stellar core collapse simulations in general relativity. Our hydrodynamics code has proved robust and accurate enough to allow for a detailed analysis of the global dynamics of the collapse. Contrary to traditional approaches based on the 3+1 formulation of the gravitational field equations, our framework uses a foliation based on a family of outgoing light cones, emanating from a regular center, and terminating at future null infinity. Such a coordinate system is well adapted to the study of interesting dynamical spacetimes in relativistic astrophysics such as stellar core collapse and neutron star formation. Perhaps most importantly this procedure allo…

GW170817: Measurements of Neutron Star Radii and Equation of State

On 17 August 2017, the LIGO and Virgo observatories made the first direct detection of gravitational waves from the coalescence of a neutron star binary system. The detection of this gravitational-wave signal, GW170817, offers a novel opportunity to directly probe the properties of matter at the extreme conditions found in the interior of these stars. The initial, minimal-assumption analysis of the LIGO and Virgo data placed constraints on the tidal effects of the coalescing bodies, which were then translated to constraints on neutron star radii. Here, we expand upon previous analyses by working under the hypothesis that both bodies were neutron stars that are described by the same equation…

Dynamics of thick discs around Schwarzschild-de Sitter black holes

We consider the effects of a cosmological constant on the dynamics of constant angular momentum discs orbiting Schwarzschild-de Sitter black holes. The motivation behind this study is to investigate whether the presence of a radial force contrasting the black hole's gravitational attraction can influence the occurrence of the runaway instability, a robust feature of the dynamics of constant angular momentum tori in Schwarzschild and Kerr spacetimes. In addition to the inner cusp near the black hole horizon through which matter can accrete onto the black hole, in fact, a positive cosmological constant introduces also an outer cusp through which matter can leave the torus without accreting on…