Observing Gravitational Waves with Advanced LIGO Laura Nuttall on - - PowerPoint PPT Presentation

observing gravitational waves with advanced ligo
SMART_READER_LITE
LIVE PREVIEW

Observing Gravitational Waves with Advanced LIGO Laura Nuttall on - - PowerPoint PPT Presentation

Observing Gravitational Waves with Advanced LIGO Laura Nuttall on behalf of the LIGO Scientific Collaboration and Virgo Collaboration Syracuse University LIGO-G1602189 LIGO Laser Interferometer Gravitational-wave Observatory LIGO-Livingston


slide-1
SLIDE 1

Observing Gravitational Waves with Advanced LIGO

Laura Nuttall on behalf of the LIGO Scientific Collaboration and Virgo Collaboration Syracuse University LIGO-G1602189

slide-2
SLIDE 2

LIGO

Laser Interferometer Gravitational-wave Observatory

LIGO-Hanford LIGO-Livingston

2

slide-3
SLIDE 3

The design

Weiss’s 1972 design study (Weiss, Electromagnetically Coupled Broadband Gravitational Antenna, 1972 Tech. Rep. MIT)

Photodetector Beam Splitter Power Recycling Laser Source 100 kW Circulating Power

b) a)

Signal Recycling T est Mass T est Mass T est Mass T est Mass Lx = 4 km 20 W H1 L1

10 ms light travel time

Ly = 4 km

Differential changes in arm length measure strain 𝜀L = Lx - Ly = hL

3

PRL 116, 131103 (2016)

slide-4
SLIDE 4

Advanced LIGO

Higher-power laser Larger mirrors Higher finesse arm cavities Signal recycling cavity Signal recycling mirror Output mode cleaner and more …

Improvements

Comprehensive upgrade of Initial LIGO instrumentation in the same vacuum system

4

slide-5
SLIDE 5

From iLIGO to aLIGO

2007

Better Seismic Isolation

I n c r e a s e d L a s e r P

  • w

e r a n d S i g n a l R e c y c l i n g Reduced Thermal Noise

2007 ~2019

https://dcc.ligo.org/LIGO-P1000103/public

5

slide-6
SLIDE 6

Sensitivity: past, present and future

6

PRL 116, 131103 (2016)

Typical range:
 BNS ~ 70 Mpc BBH ~ 580 Mpc

slide-7
SLIDE 7

In the early hours of September 14th, 2015…

7

slide-8
SLIDE 8

GW150914

  • Observed on September 14th, 2015 at 09:40:45 UTC

  • First observed in LIGO-Livingston then 7ms later at LIGO-Hanford

  • Over 0.2 seconds the signal increases in frequency and amplitude over ~8 cycles

from 35Hz to peak amplitude at 150 Hz

8

PRL 116, 061102 (2016)

~1/200th proton radius

slide-9
SLIDE 9

GW150914

  • Observed on September 14th, 2015 at 09:40:45 UTC

  • First observed in LIGO-Livingston then 7ms later at LIGO-Hanford

  • Over 0.2 seconds the signal increases in frequency and amplitude over ~8 cycles

from 35Hz to peak amplitude at 150 Hz

PRL 116, 061102 (2016)

9

~1/200th proton radius

slide-10
SLIDE 10

The big announcement…

slide-11
SLIDE 11

GW151226

PRL 116, 241103 (2016)

11

slide-12
SLIDE 12

Making a detection

12

slide-13
SLIDE 13

Template space

a)

H1 L1

10 ms light travel time

10 ms + 5 ms for uncertainly in arrival time of weak signals

  • To detect signals from compact-object

binaries, we construct a bank of template waveforms and matched-filter the data

  • An event must match the same

waveform template in both detectors within the light travel time between sites

  • Events are assigned a detection-

statistic value that ranks their likelihood

  • f being a gravitational wave signal

ρ = s|h⇥ p h|h⇥

a|b⇥ = 4Re Z fhigh

flow

˜ a(f)˜ b(f) Sn(f) d f

100 101 102

m1 [M]

100 101

m2 [M]

|χ1| < 0.9895, |χ2| < 0.05 |χ1,2| < 0.05 |χ1,2| < 0.9895 GW150914 GW151226 LVT151012 (gstlal) LVT151012 (PyCBC)

13

slide-14
SLIDE 14

x x x x x xx x x x x x events H1 L1 zero lag or foreground

  • Determined by rate at which detector noise produces an event with

a detection statistic value equal to or higher than the candidate event

  • Background set of data is created from coincident data from

multiple detectors

  • Slide the timestamps of one detector’s data by many multiples of

0.1s and computing a new set of coincident events

Calculating Significance

Usman et al., arXiv: 1508.02357 (2015)

14

slide-15
SLIDE 15

x H1 L1 x x x x x x H1 L1 x x x xx x 0.1 s 0.1 s Time shifted data 0.1 s

Usman et al., arXiv: 1508.02357 (2015)

Calculating Significance

15

  • Determined by rate at which detector noise produces an event with

a detection statistic value equal to or higher than the candidate event

  • Background set of data is created from coincident data from

multiple detectors

  • Slide the timestamps of one detector’s data by many multiples of

0.1s and computing a new set of coincident events

slide-16
SLIDE 16

x x x x x H1 L1 x x x xx Time shifted data background events

Calculating Significance

x x

Usman et al., arXiv: 1508.02357 (2015)

16

  • Determined by rate at which detector noise produces an event with

a detection statistic value equal to or higher than the candidate event

  • Background set of data is created from coincident data from

multiple detectors

  • Slide the timestamps of one detector’s data by many multiples of

0.1s and computing a new set of coincident events

slide-17
SLIDE 17

Results from the first observing run 
 (12th Sept 2015 - 19th Jan 2016)

GW151226 LVT151012

2σ 3σ 4σ 5σ > 5σ 2σ 3σ 4σ 5σ > 5σ

8 10 12 14 16 18 20 22 24

Detection statistic ˆ ρc

10−8 10−7 10−6 10−5 10−4 10−3 10−2 10−1 100 101 102 103 104

Number of events

GW150914 Search Result Search Background Background excluding GW150914

GW151226 LVT151012

17

Abbott et al., Phys. Rev. X 6, 041015 (2016)

slide-18
SLIDE 18

Results from the first observing run

18

Abbott et al., Phys. Rev. X 6, 041015 (2016)

slide-19
SLIDE 19

Parameters of the BBH systems

Posterior probability densities of the masses, spins and distance to the three events

Abbott et al., Phys. Rev. X 6, 041015 (2016)

19

  • Lowest mass is the GW151226

secondary mass

  • Highest mass is GW150914 remnant
  • Mass ratios differ:

  • GW150914 near equal mass

  • GW151226 and LVT151012 have

support for unequal mass ratios

slide-20
SLIDE 20

Parameters of the BBH systems

Posterior probability densities of the masses, spins and distance to the three events

Abbott et al., Phys. Rev. X 6, 041015 (2016)

20

All 3 remnant black holes have spins ~0.7 as expected for the merger of similar mass black holes in a binary

slide-21
SLIDE 21

Parameters of the BBH systems

Posterior probability densities of the masses, spins and distance to the three events

Abbott et al., Phys. Rev. X 6, 041015 (2016)

21

χeff = χ1m1 + χ2m2 M

1,2 = c Gm2

1,2

~ S1,2 · ˆ L

  • For GW151226 at least one black hole has

spin magnitude > 0.2

  • Large spins parallel to angular momentum

are disfavoured

slide-22
SLIDE 22

Tests of General Relativity

  • Allowing deviations in post-Newtonian waveform model
  • Parameter deviations are reasonably consistent with zero
  • GW150914 - merger-ringdown regime occurred at best instrument
  • sensitivity. Only several cycles in LIGO sensitivity band.
  • GW151226 - many cycles in sensitivity band. Signal provides
  • pportunity to probe PN inspiral

Abbott et al., Phys. Rev. X 6, 041015 (2016)

22

slide-23
SLIDE 23

Rate of BBH mergers

  • Knowledge about BBH

merger rates depend

  • n the mass

distribution - which we don’t know very well yet!

  • Assume a few different

mass distributions

  • Infer the BBH merger

rate is in the range 9-240 Gpc-3yr-1

100 101 102 103

R (Gpc−3 yr−1)

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7

R p(R) Flat Event Based Power Law

Abbott et al. arXiv: 1606.04856 (2016)

23

slide-24
SLIDE 24

Searching for BNS and NS-BH systems

  • O1 90% upper limit

BNS rate compared to

  • ther published rates
  • Constrain the merger

rate of BNS systems with component masses of 1.35±0.13 M☉ to be less than 12,600 Gpc−3 yr−1

100 101 102 103 104

BNS Rate (Gpc−3yr−1)

aLIGO 2010 rate compendium Kim et al. pulsar Fong et al. GRB Siellez et al. GRB Coward et al. GRB Petrillo et al. GRB Jin et al. kilonova Vangioni et al. r-process de Mink & Belczynski pop syn Dominik et al. pop syn

O1 O2 O3

During O1 we looking for gravitational waves from binary neutron star (BNS) and neutron star - black hole (NS-BH) systems

Abbott et al., arXiv: 1607.07456 (2016)

24

slide-25
SLIDE 25

Searching for BNS and NS-BH systems

  • O1 90% upper limit NS-BH rate

compared to other published rates

  • Dark blue assumes 1.4-5 M☉ and

light blue 1.4-10 M☉

  • Constrain the merger rate of NS-

BH systems with BH at least 5 M☉ to be less than 3,600 Gpc

−3 yr −1

(assuming isotropic distribution of component spins)

  • O2 and O3 BNS ranges are

assumed to be 1-1.9 and 1.9-2.7 times larger than O1

During O1 we looking for gravitational waves from binary neutron star (BNS) and neutron star - black hole (NS-BH) systems

Abbott et al., arXiv: 1607.07456 (2016)

10−2 10−1 100 101 102 103 104

NSBH Rate (Gpc−3yr−1)

aLIGO 2010 rate compendium Fong et al. GRB Coward et al. GRB Petrillo et al. GRB Jin et al. kilonova Vangioni et al. r-process de Mink & Belczynski pop syn Dominik et al. pop syn

O1 O2 O3 25

slide-26
SLIDE 26

Future Network

Image Credit: Caltech/MIT/LIGO Lab

26

slide-27
SLIDE 27

Advanced LIGO's sensitivity was at the upper end of that predicted for the first observing run

Abbott et al. Living Reviews in Relativity 19, 1 (2016)

Future Sensitivity

27

slide-28
SLIDE 28

Future Rates of BBH Mergers

Abbott et al. arXiv: 1606.04856 (2016)

  • The second
  • bserving run is

starting in ~month


  • Plan is to run until

christmas followed by a break for the holidays

  • Continue running

until early spring when Virgo will join

28

1 10

hV Ti/hV TiO1

0% 20% 40% 60% 80% 100%

P(N > {2, 10, 40}|hV Ti) O2 O3

slide-29
SLIDE 29

LIGO Scientific Collaboration and Virgo Collaboration

www.ligo.org

1000+ members, 90 institutions, 16 countries

Slide: Gabriela González

29

slide-30
SLIDE 30

Extra Slides

slide-31
SLIDE 31

LIGO-G1601165

31

slide-32
SLIDE 32

Localisation

Sky localization depends on:


  • the location and orientation of the detectors

  • time delay between signal arrival at spatially separated sites

32

Abbott et al., Phys. Rev. X 6, 041015 (2016)

slide-33
SLIDE 33

Electromagnetic Follow-Up

Timeline of observations of GW150914, separated by band and relative to the time of the gravitational wave event

100 101 102 t −tmerger (days) Initial GW Burst Recovery Initial GCN Circular Updated GCN Circular (identified as BBH candidate) Final sky map Fermi GBM, LAT, MAXI, IPN, INTEGRAL (archival) Swift XRT Swift XRT Fermi LAT, MAXI BOOTES-3 MASTER Swift UVOT, SkyMapper, MASTER, TOROS, TAROT, VST, iPTF, Keck, Pan-STARRS1, KWFC, QUEST, DECam, LT, P200, Pi of the Sky, PESSTO, UH Pan-STARRS1 VST TOROS VISTA MWA ASKAP, LOFAR ASKAP, MWA VLA, LOFAR VLA, LOFARVLA

Abbott et al. ApJ 826, Number 1 (2016)

33

slide-34
SLIDE 34

The first observing run (O1)

Image Credit: LIGO

SNR = 23.7 SNR = 9.7 SNR = 13.0

34

slide-35
SLIDE 35

Abbott et al. Phys. Rev. Lett. 116, 131103 (2016)

slide-36
SLIDE 36

What does better low frequency sensitivity buy us?

Time in the sensitive frequency band for binary coalescences

Alex Nitz

Lowest viable searchable frequency for Advanced LIGO (at design sensitivity)

slide-37
SLIDE 37

GravitySpy

https://www.zooniverse.org/projects/zooniverse/gravity-spy/ Help us classify glitches!

slide-38
SLIDE 38

LIGO Magazine

slide-39
SLIDE 39

Independence of time shifts

  • Different time-shifted analyses give

independent realizations of a counting experiment for noise background events.

  • It's not the length of the template

(which can be < 0.1s) that matters, but rather the autocorrelation function (the width of the peak in the SNR - 1ms)

  • The number of background events

having ρc > 9 between consecutive time shifts, where Ci denotes the number of events in the ith time shift

  • 0.1 s time shifts are independent trials
  • f a Poisson process, even with non-

Gaussian transients in the data

slide-40
SLIDE 40

How do we know this was an astrophysical source and not something the detectors made up?

We performed every check we could think of…

  • Checked for correlated (solar weather, lightning

strikes…) and uncorrelated (seismic activity, traffic…) sources of noise

  • Checked every channel (>200,000 per detector) which

monitors the instrument behavior and environmental conditions

  • Checked the whereabouts of every person on site

(physically and remotely connected)

  • Checked for ‘injections’
  • Tracked the signal throughout the interferometer

Cannot find any instrumental cause - this signal can only be produced from two black holes colliding

slide-41
SLIDE 41

Blip Glitch

A blip transient in LIGO-Livingston strain data that produced a significant background trigger in the CBC analysis in orange, and the best-match template waveform (amplitude-scaled for comparison) in black, which exhibits a few more low-SNR cycles but otherwise quite similar morphology

  • 80
  • 70
  • 60
  • 50
  • 40
  • 30
  • 20
  • 10

10 20 30 Time [milliseconds] −1.5 −1.0 −0.5 0.0 0.5 1.0 1.5 Strain amplitude ×10−21

Band-limited h(t) during blip transient Best-match NSBH waveform Best-match GW150914 waveform