I am trying to rederive the results presented in the paper, in particular equation (30). That is, I am trying to compute the correction to the ground-state energy of a dipolar condensate due to beyond-mean field effects. This involves evaluating the following integral:
$$\Delta E = \frac{1}{2} V \int \frac{d^3q}{(2\pi)^3} \left[\varepsilon_{\textbf{q}} - \frac{\hbar^2\textbf{q}^2}{2M} - nV_{\text{int}}(\textbf{q}) + \frac{V_{\text{int}}(\textbf{q})^2}{q^2}\frac{2Mn^2}{\hbar^2} \right] $$
where $\varepsilon_{\textbf{q}} = \sqrt{\frac{\hbar^2\textbf{q}^2}{2M}\left(\frac{\hbar^2\textbf{q}^2}{2M} + 2gn_0\left[1 + \epsilon_{\text{dd}}(3\cos^2{\theta} - 1) \right] \right)}$ is the Boglibov spectrum and $V_{\text{int}}(\textbf{q}) = g\left[1 + \epsilon_{\text{dd}}(3\cos^2{\theta} - 1)\right]$ is the Fourier Transform of the dipole-dipole interaction potential.
The integral is the inverse Fourier transform and $\theta$ is the Fourier polar angle. Taking $n_0 = n$, the paper finds the result to be
$$\Delta E = V \frac{2\pi\hbar^2 a_s n^2}{M} \frac{128}{15} \sqrt{\frac{a_s^3n}{\pi}}\mathcal{Q}_5(\epsilon)$$
where they have reexpressed the result in terms of $a_s$ by substituting $g = 4\pi\hbar^2a_s/M$ . Here,
$$\mathcal{Q}_5(x) = (1-x)^{5/2} {}_2F_1(-\frac{l}{2},\frac{1}{2};\frac{3}{2};\frac{3x}{x-1})$$
where ${}_2F_1(\alpha,\beta;\gamma;z)$ represents the hypergeometric function.
The last term in the integral is supposed to remove the divergence, however when I try to evaluate the integral I get still get a divergent result.
I would appreciate any attempts to explain how they reach their result or how to go about dealing with the divergence.
Thank you.
N.B: The argument in the final result should be $\epsilon_\text{dd}$ but for some reason the TeX will not compile properly.