pint.derived_quantities.omdot_to_mtot
- pint.derived_quantities.omdot_to_mtot(omdot: Unit('deg / yr'), pb: Unit('d'), e: Unit(dimensionless))[source]
Determine total mass from Post-Keplerian longitude of periastron precession rate omdot, assuming general relativity.
omdot (\(\dot \omega\)) is the relativistic advance of periastron. It relates to the total system mass (assuming GR). Can handle scalar or array inputs.
- Parameters:
omdot (astropy.units.Quantity) – relativistic advance of periastron
pb (astropy.units.Quantity) – Binary orbital period
e (astropy.units.Quantity or float) – orbital eccentricity
- Returns:
mtot – In
u.Msun- Return type:
- Raises:
astropy.units.UnitsError – If the input data are not appropriate quantities
TypeError – If the input data are not quantities
Notes
Inverts
\[\dot \omega = 3T_{\odot}^{2/3} \left(\frac{P_b}{2\pi}\right)^{-5/3} \frac{1}{1-e^2}(m_p+m_c)^{2/3}\]to calculate \(m_{\rm tot} = m_p + m_c\), with \(T_\odot = GM_\odot c^{-3}\).
More details in Timing Models. Also see [9].