def PNT2TidalOcto_Pv16(XA, beta0PNT=0, beta1PNT=0):
    """ TaylorT2 1PN Octopolar Tidal Coefficient, v^16 Phasing Term. 

   XA = mass fraction of object
   beta0PNT = 0PN Octopole Tidal Flux coefficient
   beta1PNT = 1PN Octopole Tidal Flux coefficient """
    XATo2nd = XA * XA
    XATo3rd = XATo2nd * XA
    XATo4th = XATo3rd * XA
    XATo5th = XATo4th * XA
    return 5 * (267520 + 995 * beta0PNT + 168 * beta1PNT) / 1848 + 5 / 231 * (
        74010 + 119 * beta0PNT) * XA - 85 / 66 * (1825 + 2 * beta0PNT
        ) * XATo2nd + 1625 * XATo3rd / 66
