<p>If I’m understanding everything about what you’re asking and how it relates to this other thing, I believe a kind of library you’re thinking of refers to a surrogate model.</p>
<p>The idea behind a surrogate model is that while post-Newtonian approximations can capture gravitational wave signals up to the ring down, that often time that isn’t enough, and so surrogate models are used as the ‘library’ to compare gravitational wave signals to for identifying the possible source.</p>
<p>The idea behind a surrogate model is that you can solve Einstein’s equations to get a pretty good looking waveform that has features like ringdown, but that just one run of the numerical model can take up to weeks, and so it’s not practical to run at every possible set of parameter points. Instead, people will run a sample of (not evenly spaced) points in the parameter space, and used sophisticated interpolation schemes and sometimes machine learning to try to fill in the gaps, creating something of an effectual library of gravitational waveforms in parameter space that can then be referenced against incoming signals.</p>
<p>As far as which surrogate model is used or even how popular they are amongst those who actually do the signal processing, I’m not sure, but I do know it’s a very active area of research that has had many papers about improving them come out in recent years.</p>
<p>Furthermore, as to your main question, a lot of times (using whatever model it might be, surrogate, PN, etc) these libraries serve as a primary method for finding these waveforms in the noise, because my (limited) understanding is that it is much easier to check if a given signal is present in data than to find a previously ungenerated waveform, and that once these models are construed, it's pretty quick to compare a signal against one.</p>
<p>I say all of this with some uncertainty, as working with these was only a brief side project a year ago, so if anyone sees I've misrepresented anything, please let me know.</p>

CC BY-SA 4.0
2022-03-25T03:14:18.493Z
Thank you for your answer! My term "library" refers to their "dense bank of templates" so I'm really asking if the early identification of potential GW events from mergers *is done by searching template banks or not*; i.e. is this now how it's done. I didn't realize until now that the ringdown was such an important factor for near-real time flagging potential events, but including it certainly could reduce the rate of false positives.
uhoh
CC BY-SA 4.0
2022-03-25T03:17:13.673Z
btw [What does "GPU-accelerated butterfly matched filtering over dense bank of time-symmetric chirp-like templates" mean? (GW170817)](https://astronomy.stackexchange.com/q/28404/7982) currently has a bounty, and you may find [What can be learned from, or noted in this LIGO Orrery video?](https://astronomy.stackexchange.com/q/28610/7982) interesting as well.
uhoh
CC BY-SA 4.0
2022-03-26T00:51:17.560Z
@uhoh Rather conveniently, I was able to attend a guest lecture series from Dr. Nicolas Yunes on the matter today and it does seem seem that looking for a specific template is the principal method for determining a match, and he referenced a theorem for stating that this can always be done given a signal, although he did not reference the name of the theorem
Justin T
CC BY-SA 4.0
2022-03-26T00:54:35.040Z
He also mentioned that while surrogate models work great where they exist, they actually don’t really exist in high mass ratios currently, and so perturbative theory based templates like those of post-Newtonian seem to be the principal method for determining those methods, and then using a Markov Chain Monte Carlo algorithm, they find the parameters that best match the signal from the templates
Justin T
CC BY-SA 4.0
2022-03-26T01:03:16.327Z
*Excellent!* That's great to hear, thanks for update and insight and for sharing the anecdote. I don't know what the theorem is, but it might be simply something about "spanning the space" of possible functions, or related to proofs for (some) wavelet transforms being (over-)complete sets. I'm just using my "math words" without any specific expertise though. :-)
uhoh
CC BY-SA 4.0
2022-03-25T00:50:17.527Z
2022-03-25T00:50:17.527Z
Justin T
<p>LIGO/Virgo has multiple detection pipelines. Several of them (GstLAL, MBTA, PyCBC Live and SPIIR) are &quot;modelled&quot; searches which use large grids of pre-computed models to compare the signal against. There are also &quot;unmodelled&quot; searches (cWB) which simply look for something &quot;bright&quot; that occurs in multiple detectors at the same time.</p>
<p>See <a href="https://emfollow.docs.ligo.org/userguide/analysis/searches.html" rel="nofollow noreferrer">https://emfollow.docs.ligo.org/userguide/analysis/searches.html</a> for further details</p>

CC BY-SA 4.0
2022-03-25T15:52:13.180Z
Is it possible to say if it's the library/template search that identifies most of them, or is it mostly a red herring?  "Are chirped gravitational wave events *generally first identified by searching through libraries of chirps?*"
uhoh
CC BY-SA 4.0
2022-03-25T16:20:02.727Z
@uhoh I think it's generally not true. Detection depends on many things, e.g.,  false alarm rate, signal to noise of the detectors, the properties of the source itself, etc... For example, see Table IV of the third catalog which shows that many sources are detected by many methods, but some only by others. https://arxiv.org/abs/2111.03606
Daddy Kropotkin
CC BY-SA 4.0
2022-03-25T16:23:05.523Z
@DaddyKropotkin then perhaps your link provides enough information to complete an answer to my question since I'm not asking for the whys or hows. *Thanks!*
uhoh
CC BY-SA 4.0
2022-03-25T09:03:14.673Z
2022-03-25T10:18:53.877Z
2022-03-25T10:18:53.877Z
Rob
<p>Comments below <a href="https://astronomy.stackexchange.com/q/28404/7982">What does &quot;GPU-accelerated butterfly matched filtering over dense bank of time-symmetric chirp-like templates&quot; mean? (GW170817)</a> suggest that for this technique a library of strain &quot;chirps&quot; is built, then in real time a fast computer (in this case GPU-assisted) <em>constantly compares</em> the the most recent signals from a gravitational wave detecting interferometer to all of them, looking for matches above a chosen matching threshold.</p>
<p>As long as the matching rate is not too high, these can be passed on to more sophisticated analyses to better decide if a chirped gravitational wave event might have been detected, and if so, to start figuring out where it was so that alerts can be generated and other instruments (e.g. optical and X-ray telescopes) can start looking for a signature for a source.</p>
<p><strong>Question:</strong> Are chirped gravitational wave events generally first identified by searching through libraries of chirps?</p>
<p>Is this how they are first identified <em>in general</em>, or does the paper cited in the linked question describe just one possible method, and library searches are not always a part of the process?</p>

<observational-astronomy><data-analysis><gravitational-waves><time-domain-astronomy>
Are chirped gravitational wave events generally first identified by searching through libraries of chirps?
uhoh
# Are chirped gravitational wave events generally first identified by searching through libraries of chirps?

(Asked by: username_0 on 2022-03-24)

Comments below [What does "GPU-accelerated butterfly matched filtering over dense bank of time-symmetric chirp-like templates" mean? (GW170817)](https://astronomy.stackexchange.com/q/28404/7982) suggest that for this technique a library of strain "chirps" is built, then in real time a fast computer (in this case GPU-assisted) *constantly compares* the the most recent signals from a gravitational wave detecting interferometer to all of them, looking for matches above a chosen matching threshold.

As long as the matching rate is not too high, these can be passed on to more sophisticated analyses to better decide if a chirped gravitational wave event might have been detected, and if so, to start figuring out where it was so that alerts can be generated and other instruments (e.g. optical and X-ray telescopes) can start looking for a signature for a source.

**Question:** Are chirped gravitational wave events generally first identified by searching through libraries of chirps?

Is this how they are first identified *in general*, or does the paper cited in the linked question describe just one possible method, and library searches are not always a part of the process?

**Tags:** observational-astronomy, data-analysis, gravitational-waves, time-domain-astronomy


---

## 2 answers

### Answer by username_1 on 2022-03-25. Score: 2 | Accepted answer

If I’m understanding everything about what you’re asking and how it relates to this other thing, I believe a kind of library you’re thinking of refers to a surrogate model.

The idea behind a surrogate model is that while post-Newtonian approximations can capture gravitational wave signals up to the ring down, that often time that isn’t enough, and so surrogate models are used as the ‘library’ to compare gravitational wave signals to for identifying the possible source.

The idea behind a surrogate model is that you can solve Einstein’s equations to get a pretty good looking waveform that has features like ringdown, but that just one run of the numerical model can take up to weeks, and so it’s not practical to run at every possible set of parameter points. Instead, people will run a sample of (not evenly spaced) points in the parameter space, and used sophisticated interpolation schemes and sometimes machine learning to try to fill in the gaps, creating something of an effectual library of gravitational waveforms in parameter space that can then be referenced against incoming signals.

As far as which surrogate model is used or even how popular they are amongst those who actually do the signal processing, I’m not sure, but I do know it’s a very active area of research that has had many papers about improving them come out in recent years.

Furthermore, as to your main question, a lot of times (using whatever model it might be, surrogate, PN, etc) these libraries serve as a primary method for finding these waveforms in the noise, because my (limited) understanding is that it is much easier to check if a given signal is present in data than to find a previously ungenerated waveform, and that once these models are construed, it's pretty quick to compare a signal against one.

I say all of this with some uncertainty, as working with these was only a brief side project a year ago, so if anyone sees I've misrepresented anything, please let me know.

> **Comments:**
> **username_0**: Thank you for your answer! My term "library" refers to their "dense bank of templates" so I'm really asking if the early identification of potential GW events from mergers *is done by searching template banks or not*; i.e. is this now how it's done. I didn't realize until now that the ringdown was such an important factor for near-real time flagging potential events, but including it certainly could reduce the rate of false positives.  
> **username_0**: btw [What does "GPU-accelerated butterfly matched filtering over dense bank of time-symmetric chirp-like templates" mean? (GW170817)](https://astronomy.stackexchange.com/q/28404/7982) currently has a bounty, and you may find [What can be learned from, or noted in this LIGO Orrery video?](https://astronomy.stackexchange.com/q/28610/7982) interesting as well.  
> **username_1**: @username_0 Rather conveniently, I was able to attend a guest lecture series from Dr. Nicolas Yunes on the matter today and it does seem seem that looking for a specific template is the principal method for determining a match, and he referenced a theorem for stating that this can always be done given a signal, although he did not reference the name of the theorem  
> **username_1**: He also mentioned that while surrogate models work great where they exist, they actually don’t really exist in high mass ratios currently, and so perturbative theory based templates like those of post-Newtonian seem to be the principal method for determining those methods, and then using a Markov Chain Monte Carlo algorithm, they find the parameters that best match the signal from the templates  
> **username_0**: *Excellent!* That's great to hear, thanks for update and insight and for sharing the anecdote. I don't know what the theorem is, but it might be simply something about "spanning the space" of possible functions, or related to proofs for (some) wavelet transforms being (over-)complete sets. I'm just using my "math words" without any specific expertise though. :-)  

---

### Answer by username_2 on 2022-03-25. Score: 3

LIGO/Virgo has multiple detection pipelines. Several of them (GstLAL, MBTA, PyCBC Live and SPIIR) are "modelled" searches which use large grids of pre-computed models to compare the signal against. There are also "unmodelled" searches (cWB) which simply look for something "bright" that occurs in multiple detectors at the same time.

See <https://emfollow.docs.ligo.org/userguide/analysis/searches.html> for further details

> **Comments:**
> **username_0**: Is it possible to say if it's the library/template search that identifies most of them, or is it mostly a red herring?  "Are chirped gravitational wave events *generally first identified by searching through libraries of chirps?*"  
> **username_3**: @username_0 I think it's generally not true. Detection depends on many things, e.g.,  false alarm rate, signal to noise of the detectors, the properties of the source itself, etc... For example, see Table IV of the third catalog which shows that many sources are detected by many methods, but some only by others. https://arxiv.org/abs/2111.03606  
> **username_0**: @DaddyKropotkin then perhaps your link provides enough information to complete an answer to my question since I'm not asking for the whys or hows. *Thanks!*  

---

