1 probing the universe for gravitational waves barry c. barish caltech argonne national laboratory...
Post on 22-Dec-2015
214 views
TRANSCRIPT
1
Probing the Universe for Gravitational Waves
Barry C. BarishCaltech
Argonne National Laboratory16-Jan-04
"Colliding Black Holes"
Credit:National Center for Supercomputing Applications (NCSA)
LIGO-G030523-00-M
2
Einstein’s Theory of Gravitation
a necessary consequence of Special Relativity with its finite speed for information transfer
gravitational waves come from the acceleration of masses and propagate away from their sources as a space-time warpage at the speed of light
gravitational radiationbinary inspiral
of compact objects
3
Einstein’s Theory of Gravitationgravitational waves
0)1
(2
2
22
htc
• Using Minkowski metric, the information about space-time curvature is contained in the metric as an added term, h. In the weak field limit, the equation can be described with linear equations. If the choice of gauge is the transverse traceless gauge the formulation becomes a familiar wave equation
)/()/( czthczthh x
• The strain h takes the form of a plane wave propagating at the speed of light (c).
• Since gravity is spin 2, the waves have two components, but rotated by 450 instead of 900 from each other.
4
Detectionof
Gravitational Waves
Detectors in space
LISA
Gravitational Wave
Astrophysical Source
Terrestrial detectorsVirgo, LIGO, TAMA, GEO
AIGO
5
Frequency range for EM astronomy
Electromagnetic waves over ~16 orders of
magnitude Ultra Low Frequency radio
waves to high energy gamma rays
6
Frequency range for GW Astronomy
Gravitational waves over ~8 orders of
magnitude Terrestrial and space
detectors
Audio band
Space Terrestrial
7
Detecting a passing wave ….
Free masses
8
Detecting a passing wave ….
Interferometer
9
Interferometer Concept Laser used to measure
relative lengths of two orthogonal arms
As a wave passes, the arm lengths change in different ways….
…causing the interference
pattern to change at the photodiode
Arms in LIGO are 4km Measure difference in
length to one part in 1021 or 10-18 meters
SuspendedMasses
10
Simultaneous DetectionLIGO
3002 km
(L/c = 10 ms)
Hanford Observatory
Caltech
LivingstonObservatory
MIT
11
LIGO Livingston Observatory
12
LIGO Hanford Observatory
13
LIGO Facilitiesbeam tube enclosure
• minimal enclosure
• reinforced concrete
• no services
14
LIGObeam tube
LIGO beam tube under construction in January 1998
65 ft spiral welded sections
girth welded in portable clean room in the field
1.2 m diameter - 3mm stainless50 km of weld
15
Vacuum Chambersvibration isolation systems
» Reduce in-band seismic motion by 4 - 6 orders of magnitude» Compensate for microseism at 0.15 Hz by a factor of ten» Compensate (partially) for Earth tides
16
Seismic Isolation springs and masses
ConstrainedLayer
damped spring
17
LIGOvacuum equipment
18
Seismic Isolationsuspension system
• support structure is welded tubular stainless steel • suspension wire is 0.31 mm diameter steel music wire
• fundamental violin mode frequency of 340 Hz
suspension assembly for a core optic
19
LIGO Opticsfused silica
Caltech data CSIRO data
Surface uniformity < 1 nm rms Scatter < 50 ppm Absorption < 2 ppm ROC matched < 3% Internal mode Q’s > 2 x 106
20
Core Optics installation and
alignment
21
LIGO Commissioning and Science Timeline
22
Lock Acquisition
23
An earthquake occurred, starting at UTC 17:38.
From electronic logbook 2-Jan-02
Detecting Earthquakes
24
Detecting the Earth Tides
Sun and Moon
Eric MorgensonCaltech Sophomore
25
Tidal Compensation DataTidal evaluation 21-hour locked section of S1 data
Residual signal on voice coils
Predicted tides
Residual signal on laser
Feedforward
Feedback
26
Controlling angular degrees of freedom
27
Interferometer Noise Limits
Thermal (Brownian)
Noise
LASER
test mass (mirror)
Beamsplitter
Residual gas scattering
Wavelength & amplitude fluctuations photodiode
Seismic Noise
Quantum Noise
"Shot" noise
Radiation pressure
28
What Limits LIGO Sensitivity? Seismic noise limits low
frequencies
Thermal Noise limits middle frequencies
Quantum nature of light (Shot Noise) limits high frequencies
Technical issues - alignment, electronics, acoustics, etc limit us before we reach these design goals
29
LIGO Sensitivity EvolutionHanford 4km Interferometer
Dec 01
Nov 03
30
Science Runs
S2 ~ 0.9Mpc
S1 ~ 100 kpc
E8 ~ 5 kpc
NN Binary Inspiral Range
S3 ~ 3 Mpc
Design~ 18 Mpc
A Measure of Progress
Milky WayAndromedaVirgo Cluster
31
Best Performance to Date ….
Range ~ 6 Mpc
32
Astrophysical Sourcessignatures
Compact binary inspiral: “chirps”» NS-NS waveforms are well described» BH-BH need better waveforms » search technique: matched templates
Supernovae / GRBs: “bursts” » burst signals in coincidence with signals in
electromagnetic radiation » prompt alarm (~ one hour) with neutrino
detectors
Pulsars in our galaxy: “periodic”» search for observed neutron stars
(frequency, doppler shift)» all sky search (computing challenge)» r-modes
Cosmological Signal “stochastic background”
33
Compact binary collisions
» Neutron Star – Neutron Star
– waveforms are well described
» Black Hole – Black Hole – need better waveforms
» Search: matched templates
“chirps”
34
Template Bank
Covers desiredregion of massparam space
Calculatedbased on L1noise curve
Templatesplaced formax mismatchof = 0.03
2110 templatesSecond-orderpost-Newtonian
35
Optimal Filtering
Transform data to frequency domain : Generate template in frequency domain : Correlate, weighting by power spectral density of
noise:
)(~fh
)(~ fs
|)(|)(
~)(~ *
fSfhfs
h
|)(| tzFind maxima of over arrival time and phaseCharacterize these by signal-to-noise ratio (SNR) and effective distance
dfefSfhfs
tz tfi
h
2
0
*
|)(|)(
~)(~
4)(
Then inverse Fourier transform gives you the filter output
at all times:
frequency domain
36
Matched Filtering
37
Loudest Surviving Candidate Not NS/NS inspiral event 1 Sep 2002, 00:38:33 UTC S/N = 15.9, 2/dof = 2.2 (m1,m2) = (1.3, 1.1) Msun
What caused this? Appears to be due to
saturation of a photodiode
38
Sensitivity
Reach for S1 Data Inspiral sensitivity
Livingston: <D> = 176 kpc
Hanford: <D> = 36 kpc
Sensitive to inspirals in Milky Way, LMC & SMC
Star Population in our Galaxy Population includes Milky Way, LMC and SMC Neutron star masses in range 1-3 Msun LMC and SMC contribute ~12% of Milky Way
neutron binary inspirals
39
Results of Inspiral Search
Upper limit binary neutron starcoalescence rate
LIGO S1 DataR < 160 / yr / MWEG
Previous observational limits» Japanese TAMA R < 30,000 / yr / MWEG » Caltech 40m R < 4,000 / yr / MWEG
Theoretical prediction R < 2 x 10-5 / yr / MWEG
Detectable Range of S2 data will reach Andromeda!
40
Astrophysical Sourcessignatures
Compact binary inspiral: “chirps”» NS-NS waveforms are well described» BH-BH need better waveforms » search technique: matched templates
Supernovae / GRBs: “bursts” » burst signals in coincidence with signals in
electromagnetic radiation » prompt alarm (~ one hour) with neutrino
detectors
Pulsars in our galaxy: “periodic”» search for observed neutron stars
(frequency, doppler shift)» all sky search (computing challenge)» r-modes
Cosmological Signal “stochastic background”
41
Detection of Burst Sources Known sources -- Supernovae & GRBs
» Coincidence with observed electromagnetic observations.
» No close supernovae occurred during the first science run» Second science run – We are analyzing the recent very bright and close GRB030329
NO RESULT YET
Unknown phenomena » Emission of short transients of gravitational radiation of unknown waveform (e.g. black hole mergers).
42
‘Unmodeled’ Burstssearch for waveforms from sources for which we cannot currently make an accurate prediction of the waveform shape.
GOAL
METHODS
Time-Frequency Plane Search‘TFCLUSTERS’
Pure Time-Domain Search‘SLOPE’
freq
uen
cy
time
‘Raw Data’ Time-domain high pass filter
0.125s
8Hz
43
Determination of Efficiency
Efficiency measured for ‘tfclusters’ algorithm
0 10time (ms)
amp
litu
de
0
h
To measure ourefficiency, we mustpick a waveform.
1ms Gaussian burst
44
Burst Upper Limit from S1
Upper limit in strain compared to earlier (cryogenic bar) results:
• IGEC 2001 combined bar upper limit: < 2 events per day having h=1x10-20 per Hz of burst bandwidth. For a 1kHz bandwidth, limit is < 2 events/day at h=1x10-17
• Astone et al. (2002), report a 2.2 excess of one event per day at strain level of h ~ 2x10-
18
90% confidence
Result is derived using ‘TFCLUSTERS’ algorithm
1ms gaussian bursts
45
Astrophysical Sourcessignatures
Compact binary inspiral: “chirps”» NS-NS waveforms are well described» BH-BH need better waveforms » search technique: matched templates
Supernovae / GRBs: “bursts” » burst signals in coincidence with signals in
electromagnetic radiation » prompt alarm (~ one hour) with neutrino
detectors
Pulsars in our galaxy: “periodic”» search for observed neutron stars
(frequency, doppler shift)» all sky search (computing challenge)» r-modes
Cosmological Signal “stochastic background”
46
Detection of Periodic Sources
Pulsars in our galaxy: “periodic”» search for observed neutron stars » all sky search (computing challenge)» r-modes
Frequency modulation of signal due to Earth’s motion relative to the Solar System Barycenter, intrinsic frequency changes.
Amplitude modulation due to the detector’s antenna pattern.
47
Directed searches
OBSGWh0 /TfS4.11h
NO DETECTION EXPECTED
at present sensitivities
PSR J1939+2134
1283.86 Hz
Limits of detectability for rotating NS with equatorial ellipticity =I/Izz: 10-3 , 10-4 , 10-5 @ 8.5 kpc.
Crab Pulsar
48
Two Search Methods
Frequency domain
• Best suited for large parameter space searches
• Maximum likelihood detection method + Frequentist approach
Time domain
• Best suited to target known objects, even if phase evolution is complicated
Bayesian approach
First science run --- use both pipelines for the same search for cross-checking and validation
49
The Data
hS
days
hS
hS hS
days
days
days
time behavior
50
The Data
hS
hShS
hS
Hz
Hz
Hz
Hz
frequency behavior
51
Injected signal in LLO: h = 2.83 x 10-22
MeasuredF statistic
Frequency domain
• Fourier Transforms of time series
• Detection statistic: F , maximum likelihood ratio wrt unknown parameters
• use signal injections to measure F’s pdf
• use frequentist’s approach to derive upper limit
PSR J1939+2134
52
95%
h = 2.1 x 10-21
Injected signals in GEO:h=1.5, 2.0, 2.5, 3.0 x 10-21
Data
Time domain
• time series is heterodyned
• noise is estimated
• Bayesian approach in parameter estimation: express result in terms of posterior pdf for parameters of interest
PSR J1939+2134
53
Results: Periodic Sources
No evidence of continuous wave emission from PSR J1939+2134.
Summary of 95% upper limits on h:
IFO Frequentist FDS Bayesian TDS
GEO (1.940.12)x10-21 (2.1 0.1)x10-21
LLO (2.830.31)x10-22 (1.4 0.1)x10-22
LHO-2K (4.710.50)x10-22 (2.2 0.2)x10-22
LHO-4K (6.420.72)x10-22 (2.7 0.3)x10-22
• Best previous results for PSR J1939+2134: ho < 10-
20 (Glasgow, Hough et al., 1983)
54
R
R
fIh zz
20
4
2
0 c
G8
moment of inertia tensor
gravitational ellipticity of pulsar
Assumes emission is due to deviation from axisymmetry:
Upper limit on pulsar ellipticity
h0 < 3 10-22 < 3 10-4
..
(M=1.4Msun, r=10km, R=3.6kpc)
J1939+2134
55
Multi-detector upper limits
S2 Data Run
• Performed joint coherent analysis for 28 pulsars using data from all IFOs.
• Most stringent UL is for pulsar J1629-6902 (~333 Hz) where 95% confident that h0 < 2.3x10-24.
• 95% upper limit for Crab pulsar (~ 60 Hz) is h0 < 5.1 x 10-23.
• 95% upper limit for J1939+2134 (~ 1284 Hz) is h0 < 1.3 x 10-23.
95% upper limits
56
Upper limits on ellipticity
zz
yyxx
I
II
Equatorial ellipticity:
Pulsars J0030+0451 (230 pc), J2124-3358 (250 pc), and J1024-0719 (350 pc) are the nearest three pulsars in the set and their equatorial ellipticities are all constrained to less than 10-5.
S2 upper limits
Spin-down based upper limits
57
Approaching spin-down upper limits
For Crab pulsar (B0531+21) we are still a factor of ~35 above the spin-down upper limit in S2.
Hope to reach spin-down based upper limit in S3!
Note that not all pulsars analysed are constrained due to spin-down rates; some actually appear to be spinning-up (associated with accelerations in globular cluster).
Ratio of S2 upper limits to spin-down based upper limits
58
Astrophysical Sourcessignatures
Compact binary inspiral: “chirps”» NS-NS waveforms are well described» BH-BH need better waveforms » search technique: matched templates
Supernovae / GRBs: “bursts” » burst signals in coincidence with signals in
electromagnetic radiation » prompt alarm (~ one hour) with neutrino
detectors
Pulsars in our galaxy: “periodic”» search for observed neutron stars
(frequency, doppler shift)» all sky search (computing challenge)» r-modes
Cosmological Signal “stochastic background”
59
Signals from the Early Universe
Cosmic Microwave
background
WMAP 2003
stochastic background
60
Signals from the Early Universe
Strength specified by ratio of energy density in GWs to total energy density needed to close the universe:
Detect by cross-correlating output of two GW detectors:
First LIGO Science Data
Hanford - Livingston
d(lnf)
dρ
ρ
1(f)Ω GW
criticalGW
61
Limits: Stochastic Search
Non-negligible LHO 4km-2km (H1-H2) instrumental cross-
correlation; currently being investigated.
Previous best upper limits:
» Garching-Glasgow interferometers :
» EXPLORER-NAUTILUS (cryogenic bars): 60 (907Hz)ΩGW
61.0 hrs
62.3 hrs
Tobs
GW (40Hz - 314 Hz) < 23
GW (40Hz - 314 Hz) < 72.4
90% CL Upper Limit
LHO 2km-LLO 4km
LHO 4km-LLO 4km
Interferometer Pair
5GW 103(f)Ω
62
Gravitational Waves from the Early Universe
E7
S1
S2
LIGO
Adv LIGO
results
projected
63
Advanced LIGOimproved subsystems
Active Seismic
Multiple Suspensions
Sapphire Optics
Higher Power Laser
64
Advanced LIGOCubic Law for “Window” on the
Universe
Initial LIGO
Advanced LIGO
Improve amplitude sensitivity by a factor of 10x…
…number of sources goes up 1000x!
Virgo cluster
Today
65
Advanced LIGO
Enhanced Systems• laser• suspension• seismic isolation• test mass Rate
Improvement
~ 104
+narrow band
optical configuration
2007 +
66
LIGO Construction is complete & commissioning is well underway
New upper limits for neutron binary inspirals, a fast pulsar and stochastic backgrounds have been achieved from the first short science run
Sensitivity improvements are rapid -- second data run was 10x more sensitive and 4x duration and results are beginning to be reported ----- (e.g. improved pulsar searches)
Enhanced detectors will be installed in ~ 5 years, further increasing sensitivity
Direct detection should be achieved and gravitational-wave astronomy begun within the next decade !
67
Gravitational Wave Astronomy
LIGO will provide a new way to view the dynamics of the
Universe