Pycbc Inference: a Python-Based Parameter Estimation Toolkit for Compact Binary Coalescence Signals

Total Page:16

File Type:pdf, Size:1020Kb

Pycbc Inference: a Python-Based Parameter Estimation Toolkit for Compact Binary Coalescence Signals PyCBC Inference: A Python-based parameter estimation toolkit for compact binary coalescence signals C. M. Biwer Los Alamos National Laboratory, Los Alamos, NM 87545, USA Department of Physics, Syracuse University, Syracuse, NY 13244, USA Collin D. Capano Albert-Einstein-Institut, Max-Planck-Institut f¨urGravitationsphysik, D-30167 Hannover, Germany Soumi De Department of Physics, Syracuse University, Syracuse, NY 13244, USA Miriam Cabero Albert-Einstein-Institut, Max-Planck-Institut f¨urGravitationsphysik, D-30167 Hannover, Germany Duncan A. Brown Department of Physics, Syracuse University, Syracuse, NY 13244, USA Alexander H. Nitz Albert-Einstein-Institut, Max-Planck-Institut f¨urGravitationsphysik, D-30167 Hannover, Germany V. Raymond Albert-Einstein-Institut, Max-Planck-Institut f¨urGravitationsphysik, D-14476 Potsdam, Germany School of Physics and Astronomy, Cardiff University, Cardiff, CF243AA, Wales, UK Abstract. We introduce new modules in the open-source PyCBC gravitational- arXiv:1807.10312v1 [astro-ph.IM] 26 Jul 2018 wave astronomy toolkit that implement Bayesian inference for compact-object binary mergers. We review the Bayesian inference methods implemented and describe the structure of the modules. We demonstrate that the PyCBC Inference modules produce unbiased estimates of the parameters of a simulated population of binary black hole mergers. We show that the posterior parameter distributions obtained used our new code agree well with the published estimates for binary black holes in the first LIGO-Virgo observing run. PyCBC Inference: A Python-based parameter estimation toolkit for compact binary coalescence signals 2 1. Introduction of binary black holes and the black-hole mergers detected in the first LIGO-Virgo observing run to The observations of six binary black hole mergers [1, 2, demonstrate the use of PyCBC Inference. We provide 3, 4], and the binary neutron star merger GW170817 [5] the posterior probability densities for the events by Advanced LIGO [6] and Virgo [7] have established GW150914, GW151226, and LVT151012, and the the field of gravitational-wave astronomy. Understand- command lines and configurations to reproduce these ing the origin, evolution, and physics of gravitational- results as supplemental materials [14]. Finally, we wave sources requires accurately measuring the prop- summarize the status of the code and possible future erties of detected events. In practice, this is per- developments in Sec. 5. formed using Bayesian inference [8, 9]. Bayesian in- ference allows us to determine the signal model that 2. Bayesian Inference for Binary Mergers is best supported by observations and to obtain pos- terior probability densities for a model's parameters, In gravitational-wave astronomy, Bayesian methods hence inferring the properties of the source. In this are used to infer the properties of detected astrophys- paper, we present PyCBC Inference; a set of Python ical sources [15, 16, 17, 18]. Given the observed data modules that implement Bayesian inference in the Py- d~(t)|here this is data from a gravitational-wave de- CBC open-source toolkit for gravitational-wave astron- tector network in which a search has identified a sig- omy [10]. PyCBC Inference has been used to perform nal [19, 20, 21]|Bayes' theorem [8, 9] states that for a Bayesian inference for several astrophysical problems, hypothesis H, including: testing the black hole area increase law [11]; p(d~(t)j#;~ H)p(#~jH) combining multi-messenger obervations of GW170817 p(#~jd~(t);H) = : (1) to constrain the viewing angle of the binary [12]; and p(d~(t)jH) determining that the gravitational-wave observations In our case, hypothesis H is the model of the of GW170817 favor a model where both compact ob- gravitational-wave signal and #~ are the parameters of jects have the same equation of state, and measur- this model. Together, these describe the properties of ing the tidal deformabilities and radii of the neutron the astrophysical source of the gravitational waves. In stars [13]. Eq. (1), the prior probability density p(#~jH) describes We provide a comprehensive description of the our knowledge about the parameters before considering methods and code implemented in PyCBC Inference. the observed data d~(t), and the likelihood p(d~(t)j#;~ H) We then demonstrate that PyCBC Inference can is the probability of obtaining the observation d~(t) produce unbiased estimates of the parameters of a given the waveform model H with parameters #~. simulated population of binary black holes. We Often we are only interested in a subset of the show that PyCBC Inference can recover posterior parameters #~. To obtain a probability distribution probability distributions that are in good agreement on one or a few parameters, we marginalize the with the published measurements of the binary black posterior probability by integrating p(d~(t)j#;~ H)p(#~jH) holes detected in the first LIGO-Virgo observing over the unwanted parameters. Marginalizing over all run [1]. This paper is organized as follows: Sec. 2 parameters yields the evidence, p(d~(t)jH), which is gives an overview of the Bayesian inference methods the denominator in Eq. (1). The evidence serves as a used in gravitational-wave astronomy for compact- normalization constant of the posterior probability for object binary mergers. We provide an overview of the given model H. If we have two competing models the waveform models used; the likelihood function H and H , the evidence can be used to determine for a known signal in stationary, Gaussian noise; A B which model is favored by the data via the Bayes the sampling methods used to estimate the posterior factor [22, 23, 24], probability densities and the evidence; the selection of ~ independent samples; and the estimation of parameter p(d(t)jHA) B = : (2) ~ values from posterior probabilities. Sec. 3 describes p(d(t)jHB) the design of the PyCBC Inference software and how If B is greater than 1 then model HA is favored over the methods described in Sec. 2 are implemented HB, with the magnitude of B indicating the degree of in the code. Sec. 4 uses a simulated population belief. PyCBC Inference: A Python-based parameter estimation toolkit for compact binary coalescence signals 3 PyCBC Inference can compute Bayes factors and angular momentum), reducing the dimensionality of produce marginalized posterior probability densities the waveform parameters space. given the data from a network of gravitational-wave Given a set of parameters #~, one can obtain ~ observatories with N detectors d(t) = fdi(t); 1 < i < a model of the gravitational-wave signal from a Ng, and a model H that describes the astrophysical binary merger using a variety of different methods, source. In the remainder of this section, we review the including: post-Newtonian theory (see e.g. Ref. [33] methods used to compute these quantities. and references therein), analytic models calibrated against numerical simulations [34, 35, 36, 37, 38, 39, 2.1. Waveform Models 40], perturbation theory [41, 42], and full numerical solution of the Einstein equations (see e.g. Ref. [43] and The gravitational waves radiated in a binary merger's references therein). Obtaining posterior probabilities source frame are described by the component masses and evidences can require calculating O 109 template m1;2, the three-dimensional spin vectors ~s1;2 of the waveforms, which restricts us to models that are compact objects [25], and the binary's eccentricity computationally efficient to calculate. The cost of e [26]. A parameter φ describes the phase of the binary full numerical simulations makes them prohibitively at a fiducial reference time, although this is not usually expensive at present. Even some analytic models are of physical interest. For binaries containing neutron too costly to be used, and surrogate models have been stars, additional parameters Λ1;2 describe the star's developed that capture the features of these waveforms tidal deformabilities [27, 28], which depend on the at reduced computational cost [44, 45]. equation of state of the neutron stars. The waveform The specific choice of the waveform model H observed by the Earth-based detector network depends for an analysis depends on the physics that we wish on seven additional parameters: the signal's time of to explore, computational cost limitations, and the arrival tc, the binary's luminosity distance dL, and level of accuracy desired in the model. A variety four Euler angles that describe the transformation from of waveform models are available for use in PyCBC the binary's frame to the detector network frame [29]. Inference, either directly implemented in PyCBC or These angles are typically written as the binary's right via calls to the LIGO Algorithm Library (LAL) [46]. ascension α, declination δ, a polarization angle Ψ, and We refer to the PyCBC and LAL documentation, and the inclination angle ι (the angle between the binary's references therein, for detailed descriptions of these angular momentum axis and the line of sight). models. In this paper, we demonstrate the use of Binary mergers present a challenging problem PyCBC Inference using the IMRPhenomPv2 [47, 48] for Bayesian inference, as the dimensionality of the waveform model for binary black hole mergers. This signal parameter space is large. This is further model captures the inspiral-merger-ringdown physics complicated by correlations between the signal's of spinning, precessing binaries and parameterizing parameters. For example, at leading order the spin effects using a spin magnitude aj, an azimuthal gravitational waveform depends on the chirp mass a p angle θj , and a polar angle θj for each of the two M [30]. The mass ratio enters the waveform compact objects. Examples of using PyCBC Inference at higher orders and is more difficult to measure. with different waveform models include the analysis This results in an amplitude-dependent degeneracy of Ref. [13] that used the TaylorF2 post-Newtonian between the component masses [18]. Similarly, waveform model with tidal corrections, and Ref.
Recommended publications
  • 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]
  • 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]
  • Arxiv:1606.03939V2 [Astro-Ph.HE] 20 Sep 2016 P0(X) and P1(X) for Its Detection Statistic X
    Draft version September 21, 2016 Preprint typeset using LATEX style AASTeX6 v. 1.0 SUPPLEMENT: THE RATE OF BINARY BLACK HOLE MERGERS INFERRED FROM ADVANCED LIGO OBSERVATIONS SURROUNDING GW150914 B. P. Abbott,1 yDeceased, May 2015. zDeceased, March 2015. (LIGO Scientific Collaboration and Virgo Collaboration) (and The LIGO Scientific Collaboration and Virgo Collaboration) 1LIGO, California Institute of Technology, Pasadena, CA 91125, USA ABSTRACT Supplemental information for a Letter reporting the rate of binary black hole (BBH) coalescences in- ferred from 16 days of coincident Advanced LIGO observations surrounding the transient gravitational wave (GW) signal GW150914. In that work we reported various rate estimates whose 90% credible 3 1 intervals (CIs) fell in the range 2{600 Gpc− yr− . Here we give details of our method and com- putations, including information about our search pipelines, a derivation of our likelihood function for the analysis, a description of the astrophysical search trigger distribution expected from merging BBHs, details on our computational methods, a description of the effects and our model for calibration uncertainty, and an analytic method of estimating our detector sensitivity that is calibrated to our measurements. The first detection of a gravitational wave (GW) sig- a coincident counterpart in the other detector to form a nal from a merging binary black hole (BBH) system is multiple-detector likelihood ratio. described in Abbott et al.(2016d). Abbott et al.(2016g) Both pipelines rely on an empirical estimate of the reports on inference of the local BBH merger rate from search background, making the assumption that trig- surrounding Advanced LIGO observations. This Sup- gers of terrestrial origin occur independently in the two plement provides supporting material and methodolog- detectors.
    [Show full text]
  • Recent Observations of Gravitational Waves by LIGO and Virgo Detectors
    universe Review Recent Observations of Gravitational Waves by LIGO and Virgo Detectors Andrzej Królak 1,2,* and Paritosh Verma 2 1 Institute of Mathematics, Polish Academy of Sciences, 00-656 Warsaw, Poland 2 National Centre for Nuclear Research, 05-400 Otwock, Poland; [email protected] * Correspondence: [email protected] Abstract: In this paper we present the most recent observations of gravitational waves (GWs) by LIGO and Virgo detectors. We also discuss contributions of the recent Nobel prize winner, Sir Roger Penrose to understanding gravitational radiation and black holes (BHs). We make a short introduction to GW phenomenon in general relativity (GR) and we present main sources of detectable GW signals. We describe the laser interferometric detectors that made the first observations of GWs. We briefly discuss the first direct detection of GW signal that originated from a merger of two BHs and the first detection of GW signal form merger of two neutron stars (NSs). Finally we present in more detail the observations of GW signals made during the first half of the most recent observing run of the LIGO and Virgo projects. Finally we present prospects for future GW observations. Keywords: gravitational waves; black holes; neutron stars; laser interferometers 1. Introduction The first terrestrial direct detection of GWs on 14 September 2015, was a milestone Citation: Kro´lak, A.; Verma, P. discovery, and it opened up an entirely new window to explore the universe. The combined Recent Observations of Gravitational effort of various scientists and engineers worldwide working on the theoretical, experi- Waves by LIGO and Virgo Detectors.
    [Show full text]
  • Pycbc Live a Real-Time Search for Compact Binary Mergers in the Advanced Detector Era
    PyCBC Live A real-time search for compact binary mergers in the advanced detector era Tito Dal Canton – IJCLab 3e réunion d'animation annuelle in collaboration with GdR ondes gravitationnelles A. Nitz, M. Schäfer, B. Gadre, G. Davies, Groupe de travail “analyses de données” V. Villa, T. Dent, I. Harry, L. Xiao 30 Nov 2020 Low-latency analyses of ground-based detector data ● Astrophysics ○ Rapid discovery and characterization of transients ○ Plan and start observations in other sectors ○ Examples: GW170817, GW190521 ● Detector characterization https://emfollow.docs.ligo.org ○ See what the detector is doing on an hourly time scale ● Searches targeting generic transients ○ Identify and react to problems ○ Coherent WaveBurst as soon as possible ○ Maximize the amount of “good” data ● Searches targeting for later (archival) analyses compact binary mergers ○ Improve the detectors for future runs ○ GstLAL, MBTA, PyCBC Live, SPIIR 2 Realtime frequency-domain matched filtering Data Template waveform Signal-to-noise Threshold and Candidates ratio time series find local maxima Noise power spectral density 3 Coincidence and multi-detector significance 1. Perform matched filtering and generate single-detector candidates 2. Find two-detector coincidences, calculate FARs via time slides 3. Choose most significant coincidence 4. Followup with remaining detectors, update FAR 4 Impact and mitigation of instrumental transients SNR time series for a simulated signal in pure Gaussian noise, ● Increased rate of non-GW candidates with a glitch and after gating the glitch ○ Signal-based vetoes to down-rank the candidates ○ Astrophysical prior for the relative arrival times, phases and amplitudes at the different detectors ● Corruption of signals ○ “Gating” – replace transient with zeroes 5 Sensitivity to binary neutron star mergers 1.
    [Show full text]
  • Study of the Gravitational Waveform Created by a Binary Neutron Star (BNS) Merger
    École Normale Supérieure Paris-Saclay 'T EM LIGO Laboratory M1 PHYTEM Internship report Study of the gravitational waveform created by a Binary Neutron Star (BNS) merger Artist’s view of the merging of two neutron stars from [1] Diane INDELICATO Supervisor : Alan WEINSTEIN April 12 - July 27 2018 Summary 1 Introduction 3 1.1 Gravitational waves . .3 1.1.1 What are gravitational waves ? . .3 1.1.2 Detecting gravitational waves . .3 1.1.3 Analyzing the signal . .4 1.2 LIGO Laboratory and detectors . .5 1.3 Gravitational waves from a Binary Neutron Star system . .5 1.3.1 Differences from a Binary Black Hole system . .5 1.3.2 Identifying the equation of state of nuclear matter . .6 2 Work accomplished 6 2.1 Solving the TOV equations . .7 2.1.1 For a pure ideal neutron gas . .7 2.1.2 More complicated models . .8 2.2 Generating the waveform . .9 2.2.1 TaylorF2 . 10 2.2.2 IMRPhenomPv2_NRTidal . 10 2.3 Analyzing the waveform . 12 2.3.1 TaylorF2 . 12 2.3.2 IMRPhenomPv2_NRTidal . 13 2.4 Parameter estimation . 13 3 Teaching 14 3.1 Outreach . 14 3.2 Helping other interns . 15 4 Conclusion 16 2 1 Introduction . The study of gravitational waves can help us learn more about the universe. Their existence was conjectured by Henri Poincaré in 1905 [2], then predicted by Albert Einstein in 1916 as a consequence of his general theory of relativity [3], and they were later observed indirectly [4]. But their first direct observation was in September 2015 [5].
    [Show full text]
  • How We Searched for Merging Black Holes and Found Gw150914
    HOW WE SEARCHED FOR MERGING BLACK HOLES AND FOUND GW150914 On September 14 2015 the first direct observation of a gravitational wave was made at the LIGO Hanford and Livingston observatories. This event, named GW150914, was emitted by two black holes, each around 30 times the mass of the Sun, merging around a billion light years away. We have recently published a paper giving details about how GW150914 was found and confirmed in LIGO's data. The LIGO observatories consist of two perpendicular arms, each 4km in length, and use lasers travelling along the arms to measure tiny changes in their length. When a gravitational wave passes through a detector it stretches and squeezes space by a very small amount and it was by observing this effect that LIGO detected GW150914. But how does LIGO go from these very precise distance measurements to confirmation of a detection? FIGURES FROM THE PUBLICATION MATCHED FILTERING For more information on the meaning of The initial detection was made by a so-called 'burst' search, which these figures, see the original paper. searches for any signal that has the general pattern of a gravitational wave. We know that a gravitational-wave signal from two compact objects (black holes or neutron stars) merging should increase in frequency and strength as the objects approach each other; this is known as a 'chirp' signal. Instead of looking for any chirp-like signal we can also compare our data to carefully calculated templates. The shape of a chirp signal depends on the system's properties - the most important being the masses of the objects and how much they are spinning.
    [Show full text]
  • Arxiv:2105.09933V1 [Gr-Qc] 20 May 2021 Prepare the Gravitational Wave Data for Searches: Gat- Presented in Section IV
    Identification and removal of non-Gaussian noise transients for gravitational wave searches Benjamin Steltner,1, 2, a Maria Alessandra Papa,1, 3, 2, b and Heinz-Bernd Eggenstein1, 2, c 1Max Planck Institute for Gravitational Physics (Albert Einstein Institute), Callinstrasse 38, 30167 Hannover, Germany 2Leibniz Universit¨atHannover, D-30167 Hannover, Germany 3University of Wisconsin Milwaukee, 3135 N Maryland Ave, Milwaukee, WI 53211, USA (Dated: May 21, 2021) We present a new gating method to remove non-Gaussian noise transients in gravitational wave data. The method does not rely on any a-priori knowledge on the amplitude or duration of the transient events. In light of the character of the newly released LIGO O3a data, glitch-identification is particularly relevant for searches using this data. Our method preserves more data than previously achieved, while obtaining the same, if not higher, noise reduction. We achieve a ≈ 2-fold reduction in zeroed-out data with respect to the gates released by LIGO on the O3a data. We describe the method and characterise its performance. While developed in the context of searches for continuous signals, this method can be used to prepare gravitational wave data for any search. As the cadence of compact binary inspiral detections increases and the lower noise level of the instruments unveils new glitches, excising disturbances effectively, precisely, and in a timely manner, becomes more important. Our method does this. We release the source code associated with this new technique and the gates for the newly released O3 data. I. INTRODUCTION glitches compared with other widely used and publicly available gating methods, discarding significantly less While many loud gravitational wave signals have been data.
    [Show full text]
  • GW190521: a Binary Black Hole Merger with a Total Mass of 150 M {⊙}
    UC Berkeley UC Berkeley Previously Published Works Title GW190521: A Binary Black Hole Merger with a Total Mass of 150 M_{⊙}. Permalink https://escholarship.org/uc/item/7fb662gm Journal Physical review letters, 125(10) ISSN 0031-9007 Authors Abbott, R Abbott, TD Abraham, S et al. Publication Date 2020-09-01 DOI 10.1103/physrevlett.125.101102 Peer reviewed eScholarship.org Powered by the California Digital Library University of California PHYSICAL REVIEW LETTERS 125, 101102 (2020) Editors' Suggestion Featured in Physics GW190521: A Binary Black Hole Merger with a Total Mass of 150 M⊙ R. Abbott et al.* (LIGO Scientific Collaboration and Virgo Collaboration) (Received 30 May 2020; revised 19 June 2020; accepted 9 July 2020; published 2 September 2020) On May 21, 2019 at 03:02:29 UTC Advanced LIGO and Advanced Virgo observed a short duration gravitational-wave signal, GW190521, with a three-detector network signal-to-noise ratio of 14.7, and an estimated false-alarm rate of 1 in 4900 yr using a search sensitive to generic transients. If GW190521 is from a quasicircular binary inspiral, then the detected signal is consistent with the merger of two black þ21 þ17 holes with masses of 85−14 M⊙ and 66−18 M⊙ (90% credible intervals). We infer that the primary black hole mass lies within the gap produced by (pulsational) pair-instability supernova processes, with only a þ28 0.32% probability of being below 65 M⊙. We calculate the mass of the remnant to be 142−16 M⊙, which can be considered an intermediate mass black hole (IMBH).
    [Show full text]