Finite T topological Susceptibility with heavy Quarks

Axion cosmology needs the QCD topological susceptibility between 400 and 1100 MeV. In this range the bottom quark is inconvenient to include in lattice simulations, but not heavy enough to ignore. We estimate its effect on the susceptibility by computing the ratio of the 4-quark susceptibility and the 4+1-quark susceptibility in the caloron gas approximation. We do so by computing small-mass and large-mass expansions of the finite mass and temperature fluctuation determinant and connecting them with a Padé approximant.