We report progress in the calculation of the thermal interquark potential of bottomonium using the HAL QCD method applied to bottom quarks in the non-relativistic approximation (i.e. NRQCD). We exploit the fast Fourier transform algorithm, using a momentum space representation, to efficiently calculate NRQCD correlation functions of non-local mesonic S-wave states, and thus obtain the potential for temperatures in both the hadronic and plasma phases. This work was
performed on our anisotropic 2+1 flavour “Generation 2" FASTSUM ensembles.