def PNT2TidalOcto_Tv16(XA, beta0PNT=0, beta1PNT=0):
    """ TaylorT2 1PN Octopolar Tidal Coefficient, v^16 Timing 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 33440 / 21 + 995 * beta0PNT / 168 + beta1PNT + (24670 / 7 + 17 *
        beta0PNT / 3) * XA - 17 / 6 * (1825 + 2 * beta0PNT
        ) * XATo2nd + 325 * XATo3rd / 6
