Synthetic LISA
simulating time-delay interferometry in a model LISA
(presenting) Michele Vallisneri (in absentia) John W. Armstrong
LISA Science Office, Jet Propulsion Laboratory 12/17/2003
lisa.jpl.nasa.gov
Synthetic LISA simulating time-delay interferometry in a model LISA - - PowerPoint PPT Presentation
Synthetic LISA simulating time-delay interferometry in a model LISA (presenting) Michele Vallisneri lisa.jpl.nasa.gov (in absentia) John W. Armstrong LISA Science Office, Jet Propulsion Laboratory 12/17/2003 Why Synthetic LISA? Simulate
(presenting) Michele Vallisneri (in absentia) John W. Armstrong
LISA Science Office, Jet Propulsion Laboratory 12/17/2003
lisa.jpl.nasa.gov
GWDAW 2003: Michele Vallisneri on Synthetic LISA 2 12/17/2003
include removal of laser frequency fluctuations)
(motion of the LISA array, time- and direction-dependent armlengths, causal Doppler observables, 2nd-generation TDI observables)
GWDAW 2003: Michele Vallisneri on Synthetic LISA 3 12/17/2003
spacecraft positions → photon propagation → armlengths
laser freq. fluctuations, (optical bench), proof mass, optical path
for plane waves, work from k, h+(t), hx(t) at SSB
time-delayed combinations
laser-noise and optical- bench-noise free 3 independent observables
inter-spacecraft relative frequency fluctuations intra-spacecraft relative frequency fluctuations
GWDAW 2003: Michele Vallisneri on Synthetic LISA 4 12/17/2003
spacecraft positions → photon propagation → armlengths
laser freq. fluctuations, (optical bench), proof mass, optical path
for plane waves, work from k, h+(t), hx(t) at SSB
time-delayed combinations
laser-noise and optical- bench-noise free 3 independent observables
time-delayed combinations
laser-noise and optical- bench-noise free 3 independent observables
inter-spacecraft relative frequency fluctuations intra-spacecraft relative frequency fluctuations
GWDAW 2003: Michele Vallisneri on Synthetic LISA 5 12/17/2003
spacecraft positions → photon propagation → armlengths
laser freq. fluctuations, (optical bench), proof mass, optical path
time-delayed combinations
laser-noise and optical- bench-noise free 3 independent observables
inter-spacecraft relative frequency fluctuations intra-spacecraft relative frequency fluctuations
inter-spacecraft relative frequency fluctuations intra-spacecraft relative frequency fluctuations
for plane waves, work from k, h+(t), hx(t) at SSB
photon propagation vector GW TT tensor GW buffeting of spacecraft s at emission (t-Ll) GW buffeting of spacecraft r at reception (t)
wavefront retard.; pi are spacecraft pos. Doppler shift due to GWs (Wahlquist-Estabrook response) measured for reception at spacecraft r and emission at spacecraft s (laser travels along arm l)
GWDAW 2003: Michele Vallisneri on Synthetic LISA 6 12/17/2003
spacecraft positions → photon propagation → armlengths
for plane waves, work from k, h+(t), hx(t) at SSB
time-delayed combinations
laser-noise and optical- bench-noise free 3 independent observables
inter-spacecraft relative frequency fluctuations intra-spacecraft relative frequency fluctuations
inter-spacecraft relative frequency fluctuations intra-spacecraft relative frequency fluctuations
laser freq. fluctuations, (optical bench), proof mass, optical path shot noise at sc 1 fluctuations of laser 1* (reference) at reception (t) fluctuations of laser 3 at emission (t - L2) proof-mass 1* noise Doppler shift measured for reception at spacecraft 1 and emission at spacecraft 3 (laser travels along arm 2) Doppler shift measured between optical benches on spacecraft 1 fluctuations of lasers 1 and 1* proof-mass 1 noise
GWDAW 2003: Michele Vallisneri on Synthetic LISA 7 12/17/2003
spacecraft positions → photon propagation → armlengths
for plane waves, work from k, h+(t), hx(t) at SSB
time-delayed combinations
laser-noise and optical- bench-noise free 3 independent observables
inter-spacecraft relative frequency fluctuations intra-spacecraft relative frequency fluctuations
laser freq. fluctuations, (optical bench), proof mass, optical path theory rand+digital filter
Nyquist f: πf∆t = π/2
theory rand+digital filter
Nyquist f: πf∆t = π/2
LISA noises: 18 time series (6 proof mass + 6 optical path + 6 laser)
uncorrelated white noise produced at fixed sampling time, then interpolate
correlation) and low-pass digital filter
GWDAW 2003: Michele Vallisneri on Synthetic LISA 8 12/17/2003
spacecraft positions → photon propagation → armlengths
laser freq. fluctuations, (optical bench), proof mass, optical path
for plane waves, work from k, h+(t), hx(t) at SSB
time-delayed combinations
laser-noise and optical- bench-noise free 3 independent observables
inter-spacecraft relative frequency fluctuations intra-spacecraft relative frequency fluctuations
spacecraft positions → photon propagation → armlengths
1.One Solar orbit/yr; LISA triangle spins through 360°/orbit 2.Armlengths deviate from equilateral triangle at ~ 2% 3.Armlengths are time and direction dependent
Motion hinders noise suppression (1,2,3):
armlengths
Motion improves sensitivity to GW (1):
Motion complicates GW signals (1):
(power spread through ~9 bins)
(due to relative motion, dominates for f>10-3 Hz; bandwidth ~(ΩR/c)f)
GWDAW 2003: Michele Vallisneri on Synthetic LISA 9 12/17/2003
Implements the LISA block structure as a collection of C++ classes
Defines the LISA time-evolving geometry (positions of spacecraft, armlengths) OriginalLISA: static configuration with fixed (arbitrary) armlengths ModifiedLISA: stationary configuration, rotating with T=1yr; different cw and ccw armlengths CircularRotating: spacecraft on circular, inclined orbits; cw/ccw, time-evolving, causal armlengths EccentricInclined: spacecraft on eccentric, inclined orbits; cw/ccw, time-evolving, causal armlengths NoisyLISA (use with any LISA): adds white noise to armlengths used for TDI delays ...
Return time series of noise and GW TDI
and 2nd-generation observables) TDInoise: demonstrates laser-noise subtraction TDIsignal: causal, validated vs. LISA Simulator TDIfast: cached for multiple sources (Edlund)
Defines the position and time evolution of a GW source SimpleBinary: GW from a physical monochromatic binary SimpleMonochromatic: simpler parametrization InterpolateMemory: interpolate user-provided buffers for h+, hx ...
GWDAW 2003: Michele Vallisneri on Synthetic LISA 10 12/17/2003
Return time series of noise and GW TDI
and 2nd-generation observables) TDInoise: demonstrates laser-noise subtraction TDIsignal: causal, validated vs. LISA Simulator TDIfast: cached for multiple sources (Edlund)
Defines the position and time evolution of a GW source SimpleBinary: GW from a physical monochromatic binary SimpleMonochromatic: simpler parametrization InterpolateMemory: interpolate user-provided buffers for h+, hx ...
Defines the LISA time-evolving geometry (positions of spacecraft, armlengths) OriginalLISA: static configuration with fixed (arbitrary) armlengths ModifiedLISA: stationary configuration, rotating with T=1yr; different cw and ccw armlengths CircularRotating: spacecraft on circular, inclined orbits; cw/ccw, time-evolving, causal armlengths EccentricInclined: spacecraft on eccentric, inclined orbits; cw/ccw, time-evolving, causal armlengths NoisyLISA (use with any LISA): adds white noise to armlengths used for TDI delays ...
...things to do with it right now!
GWDAW 2003: Michele Vallisneri on Synthetic LISA 11 12/17/2003
Return time series of noise and GW TDI
and 2nd-generation observables) TDInoise: demonstrates laser-noise subtraction TDIsignal: causal, validated vs. LISA Simulator TDIfast: cached for multiple sources (Edlund)
Defines the position and time evolution of a GW source SimpleBinary: GW from a physical monochromatic binary SimpleMonochromatic: simpler parametrization InterpolateMemory: interpolate user-provided buffers for h+, hx ...
Defines the LISA time-evolving geometry (positions of spacecraft, armlengths) OriginalLISA: static configuration with fixed (arbitrary) armlengths ModifiedLISA: stationary configuration, rotating with T=1yr; different cw and ccw armlengths CircularRotating: spacecraft on circular, inclined orbits; cw/ccw, time-evolving, causal armlengths EccentricInclined: spacecraft on eccentric, inclined orbits; cw/ccw, time-evolving, causal armlengths NoisyLISA (use with any LISA): adds white noise to armlengths used for TDI delays ...
...things to do with it right now!
imperfect knowledge of arms
GWDAW 2003: Michele Vallisneri on Synthetic LISA 12 12/17/2003
Return time series of noise and GW TDI
and 2nd-generation observables) TDInoise: demonstrates laser-noise subtraction TDIsignal: causal, validated vs. LISA Simulator TDIfast: cached for multiple sources (Edlund)
Defines the position and time evolution of a GW source SimpleBinary: GW from a physical monochromatic binary SimpleMonochromatic: simpler parametrization InterpolateMemory: interpolate user-provided buffers for h+, hx ...
Defines the LISA time-evolving geometry (positions of spacecraft, armlengths) OriginalLISA: static configuration with fixed (arbitrary) armlengths ModifiedLISA: stationary configuration, rotating with T=1yr; different cw and ccw armlengths CircularRotating: spacecraft on circular, inclined orbits; cw/ccw, time-evolving, causal armlengths EccentricInclined: spacecraft on eccentric, inclined orbits; cw/ccw, time-evolving, causal armlengths NoisyLISA (use with any LISA): adds white noise to armlengths used for TDI delays ...
...things to do with it right now!
GWDAW 2003: Michele Vallisneri on Synthetic LISA 13 12/17/2003
This is a Python script! Import the Synthetic LISA library (lisaswig.py, _lisaswig.so) so we can use it Create a LISA (geometry) object; use static LISA, with equal arms Armlengths (s) Create a TDI object based on our chosen LISA Noise sampling time (s) Proof mass Sn × f2 (Hz-
1)
Laser correlation (s) Print X TDI noise to disk! File name # samples requested, sampling time TDI variables to print
#!/usr/bin/python import lisaswig; unequalarmlisa = lisaswig.ModifiedLISA(15.0,16.0,17.0); unequalarmnoise = lisaswig.TDInoise(unequalarmlisa, 1.0,2.5e-48,1.0,1.8e-37,1.0,1.1e-26,1.0); lisaswig.printtdi("noise-X.txt",unequalarmnoise,1048576,1.0,"X");
GWDAW 2003: Michele Vallisneri on Synthetic LISA 14 12/17/2003
... lisaswig.printtdi("noise-a.txt",unequalarmnoise,1048576,1.0,"a"); lisaswig.printtdi("noise-z.txt",unequalarmnoise,1048576,1.0,"z"); lisaswig.printtdi("noise-E.txt",unequalarmnoise,1048576,1.0,"E");
Note laser noise subtraction! 10-
25
GWDAW 2003: Michele Vallisneri on Synthetic LISA 15 12/17/2003
noisylisa = lisaswig.NoisyLISA(originallisa,1.0,measurement noise)
1.0,2.5e-48,1.0,1.8e-37,1.0,1.1e-26,0.1) noisynoise = lisaswig.TDInoise(noisylisa,originallisa, 1.0,2.5e-48,1.0,1.8e-37,1.0,1.1e-26,0.1)
measurement noise Sn (s2 Hz-
1)
Use different LISA for noise and TDI delays
GWDAW 2003: Michele Vallisneri on Synthetic LISA 16 12/17/2003
mylisa = lisaswig.CircularRotating(0.0,0.0,1.0) mybinary = lisaswig.SimpleBinary(frequency,initial phase,inclination,amplitude, ecliptic latitude,ecliptic longitude,polarization angle) mysignal = lisaswig.TDIsignal(mylisa,mybinary) lisaswig.printtdi("signal-X.txt",mysignal,secondsperyear/16.0,16.0,"X")
ecliptic lat. = π/2 ecliptic long. = 0
f = 2 mHz T = 1 yr LISA array parameters # samples requested, sampling time
GWDAW 2003: Michele Vallisneri on Synthetic LISA 17 12/17/2003
Synthetic LISA LISA Simulator
TDI X (no noise), T = 1 yr
f = 1.94 mHz inc = 1.60 ecliptic lat. ≈ 0, long. = 0
GWDAW 2003: Michele Vallisneri on Synthetic LISA 18 12/17/2003
Hughes- Glampedakis- Kennefick integrator (C++): output h+, hx (Python ) Synthetic LISA: generate A, E, T, X GW & noise time series Matlab: compute S/Ns
GWDAW 2003: Michele Vallisneri on Synthetic LISA 19 12/17/2003
Jet Propulsion Laboratory 12/12/2003
lisa.jpl.nasa.gov
GWDAW 2003: Michele Vallisneri on Synthetic LISA 21 12/17/2003
spacecraft positions → photon propagation → armlengths
laser freq. fluctuations, (optical bench), proof mass, optical path
for plane waves, work from k, h+(t), hx(t) at SSB
time-delayed combinations
laser-noise and optical- bench-noise free 3 independent observables
inter-spacecraft relative frequency fluctuations intra-spacecraft relative frequency fluctuations
spacecraft positions → photon propagation → armlengths
dominates for f>10-3 Hz; bandwidth ~(ΩR/c)f)
GWDAW 2003: Michele Vallisneri on Synthetic LISA 22 12/17/2003
spacecraft positions → photon propagation → armlengths
for plane waves, work from k, h+(t), hx(t) at SSB
time-delayed combinations
laser-noise and optical- bench-noise free 3 independent observables
inter-spacecraft relative frequency fluctuations intra-spacecraft relative frequency fluctuations
laser freq. fluctuations, (optical bench), proof mass, optical path Proof-mass ∆f/f noise: six time series
interval ∆t
DC) theory rand+digital filter
Nyquist f: πf∆t = π/2
GWDAW 2003: Michele Vallisneri on Synthetic LISA 23 12/17/2003
spacecraft positions → photon propagation → armlengths
for plane waves, work from k, h+(t), hx(t) at SSB
time-delayed combinations
laser-noise and optical- bench-noise free 3 independent observables
inter-spacecraft relative frequency fluctuations intra-spacecraft relative frequency fluctuations
laser freq. fluctuations, (optical bench), proof mass, optical path Optical path ∆f/f noise: six time series
sampling interval ∆t
theory rand+digital filter
Nyquist f: πf∆t = π/2
GWDAW 2003: Michele Vallisneri on Synthetic LISA 24 12/17/2003
spacecraft positions → photon propagation → armlengths
for plane waves, work from k, h+(t), hx(t) at SSB
time-delayed combinations
laser-noise and optical- bench-noise free 3 independent observables
inter-spacecraft relative frequency fluctuations intra-spacecraft relative frequency fluctuations
laser freq. fluctuations, (optical bench), proof mass, optical path
Noise interpolation:
Fourier series resummation (but this requires the entire data train!)
the effective Nyquist frequency (of noise generation)
interest
theory rand+digital filter (sampling time = 1 s)
GWDAW 2003: Michele Vallisneri on Synthetic LISA 25 12/17/2003
spacecraft positions → photon propagation → armlengths
for plane waves, work from k, h+(t), hx(t) at SSB
time-delayed combinations
laser-noise and optical- bench-noise free 3 independent observables
inter-spacecraft relative frequency fluctuations intra-spacecraft relative frequency fluctuations
laser freq. fluctuations, (optical bench), proof mass, optical path
Laser ∆f/f noise: six time series
between 1 Hz and 10 Hz, Sn ≈ 1.1 10-26 Hz-1
subtraction, it is crucial to model correctly the short-time correlation structure of the noise: residual n(t) ≈ n(t + L est. error.) - n(t)
interval and then interpolating overestimates this correlation (imposing lax requirements on armlength-measurement error)
correlated noise at arbitrary times using an unequal-timestep Markov process (Ornstein- Uhlenbeck process); this underestimates the real laser-noise correlation (imposing exacting requirements on armlength-measurement error)
producing noise with a Markov chain, followed by a digital filter
GWDAW 2003: Michele Vallisneri on Synthetic LISA 26 12/17/2003
spacecraft positions → photon propagation → armlengths
laser freq. fluctuations, (optical bench), proof mass, optical path
for plane waves, work from k, h+(t), hx(t) at SSB
time-delayed combinations
laser-noise and optical- bench-noise free 3 independent observables
inter-spacecraft relative frequency fluctuations intra-spacecraft relative frequency fluctuations
for plane waves, work from k, h+(t), hx(t) at SSB For the purpose of LISA detection, plane gravitational waves are completely specified by their ecliptic coordinates (λ,β) and by their h+(t) and hx(t) time series at the solar system baricenter
structure
GWDAW 2003: Michele Vallisneri on Synthetic LISA 27 12/17/2003
Return time series of noise and GW TDI
and 2nd-generation observables) TDInoise: demonstrates laser-noise subtraction TDIsignal: causal, validated vs. LISA Simulator TDIfast: cached for multiple GW sources (Jeff)
Defines the position and time evolution of a GW source SimpleBinary: GW from a physical monochromatic binary SimpleMonochromatic: simpler parametrization InterpolateMemory: interpolate user-provided buffers for h+, hx ...
Defines the LISA time-evolving geometry (positions of spacecraft, armlengths) OriginalLISA: static configuration with fixed (arbitrary) armlengths ModifiedLISA: stationary configuration, rotating with T=1yr; different cw and ccw armlengths CircularRotating: spacecraft on circular, inclined orbits; cw/ccw, time-evolving, causal armlengths EccentricInclined: spacecraft on eccentric, inclined orbits; cw/ccw, time-evolving, causal armlengths NoisyLISA (use with any LISA): adds white noise to armlengths used for TDI delays ...
...things to do with it right now!
GWDAW 2003: Michele Vallisneri on Synthetic LISA 28 12/17/2003
equalarmlisa = lisaswig.OriginalLISA(16.6782,16.6782,16.6782); equalarmnoise = lisaswig.TDInoise(equalarmlisa, 1.0,2.5e-48,1.0,1.8e-37,1.0,1.1e-26,1.0); lisaswig.printtdi("noise-X.txt",equalarmnoise,1048576,1.0,"X");
GWDAW 2003: Michele Vallisneri on Synthetic LISA 29 12/17/2003
... lisaswig.printtdi("noise-a.txt",equalarmnoise,1048576,1.0,"z"); lisaswig.printtdi("noise-z.txt",equalarmnoise,1048576,1.0,"z"); lisaswig.printtdi("noise-E.txt",equalarmnoise,1048576,1.0,"E");
GWDAW 2003: Michele Vallisneri on Synthetic LISA 30 12/17/2003
modifiedlisa = lisaswig.ModifiedLISA(16.6782,16.6782,16.6782) modifiednoise = lisaswig.TDInoise(equalarmlisa,modifiedlisa, 1.0,2.5e-48,1.0,1.8e-37,1.0,1.1e-26,1.0e-6) lisaswig.printtdi("noise-Xm.txt",modifiednoise,samples,1.0,"X"); correctednoise = lisaswig.TDInoise(modifiedlisa, 1.0,2.5e-48,1.0,1.8e-37,1.0,1.1e-26,1.0e-6) lisaswig.printtdi("noise-Xmc.txt",correctednoise,samples,1.0,"Xm");
Use different LISA for noise and TDI delays modified TDI obs
GWDAW 2003: Michele Vallisneri on Synthetic LISA 31 12/17/2003
For 1 yr of integration, including galactic-WD confusion noise
“short LISA” (L = 1.66x106 km) baseline LISA (L = 1.66x106 km)