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.interp2d that computes g(n,e) from Peters (1964). The code assumes that the function returns the output sorted as with the interp2d returned functions (and thus unsorts). Default is None and uses exact g(n,e) in this case.

Returns

h_0float/array

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