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…
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…
On the convexity of Relativistic Hydrodynamics
The relativistic hydrodynamic system of equations for a perfect fluid obeying a causal equation of state is hyperbolic (Anile 1989 {\it Relativistic Fluids and Magneto-Fluids} (Cambridge: Cambridge University Press)). In this report, we derive the conditions for this system to be convex in terms of the fundamental derivative of the equation of state (Menikoff and Plohr 1989 {\it Rev. Mod. Phys.} {\bf 61} 75). The classical limit is recovered.
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…
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…
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…
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…
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…
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…
Partially Implicit Runge-Kutta Methods for Wave-Like Equations
Runge-Kutta methods are used to integrate in time systems of differential equations. Implicit methods are designed to overcome numerical instabilities appearing during the evolution of a system of equations. We will present partially implicit Runge-Kutta methods for a particular structure of equations, generalization of a wave equation; the partially implicit term refers to this structure, where the implicit term appears only in a subset of the system of equations. These methods do not require any inversion of operators and the computational costs are similar to those of explicit Runge-Kutta methods. Partially implicit Runge-Kutta methods are derived up to third-order of convergence. We ana…
Fully relativistic non-linear cosmological evolution in spherical symmetry using the BSSN formalism
We present a fully relativistic numerical method for the study of cosmological problems using the Baumgarte-Shapiro-Shibata-Nakamura formalism on a dynamical Friedmann-Lema\^itre-Robertson-Walker background. This has many potential applications including the study of the growth of structures beyond the linear regime. We present one such application by reproducing the Lema\^itre-Tolman-Bondi solution for the collapse of pressureless matter with arbitrary lapse function. The regular and smooth numerical solution at the center of coordinates proceeds in a natural way by relying on the Partially Implicit Runge-Kutta algorithm described in Montero and Cordero-Carri\'on [arXiv:1211.5930]. We gene…
Improved constrained scheme for the Einstein equations: An approach to the uniqueness issue
Uniqueness problems in the elliptic sector of constrained formulations of Einstein equations have a dramatic effect on the physical validity of some numerical solutions, for instance when calculating the spacetime of very compact stars or nascent black holes. The fully constrained formulation (FCF) proposed by Bonazzola, Gourgoulhon, Grandcl\'ement, and Novak is one of these formulations. It contains, as a particular case, the approximation of the conformal flatness condition (CFC) which, in the last ten years, has been used in many astrophysical applications. The elliptic part of the FCF basically shares the same differential operators as the elliptic equations in CFC scheme. We present he…
Gravitational waves in Fully Constrained Formulation in a dynamical spacetime with matter content
We analyze numerically the behaviour of the hyperbolic sector of the Fully Constrained Formulation (FCF) (Bonazzola et al. 2004). The numerical experiments allow us to be confident in the performances of the upgraded version of the CoCoNuT code (Dimmelmeier et al. 2005) by replacing the Conformally Flat Condition (CFC), an approximation of Einstein equations, by FCF. First gravitational waves in FCF in a dynamical spacetime with matter content will be shown.
On numerical relativistic hydrodynamics and barotropic equations of state
The characteristic formulation of the relativistic hydrodynamic equations (Donat et al 1998 J. Comput. Phys. 146 58), which has been implemented in many relativistic hydro-codes that make use of Godunov-type methods, has to be slightly modified in the case of evolving barotropic flows. For a barotropic equation of state, a removable singularity appears in one of the eigenvectors. The singularity can be avoided by means of a simple renormalization which makes the system of eigenvectors well defined and complete. An alternative strategy for the particular case of barotropic flows is discussed.
Minimally implicit Runge-Kutta methods for Resistive Relativistic MHD
The Relativistic Resistive Magnetohydrodynamic (RRMHD) equations are a hyperbolic system of partial differential equations used to describe the dynamics of relativistic magnetized fluids with a finite conductivity. Close to the ideal magnetohydrodynamic regime, the source term proportional to the conductivity becomes potentially stiff and cannot be handled with standard explicit time integration methods. We propose a new class of methods to deal with the stiffness fo the system, which we name Minimally Implicit Runge-Kutta methods. These methods avoid the development of numerical instabilities without increasing the computational costs in comparison with explicit methods, need no iterative …
Master Majorana neutrino mass parametrization
After introducing a master formula for the Majorana neutrino mass matrix, we present a master parametrization for the Yukawa matrices automatically in agreement with neutrino oscillation data. This parametrization can be used for any model that induces Majorana neutrino masses. The application of the master parametrization is also illustrated in an example model, with special focus on its lepton flavor violating phenomenology.
GW190521: A Binary Black Hole Merger with a Total Mass of 150 M⊙
LIGO Scientific Collaboration and Virgo Collaboration: et al.
Mathematical Issues in a Fully-Constrained Formulation of Einstein Equations
Bonazzola, Gourgoulhon, Grandcl\'ement, and Novak [Phys. Rev. D {\bf 70}, 104007 (2004)] proposed a new formulation for 3+1 numerical relativity. Einstein equations result, according to that formalism, in a coupled elliptic-hyperbolic system. We have carried out a preliminary analysis of the mathematical structure of that system, in particular focusing on the equations governing the evolution for the deviation of a conformal metric from a flat fiducial one. The choice of a Dirac's gauge for the spatial coordinates guarantees the mathematical characterization of that system as a (strongly) hyperbolic system of conservation laws. In the presence of boundaries, this characterization also depen…
Characteristic structure of the resistive relativistic magnetohydrodynamic equations
We present the analysis of the characteristic structure of the resistive (non-ideal) relativistic magnetohydrodynamics system of equations. This is a necessary step to develop high-resolution shock-capturing schemes that use the full characteristic information (Godunov-type methods), and it is convenient to establish proper boundary conditions.
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…
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 …
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…
Dynamical spacetimes and gravitational radiation in a Fully Constrained Formulation
This contribution summarizes the recent work carried out to analyze the behavior of the hyperbolic sector of the Fully Constrained Formulation (FCF) derived in Bonazzola et al. 2004. The numerical experiments presented here allows one to be confident in the performances of the upgraded version of CoCoNuT's code by replacing the Conformally Flat Condition (CFC) approximation of the Einstein equations by the FCF.
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=…
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 …
On the equivalence between the Scheduled Relaxation Jacobi method and Richardson's non-stationary method
The Scheduled Relaxation Jacobi (SRJ) method is an extension of the classical Jacobi iterative method to solve linear systems of equations ($Au=b$) associated with elliptic problems. It inherits its robustness and accelerates its convergence rate computing a set of $P$ relaxation factors that result from a minimization problem. In a typical SRJ scheme, the former set of factors is employed in cycles of $M$ consecutive iterations until a prescribed tolerance is reached. We present the analytic form for the optimal set of relaxation factors for the case in which all of them are different, and find that the resulting algorithm is equivalent to a non-stationary generalized Richardson's method. …
Gravitational waves in dynamical spacetimes with matter content in the fully constrained formulation
The Fully Constrained Formulation (FCF) of General Relativity is a novel framework introduced as an alternative to the hyperbolic formulations traditionally used in numerical relativity. The FCF equations form a hybrid elliptic-hyperbolic system of equations including explicitly the constraints. We present an implicit-explicit numerical algorithm to solve the hyperbolic part, whereas the elliptic sector shares the form and properties with the well known Conformally Flat Condition (CFC) approximation. We show the stability andconvergence properties of the numerical scheme with numerical simulations of vacuum solutions. We have performed the first numerical evolutions of the coupled system of…
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 …
On the local existence of maximal slicings in spherically symmetric spacetimes
In this talk we show that any spherically symmetric spacetime admits locally a maximal spacelike slicing. The above condition is reduced to solve a decoupled system of first order quasi-linear partial differential equations. The solution may be accomplished analytical or numerically. We provide a general procedure to construct such maximal slicings.
Analysis of the Characteristics in the Meudon Constrained Evolution Scheme
A first analysis of the characteristics associated with the evolving modes in the constraint evolution scheme proposed by the Meudon group in 2004 is presented. The system is written in a first-order hyperbolic form and a so-called generalized Dirac gauge is considered. Applications to inner boundary conditions in an excised approach to black hole evolutions are discussed.
Scheduled Relaxation Jacobi method: improvements and applications
Elliptic partial differential equations (ePDEs) appear in a wide variety of areas of mathematics, physics and engineering. Typically, ePDEs must be solved numerically, which sets an ever growing demand for efficient and highly parallel algorithms to tackle their computational solution. The Scheduled Relaxation Jacobi (SRJ) is a promising class of methods, atypical for combining simplicity and efficiency, that has been recently introduced for solving linear Poisson-like ePDEs. The SRJ methodology relies on computing the appropriate parameters of a multilevel approach with the goal of minimizing the number of iterations needed to cut down the residuals below specified tolerances. The efficien…
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…
Nonlinear cosmological spherical collapse of quintessence
We present a study of the fully relativistic spherical collapse in the presence of quintessence using on numerical relativity, following the method proposed by the authors in a previous article [Phys. Rev. D 91, 024025 (2015)]. We ascertain the validity of the method by studying the evolution of a spherically symmetric quintessence inhomogeneity on a de Sitter background and we find that it has an impact on the local expansion around the center of coordinates. We then proceed to compare the results of our method to those of the more largely adopted top-hat model. We find that quintessence inhomogeneities do build up under the effect that matter inhomogeneities have on the local space-time, …
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…
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 …
General parametrization of Majorana neutrino mass models
We discuss a general formula which allows to automatically reproduce experimental data for Majorana neutrino mass models, while keeping the complete set of the remaining model parameters free for general scans, as necessary in order to provide reliable predictions for observables outside the neutrino sector. We provide a proof of this master parametrization and show how to apply it for several well-known neutrino mass models from the literature. We also discuss a list of special cases, in which the Yukawa couplings have to fulfill some particular additional conditions.
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…
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…
Trapping Horizons as inner boundary conditions for black hole spacetimes
We present a set of inner boundary conditions for the numerical construction of dynamical black hole space-times, when employing a 3+1 constrained evolution scheme and an excision technique. These inner boundary conditions are heuristically motivated by the dynamical trapping horizon framework and are enforced in an elliptic subsystem of the full Einstein equation. In the stationary limit they reduce to existing isolated horizon boundary conditions. A characteristic analysis completes the discussion of inner boundary conditions for the radiative modes.
Excision technique in constrained formulations of Einstein equations: collapse scenario
We present a new excision technique used in constrained formulations of Einstein equations to deal with black hole in numerical simulations. We show the applicability of this scheme in several scenarios. In particular, we present the dynamical evolution of the collapse of a neutron star to a black hole, using the CoCoNuT code and this excision technique.
On the convexity of relativistic ideal magnetohydrodynamics
We analyze the influence of the magnetic field in the convexity properties of the relativistic magnetohydrodynamics system of equations. To this purpose we use the approach of Lax, based on the analysis of the linearly degenerate/genuinely non-linear nature of the characteristic fields. Degenerate and non-degenerate states are discussed separately and the non-relativistic, unmagnetized limits are properly recovered. The characteristic fields corresponding to the material and Alfv\'en waves are linearly degenerate and, then, not affected by the convexity issue. The analysis of the characteristic fields associated with the magnetosonic waves reveals, however, a dependence of the convexity con…
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…
Maximal slicings in spherical symmetry: Local existence and construction
We show that any spherically symmetric spacetime locally admits a maximal spacelike slicing and we give a procedure allowing its construction. The construction procedure that we have designed is based on purely geometrical arguments and, in practice, leads to solve a decoupled system of first order quasi-linear partial differential equations. We have explicitly built up maximal foliations in Minkowski and Friedmann spacetimes. Our approach admits further generalizations and efficient computational implementation. As by product, we suggest some applications of our work in the task of calibrating Numerical Relativity complex codes, usually written in Cartesian coordinates.