Surrogate gravitational waveform models
|
|
- Spencer Terry
- 5 years ago
- Views:
Transcription
1 Department of Physics, University of Maryland Theoretical astrophysics seminar Cornell University, November 30, 2013
2 Previous works - Field, Galley, Herrmann, Hesthaven, Ochsner, Tiglio (Phys. Rev. Lett. 2011) - Caudill, Field, Galley, Herrmann, Tiglio (Class. Quant. Grav. 2012) - Field, Galley, Ochsner (Phys. Rev. D 2012) - Jason Kaye (Dissertation, Brown University, 2012) - Antil, Field, Herrmann, Nochetto, Tiglio (Journal Sci. Comp. 2013) This talk (arxiv: , arxiv: ) Priscilla Canizares (Cambridge), Jonathan Gair (Cambridge), Chad Galley (Caltech), Jan Hesthaven (Brown), Jason Kaye (Brown), Manuel Tiglio (UMD)
3 Outline Introduction Likelihood computations EOB surrogates
4 Outline Introduction Likelihood computations EOB surrogates
5 Motivation Modeling of gravitational waves from compact binary coalescences and/or analysis of data represents a high dimensional challenge... Large number of intrinsic/extrinsic parametric dimensions Waveform s physical dimension Long durations with large number of cycles
6 Motivation Modeling of gravitational waves from compact binary coalescences and/or analysis of data represents a high dimensional challenge... Large number of intrinsic/extrinsic parametric dimensions Waveform s physical dimension Long durations with large number of cycles The high dimensionality of the problem is a bottleneck for most tasks... Waveform generation through solving ODEs or PDEs Parameter estimation using effective models Template based detection algorithms
7 Even easy problems are hard Consider the simple TaylorF2 frequency-domain inspiral waveform where µ labels the parameters Timings h(f ; µ) = A(µ)f 7/6 iψ(f ;µ) e Evaluation at a single parameter and frequency value takes 10 7 s Typical BNS waveform starting at 40Hz s LIGO parameter estimation study with a BNS signal days ALIGO parameter estimation study at 10Hz weeks
8 Strategy for parameterized problems Parameterized problems can be split into two phases Offline Before study/analysis begins data unknown Extra computational and human resources Train a fast to evaluate surrogate GW model Online Data is known Evaluate the model at many (data dependent) parameter values Surrogate must be accurate and fast to justify offline efforts
9 Outline Introduction Likelihood computations EOB surrogates
10 have been employed in other fields such as optimization of airplane design, Radar detection, blood flow Active area of research. Annual SIAM meeting (2013) featured 22 minisymposia with over 80 talks First GR workshop on this topic: Surrogate and reduced order models. What is it? When could it work?
11 What is a surrogate model? (I) Surrogate (Merriam-webster) : one that serves as a substitute Surrogate (This talk) : Easy-to-compute model that mimics behavior of the full, underlying model for a fixed range of the parameter/physical variable
12 What is a surrogate model? (II) Features NOT reduced physics, but reduced representations of underlying model Surrogate will converge to underlying model as representation is improved Only reproduces GW, not other quantities such as objects motion Decisions Where to sample the underlying model? How to tie together these samples? Examples Fits/interpolation Machine learning Reduced order modeling
13 What is a reduced order model? Seek a representation of the gravitational wave model h µ (t) or h µ (f ) where µ labels the parameterization, such that m h µ c i (µ)e i for as small an m as possible Leverage representation to accelerate a computation of interest i=1 Whats special about e i Application-specific basis Numerical problem s degrees of freedom = # of basis Fewer basis faster computations
14 Reduced basis (RB)-greedy algorithm Algorithm to generate the e i Definitions Kolmogorov n-width problem: From all possible basis find the minimum number to achieve accurate approximations h µ m i=1 c i(µ)e i for all µ RB-Greedy algorithm: Approximately solve n-width problem Other approaches for basis generation are possible (SVD) Key features Near-optimal basis selection (Binev 2011, DeVore 2012) Basis elements are evaluations of the physical model Parameter and basis selections carried out simultaneously
15 µ 2 µ 1 µ e 1 = h( µ 1 ) C 1 = {h µ1 } h µ P 1 (h µ ), P 1 (h µ ) = e 1 e 1, h µ e 2 e 1 C 2 = {h µ1, h µ2 }, C 1 C 2 Slide courtesy of Chad Galley
16 Example: Parameterized Heaviside (toy IMR model) Continuum: H(µ x) x [ 1, 1] µ [.2,.2] µ = -.2 µ =.2 Training set: {H(µ i x)} µ i = i i [0,..., 4000] Two representative functions
17 Example: Parameterized Heaviside (toy IMR model) µ = Select first basis (seed): H(.2 x)
18 Example: Parameterized Heaviside (toy IMR model) 1. Select first basis (seed): H(.2 x) 2. Find worst approximation: Err i = H(µ i x) ch(.2 x) µ =.2 Error
19 Example: Parameterized Heaviside (toy IMR model) 1. Select first basis (seed): H(.2 x) 2. Find worst approximation: Err i = H(µ i x) ch(.2 x) µ =.2 Error 3. Second basis: µ =.2 H(.2 x) Repeat steps 2 & 3 until an approximation threshold is achieved
20 10 0 Approximation error Greedy output (basis): µ RB = { 0.2, 0.2, , , ,... } Dimension of approximation space Greedy algorithm fails. Non-smooth w.r.t. parameter variations. If we let y(µ) = µ x then only 1 basis function H(y) needed
21 Introduction Assumption: For reduced/surrogate models to be successful there must be (smooth) structure of the solution with respect to parameter variations... M u(µ j ) u(µ) X Well chosen samples should be representative and lead to accurate approximations Figure courtesy of Jan Hesthaven
22 Checking the assumption: EOB results Example with EOBNRv2 waveforms (2,2) mode for q [1, 2], duration 12, 000M and cycles Must align at peak (Heaviside example) Fast decay of approximation error Other evidence h µ m i=1 c i(µ)e i Observed across models, regimes Observed by groups using POD/SVD Ex: Cannon et al (arxiv: ) Normalized waveform Greedy error, σ m h + h t/m Number of selected points, m
23 Waveform compression application (ex: q ) Ortho. Basis Basis # Basis # Basis # Approx: t/m h(t) c 1e 1(t) + c 2e 2(t) t/m (a) 2 term, err t/m h(t) c 1e 1(t) + c 2e 2(t) +c 3e 3(t) + c 4e 4(t) t/m (b) 4 term, err t/m h(t) c 1e 1(t) + c 2e 2(t) +c 3e 3(t) + c 4e 4(t) +c 5e 5(t) + c 6e 6(t) t/m (c) 6 term, err 10 6
24 Now what? Need a fast way to compute the coefficients c i (µ) for any parameter µ h µ (t) m i=1 c i(µ)e i (t) First we look for a convenient expression for c i (µ)...
25 Interpolation in time In principle we can find the approximation h µ (t) I m [h] = m c i (µ)e i (t) i=1 by solving an interpolation problem m c i (µ)e i (T j ) = h µ (T j ), i=1 j = 1,..., m Provided we know m good times to sample h µ (t) Naively selected points do not guarantee a solution or accuracy For application-specific basis good points are not known a-priori
26 method 1 Input: m basis {e i (t)} m i=1 Output: Nearly optimal selection of m times {T i } m i=1 These times are adapted to the problem/basis - unlike Chebyshev nodes Algorithm Sequential selection of points: {T 1 } {T 1, T 2 }... Set of points {T j } i 1 j=1 for interpolation with the first i 1 basis Extend set {T j } i 1 j=1 {T j} i j=1 to minimize the approximation error. Equivalent to selecting T i = argmax t e i (t) I i 1 [e i ](t) 1 Barrault 2004, Maday 2009, Chaturantabut 2009, Sorensen 2009
27 Example: Points for polynomial interpolation Basis are normalized Legendre polynomials the six first legendre defined polynomialson [ 1, 1] P 0 (x) = legendre polynomials 3 P 1 (x) = 2 x 5 ( P 2 (x) = 3x 2 1 ) 8. P n (x) P₀(x) P₁(x) P₂(x) P₃(x) P₄(x) P₅(x) x Q: What are the EIM points?
28 Example: Points for polynomial interpolation Basis: Residual: P 0 (x) = 1 2 P 0 (x) 0 = Point selection (no preference): x = x
29 Example: Points for polynomial interpolation Basis: P 1 (x) = 3 2 x P 1 P 1 - c 0 P 0 Residual: P 1 (x) c 0 P 0 = 3 2 x 0.5 Point selection (either ±1): x =
30 Example: Points for polynomial interpolation Basis: Residual: P 2 (x) = Point selection: 5 ( 3x 2 1 ) 8 P 2 (x) (c 0 P 0 + c 1 P 1 ) x = P 2 P 2 - c 0 P 0 - c 1 P
31 weights ω Introduction Example: 0.10 Points for polynomial interpolation Legendre ROQ 0.05 Continue 0.00 the process until # points = # basis x Chebyshev nodes EIM nodes x Distribution and approximation error properties similar to Chebyshev nodes
32 Interpolation points for EOB waveforms What are the best temporal interpolation points for an EOB-basis? Basis #2 1.5 Basis #4 1.5 Basis # t/m t/m t/m We do not know and, for example, Chebyshev nodes won t work. Identify these points by empirical interpolation method
33 Model: non-spinning EOB, q [1, 2], wave cycles (previous example) h t/m t/m 1 2 Any waveform in the above range can be recovered through its evaluation at these 5 (error 10 4 ) to 19 (error ) empirical time nodes m h µ (t) I[h µ ](t) = c i (µ)e i (t) = where c i solved the interpolation problem i m B i (t)h µ (T i ) i
34 Outline Introduction Likelihood computations EOB surrogates Introduction Likelihood computations EOB surrogates
35 Likelihood computations EOB surrogates Review Engines of surrogate and reduced order model generation RB-greedy algorithm will select points in parameter space which are most representative of the waveform family and comprises an accurate basis EIM algorithm will select most representative time samples which lead to an accurate (temporal) interpolation scheme
36 GW parameter estimation Likelihood computations EOB surrogates Time series of N data samples s(t i ) = h λ (t i ) + n(t i ), λ true parameter Ex: 32s at 4096Hz suggests N 130, 000 samples Bayesian parameter estimation to compute a posterior probability Computational cost dominated by likelihood evaluations p (s µ, H) exp ( χ 2 /2 ), χ 2 = s h µ, s h µ = f N s(f i ) h µ (f i ) 2 i=1 S n (f i ) where S n (f ) is the detector s noise curve Each particular µ 0 requires N evaluations of h µ0 (t) Overall cost scales with N
37 Likelihood computations EOB surrogates Notice χ 2 = s, s + h µ, h µ 2R s, h µ s, s computed once h µ, h µ often has simple expression (e.g. closed form) Cost dominated by s, h µ, let us focus on this piece... Plan of attack Design a custom quadrature rule tailored to our model h µ (f ) Once built the rule is reused online, e.g. when new data is available
38 Likelihood computations EOB surrogates Reduced order quadrature approximation Given data {s i } N i=1 and the empirical interpolation representation of h µ(f ) s, h µ = f N i=1 s (f i )h µ (f i ) S n (f i ) f N i=1 s (f i )I n [h µ ](f i ) S n (f i ) = m ω i h µ (F i ) i=1 where the data-specific weights ω T = E T V 1 E j := f N i=1 s (f i )e j (f i ) S n (f i ) comprise the startup cost. Here V is the interpolation matrix.
39 Likelihood computations EOB surrogates Reduced order quadrature approximation Given data {s i } N i=1 and the empirical interpolation representation of h µ(f ) s, h µ = f N i=1 s (f i )h µ (f i ) S n (f i ) f N i=1 s (f i )I n [h µ ](f i ) S n (f i ) = m ω i h µ (F i ) i=1 where the data-specific weights ω T = E T V 1 E j := f N i=1 s (f i )e j (f i ) S n (f i ) comprise the startup cost. Here V is the interpolation matrix. Properties N is a property of the experiment whereas m is a property of the model Model s approximation properties are independent of data, m N
40 Likelihood computations EOB surrogates Mock data: 4-dimensional burst signal in Gaussian noise MCMC Timing Time [s] Full ROQ About 25 times faster Anticipated speedup is (N = data samples)/(m = basis) # MCMC points Mean values accurately recovered
41 Likelihood computations EOB surrogates Roadmap to a surrogate model Goal: Fast and accurate evaluation of a parameterized underlying GW model Problem setup 1. Given a model: It is not the job of the surrogate to propose a model 2. Given a range of parameters: Surrogate will only work in this range 3. Given a ODE/PDE solver: Surrogate will not solve equations for you Attributes 1. Non-intrusive: Existing codes are complicated, we don t want to edit them
42 Test Case Introduction Likelihood computations EOB surrogates Model: Effective One Body (EOB) for inspiral-merger-ringdown of non-spinning binary black holes based on Pan et al (arxiv: ) Parameter/physical ranges: Mass ratio in the 1:2 and 9:10 ranges, for wave cycles before merger. Results will be shown for 1-2 mass ratio case Solver: Model implemented in the routine EOBNRv2 as part of the publicly available LIGO Analysis Library (LAL)
43 Likelihood computations EOB surrogates Building our surrogate model in 3 steps T 1. Greedy algorithm. Build the basis, the greedy parameter selections are here shown as red points in the horizontal axis. The rb waveforms are solved for at those parameter values. Offline t 0 q Greedy error, σ m Number of selected points, m
44 Likelihood computations EOB surrogates Building our surrogate model in 3 steps 1. Reduced basis waveforms hi RB (t), greedy parameters q i 2. Find a set of empirical time samples (blue points in the vertical axis) for accurate temporal interpolation with rb waveforms. Offline t h q t/m t/m 2
45 Likelihood computations EOB surrogates Building our surrogate model in 3 steps 1. Reduced basis waveforms hi RB (t), greedy parameters q i 2. Empirical time samples T i 3. At each empirical time build a fit for the waveform s parametric dependence only using waveforms evaluated at q i. Offline t q Polynomial errors Amplitude Phase t/m t/m Fit errors for the amplitude (relative) and phase (absolute) for the mass ratio value for which these errors are largest
46 Likelihood computations EOB surrogates Building our surrogate model in 3 steps 1. Reduced basis waveforms hi RB (t), greedy parameters q i t 2. Empirical time samples T i 3. Fits h FIT µ (T i ) at each empirical time T i 4. Online: Evaluate the surrogate at any parameter value (yellow points) by i) evaluating the fits at each empirical time which ii) specifies the full waveform via the empirical interpolant. h S µ(t) m i=1 B i (t)h FIT µ (T i ), B j (t) m i=1 hi RB (t) ( V 1) ij q where V is the interpolation matrix and {B i } are precomputed offline
47 Evaluation speedup Introduction Likelihood computations EOB surrogates Generation time (sec) Speedup factor EOB Surrogate log 2 (Sample rate) The figure of interest is a sampling rates around Hz Speedup of more than three orders of magnitude compared to solving the original EOB equations Surrogate s evaluation cost is independent of computational costs of the ODE/PDE solver
48 Likelihood computations EOB surrogates Putting the pieces together. Full surrogate model L 2 surrogate errors h q EOB Surrogate A S /A φ φ S t/m t/m Top: overlap error of the surrogate model compared to solving the original EOB equations Center: waveform for the mass ratio for which the overlap error is the largest Bottom: amplitude and phase errors for the same waveform as center plot. These errors are smaller than the errors of the model itself and those of numerical relativity simulations. Equivalent to underlying model.
49 Likelihood computations EOB surrogates What if we cannot build a dense training set? Cost of full Einstein solver prohibits a proper training set... Algorithm wish list 1. Make the most of limited number of physical model solves 2. Surrogate model should converge to physical model as solves 3. Guidance for where to solve physical model
50 Likelihood computations EOB surrogates What if we cannot build a dense training set? Cost of full Einstein solver prohibits a proper training set... Algorithm wish list 1. Make the most of limited number of physical model solves 2. Surrogate model should converge to physical model as solves 3. Guidance for where to solve physical model Status Run solver at greedy points selected for EOB model Distribution of greedy points likely more important than actual values Jonathan Blackman and Bela Szilagyi are performing runs we should know soon
51 Likelihood computations EOB surrogates Remarks Surrogate and reduced order modeling offers an exciting new approach to overcome a variety of computationally challenging problems of GW physics Proposed specific surrogate based on three offline steps. Online evaluation fast and accurate Reduced order quadratures directly offset cost of overlap inner products Future outlook Multi-mode surrogates for larger mass ratio, more cycles and with spin Potentially applicable to full numerical relativity waveforms PE accelerated by surrogates and/or reduced order quadratures Making these tools publicly available More broadly, continue to adapt tools from engineering and applied math communities for GWs
Reduced Basis in General Relativity: Select-Solve-Represent-Predict
Reduced Basis in General Relativity: Select-Solve-Represent-Predict Manuel Tiglio University of Maryland In collaboration with Scott Field, Chad Galley, Frank Herrmann, Jan Hesthaven, Evan Ochsner arxiv:1101.3765
More informationA reduced basis representation for chirp and ringdown gravitational wave templates
A reduced basis representation for chirp and ringdown gravitational wave templates 1 Sarah Caudill 2 Chad Galley 3 Frank Herrmann 1 Jan Hesthaven 4 Evan Ochsner 5 Manuel Tiglio 1 1 University of Maryland,
More informationGreedy algorithm for building a reduced basis of gravitational wave templates
Greedy algorithm for building a reduced basis of gravitational wave templates 1 Chad Galley 2 Frank Herrmann 3 Jan Hesthaven (Advisor) 4 Evan Ochsner 5 Manuel Tiglio 3 1 Brown University, Department of
More informationWhat have we learned from coalescing Black Hole binary GW150914
Stas Babak ( for LIGO and VIRGO collaboration). Albert Einstein Institute (Potsdam-Golm) What have we learned from coalescing Black Hole binary GW150914 LIGO_DCC:G1600346 PRL 116, 061102 (2016) Principles
More informationWaveform modeling for LIGO parameter estimation: status & challenges for LISA Prayush Kumar Cornell University
Waveform modeling for LIGO parameter estimation: status & challenges for LISA Prayush Kumar Cornell University The Architecture of LISA Science Analysis: Imagining the Future January 16-19, 2018 1 Outline
More informationGravitational-wave Detectability of Equal-Mass Black-hole Binaries With Aligned Spins
Intro Simulations Results Gravitational-wave Detectability of Equal-Mass Black-hole Binaries With Aligned Spins Jennifer Seiler Christian Reisswig, Sascha Husa, Luciano Rezzolla, Nils Dorband, Denis Pollney
More informationAccurate Phenomenological Waveform Models for BH Coalescence in the Frequency Domain
Accurate Phenomenological Waveform Models for BH Coalescence in the Frequency Domain Goal: synthesize inspiral-merger-ringdown models of the complete WF of Compact Binary Coalescence from pn, NR, BH perturbation
More informationCalculating Accurate Waveforms for LIGO and LISA Data Analysis
Calculating Accurate Waveforms for LIGO and LISA Data Analysis Lee Lindblom Theoretical Astrophysics, Caltech HEPL-KIPAC Seminar, Stanford 17 November 2009 Results from the Caltech/Cornell Numerical Relativity
More informationGravitational-Wave Memory Waveforms: A Generalized Approach
Gravitational-Wave Memory Waveforms: A Generalized Approach Fuhui Lin July 31, 2017 Abstract Binary black hole coalescences can produce a nonlinear memory effect besides emitting oscillatory gravitational
More informationStudying the Effects of Tidal Corrections on Parameter Estimation
Studying the Effects of Tidal Corrections on Parameter Estimation Leslie Wade Jolien Creighton, Benjamin Lackey, Evan Ochsner Center for Gravitation and Cosmology 1 Outline Background: Physical description
More informationAn eccentric binary black hole inspiral-mergerringdown gravitational waveform model from post- Newtonian and numerical relativity
An eccentric binary black hole inspiral-mergerringdown gravitational waveform model from post- Newtonian and numerical relativity Ian Hinder Max Planck Institute for Gravitational Physics (Albert Einstein
More informationGravitational-Wave Data Analysis: Lecture 2
Gravitational-Wave Data Analysis: Lecture 2 Peter S. Shawhan Gravitational Wave Astronomy Summer School May 29, 2012 Outline for Today Matched filtering in the time domain Matched filtering in the frequency
More informationModeling, detection and characterization of eccentric compact binary coalescence
Modeling, detection and characterization of eccentric compact binary coalescence E.A. Huerta et al., Phys. Rev. D 95, 024038 (2017) Daniel George, E.A. Huerta, arxiv:1701.00008 E.A. Huerta et al., in preparation
More informationSearching for Intermediate Mass Black Holes mergers
Searching for Intermediate Mass Black Holes mergers G. A. Prodi, Università di Trento and INFN for the LIGO Scientific collaboration and the Virgo collaboration special credits to Giulio Mazzolo and Chris
More informationTesting relativity with gravitational waves
Testing relativity with gravitational waves Michał Bejger (CAMK PAN) ECT* workshop New perspectives on Neutron Star Interiors Trento, 10.10.17 (DCC G1701956) Gravitation: Newton vs Einstein Absolute time
More informationSupplementary material for Black Hole Spectroscopy with Coherent Mode Stacking. Details in deriving the hypothesis test
1 Supplementary material for Black Hole Spectroscopy with Coherent Mode Stacking Details in deriving the hypothesis test The Generalized Likelihood Ratio Test (GLRT) was first presented in [1] and applied
More informationGravitational waves from binary black holes
Gravitational waves from binary black holes Hiroyuki Nakano YITP, Kyoto University DECIGO workshop, October 27, 2013 Hiroyuki Nakano Gravitational waves from binary black holes Binary black holes (BBHs)
More informationLIGO s continuing search for gravitational waves
LIGO s continuing search for gravitational waves Patrick Brady University of Wisconsin-Milwaukee LIGO Scientific Collaboration LIGO Interferometers LIGO is an interferometric detector» A laser is used
More informationSearch for Gravitational Wave Transients. Florent Robinet On behalf of the LSC and Virgo Collaborations
Search for Gravitational Wave Transients On behalf of the LSC and Virgo Collaborations 1 Gravitational Waves Gravitational waves = "ripples" in space time Weak field approximation : g = h h 1 Wave equation,
More informationThe Generalized Empirical Interpolation Method: Analysis of the convergence and application to data assimilation coupled with simulation
The Generalized Empirical Interpolation Method: Analysis of the convergence and application to data assimilation coupled with simulation Y. Maday (LJLL, I.U.F., Brown Univ.) Olga Mula (CEA and LJLL) G.
More informationGravitational waveforms for data analysis of spinning binary black holes
Gravitational waveforms for data analysis of spinning binary black holes Andrea Taracchini (Max Planck Institute for Gravitational Physics, Albert Einstein Institute Potsdam, Germany) [https://dcc.ligo.org/g1700243]
More informationAn Aperçu about Gravitational Waves and Data Analysis
Ondes Gravitationnelles, Séminaire Poincaré XXII (2016) 81 86 Séminaire Poincaré An Aperçu about Gravitational Waves and Data Analysis Eric Chassande-Mottin APC AstroParticule et Cosmologie Université
More informationThe Interpolation of Gravitational Waveforms
Brown University Division of Applied Mathematics The Interpolation of Gravitational Waveforms Author: Jason Kaye Advisor: Professor Jan S. Hesthaven May 2012 This honors thesis, written by Jason Kaye,
More informationModel Waveform and Calibration Accuracy Standards for Gravitational Wave Data Analysis
Model Waveform and Calibration Accuracy Standards for Gravitational Wave Data Analysis Lee Lindblom Theoretical Astrophysics, Caltech LIGO Seminar, Caltech 10 November 2009 Collaborators: Duncan Brown
More informationmeasuring GW polarizations beyond GR recent results and future prospects
measuring GW polarizations beyond GR recent results and future prospects LIGO Laboratory California Institute of Technology Massachusetts Institute of Technology Oct 2, 2018 Einstein Symposium Harvard
More informationNumerical Simulations of Compact Binaries
Numerical Simulations of Compact Binaries Lawrence E. Kidder Cornell University CSCAMM Workshop Matter and Electromagnetic Fields in Strong Gravity 26 August 2009, University of Maryland Cornell-Caltech
More informationThe direct detection of gravitational waves: The first discovery, and what the future might bring
The direct detection of gravitational waves: The first discovery, and what the future might bring Chris Van Den Broeck Nikhef - National Institute for Subatomic Physics Amsterdam, The Netherlands Physics
More informationFirst order BSSN formulation of Einstein s field equations
David Brown 1 Peter Diener 2 3 Jan Hesthaven 4 Frank Herrmann 3 Abdul Mroué 5 Olivier Sarbach 6 Erik Schnetter 7 Manuel Tiglio 3 Michael Wagman 4 1 North Carolina State University 2 Louisiana State University
More informationLIGO Observational Results
LIGO Observational Results Patrick Brady University of Wisconsin Milwaukee on behalf of LIGO Scientific Collaboration LIGO Science Goals Direct verification of two dramatic predictions of Einstein s general
More informationNumerical Simulations of Black Hole Spacetimes
Numerical Simulations of Black Hole Spacetimes Lee Lindblom Senior Research Associate Theoretical Astrophysics Physics Research Conference California Institute of Technology 24 May 2007 Lee Lindblom (Caltech)
More informationHow do we really look for gravitational waves?
How do we really look for gravitational waves? A tour of some applied mathematical tools used within the LIGO and Virgo collaborations Ra Inta (Texas Tech University) for the LIGO Scientific Collaboration
More informationGravitational Wave Astronomy s Next Frontier in Computation
Gravitational Wave Astronomy s Next Frontier in Computation Chad Hanna - Penn State University Penn State Physics Astronomy & Astrophysics Outline 1. Motivation 2. Gravitational waves. 3. Birth of gravitational
More informationGW150914: Observation of gravitational waves from a binary black hole merger
IL NUOVO CIMENTO 39 C (2016) 310 DOI 10.1393/ncc/i2016-16310-2 Colloquia: La Thuile 2016 GW150914: Observation of gravitational waves from a binary black hole merger F. Marion on behalf of the LIGO Scientific
More informationThe Quasi-normal Modes of Black Holes Review and Recent Updates
Ringdown Inspiral, Merger Context: Quasinormal models resulting from the merger of stellar mass BHs, and learning as much as we can from post-merger (ringdown) signals The Quasi-normal Modes of Black Holes
More informationStatus and Prospects for LIGO
Status and Prospects for LIGO Crab Pulsar St Thomas, Virgin Islands Barry C. Barish Caltech 17-March-06 LIGO Livingston, Louisiana 4 km 17-March-06 Confronting Gravity - St Thomas 2 LIGO Hanford Washington
More informationarxiv: v1 [gr-qc] 26 Sep 2011
Reduced Basis representations of multi-mode black hole ringdown gravitational waves Sarah Caudill, 1 Scott E. Field, 2, 3 Chad R. Galley, 4, 5 Frank Herrmann, 3 and Manuel Tiglio 3 1 Department of Physics
More informationTesting the strong-field dynamics of general relativity with gravitional waves
Testing the strong-field dynamics of general relativity with gravitional waves Chris Van Den Broeck National Institute for Subatomic Physics GWADW, Takayama, Japan, May 2014 Statement of the problem General
More informationResults from LIGO Searches for Binary Inspiral Gravitational Waves
Results from LIGO Searches for Binary Inspiral Gravitational Waves Peter Shawhan (LIGO Laboratory / Caltech) For the LIGO Scientific Collaboration American Physical Society April Meeting May 4, 2004 Denver,
More informationStrong field tests of Gravity using Gravitational Wave observations
Strong field tests of Gravity using Gravitational Wave observations K. G. Arun Chennai Mathematical Institute Astronomy, Cosmology & Fundamental Physics with GWs, 04 March, 2015 indig K G Arun (CMI) Strong
More informationarxiv: v2 [gr-qc] 23 Jun 2014
Investigating Binary Black Hole Mergers with Principal Component Analysis J. Clark 1, L. Cadonati 1,2, J. Healy 4, I.S. Heng 3, J. Logue 3, N. Mangini 1, L. London 4, L. Pekowsky 4, D. Shoemaker 4 arxiv:1406.5426v2
More informationarxiv:gr-qc/ v1 4 Dec 2003
Testing the LIGO Inspiral Analysis with Hardware Injections arxiv:gr-qc/0312031 v1 4 Dec 2003 Duncan A. Brown 1 for the LIGO Scientific Collaboration 1 Department of Physics, University of Wisconsin Milwaukee,
More informationDistinguishing boson stars from black holes and neutron stars with tidal interactions
Distinguishing boson stars from black holes and neutron stars with tidal interactions Noah Sennett Max Planck Institute for Gravitational Physics (Albert Einstein Institute) Department of Physics, University
More informationarxiv: v2 [gr-qc] 18 Mar 2012
Reduced Basis representations of multi-mode black hole ringdown gravitational waves Sarah Caudill, 1 Scott E. Field, 2, 3 Chad R. Galley, 4, 5 Frank Herrmann, 3 and Manuel Tiglio 3 1 Department of Physics
More informationThe Impact of Gravitational Waves: Detectability and Signatures
The Impact of Gravitational Waves: Detectability and Signatures Jenna Klemkowsky Wheeling Jesuit University Summer 2013 University of Florida International Research Experience for Undergraduates University
More informationLecture 3. Alex Nielsen Max Planck Institute for Gravitational Physics Hanover, Germany. How can we detect gravitational wave signals?
Lecture 3 Alex Nielsen Max Planck Institute for Gravitational Physics Hanover, Germany How can we detect gravitational wave signals? 2015 International Summer School on Numerical Relativity and Gravitational
More informationGW Observation of Gravitational Waves from a Binary Black Hole Merger
GW150914 Observation of Gravitational Waves from a Binary Black Hole Merger F. Marion for the LIGO Scientific Collaboration and the Virgo Collaboration Seminar at CPPM, 2016 March 3 Introduction Sources
More informationData Analysis Pipeline: The Search for Gravitational Waves in Real life
Data Analysis Pipeline: The Search for Gravitational Waves in Real life Romain Gouaty LAPP - Université de Savoie - CNRS/IN2P3 On behalf of the LIGO Scientific Collaboration and the Virgo Collaboration
More informationDetecting the next Galactic supernova
Detecting the next Galactic supernova Nicolas Arnaud on behalf of the Virgo-LAL group Now fellow at LHCb-CERN Moriond Gravitation 2003 GW supernova amplitudes Burst online filter performances Comparison
More informationGravity and the Unseen Sky. Sydney J. Chamberlin
Gravity and the Unseen Sky Sydney J. Chamberlin Outline Gravitational-wave detectors The road to detecting stochastic backgrounds with pulsar timing arrays Testing general relativity with pulsar timing
More informationSearch for compact binary systems in LIGO data
Search for compact binary systems in LIGO data Thomas Cokelaer On behalf of the LIGO Scientific Collaboration Cardiff University, U.K. LIGO-G060630-00-Z Plan 1) Overview What kind of gravitational waves
More informationANALYSIS OF BURST SIGNALS IN LIGO DATA. Irena Zivkovic, Alan Weinstein
LASER INTERFEROMETER GRAVITATIONAL WAVE OBSERVATORY LIGO CALIFORNIA INSTITUTE OF TECHNOLOGY MASSACHUSETTS INSTITUTE OF TECHNOLOGY Document Type LIGO-T010157-00-R 10/15/01 ANALYSIS OF BURST SIGNALS IN LIGO
More informationSearching for Gravitational Waves from Binary Inspirals with LIGO
Searching for Gravitational Waves from Binary Inspirals with LIGO Duncan Brown University of Wisconsin-Milwaukee for the LIGO Scientific Collaboration Inspiral Working Group LIGO-G030671-00-Z S1 Binary
More informationLISA: Probing the Universe with Gravitational Waves. Tom Prince Caltech/JPL. Laser Interferometer Space Antenna LISA
: Probing the Universe with Gravitational Waves Tom Caltech/JPL Laser Interferometer Space Antenna http://lisa.nasa.gov Gravitational Wave Astronomy is Being Born LIGO, VIRGO, GEO, TAMA 4000m, 3000m, 2000m,
More informationBlack Hole Physics via Gravitational Waves
Black Hole Physics via Gravitational Waves Image: Steve Drasco, California Polytechnic State University and MIT How to use gravitational wave observations to probe astrophysical black holes In my entire
More informationThe Dynamical Strong-Field Regime of General Relativity
The Dynamical Strong-Field Regime of General Relativity Frans Pretorius Princeton University IFT Colloquium Sao Paulo, March 30, 2016 Outline General Relativity @100 the dynamical, strong-field regime
More informationModeling gravitational waves from compact-object binaries
Modeling gravitational waves from compact-object binaries Andrea Taracchini (Max Planck Institute for Gravitational Physics, Albert Einstein Institute Potsdam, Germany) [https://dcc.ligo.org/g1602133]
More informationBlack-hole binaries in Einstein-dilaton Gauss Bonnet gravity
Black-hole binaries in Einstein-dilaton Gauss Bonnet gravity Helvi Witek Theoretical Particle Physics and Cosmology Department of Physics, King s College London work in progress with L. Gualtieri, P. Pani,
More informationBBH coalescence in the small mass ratio limit: Marrying black hole perturbation theory and PN knowledge
BBH coalescence in the small mass ratio limit: Marrying black hole perturbation theory and PN knowledge Alessandro Nagar INFN (Italy) and IHES (France) Small mass limit: Nagar Damour Tartaglia 2006 Damour
More informationAstrophysical Rates of Gravitational-Wave Compact Binary Sources in O3
Astrophysical Rates of Gravitational-Wave Compact Binary Sources in O3 Tom Dent (Albert Einstein Institute, Hannover) Chris Pankow (CIERA/Northwestern) for the LIGO and Virgo Collaborations DCC: LIGO-G1800370
More informationSearching for gravitational waves from neutron stars
Searching for gravitational waves from neutron stars Ian Jones D.I.Jones@soton.ac.uk General Relativity Group, Southampton University Ian Jones Searching for gravitational waves from neutron stars 1/23
More informationStrategy for signal classification to improve data quality for Advanced Detectors gravitational-wave searches
Strategy for signal classification to improve data quality for Advanced Detectors gravitational-wave searches E. Cuoco 1, A.Torres-Forné 2,J.A. Font 2,7, J.Powell 3, R.Lynch 4, D.Trifiró 5 2 Universitat
More informationCoalescing binary black holes in the extreme mass ratio limit
Coalescing binary black holes in the extreme mass ratio limit Alessandro Nagar Relativity and Gravitation Group, Politecnico di Torino and INFN, sez. di Torino www.polito.it/relgrav/ alessandro.nagar@polito.it
More informationSearches for Gravitational waves associated with Gamma-ray bursts
Searches for Gravitational waves associated with Gamma-ray bursts Raymond Frey University of Oregon for the LIGO Scientific Collaboration and the Virgo Collaboration 1 Current network of groundbased GW
More informationSearching for Gravitational Waves from Coalescing Binary Systems
Searching for Gravitational Waves from Coalescing Binary Systems Stephen Fairhurst Cardiff University and LIGO Scientific Collaboration 1 Outline Motivation Searching for Coalescing Binaries Latest Results
More informationGRAVITATIONAL WAVE SOURCES AND RATES FOR LISA
GRAVITATIONAL WAVE SOURCES AND RATES FOR LISA W. Z. Korth, PHZ6607, Fall 2008 Outline Introduction What is LISA? Gravitational waves Characteristics Detection (LISA design) Sources Stochastic Monochromatic
More informationLaser Interferometer Gravitational-Wave Observatory (LIGO)! A Brief Overview!
Laser Interferometer Gravitational-Wave Observatory (LIGO) A Brief Overview Sharon Brunett California Institute of Technology Pacific Research Platform Workshop October 15, 2015 Credit: AEI, CCT, LSU LIGO
More informationAstrophysics to z~10 with Gravitational Waves
Astrophysics to z~10 with Gravitational Waves Robin Stebbins U.S. LISA Project Scientist University of Maryland Physics Seminar College Park, MD 1 May 2007 My Problem Gravitational wave detection is capability-driven,
More informationPrecessing NR simulations
Precessing NR simulations Harald Pfeiffer, CITA ICTS Program on Numerical Relativity ICTS/TIFR Bengaluru June 26, 2013 Abdul Mroue (CITA) Today s goal: How does one compute these? Where are we in terms
More informationTowards Reduced Order Modeling (ROM) for Gust Simulations
Towards Reduced Order Modeling (ROM) for Gust Simulations S. Görtz, M. Ripepi DLR, Institute of Aerodynamics and Flow Technology, Braunschweig, Germany Deutscher Luft und Raumfahrtkongress 2017 5. 7. September
More informationOutline. 1. Basics of gravitational wave transient signal searches. 2. Reconstruction of signal properties
Gravitational Wave Transients state-of-the-arts: detection confidence and signal reconstruction G.A.Prodi, University of Trento and INFN, for the LIGO Scientific Collaboration and the Virgo Collaboration
More informationThe Advanced LIGO detectors at the beginning of the new gravitational wave era
The Advanced LIGO detectors at the beginning of the new gravitational wave era Lisa Barsotti MIT Kavli Institute LIGO Laboratory on behalf of the LIGO Scientific Collaboration LIGO Document G1600324 LIGO
More informationGW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral
GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral Lazzaro Claudia for the LIGO Scientific Collaboration and the Virgo Collaboration 25 October 2017 GW170817 PhysRevLett.119.161101
More informationAnalytic methods in the age of numerical relativity
Analytic methods in the age of numerical relativity vs. Marc Favata Kavli Institute for Theoretical Physics University of California, Santa Barbara Motivation: Modeling the emission of gravitational waves
More informationOverview of Gravitational Wave Observations by LIGO and Virgo
Overview of Gravitational Wave Observations by LIGO and Virgo Giovanni Andrea Prodi Virgo Group at Padova-Trento, LIGO Scientific Collaboration and Virgo Collaboration Vulcano Workshop 2016, May 23 Published
More informationNonlinear Model Reduction for Uncertainty Quantification in Large-Scale Inverse Problems
Nonlinear Model Reduction for Uncertainty Quantification in Large-Scale Inverse Problems Krzysztof Fidkowski, David Galbally*, Karen Willcox* (*MIT) Computational Aerospace Sciences Seminar Aerospace Engineering
More informationInadequacies of the Fisher Matrix in Gravitational-wave Parameter Estimation
Inadequacies of the Fisher Matrix in Gravitational-wave Parameter Estimation Carl Rodriguez Ilya Mandel Ben Farr The Problem Advanced LIGO receives a signal: Past the initial detection, all science potential
More informationSearching for Big Bang relicts with LIGO. Stefan W. Ballmer For the LSC/VIRGO collaboration GW2010, University of Minnesota October 16, 2010
Searching for Big Bang relicts with LIGO Stefan W. Ballmer For the LSC/VIRGO collaboration GW2010, University of Minnesota October 16, 2010 Outline The search for a stochastic GW background» Motivation»
More informationPerformance Comparison of K-Means and Expectation Maximization with Gaussian Mixture Models for Clustering EE6540 Final Project
Performance Comparison of K-Means and Expectation Maximization with Gaussian Mixture Models for Clustering EE6540 Final Project Devin Cornell & Sushruth Sastry May 2015 1 Abstract In this article, we explore
More informationOn the detectability of post-newtonian eects. in gravitational-wave emission of a coalescing. binary 1. Institute of Mathematics
On the detectability of post-newtonian eects in gravitational-wave emission of a coalescing binary 1 ANDRZEJ KR OLAK a KOSTAS D. KOKKOTAS b GERHARD SCH AFER c PostScript processed by the SLAC/DESY Libraries
More informationIntroduction to machine learning and pattern recognition Lecture 2 Coryn Bailer-Jones
Introduction to machine learning and pattern recognition Lecture 2 Coryn Bailer-Jones http://www.mpia.de/homes/calj/mlpr_mpia2008.html 1 1 Last week... supervised and unsupervised methods need adaptive
More informationSearching for gravitational waves
Searching for gravitational waves Matteo Barsuglia (barsuglia@apc.univ-paris7.fr) CNRS - Laboratoire Astroparticule et Cosmologie 1 The gravitational waves (GW) Perturbations of the space-time metrics
More informationS. V. Dhurandhar IUCAA Pune, India
Extended hierarchical search (EHS) for inspiraling compact binaries S. V. Dhurandhar IUCAA Pune, India * A. Sengupta, A. Lazzarini, P. Shawhan LIGO-G040556-00-R Plan of the Talk The 1-step (flat) search
More informationAstrophysical Stochastic Gravitational Waves. Jonah Kanner PHYS 798G March 27, 2007
Astrophysical Stochastic Gravitational Waves Jonah Kanner PHYS 798G March 27, 2007 Introduction Gravitational Waves come from space Require acceleration of dense mass (Think black holes and neutron stars!)
More informationBinary Black Holes. Deirdre Shoemaker Center for Relativistic Astrophysics School of Physics Georgia Tech
Binary Black Holes Deirdre Shoemaker Center for Relativistic Astrophysics School of Physics Georgia Tech NR confirmed BBH GW detections LIGO-P150914-v12 Abbott et al. 2016a, PRL 116, 061102 an orbital
More informationStochastic Spectral Approaches to Bayesian Inference
Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to
More informationThe nonlinear gravitational-wave memory in binary black hole mergers
The nonlinear gravitational-wave memory in binary black hole mergers Marc Favata Kavli Institute for Theoretical Physics University of California, Santa Barbara What is memory? Generally think of GW s
More informationGravitational Wave Memory Revisited:
Gravitational Wave Memory Revisited: Memory from binary black hole mergers Marc Favata Kavli Institute for Theoretical Physics arxiv:0811.3451 [astro-ph] and arxiv:0812.0069 [gr-qc] What is the GW memory?
More informationReduced-Order Greedy Controllability of Finite Dimensional Linear Systems. Giulia Fabrini Laura Iapichino Stefan Volkwein
Universität Konstanz Reduced-Order Greedy Controllability of Finite Dimensional Linear Systems Giulia Fabrini Laura Iapichino Stefan Volkwein Konstanzer Schriften in Mathematik Nr. 364, Oktober 2017 ISSN
More informationGravitational Radiation of Binaries Coalescence into Intermediate Mass Black Holes
Commun Theor Phys 57 (22) 56 6 Vol 57 No January 5 22 Gravitational Radiation of Binaries Coalescence into Intermediate Mass Black Holes LI Jin (Ó) HONG Yuan-Hong ( ) 2 and PAN Yu ( ) 3 College of Physics
More informationMCMC Sampling for Bayesian Inference using L1-type Priors
MÜNSTER MCMC Sampling for Bayesian Inference using L1-type Priors (what I do whenever the ill-posedness of EEG/MEG is just not frustrating enough!) AG Imaging Seminar Felix Lucka 26.06.2012 , MÜNSTER Sampling
More informationGravitational Waves: From Einstein to a New Science
Gravitational Waves: From Einstein to a New Science LIGO-G1602199 Barry C Barish Caltech - LIGO 1.3 Billion Years Ago 2 Black Holes Regions of space created by super dense matter from where nothing can
More informationBayesian methods in the search for gravitational waves
Bayesian methods in the search for gravitational waves Reinhard Prix Albert-Einstein-Institut Hannover Bayes forum Garching, Oct 7 2016 Statistics as applied Probability Theory Probability Theory: extends
More informationSearching for Binary Coalescences with Inspiral Templates: Detection and Parameter Estimation
Rochester Institute of Technology RIT Scholar Works Presentations and other scholarship 5-9-2008 Searching for Binary Coalescences with Inspiral Templates: Detection and Parameter Estimation Benjamin Farr
More informationGravitational-Wave Astronomy - a Long Time Coming Livia Conti, for the Virgo Collaboration Fred Raab, for the LIGO Scientific Collaboration
Gravitational-Wave Astronomy - a Long Time Coming Livia Conti, for the Virgo Collaboration Fred Raab, for the LIGO Scientific Collaboration LIGO Hanford, WA LIGO Livingston, LA Virgo (Cascina, Italy) What
More informationKey ideas on how inspiral-merger-ringdown waveforms are built within the effective-one-body formalism
Key ideas on how inspiral-merger-ringdown waveforms are built within the effective-one-body formalism Alessandra Buonanno Maryland Center for Fundamental Physics & Joint Space-Science Institute Department
More informationGravitational Wave Memory Revisited:
Gravitational Wave Memory Revisited: Memories from the merger and recoil Marc Favata Kavli Institute for Theoretical Physics Metals have memory too What is the GW memory? Generally think of GW s as oscillating
More informationClassification CE-717: Machine Learning Sharif University of Technology. M. Soleymani Fall 2012
Classification CE-717: Machine Learning Sharif University of Technology M. Soleymani Fall 2012 Topics Discriminant functions Logistic regression Perceptron Generative models Generative vs. discriminative
More informationKicked Waveforms Observing Black Hole Recoils in Gravitational Wave Signals
Kicked Waveforms Observing Black Hole Recoils in Gravitational Wave Signals Christopher Moore, DAMTP, Cambridge, UK StronG BaD, Mississippi 1st March 2017 Work done in collaboration with Davide Gerosa
More informationGRAVITATIONAL-WAVE ASTRONOMY WITH COALESCING COMPACT BINARIES: DETECTION AND PARAMETER ESTIMATION WITH ADVANCED DETECTORS
GRAVITATIONAL-WAVE ASTRONOMY WITH COALESCING COMPACT BINARIES: DETECTION AND PARAMETER ESTIMATION WITH ADVANCED DETECTORS RORY JAMES EDWIN SMITH A thesis submitted to the University of Birmingham for the
More informationCS168: The Modern Algorithmic Toolbox Lecture #6: Regularization
CS168: The Modern Algorithmic Toolbox Lecture #6: Regularization Tim Roughgarden & Gregory Valiant April 18, 2018 1 The Context and Intuition behind Regularization Given a dataset, and some class of models
More information