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

Total Page:16

File Type:pdf, Size:1020Kb

GW190521: a Binary Black Hole Merger with a Total Mass of 150M ⊙ GW190521: A Binary Black Hole Merger with a Total Mass of 150M - Supplemental Material The LIGO Scientific Collaboration and the Virgo Collaboration SUPPLEMENTAL MATERIAL The CWB searches are performed with two pipeline configurations targeting detection of low mass (binary Coherent Wave Burst search BH configuration) and high mass (IMBH configuration) binary systems. They use different selection cuts de- fined a priori to alleviate the large variability of the non- Coherent Wave Burst (CWB) is a search algorithm for stationary noise in the detectors. The GW190521 event detection and reconstruction of gravitational wave (GW) was detected by CWB in its IMBH configuration. To fur- transient signals operating without a specific waveform ther improve the background the CWB analysis employs model. CWB searches for coincident signal power in additional selection cuts based on the generic properties multiple detectors with duration up to a few seconds. of IMBH binary signals: (a) the low central frequency The analysis is performed in the wavelet domain [1] on f of the GW signal and (b) their large chirp mass. In data normalized by the amplitude spectral density of the c CWB, the redshifted chirp mass M of a binary system detector noise. CWB selects wavelets with amplitudes is estimated by using the time-frequency evolution of the above the fluctuations of the detector noise and groups GW signal [4]. To select IMBH systems we require that them into clusters. For clusters correlated in multiple de- jMj > 10M and the signal central frequency f is below tectors, CWB identifies coherent events and reconstructs c 80 Hz. Both parameters M = 22:6M and f = 58 Hz the source sky location and signal waveforms with the c reconstructed for GW190521 are well within the limits constrained maximum likelihood method [2]. defined by the analysis selection cuts. The CWB detection statistic is based on the coherent To estimate the statistical significance of selected can- network energy E obtained by cross-correlating data in c p didates, each event was ranked against a set of back- different detectors; E is the estimator of the coherent c ground triggers obtained by repeating the analysis on network signal-to-noise ration (SNR) [2]. The agreement time-shifted data. The time shifts were selected to be between a GW signal and CWB reconstruction is charac- much longer (≥ 1 s) than the expected signal time delay terized by a reduced chi-squared statistic χ2 = E =N , n df between the detectors, in order to exclude astrophysical where E is the residual energy estimated after the re- n events from the background set. By using multiple time constructed waveforms are subtracted from the whitened shifts, a set of background events equivalent to 9800 yr of data, and N is the number of independent wavelet am- df null observation was obtained. The candidate events that plitudes describing the event. GW190521 is identified survived the selection criteria were then assigned a false- in LIGO data with χ2 = 0:68, showing no evidence for alarm rate given by the rate of background events with residual energy inconsistent with Gaussian noise. the CWB detection statistics greater than the detection To improve the robustness of the algorithm against statistic of the candidate event. non-stationary detector noise (glitches) and reduce the rate of false alarms, CWB uses signal-independent ve- toes. The primary veto cuts are on the network correla- 2 GstLAL and PyCBC searches tion coefficient cc = Ec=(Ec +En) and the χ . For a GW signal the expected values of both statistics are close to 2 unity and candidate events with cc < 0:7 and χ > 2:5 GstLAL [5{10] and PyCBC [11{16] independently are rejected as potential glitches. The blip glitches, de- implement search algorithms targeting compact binary scribed in the Data section of the letter, are decimated mergers with banks of template waveforms and matched by a selection cut based on a distinct shape of the glitch filtering. Both algorithms are used for low-latency as well waveform. It is best described by a cosine wave in a Gaus- as offline analyses. sian envelope with the standard deviation less than half All analyses based on GstLAL and PyCBC in O3 to a cycle. This selection cut starts to affect the detection date rely on SEOBNRv4 waveforms [17] as templates for of intermediate mass black hole (IMBH) binary systems the heavier-mass BH binary systems. These templates with total redshifted mass M > 300M , but the loss of neglect several physical effects. First, the spins of the detection efficiency is marginal (< 20%) as was tested compact objects are assumed to be aligned with the or- with simulated numerical relativity waveforms. In addi- bital angular momentum, such that the orbital plane does tion, the periodic variations of the noise in the frequency not precess. Second, the orbit is assumed to be quasicir- band 16-48 Hz, due to scattered light, are corrected with cular: its eccentricity is assumed to be negligible when the linear prediction error filter [3]. the signal becomes visible in the interferometers. Third, 2 only the quadrupolar radiation multipole moment is in- In order to assign a false-alarm rate to each candidate cluded. Fourth, finite-size effects (such as tidal deforma- coincidence, both the online and offline PyCBC analy- bility) are neglected. Simulation studies, as well as the ses use time shifts to construct background samples and published LIGO and Virgo discoveries, demonstrate that measure the cumulative background rate as a function of the above simplifications reduce the computational cost the ranking statistic. The two pipelines differ in how this of the search by at least an order of magnitude without process is implemented. The online analysis includes the compromising its sensitivity to most detectable signals data from the third interferometer, if available, in the [18, 19]. Nevertheless, there can be extreme systems for calculation of the false-alarm rate of a coincident can- which one or more of the above assumptions are not ap- didate. The offline pipeline calculates the background plicable, implying that the search may detect the system over an amount of data larger by more than an order with a reduced significance, or miss it altogether [19{21]. of magnitude, and can therefore claim correspondingly The template banks also necessarily cover a limited range lower false-alarm rates for events louder than the avail- of redshifted masses. In O3, the range of redshifted total able background. masses is [2; 758] M for GstLAL and [2; 500] M for Contrary to previously discovered compact binary 2 PyCBC [22, 23]; however, the effective high-mass bound- mergers, PyCBC's χr signal-based veto [24] reported a ary depends on mass ratio and spins. In the bank used disagreement between the event and the template wave- by PyCBC, templates with nearly equal mass and small form producing the highest significance in the search 2 effective spin have a maximum redshifted total mass just [17, 25]. Specifically, χr, expected to be ∼ 1 for sig- above 100 M , while this limit is ≈ 600M for the bank nals well described by the search templates, resulted in used by GstLAL. 1:1 (LHO) and 3:2 (LLO) in the offline analysis; the on- GstLAL and PyCBC differ in several important tech- line analysis showed a similar behavior. The large LLO 2 nical choices, notably the signal-based vetoes used to re- χr reflects a deficit of signal power between ∼ 23 Hz > ject instrumental transients, the ranking of GW candi- and ∼ 43 Hz, and similarly for ∼ 68 Hz, with respect dates resulting from matched filtering, and the method to the template waveform. This effect can be due to a used to assign false-alarm rates to the ranked GW can- combination of high redshifted total mass, which places didates. These differences can lead to vastly different GW190521 at the limit of the mass space covered by Py- significances for the same astrophysical signal. CBC's configuration in O3 [23], and possibly waveform GstLAL uses a likelihood-ratio statistic to rank GW effects not modeled by the templates, such as precession, candidates by (i) the matched-filter SNRs, (ii) a sig- higher-order multipole moments, or eccentricity. 2 nal consistency check based on waveform autocorrelation, In order to investigate this point, we recalculated χr at (iii) the sensitivities of the detectors, (iv) the time and LLO using different quasicircular waveform models with phase delays observed across the detector network, and parameters that maximize the network SNR, obtaining (v) a template-dependent factor that describes the prob- values between 0:8 and 1:8, which no longer indicate a ability that a signal from the observed source population disagreement. The LHO{LLO network matched-filter is identified by a specific template in our bank. Informa- SNR correspondingly increases to 14:5{15:0, consistent tion from auxiliary channels that monitor the detectors with the values obtained by CWB and GstLAL. As an and their environment is also included in the ranking. additional check, we then added the maximum-likelihood GstLAL estimates the search background from trig- SEOBNRv4 template to the bank used by PyCBC and gers not found in coincidence at times when both LIGO repeated the offline analysis. The new template was pre- interferometers are operating. The background triggers ferred by the ranking over the previous template, increas- are used to estimate the probability density of noise-like ing the LHO{LLO network SNR to 14:5 and decreasing 2 events, and Monte Carlo methods are used to sample the the LLO χr to 1:7. These values are consistent with the noise probability density and produce a mapping from other pipelines.
Recommended publications
  • Numerical Relativity for Vacuum Binary Black Hole Luisa T
    Numerical Relativity for Vacuum Binary Black Hole Luisa T. Buchman UW Bothell SXS collaboration (Caltech) Max Planck Institute for Gravitational Physics Albert Einstein Institute GWA NW @ LIGO Hanford Observatory June 25, 2019 Outline • What is numerical relativity? • What is its role in the detection and interpretation of gravitational waves? What is numerical relativity for binary black holes? • 3 stages for binary black hole coalescence: • inspiral • merger • ringdown • Inspiral waveform: Post-Newtonian approximations • Ringdown waveform: perturbation theory • Merger of 2 black holes: • extremely energetic, nonlinear, dynamical and violent • the strongest GW signal • warpage of spacetime • strong gravity regime The Einstein equations: • Gμν = 8πTμν • encodes all of gravity • 10 equations which would take up ~100 pages if no abbreviations or simplifications used (and if written in terms of the spacetime metric alone) • pencil and paper solutions possible only with spherical (Schwarzschild) or axisymmetric (Kerr) symmetries. • Vacuum: G = 0 μν (empty space with no matter) Numerical Relativity is: • Directly solving the full dynamical Einstein field equations using high-performance computing. • A very difficult problem: • the first attempts at numerical simulations of binary black holes was in the 1960s (Hahn and Lindquist) and mid-1970s (Smarr and others). • the full 3D problem remained unsolved until 2005 (Frans Pretorius, followed quickly by 2 other groups- Campanelli et al., Baker et al.). Solving Einstein’s Equations on a Computer •
    [Show full text]
  • Towards a Gauge Polyvalent Numerical Relativity Code: Numerical Methods, Boundary Conditions and Different Formulations
    UNIVERSITAT DE LES ILLES BALEARS Towards a gauge polyvalent numerical relativity code: numerical methods, boundary conditions and different formulations per Carles Bona-Casas Tesi presentada per a l'obtenci´o del t´ıtolde doctor a la Facultat de Ci`encies Departament de F´ısica Dirigida per: Prof. Carles Bona Garcia i Dr. Joan Mass´oBenn`assar "The fact that we live at the bottom of a deep gravity well, on the surface of a gas covered planet going around a nuclear fireball 90 million miles away and think this to be normal is obviously some indication of how skewed our perspective tends to be." Douglas Adams. Acknowledgements I would like to acknowledge everyone who ever taught me something. Specially my supervisors who, at least in the beginning, had to suffer my eyes, puzzling at their faces, saying that I had the impression that what they were explaining to me was some unreachable knowledge. It turns out that there is not such a thing and they have finally managed to make me gain some insight in this world of relativity and computing. They have also infected me the disease of worrying about calculations and stuff that most people wouldn't care about and my hair doesn't seem very happy with it, but I still love them for that. Many thanks to everyone in the UIB relativity group and to all the PhD. students who have shared some time with me. Work is less hard in their company. Special thanks to Carlos, Denis, Sasha and Dana for all their very useful comments and for letting me play with their work and codes.
    [Show full text]
  • Numerical Relativity
    Paper presented at the 13th Int. Conf on General Relativity and Gravitation 373 Cordoba, Argentina, 1992: Part 2, Workshop Summaries Numerical relativity Takashi Nakamura Yulcawa Institute for Theoretical Physics, Kyoto University, Kyoto 606, JAPAN In GR13 we heard many reports on recent. progress as well as future plans of detection of gravitational waves. According to these reports (see the report of the workshop on the detection of gravitational waves by Paik in this volume), it is highly probable that the sensitivity of detectors such as laser interferometers and ultra low temperature resonant bars will reach the level of h ~ 10—21 by 1998. in this level we may expect the detection of the gravitational waves from astrophysical sources such as coalescing binary neutron stars once a year or so. Therefore the progress in numerical relativity is urgently required to predict the wave pattern and amplitude of the gravitational waves from realistic astrophysical sources. The time left for numerical relativists is only six years or so although there are so many difficulties in principle as well as in practice. Apart from detection of gravitational waves, numerical relativity itself has a final goal: Solve the Einstein equations numerically for (my initial data as accurately as possible and clarify physics in strong gravity. in GRIIS there were six oral presentations and ll poster papers on recent progress in numerical relativity. i will make a brief review of six oral presenta— tions. The Regge calculus is one of methods to investigate spacetimes numerically. Brewin from Monash University Australia presented a paper Particle Paths in a Schwarzshild Spacetime via.
    [Show full text]
  • Arxiv:1812.04350V2 [Gr-Qc] 16 Jan 2020 PACS Numbers: 04.70.Bw, 04.80.Nn, 95.10.Fh
    International Journal of Modern Physics D Testing dispersion of gravitational waves from eccentric extreme-mass-ratio inspirals Shu-Cheng Yang*,y , Wen-Biao Han*,y,x, Shuo Xinz and Chen Zhang*,y *Shanghai Astronomical Observatory, Chinese Academy of Sciences, Shanghai, 200030, P. R. China ySchool of Astronomy and Space Science, University of Chinese Academy of Sciences, Beijing, 100049, P. R. China zSchool of Physics Sciences and Engineering, Tongji University, Shanghai 200092, P. R. China [email protected] In general relativity, there is no dispersion in gravitational waves, while some modified gravity theories predict dispersion phenomena in the propagation of gravitational waves. In this paper, we demonstrate that this dispersion will induce an observable deviation of waveforms if the orbits have large eccentricities. The mechanism is that the waveform modes with different frequencies will be emitted at the same time due to the existence of eccentricity. During the propagation, because of the dispersion, the arrival time of different modes will be different, then produce the deviation and dephasing of wave- forms compared with general relativity. This kind of dispersion phenomena related with extreme-mass-ratio inspirals could be observed by space-borne detectors, and the con- straint on the graviton mass could be improved . Moreover, we find that the dispersion effect may also be constrained by ground detectors better than the current result if a highly eccentric intermediate-mass-ratio inspirals be observed. Keywords: Gravitational waves; extreme-mass-ratio inspirals; dispersion. arXiv:1812.04350v2 [gr-qc] 16 Jan 2020 PACS numbers: 04.70.Bw, 04.80.Nn, 95.10.Fh 1.
    [Show full text]
  • GW150914: First Results from the Search for Binary Black Hole Coalescence with Advanced LIGO
    GW150914: First results from the search for binary black hole coalescence with Advanced LIGO The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation Abbott, B.P.; Abbott, R.; Abbott, T.D.; Abernathy, M.R.; Acernese, F.; Ackley, K.; Adams, C.; Adams, T.; Addesso, P.; Adhikari, R.X.; Adya, V.B. et al. "GW150914: First results from the search for binary black hole coalescence with Advanced LIGO." Physical Review D 93, 122003 (June 2016): 1-21 © 2016 American Physical Society As Published http://dx.doi.org/10.1103/PhysRevD.93.122003 Publisher American Physical Society Version Final published version Citable link http://hdl.handle.net/1721.1/110492 Terms of Use Article is made available in accordance with the publisher's policy and may be subject to US copyright law. Please refer to the publisher's site for terms of use. PHYSICAL REVIEW D 93, 122003 (2016) GW150914: First results from the search for binary black hole coalescence with Advanced LIGO B. P. Abbott et al.* (LIGO Scientific Collaboration and Virgo Collaboration) (Received 9 March 2016; published 7 June 2016) On September 14, 2015, at 09∶50:45 UTC the two detectors of the Laser Interferometer Gravitational- Wave Observatory (LIGO) simultaneously observed the binary black hole merger GW150914. We report the results of a matched-filter search using relativistic models of compact-object binaries that recovered GW150914 as the most significant event during the coincident observations between the two LIGO detectors from September 12 to October 20, 2015 GW150914 was observed with a matched-filter signal-to- noise ratio of 24 and a false alarm rate estimated to be less than 1 event per 203000 years, equivalent to a significance greater than 5.1 σ.
    [Show full text]
  • Stringent Constraints on Neutron-Star Radii from Multimessenger Observations and Nuclear Theory
    Stringent constraints on neutron-star radii from multimessenger observations and nuclear theory Collin D. Capano1;2∗, Ingo Tews3, Stephanie M. Brown1;2, Ben Margalit4;5;6, Soumi De6;7, Sumit Kumar1;2, Duncan A. Brown6;7, Badri Krishnan1;2, Sanjay Reddy8;9 1 Albert-Einstein-Institut, Max-Planck-Institut fur¨ Gravitationsphysik, Callinstraße 38, 30167 Hannover, Germany, 2 Leibniz Universitat¨ Hannover, 30167 Hannover, Germany, 3 Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87545, USA, 4 Department of Astronomy and Theoretical Astrophysics Center, University of California, Berkeley, CA 94720, USA, 5 NASA Einstein Fellow 6 Kavli Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106, USA, 7 Department of Physics, Syracuse University, Syracuse NY 13244, USA, 8 Institute for Nuclear Theory, University of Washington, Seattle, WA 98195-1550, USA, 9 JINA-CEE, Michigan State University, East Lansing, MI, 48823, USA The properties of neutron stars are determined by the nature of the matter that they contain. These properties can be constrained by measurements of the star’s size. We obtain stringent constraints on neutron-star radii by combining multimessenger observations of the binary neutron-star merger GW170817 arXiv:1908.10352v3 [astro-ph.HE] 24 Mar 2020 with nuclear theory that best accounts for density-dependent uncertainties in the equation of state. We construct equations of state constrained by chi- ral effective field theory and marginalize over these using the gravitational- wave observations. Combining this with the electromagnetic observations of the merger remnant that imply the presence of a short-lived hyper-massive 1 neutron star, we find that the radius of a 1:4 M neutron star is R1:4 M = +0:9 11:0−0:6 km (90% credible interval).
    [Show full text]
  • The Initial Data Problem for 3+1 Numerical Relativity Part 2
    The initial data problem for 3+1 numerical relativity Part 2 Eric Gourgoulhon Laboratoire Univers et Th´eories (LUTH) CNRS / Observatoire de Paris / Universit´eParis Diderot F-92190 Meudon, France [email protected] http://www.luth.obspm.fr/∼luthier/gourgoulhon/ 2008 International Summer School on Computational Methods in Gravitation and Astrophysics Asia Pacific Center for Theoretical Physics, Pohang, Korea 28 July - 1 August 2008 Eric Gourgoulhon (LUTH) Initial data problem 2 / 2 APCTP School, 31 July 2008 1 / 41 Plan 1 Helical symmetry for binary systems 2 Initial data for orbiting binary black holes 3 Initial data for orbiting binary neutron stars 4 Initial data for orbiting black hole - neutron star systems 5 References for lectures 1-3 Eric Gourgoulhon (LUTH) Initial data problem 2 / 2 APCTP School, 31 July 2008 2 / 41 Helical symmetry for binary systems Outline 1 Helical symmetry for binary systems 2 Initial data for orbiting binary black holes 3 Initial data for orbiting binary neutron stars 4 Initial data for orbiting black hole - neutron star systems 5 References for lectures 1-3 Eric Gourgoulhon (LUTH) Initial data problem 2 / 2 APCTP School, 31 July 2008 3 / 41 Helical symmetry for binary systems Helical symmetry for binary systems Physical assumption: when the two objects are sufficiently far apart, the radiation reaction can be neglected ⇒ closed orbits Gravitational radiation reaction circularizes the orbits ⇒ circular orbits Geometrical translation: spacetime possesses some helical symmetry Helical Killing vector ξ:
    [Show full text]
  • 3+1 Formalism and Bases of Numerical Relativity
    3+1 Formalism and Bases of Numerical Relativity Lecture notes Eric´ Gourgoulhon Laboratoire Univers et Th´eories, UMR 8102 du C.N.R.S., Observatoire de Paris, Universit´eParis 7 arXiv:gr-qc/0703035v1 6 Mar 2007 F-92195 Meudon Cedex, France [email protected] 6 March 2007 2 Contents 1 Introduction 11 2 Geometry of hypersurfaces 15 2.1 Introduction.................................... 15 2.2 Frameworkandnotations . .... 15 2.2.1 Spacetimeandtensorfields . 15 2.2.2 Scalar products and metric duality . ...... 16 2.2.3 Curvaturetensor ............................... 18 2.3 Hypersurfaceembeddedinspacetime . ........ 19 2.3.1 Definition .................................... 19 2.3.2 Normalvector ................................. 21 2.3.3 Intrinsiccurvature . 22 2.3.4 Extrinsiccurvature. 23 2.3.5 Examples: surfaces embedded in the Euclidean space R3 .......... 24 2.4 Spacelikehypersurface . ...... 28 2.4.1 Theorthogonalprojector . 29 2.4.2 Relation between K and n ......................... 31 ∇ 2.4.3 Links between the and D connections. .. .. .. .. .. 32 ∇ 2.5 Gauss-Codazzirelations . ...... 34 2.5.1 Gaussrelation ................................. 34 2.5.2 Codazzirelation ............................... 36 3 Geometry of foliations 39 3.1 Introduction.................................... 39 3.2 Globally hyperbolic spacetimes and foliations . ............. 39 3.2.1 Globally hyperbolic spacetimes . ...... 39 3.2.2 Definition of a foliation . 40 3.3 Foliationkinematics .. .. .. .. .. .. .. .. ..... 41 3.3.1 Lapsefunction ................................. 41 3.3.2 Normal evolution vector . 42 3.3.3 Eulerianobservers ............................. 42 3.3.4 Gradients of n and m ............................. 44 3.3.5 Evolution of the 3-metric . 45 4 CONTENTS 3.3.6 Evolution of the orthogonal projector . ....... 46 3.4 Last part of the 3+1 decomposition of the Riemann tensor .
    [Show full text]
  • Advanced Virgo: Status of the Detector, Latest Results and Future Prospects
    universe Review Advanced Virgo: Status of the Detector, Latest Results and Future Prospects Diego Bersanetti 1,* , Barbara Patricelli 2,3 , Ornella Juliana Piccinni 4 , Francesco Piergiovanni 5,6 , Francesco Salemi 7,8 and Valeria Sequino 9,10 1 INFN, Sezione di Genova, I-16146 Genova, Italy 2 European Gravitational Observatory (EGO), Cascina, I-56021 Pisa, Italy; [email protected] 3 INFN, Sezione di Pisa, I-56127 Pisa, Italy 4 INFN, Sezione di Roma, I-00185 Roma, Italy; [email protected] 5 Dipartimento di Scienze Pure e Applicate, Università di Urbino, I-61029 Urbino, Italy; [email protected] 6 INFN, Sezione di Firenze, I-50019 Sesto Fiorentino, Italy 7 Dipartimento di Fisica, Università di Trento, Povo, I-38123 Trento, Italy; [email protected] 8 INFN, TIFPA, Povo, I-38123 Trento, Italy 9 Dipartimento di Fisica “E. Pancini”, Università di Napoli “Federico II”, Complesso Universitario di Monte S. Angelo, I-80126 Napoli, Italy; [email protected] 10 INFN, Sezione di Napoli, Complesso Universitario di Monte S. Angelo, I-80126 Napoli, Italy * Correspondence: [email protected] Abstract: The Virgo detector, based at the EGO (European Gravitational Observatory) and located in Cascina (Pisa), played a significant role in the development of the gravitational-wave astronomy. From its first scientific run in 2007, the Virgo detector has constantly been upgraded over the years; since 2017, with the Advanced Virgo project, the detector reached a high sensitivity that allowed the detection of several classes of sources and to investigate new physics. This work reports the Citation: Bersanetti, D.; Patricelli, B.; main hardware upgrades of the detector and the main astrophysical results from the latest five years; Piccinni, O.J.; Piergiovanni, F.; future prospects for the Virgo detector are also presented.
    [Show full text]
  • Search for Gravitational Waves from High-Mass-Ratio Compact-Binary Mergers of Stellar Mass and Subsolar Mass Black Holes
    PHYSICAL REVIEW LETTERS 126, 021103 (2021) Search for Gravitational Waves from High-Mass-Ratio Compact-Binary Mergers of Stellar Mass and Subsolar Mass Black Holes Alexander H. Nitz * and Yi-Fan Wang (王一帆) Max-Planck-Institut für Gravitationsphysik (Albert-Einstein-Institut), D-30167 Hannover, Germany and Leibniz Universität Hannover, D-30167 Hannover, Germany (Received 8 July 2020; accepted 11 December 2020; published 12 January 2021) We present the first search for gravitational waves from the coalescence of stellar mass and subsolar mass 3 black holes with masses between 20–100 M⊙ and 0.01–1 M⊙ð10–10 MJÞ, respectively. The observation of a single subsolar mass black hole would establish the existence of primordial black holes and a possible component of dark matter. We search the ∼164 day of public LIGO data from 2015–2017 when LIGO- Hanford and LIGO-Livingston were simultaneously observing. We find no significant candidate gravitational-wave signals. Using this nondetection, we place a 90% upper limit on the rate of 6 4 −3 −1 30–0.01 M⊙ and 30–0.1 M⊙ mergers at < 1.2 × 10 and < 1.6 × 10 Gpc yr , respectively. If we consider binary formation through direct gravitational-wave braking, this kind of merger would be exceedingly rare if only the lighter black hole were primordial in origin (< 10−4 Gpc−3 yr−1). If both black holes are primordial in origin, we constrain the contribution of 1ð0.1ÞM⊙ black holes to dark matter to < 0.3ð3Þ%. DOI: 10.1103/PhysRevLett.126.021103 Introduction.—The first gravitational wave observation known model which can produce subsolar mass black holes from the merger of black holes was detected on September by conventional formation mechanisms; the observation of 14, 2015 [1].
    [Show full text]
  • Probing the Nuclear Equation of State from the Existence of a 2.6 M Neutron Star: the GW190814 Puzzle
    S S symmetry Article Probing the Nuclear Equation of State from the Existence of a ∼2.6 M Neutron Star: The GW190814 Puzzle Alkiviadis Kanakis-Pegios †,‡ , Polychronis S. Koliogiannis ∗,†,‡ and Charalampos C. Moustakidis †,‡ Department of Theoretical Physics, Aristotle University of Thessaloniki, 54124 Thessaloniki, Greece; [email protected] (A.K.P.); [email protected] (C.C.M.) * Correspondence: [email protected] † Current address: Aristotle University of Thessaloniki, 54124 Thessaloniki, Greece. ‡ These authors contributed equally to this work. Abstract: On 14 August 2019, the LIGO/Virgo collaboration observed a compact object with mass +0.08 2.59 0.09 M , as a component of a system where the main companion was a black hole with mass ∼ − 23 M . A scientific debate initiated concerning the identification of the low mass component, as ∼ it falls into the neutron star–black hole mass gap. The understanding of the nature of GW190814 event will offer rich information concerning open issues, the speed of sound and the possible phase transition into other degrees of freedom. In the present work, we made an effort to probe the nuclear equation of state along with the GW190814 event. Firstly, we examine possible constraints on the nuclear equation of state inferred from the consideration that the low mass companion is a slow or rapidly rotating neutron star. In this case, the role of the upper bounds on the speed of sound is revealed, in connection with the dense nuclear matter properties. Secondly, we systematically study the tidal deformability of a possible high mass candidate existing as an individual star or as a component one in a binary neutron star system.
    [Show full text]
  • Introduction to Numerical Relativity Erik SchneEr Perimeter InsTute for TheoreCal Physics CGWAS 2013, Caltech, 2013-07-26 What Is Numerical Relativity?
    Introduction to Numerical Relativity Erik Schneer Perimeter Ins1tute for Theore1cal Physics CGWAS 2013, Caltech, 2013-07-26 What is Numerical Relativity? • Solving Einstein equaons numerically • Can handle arbitrarily complex systems • Sub-field of Computaonal Astrophysics • Also beginning to be relevant in cosmology • Einstein equaons relevant only for dense, compact objects: • Black holes • Neutron stars Overview • General Relavity: Geometry, Coordinates • Solving the Einstein equaons • Relavis1c Hydrodynamics • Analyzing Space1mes: horizons, gravitaonal waves Some relevant concepts GENERAL RELATIVITY Einstein Equations • Gab = 8π Tab • Gab: Einstein tensor, one measure of curvature • Tab: stress-energy tensor, describes mass/energy/momentum/ pressure/stress densi1es • Gab and Tab are symmetric: 10 independent components in 4D • Loose reading: (some part of) the space-1me curvature equals (“is generated by”) its maer content Spacetime Curvature • Difference between special and general relavity: in GR, space1me is curved • 2D example of a curved manifold: earth’s surface • Can’t use a straight coordinate system for a curved manifold! • E.g. Cartesian coordinate system doesn’t “fit” earth’s surface • In GR, one needs to use curvilinear, 1me-dependent coordinate systems • In fact, if one knows how to use arbitrary coordinate systems for a theory (e.g. hydrodynamics or electrodynamics), then working with this theory in GR becomes trivial Riemann, Ricci, Weyl • Curvature is measured by Riemann tensor Rabcd; has 20 independent components in 4D
    [Show full text]