Snapshots of Extreme Mass Ratio Inspirals
Steve Drasco
and
Scott Hughes
Waveform Generation Program
Here is a short program that can be run in Matlab (or Octave, which is free) to turn data files into waveforms
and fluxes:
For instructions on how to use the program, do "help snapshot" when you are in the same directory
as the file snapshot.m. Here is what the help file says:
Data Files
You will need data files before you can make waveforms. The data files below were generated
with a code described in
Phys. Rev. D 73 024027 (2006), or equivalently
gr-qc/0509101.
In the table of links below, you will find all the data used to produce tables III-VII in the Phys. Rev. D version of our paper.
For these links, e is eccentricity, and theta_inc is the inclination angle.
|
IMPORTANT: Use of this data set requires the following acknowledgement:
"The supercomputers used in this investigation were provided by funding
from the JPL Office of the Chief Information Officer."
|
| Prograde orbits with spin a = 0.9M and semilatus rectum p = 6.
|
| e = 0.1, theta_inc = 20 deg
| e = 0.1, theta_inc = 40 deg
| e = 0.1, theta_inc = 60 deg
| e = 0.1, theta_inc = 80 deg
|
| e = 0.3, theta_inc = 20 deg
| e = 0.3, theta_inc = 40 deg
| e = 0.3, theta_inc = 60 deg
| e = 0.3, theta_inc = 80 deg
|
| e = 0.5, theta_inc = 20 deg
| e = 0.5, theta_inc = 40 deg
| e = 0.5, theta_inc = 60 deg
| e = 0.5, theta_inc = 80 deg
|
| e = 0.7, theta_inc = 20 deg
| e = 0.7, theta_inc = 40 deg
| e = 0.7, theta_inc = 60 deg
| e = 0.7, theta_inc = 80 deg
|
| Retrograde orbits with spin a = 0.9M and semilatus rectum p = 12.
|
| e = 0.1, theta_inc = 100 deg
| e = 0.1, theta_inc = 120 deg
| e = 0.1, theta_inc = 140 deg
| e = 0.1, theta_inc = 160 deg
|
| e = 0.3, theta_inc = 100 deg
| e = 0.3, theta_inc = 120 deg
| e = 0.3, theta_inc = 140 deg
| e = 0.3, theta_inc = 160 deg
|
| e = 0.5, theta_inc = 100 deg
| e = 0.5, theta_inc = 120 deg
| e = 0.5, theta_inc = 140 deg
| e = 0.5, theta_inc = 160 deg
|
| e = 0.7, theta_inc = 100 deg
| e = 0.7, theta_inc = 120 deg
| e = 0.7, theta_inc = 140 deg
| e = 0.7, theta_inc = 160 deg
|
Here are some other data sets that people have requested.
-
Circular-equatorial Kerr (CircEq.tar.gz)
[63 snapshots, 1 MB gzipped, about 36 hours CPU-time, requested by Chao Li (Caltech), posted 8 March 2006]
For spins of a/M = 0.891, 0.9, and 0.909, and for 21 radii ranging from r/r_ISCO
of 1.001 to 10, where r_ISCO is the radius of the Innermost Stable Circular Orbit.
See file headers and the included README for details.
-
Circular Schwarzschild (CircSchw.tar.gz)
[20 snapshots, 154 KB gzipped, about 2 hours CPU-time, requested by Curt Cutler (JPL), posted 8 March 2006]
Twenty snapshots with radii ranging from r = 6.1M to r = 15M.
See file headers and the included README for details. [29 Sept 2006, added Ain_{lmkn} to each line]
-
generic (generic_0.8.tar.gz)
[1,742 snapshots, 1.4 GB gzipped, about 2,900 hours CPU-time, requested by Jonathan Gair (Cambridge, UK), posted 19 June 2006]
Two thousand snapshots with black hole spin a/M = 0.8, eccentricity from 0.1 to 0.5, inclination from 20 degrees to 70 degrees,
semilatus rectum from 1.2 to 3 times the value of that for the innermost stable circular orbit, and a requested total flux accuracy of 1%.
See file headers and the included README for details.
IMPORTANT: Use of this data set requires the following acknowledgement:
"The supercomputers used in this investigation were provided by funding
from the JPL Office of the Chief Information Officer."
-
Table of fluxes for generic orbits (FluxSummary_0.8)
[1,742 line table, 190 kb, trivial CPU time used to post-process generic_0.8.tar.gz data set, requested by Jonathan Gair (Cambridge, UK), posted 26 June 2006]
A single file that tabulates orbital parameters and total radiative fluxes for the waveforms from the preceeding data set (generic_0.8.tar.gz).
The fluxes should be accurate to 1% or better.
See description of generic_0.8.tar for the range of parameters used, etc. See 3-line header of file for explanation of format.
IMPORTANT: Use of this data set requires the following acknowledgement:
"The supercomputers used in this investigation were provided by funding
from the JPL Office of the Chief Information Officer."
-
equatorial (equatorial.tar.gz)
[52 snapshots, 21 MB gzipped, about 126 hours CPU-time, requested by Pranesh Sundararajan (MIT), posted 30 Nov 2006]
Equatorial snapshots with black hole spin a/M = [0.1 0.3 0.5 0.7], eccentricity = [0.1 0.3 0.5 0.7], inclination 0 degrees,
semilatus rectum = [1.1 1.3 1.5 1.7] times the value of that for the innermost stable circular orbit (note value of 1.1 used only for case of eccentricity = 0.1), and a requested total flux accuracy of 1e-6.
See file headers and the included README for details.
IMPORTANT: Use of this data set requires the following acknowledgement:
"The supercomputers used in this investigation were provided by funding
from the JPL Office of the Chief Information Officer."
- Horizon Skimming (HorizonSkimming.tar.gz)
[130 snapshots, 337 MB gzipped, about 1,300 hours CPU-time, requested by Enrico Barausse (SISSA, Trieste), posted 15 Oct 2006]
Snapshots of generalized versions of Hughes' nearly horizon skimming orbits described in
Phys. Rev. D 63, 064016 (2001)
(gr-qc/0101023). The orbital parameters used are:
a/M = 0.998,
eccentricity = [0.1 0.2 0.3 0.4 0.5],
semilatus rectum = [1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.0],
inclination = [0 8 16 24 32 40] degrees,
and a requested flux accuracy is 1e-3.
See file headers and the included README for details.
IMPORTANT: Use of this data set requires the following acknowledgement:
"The supercomputers used in this investigation were provided by funding
from the JPL Office of the Chief Information Officer."
- for time domain comparisons (FDTD.tar.gz)
[24 snapshots, 23 MB gzipped, about 75 hours CPU-time, requested by Pranesh Sundararajan (MIT), posted 9 Aug 2007]
Miscellaneous collection of snapshots (16 circular, and 8 generic) to be used in comparions with a time domain code.
The orbital parameters used are:
a/M = [0.5 0.9],
eccentricity = [1e-4 0.3 0.7],
semilatus rectum = [6 8],
inclination = [40 45 60] degrees,
and a requested flux accuracy is [1e-3 1e-4].
See file headers for details.
IMPORTANT: Use of this data set requires the following acknowledgement:
"The supercomputers used in this investigation were provided by funding
from the JPL Office of the Chief Information Officer."
If you would like to request your own custom data, write the orbital parameters
(a, eccentricity, inclination, semilatus rectum, desired flux accuracy) on the back of
a twenty dollar bill, and send it to Steve.
This site was last updated on 26 Feb 2008, by Steve Drasco.