For two distinguishable particles in a haromonic confinining potential in three dimensions, the imaginary time propagator is
(1)
The partition function is the trace of the propagator for
,
(2)![Z &= \operatorname{tr} K \\
&= \int d\mathbf{r}_1^3 \int d\mathbf{r}_1^3 \,
K(\mathbf{r}_1, \mathbf{r}_2, \mathbf{r}_1, \mathbf{r}_2; \beta\hbar) \\
&= \left(\frac{m\omega}{2\pi\hbar\sinh\beta\hbar\omega}\right)^3
\int d\mathbf{r}_1^3 \int d\mathbf{r}_1^3
\exp\left(-\frac{2m\omega(\cosh\beta\hbar\omega -1)(r_1^2 +r_2^2)}
{2\hbar\sinh\omega\tau}\right)\\
&= \left(2(\cosh\beta\hbar\omega-1)\right)^{-3} \\
&= \left[2\sinh\left(\frac{\hbar\omega}{2k_BT}\right)\right]^{-6}](_images/math/3c058e1b1d73cd00dc22e31d2c73d93f6eb7df95.png)
For identical particles, we need to symmetrize the states for fermions, or antisymmetrize the states for bosons. The trace of the permuted propagator is,

Next we change coordinates to
and
.
Then
and
,
and we find,
(3)![Z_P = \left[2\cosh\left(\frac{\hbar\omega}{2k_BT}\right)
\sinh\left(\frac{\hbar\omega}{2k_BT}\right)\right]^{-3}](_images/math/e49293d78c1fc5b3764c6912c6e893f575b7e577.png)