Let $K({\bf r}_1,{\bf r}_2)= e^{-\frac12|{\bf r}_1-{\bf r}_2|^2}$ be the translation and rotation invariant Gaussian kernel. I want to compute the trace Tr$\hspace{1pt}K^3$ over the $n$-dimensional ball $B_n$ ($n=2$ is the disc) of radius $R$, that is $${\rm Tr}\hspace{1pt}K^3 = \int_{B_n^3} d^n{\bf r}_1 d^n{\bf r}_2 d^n{\bf r}_3\, K({\bf r}_1,{\bf r}_2)K({\bf r}_2,{\bf r}_3)K({\bf r}_3,{\bf r}_1)\,.$$ In particular, I want to obtain the expansion at $R\rightarrow\infty$ up to the $O(1)$ term included.
I was wondering if a closed formula for the Gaussian kernel existed? I can perform one radial integral but after that it seems difficult, especially in $n>2$ dimensions. Symmetry allows to reduce the number of angular integrals, but it still remain quite complicated.