WaSABI

A package to generate gravitational waveforms from SF theory and SF-PN Hybrids.

Currently 0PA and 1PAT1 quasicircular, non-spinning models are available. As an example, the gravitational waveform for the $l=2,m=2$ mode is easily computed using:

inspiral = BinaryInspiral[<|"r0" -> 10, "M" -> 1, "\[Nu]" -> 1/10, "\[Phi]" -> 0|>, "Model" -> "1PAT1"];
tmax = inspiral["Duration"];
Plot[Evaluate[ReIm[inspiral["Waveform"][2, 2][t]]], {t, 0, tmax}]

We can also plot the orbital trajectory:

r0 = inspiral["Trajectory"]["r0"];
\[Phi] = inspiral["Trajectory"]["\[Phi]"];
ParametricPlot[{r0[t] Cos[\[Phi][t]], r0[t] Sin[\[Phi][t]]}, {t, 0, tmax}]

Citing

In addition to acknowledging the Black Hole Perturbation Toolkit as suggested on the front page we also recommend citing the specific package version you use via the citation information on the package’s Zenodo page linked from the above DOI.

Authors and contributors

Barry Wardell, Josh Mathews, Loïc Honet