LALSimulation  5.4.0.1-fe68b98
LALSimIMRPhenomD.h File Reference

Tabulated Quasi-Normal Mode Information for Ringdown. More...

Prototypes

int IMRPhenomDPhaseFrequencySequence (REAL8Sequence *phases, REAL8Sequence *freqs, size_t ind_min, size_t ind_max, REAL8 m1, REAL8 m2, REAL8 chi1x, REAL8 chi1y, REAL8 chi1z, REAL8 chi2x, REAL8 chi2y, REAL8 chi2z, REAL8 Rholm, REAL8 Taulm, LALDict *extraParams)
 Helper function used in PhenomHM and PhenomPv3HM Returns the phenomD phase, with modified QNM. More...
 
int IMRPhenomDAmpFrequencySequence (REAL8Sequence *phases, REAL8Sequence *freqs, size_t ind_min, size_t ind_max, REAL8 m1, REAL8 m2, REAL8 chi1x, REAL8 chi1y, REAL8 chi1z, REAL8 chi2x, REAL8 chi2y, REAL8 chi2z)
 Helper function used in PhenomHM and PhenomPv3HM Returns the phenomD amplitude. More...
 
REAL8 IMRPhenomDComputet0 (REAL8 eta, REAL8 chi1z, REAL8 chi2z, REAL8 finspin, LALDict *extraParams)
 computes the time shift as the approximate time of the peak of the 22 mode. More...
 

Detailed Description

Tabulated Quasi-Normal Mode Information for Ringdown.

This file contains listed final spin values, and corresponding Quasi-Normal-Mode(QNM) central frequencies and decay rates. The spins are on [-1,1] where values are the dimensionless Kerr parameters S/M^2. The frequencies and decay rates are in units of 1/(s*pi*M), where M is 1. To compare values with tables given at, for example, www.phy.olemiss.edu/~berti/ringdown/, multiply by 2*pi, and note that the decay rate corresponds to the imaginary part of the QNM frequency.

Values for spins between -0.994 and 0.994 were sourced from the above website (2014), while qualitatively accurate values for the remaining spins where sourced from the implementation used in arxiv:1404:3197. Both references use the analytic representation of QNMs developed by Leaver in 1986, and for spin values within +-0.994, have identical values within numerical error.

Definition in file LALSimIMRPhenomD.h.

Go to the source code of this file.

Macros

#define f_CUT   0.2
 Dimensionless frequency (Mf) at which define the end of the waveform. More...
 
#define AMP_fJoin_INS   0.014
 Dimensionless frequency (Mf) at which the inspiral amplitude switches to the intermediate amplitude. More...
 
#define PHI_fJoin_INS   0.018
 Dimensionless frequency (Mf) at which the inspiral phase switches to the intermediate phase. More...
 
#define MIN_FINAL_SPIN   -0.717
 Minimal final spin value below which the waveform might behave pathological because the ISCO frequency is too low. More...
 
#define MAX_ALLOWED_MASS_RATIO   5000
 A large mass ratio causes memory over-runs. More...
 

Variables

static const int QNMData_length = 1003
 
static const double QNMData_a []
 
static const double QNMData_fring []
 
static const double QNMData_fdamp []
 

Macro Definition Documentation

◆ f_CUT

#define f_CUT   0.2

Dimensionless frequency (Mf) at which define the end of the waveform.

Definition at line 20 of file LALSimIMRPhenomD.h.

◆ AMP_fJoin_INS

#define AMP_fJoin_INS   0.014

Dimensionless frequency (Mf) at which the inspiral amplitude switches to the intermediate amplitude.

Definition at line 26 of file LALSimIMRPhenomD.h.

◆ PHI_fJoin_INS

#define PHI_fJoin_INS   0.018

Dimensionless frequency (Mf) at which the inspiral phase switches to the intermediate phase.

Definition at line 32 of file LALSimIMRPhenomD.h.

◆ MIN_FINAL_SPIN

#define MIN_FINAL_SPIN   -0.717

Minimal final spin value below which the waveform might behave pathological because the ISCO frequency is too low.

For more details, see the review wiki page https://www.lsc-group.phys.uwm.edu/ligovirgo/cbcnote/WaveformsReview/IMRPhenomDCodeReview/PhenD_LargeNegativeSpins

Definition at line 39 of file LALSimIMRPhenomD.h.

◆ MAX_ALLOWED_MASS_RATIO

#define MAX_ALLOWED_MASS_RATIO   5000

A large mass ratio causes memory over-runs.

We test and put the limit an order of magnitude above that of previous waveform models (which were around q=100).

Definition at line 45 of file LALSimIMRPhenomD.h.

Function Documentation

◆ IMRPhenomDPhaseFrequencySequence()

int IMRPhenomDPhaseFrequencySequence ( REAL8Sequence phases,
REAL8Sequence freqs,
size_t  ind_min,
size_t  ind_max,
REAL8  m1,
REAL8  m2,
REAL8  chi1x,
REAL8  chi1y,
REAL8  chi1z,
REAL8  chi2x,
REAL8  chi2y,
REAL8  chi2z,
REAL8  Rholm,
REAL8  Taulm,
LALDict *  extraParams 
)

Helper function used in PhenomHM and PhenomPv3HM Returns the phenomD phase, with modified QNM.

Parameters
[out]phasesphase evaluated at input freqs
freqsSequency of Geometric frequencies
ind_minstart index for frequency loop
ind_maxend index for frequency loop
m1mass of primary in solar masses
m2mass of secondary in solar masses
chi1xx-component of the dimensionless spin of object 1 w.r.t. Lhat = (0,0,1)
chi1yy-component of the dimensionless spin of object 1 w.r.t. Lhat = (0,0,1)
chi1zz-component of the dimensionless spin of object 1 w.r.t. Lhat = (0,0,1)
chi2xx-component of the dimensionless spin of object 2 w.r.t. Lhat = (0,0,1)
chi2yy-component of the dimensionless spin of object 2 w.r.t. Lhat = (0,0,1)
chi2zz-component of the dimensionless spin of object 2 w.r.t. Lhat = (0,0,1)
Rholmratio of ringdown frequencies f_RD_22/f_RD_lm
Taulmratio of ringdown damping times f_RM_22/f_RM_lm
extraParamslinked list containing the extra testing GR parameters

Definition at line 758 of file LALSimIMRPhenomD.c.

◆ IMRPhenomDAmpFrequencySequence()

int IMRPhenomDAmpFrequencySequence ( REAL8Sequence phases,
REAL8Sequence freqs,
size_t  ind_min,
size_t  ind_max,
REAL8  m1,
REAL8  m2,
REAL8  chi1x,
REAL8  chi1y,
REAL8  chi1z,
REAL8  chi2x,
REAL8  chi2y,
REAL8  chi2z 
)

Helper function used in PhenomHM and PhenomPv3HM Returns the phenomD amplitude.

Parameters
[out]phasesphase evaluated at input freqs
freqsSequency of Geometric frequencies
ind_minstart index for frequency loop
ind_maxend index for frequency loop
m1mass of primary in solar masses
m2mass of secondary in solar masses
chi1xx-component of the dimensionless spin of object 1 w.r.t. Lhat = (0,0,1)
chi1yy-component of the dimensionless spin of object 1 w.r.t. Lhat = (0,0,1)
chi1zz-component of the dimensionless spin of object 1 w.r.t. Lhat = (0,0,1)
chi2xx-component of the dimensionless spin of object 2 w.r.t. Lhat = (0,0,1)
chi2yy-component of the dimensionless spin of object 2 w.r.t. Lhat = (0,0,1)
chi2zz-component of the dimensionless spin of object 2 w.r.t. Lhat = (0,0,1)

Definition at line 821 of file LALSimIMRPhenomD.c.

◆ IMRPhenomDComputet0()

REAL8 IMRPhenomDComputet0 ( REAL8  eta,
REAL8  chi1z,
REAL8  chi2z,
REAL8  finspin,
LALDict *  extraParams 
)

computes the time shift as the approximate time of the peak of the 22 mode.

Parameters
etasymmetric mass-ratio
chi1zdimensionless aligned-spin of primary
chi2zdimensionless aligned-spin of secondary
finspinfinal spin
extraParamslinked list containing the extra testing GR parameters

Definition at line 899 of file LALSimIMRPhenomD.c.

Variable Documentation

◆ QNMData_length

const int QNMData_length = 1003
static

Definition at line 50 of file LALSimIMRPhenomD.h.

◆ QNMData_a

const double QNMData_a[]
static

Definition at line 51 of file LALSimIMRPhenomD.h.

◆ QNMData_fring

const double QNMData_fring[]
static

Definition at line 159 of file LALSimIMRPhenomD.h.

◆ QNMData_fdamp

const double QNMData_fdamp[]
static

Definition at line 325 of file LALSimIMRPhenomD.h.