h_0_n

legwork.strain.h_0_n(m_c, f_orb, ecc, n, dist, position=None, polarisation=None, inclination=None, interpolated_g=None)

Computes strain amplitude

Computes the dimensionless power of a general binary radiating gravitational waves in the quadrupole approximation at n harmonics of the orbital frequency

In the docs below, x refers to the number of sources, y to the number of timesteps and z to the number of harmonics.

Parameters:

m_cfloat/array

Chirp mass of each binary. Shape should be (x,).

f_orbfloat/array

Orbital frequency of each binary at each timestep. Shape should be (x, y), or (x,) if only one timestep.

eccfloat/array

Eccentricity of each binary at each timestep. Shape should be (x, y), or (x,) if only one timestep.

nint/array

Harmonic(s) at which to calculate the strain. Either a single int or shape should be (z,)

distfloat/array

Distance to each binary. Shape should be (x,)

positionSkyCoord/array, optional

Sky position of source. Must be specified using Astropy’s astropy.coordinates.SkyCoord class.

polarisationfloat/array, optional

GW polarisation angle of the source. Must have astropy angular units.

inclinationfloat/array, optional

Inclination of the source. Must have astropy angular units.

interpolated_gfunction

A function returned by scipy.interpolate.RectBivariateSpline that computes g(n,e) from Peters (1964). Default is None and uses exact g(n,e) in this case.

Returns:

h_0float/array

Strain amplitude. Shape is (x, y, z).