Measurement of sub-dominant modes in a BBH population

Similar documents
Status of the untriggered burst search in S3 LIGO data

Model Waveform Accuracy Standards for Gravitational Wave Data Analysis

Data-analysis driven comparison of analytic and numerical coalescing binary waveforms: Nonspinning case

Application of Non-Harmonic Analysis for Gravitational Wave Detection

Galactic binary foregrounds Resolving, identifying and subtracting binary stars

The Simulation of Lisa and Data Analysis. E.Plagnol for the LISA_APC group

Parameter Estimation Techniques for Ultrasound Phase Reconstruction. Fatemeh Vakhshiteh Sept. 16, 2010

The NINJA2 Waveform Model Data Set

The International Pulsar Timing Array. Maura McLaughlin West Virginia University June

arxiv: v2 [gr-qc] 25 Aug 2007

The short FFT database and the peak map for the hierarchical search of periodic sources

LabVIEW Basics Peter Avitabile,Jeffrey Hodgkins Mechanical Engineering Department University of Massachusetts Lowell

CHAPTER 3 Noise in Amplitude Modulation Systems

Advanced Digital Signal Processing Part 2: Digital Processing of Continuous-Time Signals

THE FUTURE OF VIRGO BEYOND ADVANCED DETECTORS. Gianluca Gemme INFN Genova for the Virgo Collaboration

Signal segmentation and waveform characterization. Biosignal processing, S Autumn 2012

Digital Image Processing

UNIT-1. Basic signal processing operations in digital communication

Cognitive Ultra Wideband Radio

Key Questions. What is an LED and how does it work? How does a laser work? How does a semiconductor laser work? ECE 340 Lecture 29 : LEDs and Lasers

Design IV. E232 Spring 07

ECE 340 Lecture 29 : LEDs and Lasers Class Outline:

Progress on Burst Upper Limits simulation development. AJW, 4/3/01. Here s my overview of the task at hand:

Modeling Photonic Crystal Light Emitting Diode (PhCLED) Using APSYS. Copyright 2007 Crosslight Software Inc.

145M Final Exam Solutions page 1 May 11, 2010 S. Derenzo R/2. Vref. Address encoder logic. Exclusive OR. Digital output (8 bits) V 1 2 R/2

Filters. Motivating Example. Tracking a fly, oh my! Moving Weighted Average Filter. General Picture

Enhanced Blind Reception of WiGig ad Multicarrier PHY using MIMO Beam Analysis

Extraction of Musical Pitches from Recorded Music. Mark Palenik

Lecture 3: Wireless Physical Layer: Modulation Techniques. Mythili Vutukuru CS 653 Spring 2014 Jan 13, Monday

Noise simulation from data

Module 3 Greedy Strategy

Chapter-1: Introduction

Empirical Mode Decomposition (EMD) of Turner Valley Airborne Gravity Data in the Foothills of Alberta, Canada

Orthogonal Frequency Division Multiplexing (OFDM)

1. General Outline Project Proposal April 9, 2014 Kayla Esquivel and Jason Yang

Chapter 2 Direct-Sequence Systems

Presentation Outline. Advisors: Dr. In Soo Ahn Dr. Thomas L. Stewart. Team Members: Luke Vercimak Karl Weyeneth

Electromagnetic-based Side Channel Attacks

Detection of Targets in Noise and Pulse Compression Techniques

Presentation Outline. Advisors: Dr. In Soo Ahn Dr. Thomas L. Stewart. Team Members: Luke Vercimak Karl Weyeneth. Karl. Luke

GEARBOX FAULT DETECTION BY MOTOR CURRENT SIGNATURE ANALYSIS. A. R. Mohanty

THE STATISTICAL ANALYSIS OF AUDIO WATERMARKING USING THE DISCRETE WAVELETS TRANSFORM AND SINGULAR VALUE DECOMPOSITION

Carrier Phase Recovery. EE3723 : Digital Communications. Synchronization. Carrier Phase Recovery. Carrier Phase Synchronization Techniques.

Experiment begins this autumn. This Talk: Motivation for TAMA300 recycling Length sensing/control system Lock acquisition

Sound pressure level calculation methodology investigation of corona noise in AC substations

Communications I (ELCN 306)

Sampling and Reconstruction of Analog Signals

Detectors. RIT Course Number Lecture Noise

Empirical Mode Decomposition (EMD) of Turner Valley Airborne Gravity Data in the Foothills of Alberta, Canada*

ASD and Speckle Interferometry. Dave Rowe, CTO, PlaneWave Instruments

Target Echo Information Extraction

3D Distortion Measurement (DIS)

CHAPTER 4 IMPLEMENTATION OF ADALINE IN MATLAB

Prof. P. Subbarao 1, Veeravalli Balaji 2

Computational Complexity of Multiuser. Receivers in DS-CDMA Systems. Syed Rizvi. Department of Electrical & Computer Engineering

JGW-G ikagra calibration offline h(t) of ikagra

Ambipolar electronics

United States Patent [19] Adelson

Objectives. Abstract. This PRO Lesson will examine the Fast Fourier Transformation (FFT) as follows:

Cooperative Compressed Sensing for Decentralized Networks

ECE 630: Statistical Communication Theory

Compact Torus Injection for Fuelling* C. Xiao, A. Hirose, STOR-M team Plasma Physics Laboratory University of Saskatchewan

YEDITEPE UNIVERSITY ENGINEERING FACULTY COMMUNICATION SYSTEMS LABORATORY EE 354 COMMUNICATION SYSTEMS

Wavelet Speech Enhancement based on the Teager Energy Operator

Digital Signal Processing

MODULATION AND MULTIPLE ACCESS TECHNIQUES

Deep Learning Overview

Using Image Processing to Enhance Vehicle Safety

Audio Similarity. Mark Zadel MUMT 611 March 8, Audio Similarity p.1/23

Universitas Sumatera Utara

Introduction to Wavelets. For sensor data processing

Nonuniform multi level crossing for signal reconstruction

APPLICATION NOTE 3942 Optimize the Buffer Amplifier/ADC Connection

The KNIME Image Processing Extension User Manual (DRAFT )

EXPERIMENT 4 SIGNAL RECOVERY

Imaging with Wireless Sensor Networks

14 fasttest. Multitone Audio Analyzer. Multitone and Synchronous FFT Concepts

A gravitational wave is a differential strain in spacetime. Equivalently, it is a differential tidal force that can be sensed by multiple test masses.

Sub-daily signals in GPS. at semi-annual and annual periods

Prognostic Health Monitoring for Wind Turbines

Introducing COVAREP: A collaborative voice analysis repository for speech technologies

EC209 - Improving Signal-To-Noise Ratio (SNR) for Optimizing Repeatable Auditory Brainstem Responses

On-off keying, which consists of keying a sinusoidal carrier on and off with a unipolar binary signal

Lecture 9, ANIK. Data converters 1

Communication Engineering Prof. Surendra Prasad Department of Electrical Engineering Indian Institute of Technology, Delhi

Virgo and the quest for low frequency sensitivity in GW detectors. Adalberto Giazotto INFN Pisa

The quality of the transmission signal The characteristics of the transmission medium. Some type of transmission medium is required for transmission:

Optimization Techniques for Alphabet-Constrained Signal Design

Analysis Of Induction Motor With Broken Rotor Bars Using Discrete Wavelet Transform Princy P 1 and Gayathri Vijayachandran 2

FPGA implementation of DWT for Audio Watermarking Application

DESIGN AND IMPLEMENTATION OF AN ALGORITHM FOR MODULATION IDENTIFICATION OF ANALOG AND DIGITAL SIGNALS

The Fundamental Triad System

Diversity. Presented by ENG.: Ahmed Hamza Supervisor: Dr. Mohab Mangoud

A 330 GHz active terahertz imaging system for hidden objects detection

Department of Electronics and Communication Engineering 1

Diversity Techniques

Strong Noise Removal and Replacement on Seismic Data

An integrated approach of signature hole vibration monitoring and modeling for quarry vibration control

Forced Oscillation Detection Fundamentals Fundamentals of Forced Oscillation Detection

Transfer Function (TRF)

Transcription:

Measurement of sub-dominant modes in a BBH population Brendan O Brien, Filipe Da Silva, Sergey Klimenko University of Florida 1

Outline Introduction: motivation, previous work

Outline Introduction: motivation, previous work Synchronization Method: coherent summation of binary black hole (BBH) events

Outline Introduction: motivation, previous work Synchronization Method: coherent summation of binary black hole (BBH) events Results: method applied to simulated events, figures of merit

Outline Introduction: motivation, previous work Synchronization Method: coherent summation of binary black hole (BBH) events Results: method applied to simulated events, figures of merit Conclusion: upshot, future work

Outline Introduction: motivation, previous work Synchronization Method: coherent summation of binary black hole (BBH) events Results: method applied to simulated events, figures of merit Conclusion: upshot, future work

Introduction: Background and motivation GW emission from BBH systems is best understood in spherical harmonics: q=2 (SXS:BBH:0169) 2

Introduction: Background and motivation GW emission from BBH systems is best understood in spherical harmonics: (2,2) dominates for most systems q=2 (SXS:BBH:0169) 2

Introduction: Background and motivation GW emission from BBH systems is best understood in spherical harmonics: (2,2) dominates for most systems all other modes are known as subdominant modes (SDMs) q=2 (SXS:BBH:0169) 2

Introduction: Background and motivation GW emission from BBH systems is best understood in spherical harmonics: (2,2) dominates for most systems Ringdown: quasi- normal modes (QNMs) all other modes are known as subdominant modes (SDMs) q=2 (SXS:BBH:0169) 2

Introduction: Background and motivation GW emission from BBH systems is best understood in spherical harmonics: (2,2) dominates for most systems Ringdown: quasi- normal modes (QNMs) all other modes are known as subdominant modes (SDMs) study of QNMs is known as black hole spectroscopy q=2 (SXS:BBH:0169) 2

Introduction: Background and motivation GW emission from BBH systems is best understood in spherical harmonics: (2,2) dominates for most systems Ringdown: quasi- normal modes (QNMs) all other modes are known as subdominant modes (SDMs) study of QNMs is known as black hole spectroscopy SDMs and QNMs have yet to be observed q=2 (SXS:BBH:0169) 2

Introduction: Measuring SDMs/QNMs GW150914 Requires a BBH event with high signal-to-noise ratio 3

Introduction: Measuring SDMs/QNMs GW150914 Requires a BBH event with high signal-to-noise ratio During O3: possible we could detect sub-dominant modes, unlikely to measure QNMs 3

Introduction: Measuring SDMs/QNMs GW150914 Requires a BBH event with high signal-to-noise ratio During O3: possible we could detect sub-dominant modes, unlikely to measure QNMs Solution: coherent summation of multiple signals + + + = 3

Introduction: Previous work and main concept Previous work: 4

Introduction: Previous work and main concept Previous work: Synchronizing QNMs using a signal model (Yang et al. 2017, C. Da Silva et al. 2018) 4

Introduction: Previous work and main concept Previous work: Synchronizing QNMs using a signal model (Yang et al. 2017, C. Da Silva et al. 2018) This study is a simple approach to an ill-defined problem: 4

Introduction: Previous work and main concept Previous work: Synchronizing QNMs using a signal model (Yang et al. 2017, C. Da Silva et al. 2018) This study is a simple approach to an ill-defined problem: Synchronize without the use of a priori knowledge of the signal model 4

Introduction: Previous work and main concept Previous work: Synchronizing QNMs using a signal model (Yang et al. 2017, C. Da Silva et al. 2018) This study is a simple approach to an ill-defined problem: Synchronize without the use of a priori knowledge of the signal model Measure sub-dominant modes during merger and ringdown of combined BBH event 4

Introduction: Previous work and main concept Previous work: Synchronizing QNMs using a signal model (Yang et al. 2017, C. Da Silva et al. 2018) This study is a simple approach to an ill-defined problem: Synchronize without the use of a priori knowledge of the signal model Measure sub-dominant modes during merger and ringdown of combined BBH event Maximize the overlap between multiple BBH signals by applying a transformation set 4

Outline Introduction: motivation, previous work Synchronization Method: coherent summation of binary black hole (BBH) events Results: method applied to simulated events, figures of merit Conclusion: upshot, future work

Synchronization Method: Simulated waveforms Numerical relativity waveforms (SXS) including harmonic modes up to (ℓ=5,m=5) 5

Synchronization Method: Simulated waveforms Numerical relativity waveforms (SXS) including harmonic modes up to (ℓ=5,m=5) Parameter space chosen to emulate BBH events which have been detected: 5

Synchronization Method: Simulated waveforms Numerical relativity waveforms (SXS) including harmonic modes up to (ℓ=5,m=5) Parameter space chosen to emulate BBH events which have been detected: Mass ratio, total mass zero initial spin, zero eccentricity, chirp mass 5

Synchronization Method: Simulated waveforms Numerical relativity waveforms (SXS) including harmonic modes up to (ℓ=5,m=5) Parameter space chosen to emulate BBH events which have been detected: Mass ratio, total mass, chirp mass zero initial spin, zero eccentricity Randomly selected source distance, sky location, orientation confined so recovered SNR is less than SNRGW150914 5

Synchronization Method: Coherent WaveBurst (cwb) Use cwb for unmodeled signal detection and reconstruction GW150914 6

Synchronization Method: Coherent WaveBurst (cwb) Use cwb for unmodeled signal detection and reconstruction GW150914 strain data 6

Synchronization Method: Coherent WaveBurst (cwb) Use cwb for unmodeled signal detection and reconstruction GW150914 strain data wavelet transform/ pixel selection 6

Synchronization Method: Coherent WaveBurst (cwb) Use cwb for unmodeled signal detection and reconstruction GW150914 reconstructed waveform strain data wavelet transform/ pixel selection 6

Synchronization Method: Transformations Exclusively use the waveforms reconstructed by cwb to synchronize BBH signals Time-frequency representation of waveform reconstructed by cwb 7

Synchronization Method: Transformations Exclusively use the waveforms reconstructed by cwb to synchronize BBH signals Apply a transformation set: time rescale, time shift, frequency shift, phase shift Time-frequency representation of waveform reconstructed by cwb 7

Synchronization Method: Transformations Exclusively use the waveforms reconstructed by cwb to synchronize BBH signals Apply a transformation set: time rescale, time shift, frequency shift, phase shift rescale Time-frequency representation of waveform reconstructed by cwb 7

Synchronization Method: Transformations Exclusively use the waveforms reconstructed by cwb to synchronize BBH signals Apply a transformation set: time rescale, time shift, frequency shift, phase shift time shift Time-frequency representation of waveform reconstructed by cwb 7

Synchronization Method: Transformations Exclusively use the waveforms reconstructed by cwb to synchronize BBH signals Apply a transformation set: time rescale, time shift, frequency shift, phase shift frequency shift Time-frequency representation of waveform reconstructed by cwb 7

Synchronization Method: Transformations Exclusively use the waveforms reconstructed by cwb to synchronize BBH signals Apply a transformation set: time rescale, time shift, frequency shift, phase shift Time-frequency representation of waveform reconstructed by cwb 7

Synchronization Method: Transformations Exclusively use the waveforms reconstructed by cwb to synchronize BBH signals Apply a transformation set: time rescale, time shift, frequency shift, phase shift Optimize a functional using MINUIT: Time-frequency representation of waveform reconstructed by cwb 7

Synchronization Method: Maximizing waveform overlap Synchronization method applied to two reconstructed BBH events (black: q = 1.0, Mtot = 75.7M, Mc= 37.0; red: q = 1.5, Mtot = 76.5M, Mc = 42.0) 8

Synchronization Method: Maximizing waveform overlap Synchronization method applied to two reconstructed BBH events (black: q = 1.0, Mtot = 75.7M, Mc= 37.0; red: q = 1.5, Mtot = 76.5M, Mc = 42.0) 8

Synchronization Method: Maximizing waveform overlap optimization time window Synchronization method applied to two reconstructed BBH events (black: q = 1.0, Mtot = 75.7M, Mc= 37.0; red: q = 1.5, Mtot = 76.5M, Mc = 42.0) 8

Outline Introduction: motivation, previous work Synchronization Method: coherent summation of binary black hole (BBH) events Results: method applied to simulated events, figures of merit Conclusion: upshot, future work

Results: Summation of simulated waveforms Coherent stacking method applied to 16 BBH simulated waveforms 9

Results: Summation of simulated waveforms one simulated BBH waveform Coherent stacking method applied to 16 BBH simulated waveforms 9

Results: Summation of simulated waveforms one simulated BBH waveform 16 waveforms stacked Coherent stacking method applied to 16 BBH simulated waveforms 9

Results: Summation of simulated waveforms one simulated BBH waveform 16 waveforms stacked Coherent stacking method applied to 16 BBH simulated waveforms 9

Results: Summation of simulated waveforms one simulated BBH waveform 16 waveforms stacked x16 Coherent stacking method applied to 16 BBH simulated waveforms 9

Results: Summation of simulated waveforms one simulated BBH waveform 16 waveforms stacked FFT of stacked waveform Coherent stacking method applied to 16 BBH simulated waveforms 9

Results: Summation of simulated waveforms one simulated BBH waveform 16 waveforms stacked FFT of stacked waveform Coherent stacking method applied to 16 BBH simulated waveforms 9

Results: Summation of simulated waveforms one simulated BBH waveform 16 waveforms stacked FFT of stacked waveform Coherent stacking method applied to 16 BBH simulated waveforms 9

Results: Extracting sub-dominant modes 10

Results: Extracting sub-dominant modes 10

Results: Extracting sub-dominant modes cwb 10

Results: Extracting sub-dominant modes cwb x16 x16 10

Results: Extracting sub-dominant modes cwb x16 x16 synchronization method 10

Results: Extracting sub-dominant modes cwb x16 x16 synchronization method cwb 10

Results: Extracting sub-dominant modes cwb x16 x16 synchronization method cwb an aggregate waveform is reconstructed without a signal model 10

Results: Extracting sub-dominant modes cwb x16 x16 synchronization method cwb an aggregate waveform is reconstructed without a signal model 10

Results: Extracting sub-dominant modes Coherent stacking method applied to strain data with injected numerical relativity waveforms, recovered and reconstructed with cwb 11

Results: Extracting sub-dominant modes single reconstructed BBH waveform Coherent stacking method applied to strain data with injected numerical relativity waveforms, recovered and reconstructed with cwb 11

Results: Extracting sub-dominant modes single reconstructed BBH waveform reconstructed aggregate BBH waveform Coherent stacking method applied to strain data with injected numerical relativity waveforms, recovered and reconstructed with cwb 11

Results: Extracting sub-dominant modes single reconstructed BBH waveform reconstructed aggregate BBH waveform x4 Coherent stacking method applied to strain data with injected numerical relativity waveforms, recovered and reconstructed with cwb 11

Results: Extracting sub-dominant modes single reconstructed BBH waveform reconstructed aggregate BBH waveform FFT of aggregate BBH waveform Coherent stacking method applied to strain data with injected numerical relativity waveforms, recovered and reconstructed with cwb 11

Results: Extracting sub-dominant modes single reconstructed BBH waveform reconstructed aggregate BBH waveform FFT of aggregate BBH waveform Coherent stacking method applied to strain data with injected numerical relativity waveforms, recovered and reconstructed with cwb 11

Conclusion: We introduce a method for coherently stacking multiple BBH events 12

Conclusion: We introduce a method for coherently stacking multiple BBH events We demonstrate the possibility to extract sub-dominant modes from noise without using signal model 12

Conclusion: We introduce a method for coherently stacking multiple BBH events We demonstrate the possibility to extract sub-dominant modes from noise without using signal model Future work: applying this method to BBH events detected by LIGO 12

Conclusion: We introduce a method for coherently stacking multiple BBH events We demonstrate the possibility to extract sub-dominant modes from noise without using signal model Future work: applying this method to BBH events detected by LIGO Thank you! 12

Extra slides Signals with similar astrophysical parameters are easier to synchronize We use chirp mass as a measure of signal morphology Binary tree for coherent stacking of 16 BBH events

Extra slides Need a way to estimate the performance of summation procedure Introduce root-sum-squared amplitude of a waveform: We define the efficiency of stacking as: Can be calculated for entire simulated signal or individual modes

Extra slides We stack 16 BBH events by maximizing overlap of cwb reconstructed waveforms We calculate the efficiency for several harmonic modes: Harmonic Mode H1 Efficiency (%) L1 Efficiency (2,2) 97.6 97.7 (3,3) 51.1 55.0 (4,4) 45.0 47.2 Sub-dominant 47.3 50.8 (%)