Skip to content

TaylorF2

This file implements the TaylorF2 waveform, as described in the LALSuite library.

Functions:

Name Description
gen_TaylorF2

Generate TaylorF2 frequency domain waveform

gen_TaylorF2_hphc

Generate PhenomD frequency domain waveform following 1508.07253.

get_PNPhasing_F2

Gets dictionaries giving the phasing coefficients to be used in the approximant.

gen_TaylorF2(f: Float[Array, ' n_freq'], params: Float[Array, 9], f_ref: float, use_lambda_tildes: bool = True) -> Complex[Array, ' n_freq'] ¤

Generate TaylorF2 frequency domain waveform

vars array contains both intrinsic and extrinsic variables theta = [Mchirp, eta, chi1, chi2, lambda1, lambda1, D, tc, phic] Mchirp: Chirp mass of the system [solar masses] eta: Symmetric mass ratio [between 0.0 and 0.25] chi1: Dimensionless aligned spin of the primary object [between -1 and 1] chi2: Dimensionless aligned spin of the secondary object [between -1 and 1] lambda tilde: Dimensionless tidal deformability first object [between 0 and 5000] delta lamda tilde: Dimensionless tidal deformability second object [between 0 and 5000] D: Luminosity distance to source [Mpc] tc: Time of coalesence. This only appears as an overall linear in f contribution to the phase phic: Phase of coalesence

f_ref: Reference frequency for the waveform

Returns:

Name Type Description
h0 array

Strain

gen_TaylorF2_hphc(f: Float[Array, ' n_freq'], params: Float[Array, 10], f_ref: float, use_lambda_tildes: bool = True) -> tuple[Complex[Array, ' n_freq'], Complex[Array, ' n_freq']] ¤

Generate PhenomD frequency domain waveform following 1508.07253. vars array contains both intrinsic and extrinsic variables theta = [Mchirp, eta, chi1, chi2, D, tc, phic] Mchirp: Chirp mass of the system [solar masses] eta: Symmetric mass ratio [between 0.0 and 0.25] chi1: Dimensionless aligned spin of the primary object [between -1 and 1] chi2: Dimensionless aligned spin of the secondary object [between -1 and 1] lambda1: Dimensionless tidal deformability of the primary object [between 0 and 5000] lambda2: Dimensionless tidal deformability of the secondary object [between 0 and 5000] D: Luminosity distance to source [Mpc] tc: Time of coalesence. This only appears as an overall linear in f contribution to the phase phic: Phase of coalesence inclination: Inclination angle of the binary [between 0 and PI]

f_ref: Reference frequency for the waveform

Returns:

Name Type Description
hp array

Strain of the plus polarization

hc array

Strain of the cross polarization

get_PNPhasing_F2(m1: float, m2: float, S1z: float, S2z: float, lambda1: float, lambda2: float) -> tuple[dict[str, float], dict[str, float]] ¤

Gets dictionaries giving the phasing coefficients to be used in the approximant. Keys are the different PN orders, with values being the corresponding coefficient. This follows the implementation of XLALSimInspiralPNPhasing_F2 from lalsuite.

Parameters:

Name Type Description Default
m1 float

Mass of first (heavier) object

required
m2 float

Mass of second (lighter) object

required
S1z float

z-component of spin of first object

required
S2z float

z-component of spin of second object

required
lambda1 float

Tidal deformability first object

required
lambda2 float

Tidal deformability first object

required

Returns:

Type Description
tuple[dict[str, float], dict[str, float]]

tuple[dict, dict]: phasing_coeffs, phasing_log_coeffs as defined in the LAL source code, coefficients for various PN orders.