DOCSLIB.ORG

Arxiv:2106.13785V1 [Astro-Ph.IM] 25 Jun 2021

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Arxiv:2106.13785V1 [Astro-Ph.IM] 25 Jun 2021 Inference with finite time series: Observing the gravitational Universe through windows Colm Talbot,1, 2, 3,a Eric Thrane,4, 5 Sylvia Biscoveanu,2 and Rory Smith4, 5 1LIGO Laboratory, California Institute of Technology, Pasadena, CA 91125, USA 2LIGO Laboratory, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA 3Kavli Institute for Astrophysics and Space Research, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA 4School of Physics and Astronomy, Monash University, VIC 3800, Australia 5OzGrav: The ARC Centre of Excellence for Gravitational-Wave Discovery, Clayton, VIC 3800, Australia Time series analysis is ubiquitous in many fields of science including gravitational-wave astronomy, where strain time series are analyzed to infer the nature of gravitational-wave sources, e.g., black holes and neutron stars. It is common in gravitational-wave transient studies to apply a tapered window function to reduce the effects of spectral artifacts from the sharp edges of data segments. We show that the conventional analysis of tapered data fails to take into account covariance between frequency bins, which arises for all finite time series|no matter the choice of window function. We discuss the origin of this covariance and derive a framework that models the correlation induced by the window function. We demonstrate this solution using both simulated Gaussian noise and real Advanced LIGO/Advanced Virgo data. We show that the effect of these correlations is similar in scale to widely studied systematic errors, e.g., uncertainty in detector calibration and power spectral density estimation. I. INTRODUCTION for the correlations between neighboring frequencies in- troduced by the window function applied to obtain finite Time-series analysis underpins recent advances in time series. We show how correlations between frequency gravitational-wave astronomy. The vast majority of bins arise from the fact that quasi-stationary Gaussian gravitational-wave data analysis relies on windowing, a noise processes are fundamentally described in the fre- procedure that multiplies the time-domain data segment quency domain by continuous functions, which imply by a window function that tapers off at the beginning infinite-duration time series. We derive a simple expres- and end of the segment. Analysts apply tapered win- sion for the “finite-duration” covariance matrix, which dows to mitigate two effects: (1) spectral artifacts aris- encodes the correlations naturally present in all finite ing from the Fourier transform of the data segment edges time series and identify our result as a specific basis for (Gibbs phenomena) and (2) correlations between neigh- a Karhunen-Lo`eve transformation (KLT) (see, e.g., [11]). boring frequency bins. While correlations between neigh- We show that there are practical applications where the boring frequency bins can be reduced, they are never current conventional approach of windowing data incurs eliminated. systematic errors, which though small, produce invalid Choosing a suitable window requires balancing vari- inferences when data are combined in large sets or when ous considerations including the spectral leakage from we analyze gravitational-wave events with high signal-to- instrumental lines and the low-frequency \seismic wall," noise ratio (SNR). effectiveness mitigating the Gibbs phenomenon, and the loss of signal. For a systematic study of the properties of different windows, we refer the reader to, e.g., [1,2] for theoretical introductions, and [3] for a specific discus- The remainder of this paper is organized as follows. sion in the context of gravitational-wave data analysis. In Section II, we present the formalism underlying our Once the data are windowed, they are typically analyzed framework. We derive the finite-duration covariance ma- in the frequency domain where the noise is described by trix for the analysis of finite time series. In Section III, arXiv:2106.13785v2 [astro-ph.IM] 29 Sep 2021 a power spectral density (PSD), and it is assumed that we perform a demonstration, applying our method to each frequency bin is statistically independent. However, the binary black hole merger events, GW150914 [12], this assumption is not true for a finite stretch of a longer GW170814 [13], and GW190521 [14] and contrast with noise process. results neglecting covariances. We show how the cur- The assumptions underpinning the Whittle approxi- rent conventional windowing procedure can lead to faulty mation have been thoroughly studied and many refine- inferences when many gravitational-wave measurements ments have been proposed (e.g., [4{10]). In this paper, are combined. While our demonstration uses data from we derive from first principles a formalism which accounts gravitational-wave astronomy, the framework we put for- ward is broadly applicable to all time-domain analysis. We show that the problem is fixed by using the finite- duration covariance matrix. We provide closing thoughts a [email protected] in SectionIV. 2 II. FORMALISM We note that, even if the true covariance matrix is di- agonal, a naive empirical estimate of this quantity neces- In this section, we derive a likelihood that enables sarily has some statistical uncertainty and will not generi- us to analyze time-series data characterized by station- cally be diagonal, even for a stationary Gaussian process. ary, Gaussian noise in a way that correctly takes into We emphasize that in this section we seek to derive ex- account covariance between neighboring frequency bins pressions for the true noise covariance matrix assuming µν is known. and neglect the impact of empirical esti- that arises for all finite time series. C mates. In AppendixB, we describe how we estimate µν in practice. C A. Basic notation We consider time-series data d(t) consisting of signal s(t) and noise n(t) C. Non-continuous, finite-duration noise d(t) = s(t) + n(t): (1) In practice, we only consider finite stretches of data. In In gravitational-wave observatories like LIGO [15] and this subsection, we derive the properties of finite stretches Virgo [16], n(t) is a time series of dimensionless of continuous noise. Let us consider data measured with strain (change in length per unit length). Transient sampling rate fs over data segment duration T . There gravitational-wave signals from merging binaries are are characterized by comparing the data to gravitational waveform templates h(t). N = fsT (5) B. Continuous, infinite-duration noise independent measurements. We assume that the noise has no content above half the sampling rate and so we To start, we focus on noise in the absence of signals. can probe every frequency without aliasing. In practice, The noise can be expressed in the frequency domain as applying an aggressive low-pass filter removes this high- frequency content. Z 1 n~(f) = dt e−2πiftn(t): (2) These data can be represented either in the time- −∞ domain as a real N-component time series with spacing 1=fs or in the frequency domain as a complex frequency The noise is best described as continuous because it is series d~i with fs=2 f fs=2 with spacing 1=T where defined for an arbitrary choice of frequency: with a suf- the endpoints− and zero-frequency≤ ≤ component are required ficiently long measurement, it is possible in principle to to be real. These two domains are related via the discrete achieve sufficient resolution to measuren ~(f) for any value Fourier transform [17] of f. If we assume the noise is Gaussian, the likelihood of N−1 observing a specific noise realization is characterized by 1 X d~0 = d0 e−2πijk=N : (6) a covariance matrix k f j s j=0 1 ∗ µν = n~(fµ)~n (fν ) (3) C 2 h i The frequency-domain covariance matrix for finite- the diagonal of which is equal to the noise PSD duration, non-continuous noise is: S µ = diag ( µν ) : (4) 1 D 0 ∗0E S C Cij = d~ d~ : (7) 2 i j We refer to µν as the “infinite-duration” covariance matrix. It is definedC continuously for arbitrary values Here, the angled brackets denote ensemble averages over of fµ and fν and, in the time domain, it is defined for all times: ( ; + ). Throughout, repeated indices are noise realizations. The widely-used Whittle approxima- summed over−∞ unless1 otherwise specified. In the next sub- tion assumes that data at each of the analyzed frequen- cies are independent, i.e., C is a diagonal matrix. This is section, we contrast µν (calligraphic script and greek ij indices) with the finite-durationC covariance matrix, de- generally a good approximation. However, as we show in this paper, the assumption of independence is not strictly noted Cij (no calligraphic script and roman indices), valid when analyzing a finite stretch of data, especially which is defined only for discrete frequency bins fi and when using a tapered window. fj (or equivalently, for a finite duration). If we further assume that the noise is stationary (the PSD does not We begin by defining our window function w, which vary in time), µν is diagonal. describes how we measure some segment of noise from C 3 what is, in theory, an infinite-duration noise process: 42 (a) − 1 ~0 1 X −2πi k=N dk = d w e (8) 520 44 fs − =−∞ N 1 X −2πijk=N 46 = djwje (9) 510 − | fs ij C j=0 | ~ 48 10 = (d w~)k: (10) log ∗ 500 − Frequency [Hz] w Here, j is a time-domain window function and the 50 − frequency-domain noise is now the convolution of the 490 original frequency-domain noise with the Fourier trans- formed window function. The prime denotes quantities 52 − associated with the windowed data. 480 We stress that this window function is always present 480 490 500 510 520 Frequency [Hz] in gravitational-wave data analysis problems and is de- 42 fined for all times, not just the analysis segment.
Load more

Recommended publications
  • Windowing Techniques, the Welch Method for Improvement of Power Spectrum Estimation
    Computers, Materials & Continua Tech Science Press DOI:10.32604/cmc.2021.014752 Article Windowing Techniques, the Welch Method for Improvement of Power Spectrum Estimation Dah-Jing Jwo1, *, Wei-Yeh Chang1 and I-Hua Wu2 1Department of Communications, Navigation and Control Engineering, National Taiwan Ocean University, Keelung, 202-24, Taiwan 2Innovative Navigation Technology Ltd., Kaohsiung, 801, Taiwan *Corresponding Author: Dah-Jing Jwo. Email: [email protected] Received: 01 October 2020; Accepted: 08 November 2020 Abstract: This paper revisits the characteristics of windowing techniques with various window functions involved, and successively investigates spectral leak- age mitigation utilizing the Welch method. The discrete Fourier transform (DFT) is ubiquitous in digital signal processing (DSP) for the spectrum anal- ysis and can be efciently realized by the fast Fourier transform (FFT). The sampling signal will result in distortion and thus may cause unpredictable spectral leakage in discrete spectrum when the DFT is employed. Windowing is implemented by multiplying the input signal with a window function and windowing amplitude modulates the input signal so that the spectral leakage is evened out. Therefore, windowing processing reduces the amplitude of the samples at the beginning and end of the window. In addition to selecting appropriate window functions, a pretreatment method, such as the Welch method, is effective to mitigate the spectral leakage. Due to the noise caused by imperfect, nite data, the noise reduction from Welch’s method is a desired treatment. The nonparametric Welch method is an improvement on the peri- odogram spectrum estimation method where the signal-to-noise ratio (SNR) is high and mitigates noise in the estimated power spectra in exchange for frequency resolution reduction.
    [Show full text]
  • Windowing Design and Performance Assessment for Mitigation of Spectrum Leakage
    E3S Web of Conferences 94, 03001 (2019) https://doi.org/10.1051/e3sconf/20199403001 ISGNSS 2018 Windowing Design and Performance Assessment for Mitigation of Spectrum Leakage Dah-Jing Jwo 1,*, I-Hua Wu 2 and Yi Chang1 1 Department of Communications, Navigation and Control Engineering, National Taiwan Ocean University 2 Pei-Ning Rd., Keelung 202, Taiwan 2 Bison Electronics, Inc., Neihu District, Taipei 114, Taiwan Abstract. This paper investigates the windowing design and performance assessment for mitigation of spectral leakage. A pretreatment method to reduce the spectral leakage is developed. In addition to selecting appropriate window functions, the Welch method is introduced. Windowing is implemented by multiplying the input signal with a windowing function. The periodogram technique based on Welch method is capable of providing good resolution if data length samples are selected optimally. Windowing amplitude modulates the input signal so that the spectral leakage is evened out. Thus, windowing reduces the amplitude of the samples at the beginning and end of the window, altering leakage. The influence of various window functions on the Fourier transform spectrum of the signals was discussed, and the characteristics and functions of various window functions were explained. In addition, we compared the differences in the influence of different data lengths on spectral resolution and noise levels caused by the traditional power spectrum estimation and various window-function-based Welch power spectrum estimations. 1 Introduction knowledge. Due to computer capacity and limitations on processing time, only limited time samples can actually The spectrum analysis [1-4] is a critical concern in many be processed during spectrum analysis and data military and civilian applications, such as performance processing of signals, i.e., the original signals need to be test of electric power systems and communication cut off, which induces leakage errors.
    [Show full text]
  • Stochastic Gravitational Wave Backgrounds
    Stochastic Gravitational Wave Backgrounds Nelson Christensen1;2 z 1ARTEMIS, Universit´eC^oted'Azur, Observatoire C^oted'Azur, CNRS, 06304 Nice, France 2Physics and Astronomy, Carleton College, Northfield, MN 55057, USA Abstract. A stochastic background of gravitational waves can be created by the superposition of a large number of independent sources. The physical processes occurring at the earliest moments of the universe certainly created a stochastic background that exists, at some level, today. This is analogous to the cosmic microwave background, which is an electromagnetic record of the early universe. The recent observations of gravitational waves by the Advanced LIGO and Advanced Virgo detectors imply that there is also a stochastic background that has been created by binary black hole and binary neutron star mergers over the history of the universe. Whether the stochastic background is observed directly, or upper limits placed on it in specific frequency bands, important astrophysical and cosmological statements about it can be made. This review will summarize the current state of research of the stochastic background, from the sources of these gravitational waves, to the current methods used to observe them. Keywords: stochastic gravitational wave background, cosmology, gravitational waves 1. Introduction Gravitational waves are a prediction of Albert Einstein from 1916 [1,2], a consequence of general relativity [3]. Just as an accelerated electric charge will create electromagnetic waves (light), accelerating mass will create gravitational waves. And almost exactly arXiv:1811.08797v1 [gr-qc] 21 Nov 2018 a century after their prediction, gravitational waves were directly observed [4] for the first time by Advanced LIGO [5, 6].
    [Show full text]
  • Effect of Data Gaps: Comparison of Different Spectral Analysis Methods
    Ann. Geophys., 34, 437–449, 2016 www.ann-geophys.net/34/437/2016/ doi:10.5194/angeo-34-437-2016 © Author(s) 2016. CC Attribution 3.0 License. Effect of data gaps: comparison of different spectral analysis methods Costel Munteanu1,2,3, Catalin Negrea1,4,5, Marius Echim1,6, and Kalevi Mursula2 1Institute of Space Science, Magurele, Romania 2Astronomy and Space Physics Research Unit, University of Oulu, Oulu, Finland 3Department of Physics, University of Bucharest, Magurele, Romania 4Cooperative Institute for Research in Environmental Sciences, Univ. of Colorado, Boulder, Colorado, USA 5Space Weather Prediction Center, NOAA, Boulder, Colorado, USA 6Belgian Institute of Space Aeronomy, Brussels, Belgium Correspondence to: Costel Munteanu ([email protected]) Received: 31 August 2015 – Revised: 14 March 2016 – Accepted: 28 March 2016 – Published: 19 April 2016 Abstract. In this paper we investigate quantitatively the ef- 1 Introduction fect of data gaps for four methods of estimating the ampli- tude spectrum of a time series: fast Fourier transform (FFT), Spectral analysis is a widely used tool in data analysis and discrete Fourier transform (DFT), Z transform (ZTR) and processing in most fields of science. The technique became the Lomb–Scargle algorithm (LST). We devise two tests: the very popular with the introduction of the fast Fourier trans- single-large-gap test, which can probe the effect of a sin- form algorithm which allowed for an extremely rapid com- gle data gap of varying size and the multiple-small-gaps test, putation of the Fourier Transform. In the absence of modern used to study the effect of numerous small gaps of variable supercomputers, this was not just useful, but also the only size distributed within the time series.
    [Show full text]
  • Chapter 2: the Fourier Transform Chapter 2
    “EEE305”, “EEE801 Part A”: Digital Signal Processing Chapter 2: The Fourier Transform Chapter 2 The Fourier Transform 2.1 Introduction The sampled Fourier transform of a periodic, discrete-time signal is known as the discrete Fourier transform (DFT). The DFT of such a signal allows an interpretation of the frequency domain descriptions of the given signal. 2.2 Derivation of the DFT The Fourier transform pair of a continuous-time signal is given by Equation 2.1 and Equation 2.2. ∞ − j2πft X ( f ) = ∫ x(t) ⋅e dt (2.1) −∞ ∞ j2πft x(t) = ∫ X ( f ) ⋅e df (2.2) −∞ Now consider the signal x(t) to be periodic i.e. repeating. The Fourier transform pair from above can now be represented by Equation 2.3 and Equation 2.4 respectively. T0 1 − π = ⋅ j2 kf0t X (kf 0 ) ∫ x(t) e dt (2.3) T0 0 ∞ π = ⋅ j 2 kf0t x(t) ∑ X (kf0 ) e (2.4) k=−∞ = 1 where f 0 . Now suppose that the signal x(t) is sampled N times per period, as illustrated in Figure 2.1, with a T0 sampling period of T seconds. The discrete signal can now be represented by Equation 2.5, where δ is the dirac delta impulse function and has a unit area of one. ∞ x[nT ] = ∑ x(t) ⋅δ (t − nT ) (2.5) n=−∞ Now replace x(t) in the Fourier integral of Equation 2.3, with x[nT] of Equation 2.5. The Equation can now be re- written as: ∞ T0 1 − π = ⋅δ − ⋅ jk 2 f0t X (kf0 ) ∑ ∫ x(t) (t nT ) e dt (2.6) T0 n=−∞ 0 1 for t = nT Since the dirac delta function δ ()t − nT = , Equation 2.6 can be modified even further to become 0 otherwise Equation 2.7 below: N −1 1 − π = ⋅ jk 2 f0nT X[kf 0 ] ∑ x[nT] e (2.7) N n=0 University of Newcastle upon Tyne Page 2.1 “EEE305”, “EEE801 Part A”: Digital Signal Processing Chapter 2: The Fourier Transform Sampled sequence δ(t-nT) litude p Am nT Time (t) Figure 2.1: Sampled sequence for the DFT.
    [Show full text]
  • New Physics and the Black Hole Mass Gap
    New physics and the Black Hole Mass Gap Djuna Lize Croon (TRIUMF) University of Michigan, September 2020 [email protected] | djunacroon.com GW190521 LIGO/Virgo’s biggest discovery yet: the impossible black holes Phys. Rev. Le?. 125, 101102 (2020). From R. Abbo? et al. (LIGO ScienCfic CollaboraCon and Virgo CollaboraCon), GW190521 LIGO/Virgo’s biggest discovery yet: the impossible black holes Phys. Rev. Le?. 125, 101102 (2020). From R. Abbo? et al. (LIGO ScienCfic CollaboraCon and Virgo CollaboraCon), GW190521 LIGO/Virgo’s biggest discovery yet: the impossible black holes … let’s wind back a bit “The Stellar Binary mergers in LIGO/Virgo O1+O2 Graveyard” Adapted from LIGO-Virgo, Frank Elavsky, Aaron Geller “The Stellar Binary mergers in LIGO/Virgo O1+O2 Graveyard” Mass Gap? “Mass Gap” Adapted from LIGO-Virgo, Frank Elavsky, Aaron Geller What populates the stellar graveyard? • In the LIGO/Virgo mass range: remnants of heavy, low-metallicity population-III stars • Primarily made of hydrogen (H) and helium (He) • Would have existed for z ≳ 6, M ∼ 20 − 130 M⊙ • Have not been directly observed yet (JWST target) • Collapsed into black holes in core-collapse supernova explosions. (Or did they?) • We study their evolution from the Zero-Age Helium Branch (ZAHB) <latexit sha1_base64="PYrKTlHPDF1/X3ZrAR6z/bngQPE=">AAACBHicdVDLSgMxFM34rPU16rKbYBFcDTN1StuFUHTjplDBPqAtQyZN29BMMiQZoQxduPFX3LhQxK0f4c6/MX0IKnrgwuGce7n3njBmVGnX/bBWVtfWNzYzW9ntnd29ffvgsKlEIjFpYMGEbIdIEUY5aWiqGWnHkqAoZKQVji9nfuuWSEUFv9GTmPQiNOR0QDHSRgrsXC1IuzKClE/PvYI757Vp0BV9oQM77zquXywWPeg6Z5WS75cNqZS8sleBnuPOkQdL1AP7vdsXOIkI15ghpTqeG+teiqSmmJFptpsoEiM8RkPSMZSjiKheOn9iCk+M0ocDIU1xDefq94kURUpNotB0RkiP1G9vJv7ldRI9KPdSyuNEE44XiwYJg1rAWSKwTyXBmk0MQVhScyvEIyQR1ia3rAnh61P4P2kWHM93Ktd+vnqxjCMDcuAYnAIPlEAVXIE6aAAM7sADeALP1r31aL1Yr4vWFWs5cwR+wHr7BCb+l9Y=</latexit>
    [Show full text]
  • Fourier Analysis
    FOURIER ANALYSIS Lucas Illing 2008 Contents 1 Fourier Series 2 1.1 General Introduction . 2 1.2 Discontinuous Functions . 5 1.3 Complex Fourier Series . 7 2 Fourier Transform 8 2.1 Definition . 8 2.2 The issue of convention . 11 2.3 Convolution Theorem . 12 2.4 Spectral Leakage . 13 3 Discrete Time 17 3.1 Discrete Time Fourier Transform . 17 3.2 Discrete Fourier Transform (and FFT) . 19 4 Executive Summary 20 1 1. Fourier Series 1 Fourier Series 1.1 General Introduction Consider a function f(τ) that is periodic with period T . f(τ + T ) = f(τ) (1) We may always rescale τ to make the function 2π periodic. To do so, define 2π a new independent variable t = T τ, so that f(t + 2π) = f(t) (2) So let us consider the set of all sufficiently nice functions f(t) of a real variable t that are periodic, with period 2π. Since the function is periodic we only need to consider its behavior on one interval of length 2π, e.g. on the interval (−π; π). The idea is to decompose any such function f(t) into an infinite sum, or series, of simpler functions. Following Joseph Fourier (1768-1830) consider the infinite sum of sine and cosine functions 1 a0 X f(t) = + [a cos(nt) + b sin(nt)] (3) 2 n n n=1 where the constant coefficients an and bn are called the Fourier coefficients of f. The first question one would like to answer is how to find those coefficients.
    [Show full text]
  • LIGO Listens for Gravitational Waves
    LIGO-Virgo Gravitational-Wave Findings So Far, and Current Events Peter S. Shawhan University of Maryland Physics Department and Joint Space-Science Institute University of Maryland APS Mid-Atlantic Senior Physicists Group February 19, 2020 1 LIGO = Laser Interferometer Gravitational-wave Observatory LIGO Hanford Observatory Two LIGO Observatories LIGO Hanford LIGO Livingston Observatory Plus the Virgo Observatory in Italy Having three detectors significantly improves our ability to confidently detect weak sources and to locate them in the sky 4 Improving detectors over time to reduce the noise level Initial LIGO (2002-2010) Advanced LIGO (as of 2015) 10−23 amplitude spectral density! LIGO-G1501223-v3 5 Advanced Detector Observing Runs O1 — September 12, 2015 to January 19, 2016 LIGO Hanford and Livingston O2 — November 30, 2016 to August 25, 2017 Initially, just the two LIGO observatories Virgo joined on August 1, 2017 O3 — Began April 1, 2019 Both LIGO observatories plus Virgo 6 Ten binary black hole mergers detected in O1+O2! Browse events at https://ligo.northwestern.edu/media/mass-plot/index.html 7 Highlights: Binary Black Hole Mergers Exploring the Properties of GW Events Bayesian parameter estimation: Adjust physical parameters of waveform model to see what fits the data from all detectors well Illustration by N. Cornish and T. Littenberg Get ranges of likely (“credible”) parameter values 9 Population of Merging BBHs: Masses Mass ratio (풒) consistent with 1 for all these events, but with significant uncertainty The data determines “chirp mass” best for low-mass BBHs, and total mass best for high-mass BBHs [Abbott et al.
    [Show full text]
  • Properties of the Binary Neutron Star Merger GW170817
    Properties of the Binary Neutron Star Merger GW170817 The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation Abbott, B.P. et al. “Properties of the Binary Neutron Star Merger GW170817.” Physical Review X 9, 1 (January 2019): 11001 © 2019 The Authors As Published http://dx.doi.org/10.1103/PhysRevX.9.011001 Publisher American Physical Society Version Final published version Citable link https://hdl.handle.net/1721.1/121211 Detailed Terms https://creativecommons.org/licenses/by/4.0/ PHYSICAL REVIEW X 9, 011001 (2019) Properties of the Binary Neutron Star Merger GW170817 B. P. Abbott et al.* (LIGO Scientific Collaboration and Virgo Collaboration) (Received 6 June 2018; revised manuscript received 20 September 2018; published 2 January 2019) On August 17, 2017, the Advanced LIGO and Advanced Virgo gravitational-wave detectors observed a low-mass compact binary inspiral. The initial sky localization of the source of the gravitational-wave signal, GW170817, allowed electromagnetic observatories to identify NGC 4993 as the host galaxy. In this work, we improve initial estimates of the binary’s properties, including component masses, spins, and tidal parameters, using the known source location, improved modeling, and recalibrated Virgo data. We extend the range of gravitational-wave frequencies considered down to 23 Hz, compared to 30 Hz in the initial analysis. We also compare results inferred using several signal models, which are more accurate and incorporate additional physical effects as compared to the initial analysis. We improve the localization of the gravitational-wave source to a 90% credible region of 16 deg2.
    [Show full text]
  • GW170814 : a Three-Detector Observation of Gravitational Waves from a Binary Black Hole Coalescence
    GW170814 : A three-detector observation of gravitational waves from a binary black hole coalescence The LIGO Scientific Collaboration and The Virgo Collaboration 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 +5:7 +2:8 holes are 30:5 3:0 M and 25:3 4:2 M (at the 90% credible level). The luminosity distance of the source +130 − − +0:03 is 540 210 Mpc, corresponding to a redshift of z =0:11 0:04. A network of three detectors improves the sky localization− of the source, reducing the area of the− 90% credible region from 1160 deg2 using only the two LIGO detectors to 60 deg2 using all three detectors. For the first time, we can test the nature of gravitational wave polarizations from the antenna response of the LIGO-Virgo network, thus enabling a new class of phenomenological tests of gravity. PACS numbers: INTRODUCTION in Virgo and a BBH signal only in the LIGO detectors: the three detector BBH signal model is preferred with a Bayes The era of gravitational wave (GW) astronomy began factor of more than 1600. with the detection of binary black hole (BBH) mergers, Until Advanced Virgo became operational, typical GW by the Advanced Laser Interferometer Gravitational Wave position estimates were highly uncertain compared to the Observatory (LIGO) detectors [1], during the first of the fields of view of most telescopes.
    [Show full text]
  • Deconvolution of the Fourier Spectrum
    Koninklijk Meteorologisch Instituut van Belgi¨e Institut Royal M´et´eorologique de Belgique Deconvolution of the Fourier spectrum F. De Meyer 2003 Wetenschappelijke en Publication scientifique technische publicatie et technique Nr 33 No 33 Uitgegeven door het Edit´epar KONINKLIJK METEOROLOGISCH l’INSTITUT ROYAL INSTITUUT VAN BELGIE METEOROLOGIQUE DE BELGIQUE Ringlaan 3, B -1180 Brussel Avenue Circulaire 3, B -1180 Bruxelles Verantwoordelijke uitgever: Dr. H. Malcorps Editeur responsable: Dr. H. Malcorps Koninklijk Meteorologisch Instituut van Belgi¨e Institut Royal M´et´eorologique de Belgique Deconvolution of the Fourier spectrum F. De Meyer 2003 Wetenschappelijke en Publication scientifique technische publicatie et technique Nr 33 No 33 Uitgegeven door het Edit´epar KONINKLIJK METEOROLOGISCH l’INSTITUT ROYAL INSTITUUT VAN BELGIE METEOROLOGIQUE DE BELGIQUE Ringlaan 3, B -1180 Brussel Avenue Circulaire 3, B -1180 Bruxelles Verantwoordelijke uitgever: Dr. H. Malcorps Editeur responsable: Dr. H. Malcorps 1 Abstract The problem of estimating the frequency spectrum of a real continuous function, which is measured only at a finite number of discrete times, is discussed in this tutorial review. Based on the deconvolu- tion method, the complex version of the one-dimensional CLEAN algorithm provides a simple way to remove the artefacts introduced by the sampling and the effects of missing data from the computed Fourier spectrum. The technique is very appropriate in the case of equally spaced data, as well as for data samples randomly distributed in time. An example of the application of CLEAN to a synthetic spectrum of a small number of harmonic components at discrete frequencies is shown. The case of a low signal-to-noise time sequence is also presented by illustrating how CLEAN recovers the spectrum of the declination component of the geomagnetic field, measured in the magnetic observatory of Dourbes.
    [Show full text]
  • Polarization Tests of GW170814 and GW170817 Using Waveforms Consistent with Alternative Theories of Gravity
    7th KAGRA International Workshop 2020. 12/19 Polarization tests of GW170814 and GW170817 using waveforms consistent with alternative theories of gravity Hiroki Takeda (University of Tokyo), Atsushi Nishizawa, Soichiro Morisaki Polarization Generic metric theory allows 6 polarizations. A = + , × , x, y, b, l Tensor Plus Cross Vector Vector x Plus mode Vector x mode Breathing mode Vector y Scalar Breathing Cross mode Vector y mode Longitudinal mode Longitudinal 2 Tests of GR by polarization Possible polarization modes in a specific theory. Theory Plus Cross Vector x vector y Breathing Longitudinal General Relativity Kaluza-Klein theory Brans-Dicke theory f(R) theory Bimetric theory Separating the polarization modes from detector signals. ->We can test GR by the polarization modes of the gravitational waves. 3 Detector Signal GW amplitude for polarization mode A Detector signal ̂ A ̂ hI(t, Ω) = ∑ FI (Ω)hA A Antenna pattern functions Plus Cross Vector x represent detector response. Vector y Breathing Longitudinal Detector tensor 4 Reconstruction In principle, (The number of polarizations) = (The number of detectors) % × ( ℎ" = $" ℎ% + $" ℎ× + $" ℎ( % × ( e.g. Detector =3, modes = ( + , × , b) ℎ) = $) ℎ% + $) ℎ× + $) ℎ( % × ( ℎ* = $* ℎ% + $* ℎ× + $* ℎ( Reconstruction(Inverse problem) Detector network expansion → More polarizations can be probed. 5 Motivation Bayesian model selection between GR and the theory allowing only scalar or vector by simple substitution of the antenna pattern functions. [LVC(2017)PRL, LVC(2019)PRL.] Tensor vs Vector T V log BTV > 3 (GW170814) hI = FI hT vs hI = FI hT {log BTV = 20.81 (GW170817) Tensor vs Scalar log B > 2.3 (GW170814) T vs S TS hI = FI hT hI = FI hT {log BTS = 23.09 (GW170817) 1.
    [Show full text]