Skip to content

Commit

Permalink
More paper edits
Browse files Browse the repository at this point in the history
  • Loading branch information
sfield17 committed Apr 2, 2024
1 parent 180141b commit f7932a5
Showing 1 changed file with 11 additions and 3 deletions.
14 changes: 11 additions & 3 deletions tutorial/paper/paper.md
Original file line number Diff line number Diff line change
Expand Up @@ -34,14 +34,22 @@ bibliography: paper.bib

---

![Example gravitational wave prediction from a surrogate model compared with numerical relativity for a precessing binary black hole. This particular numerical relativity simulation took 70,881 CPU-hours (about 1.75 months using 56 cores on the supercomputer Frontera), while the surrogate model can be evaluated in about 100 milliseconds. \label{fig:gws}](gwsurrogate.png)
![Example gravitational wave prediction from a surrogate model compared with numerical relativity for a precessing binary black hole system. This particular numerical relativity simulation took 70,881 CPU-hours (about 1.75 months using 56 cores on the supercomputer Frontera), while the surrogate model can be evaluated in about 100 milliseconds. \label{fig:gws}](gwsurrogate.png)

# Summary


Gravitational waves are ripples in space-time caused by the motion of massive objects. One of the most astrophysically important sources of gravitational radiation is caused by two orbiting compact objects, such as black holes and neutron stars, that slowly inspiral and merge. The motion of these massive objects generates gravitational waves that radiate to the far field where gravitational-wave detectors can observe them. Complicated partial or ordinary differential equations govern the entire process. Traditionally, the dynamics of compact binary systems and the emitted gravitational waves have been computed by expensive simulation codes that can take days to months to run. A key simulation output is the gravitational wave signal for a particular set of parameter values describing the system, such as the black holes' masses and spins. The computed signal is required for a diverse range of multiple-query applications, such as template bank generation for searches, parameter estimation, mock data analysis, studies of model bias, and tests of general relativity, to name a few. In such settings, the high-fidelity signal computed from differential equations is often too slow to be directly used.
Gravitational waves are ripples in space-time caused by the motion of massive objects. One of the most astrophysically important sources of gravitational radiation is caused by two orbiting compact objects, such as black holes and neutron stars, that slowly inspiral and merge. The motion of these massive objects generates gravitational waves that radiate to the far field where gravitational-wave detectors can observe them. Complicated partial or ordinary differential equations govern the entire process.

Surrogate models offer a practical way to dramatically accelerate model evaluation while retaining the high-fidelity accuracy of the expensive simulation code; an example is shown in Fig.\ref{fig:gws}. Surrogate models can be constructed in various ways, but what separates these models from others is that they are primarily data-driven. Given a training set of gravitational waveform data sampling the parameter space, a model is built by following three steps: (i) Feature extraction: the waveform is decomposed into *data pieces* that are simple to model, (ii) Dimensionality reduction: each data piece is approximated by a low-dimensional vector space, which reduces the degrees of freedom we need to model, and (iii) regression techniques are applied to the low-dimensional representation of each data piece over the parameter space defining the model. These model-building steps result in an HDF5 file defining the surrogate model's data and structure, which is stored on Zenodo. The GWSurrogate code provides access to these HDF5 files through its catalog interface, and all available models and their HDF5 files can be found in `gwsurrogate.catalog.list()`. For a recent overview of surrogate modeling as used in gravitational wave astrophysics, please see Section 5 of @LISAConsortiumWaveformWorkingGroup:2023arg.
Traditionally, the dynamics of compact binary systems and the emitted gravitational waves have been computed by expensive simulation codes that can take days to months to run. A key simulation output is the gravitational wave signal for a particular set of parameter values describing the system, such as the black holes' masses and spins. The computed signal is required for a diverse range of multiple-query applications, such as template bank generation for searches, parameter estimation, mock data analysis, studies of model bias, and tests of general relativity, to name a few. In such settings, the high-fidelity signal computed from differential equations is often too slow to be directly used.

Surrogate models offer a practical way to dramatically accelerate model evaluation while retaining the high-fidelity accuracy of the expensive simulation code; an example is shown in Fig.\ref{fig:gws}. Surrogate models can be constructed in various ways, but what separates these models from other modeling frameworks is that they are primarily data-driven. Given a training set of gravitational waveform data sampling the parameter space, a model is built by following three steps:

1. Feature extraction: the waveform is decomposed into *data pieces* that are simple to model,
2. Dimensionality reduction: each data piece is approximated by a low-dimensional vector space, which reduces the degrees of freedom we need to model, and
3. Regression: fitting and regression techniques are applied to the low-dimensional representation of each data piece over the parameter space defining the model.

These model-building steps result in an HDF5 file defining the surrogate model's data and structure, which is stored on Zenodo. The GWSurrogate package provides access to these HDF5 files through its catalog interface, and all available models and their HDF5 files can be found in `gwsurrogate.catalog.list()`. For a recent overview of surrogate modeling as used in gravitational wave astrophysics, please see Section 5 of @LISAConsortiumWaveformWorkingGroup:2023arg.

The development of ``GWSurrogate`` is hosted on [GitHub](https://github.com/sxs-collaboration/gwsurrogate) and distributed through both [PyPI](https://pypi.org/project/gwsurrogate/) and [Conda](https://anaconda.org/conda-forge/gwsurrogate/). Quick start guides are found on the project's [homepage](https://github.com/sxs-collaboration/gwsurrogate) while model-specific documentation is described through a collection of model-specific [Jupyter notebooks](https://github.com/sxs-collaboration/gwsurrogate/tree/master/tutorial). Automated testing is run on [GitHub Actions](https://github.com/sxs-collaboration/gwsurrogate/actions).

Expand Down

0 comments on commit f7932a5

Please sign in to comment.