h_c_n

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

Computes characteristic strain amplitude

Computes the dimensionless characteristic 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_cfloat/array

Characteristic strain. Shape is (x, y, z).